Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational physics”

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,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Application of Quantum Computing to Biochemical Systems: A Look to the Future

Chemistry is considered as one of the more promising applications to science of near-term quantum computing. Recent work in transitioning classical algorithms to a quantum computer has led to great strides in improving quantum algorithms and illustrating their quantum advantage. Because of the limitations of near-term quantum computers, the most effective strategies split the work over classical and quantum computers. There is a proven set of methods in computational chemistry and materials physics that has used this same idea of splitting a complex physical system into parts that are treated at different levels of theory to obtain solutions for the complete physical system for which a brute force solution with a single method is not feasible. These methods are variously known as embedding, multi-scale, and fragment techniques and methods. We review these methods and then propose the embedding approach as a method for describing complex biochemical systems, with the parts not only treated with different levels of theory, but computed with hybrid classical and quantum algorithms. Such strategies are critical if one wants to expand the focus to biochemical molecules that contain active regions that cannot be properly explained with traditional algorithms on classical computers. While we do not solve this problem here, we provide an overview of where the field is going to enable such problems to be tackled in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

PPPL Laboratory Directed Research and Development (Project Final Reports, FY2018 - FY2020)

The U.S. Department of Energy’s (DOE) Princeton Plasma Physics Laboratory (PPPL) is a collaborative national center for fusion energy science, basic sciences, and advanced technology. The Laboratory has three major missions: (1) to develop the scientific knowledge and advanced engineering to enable fusion to power the U.S. and the world; (2) to advance the science of nanoscale fabrication for technologies of tomorrow; and (3) to further the development of the scientific understanding of the plasma universe from laboratory to astrophysical scales. PPPL’s Laboratory Directed Research and Development (LDRD) program supports and encourages creativity and innovation and contributes to its long-term viability. New scientific and technical research areas emerge and are nurtured through the program. Furthermore, new capabilities are developed to enable the Laboratory to meet its and DOE’s missions. The program is used to systematically diversify the Laboratory’s programs and mission. In the last few years, the program has started projects in nanomaterial synthesis, microelectronics, advanced x-ray spectroscopy, high-energy-density physics, superconducting magnet technology, machine learning and artificial intelligence, 3D magnetic fields to optimize fusion plasmas, integration of permanent magnets with simple high-field magnets to reduce the cost of producing complex 3D magnetic fields, advanced computational methods for predictive understanding and control of fusion plasma, development of quantum computing algorithms for plasma physics, liquid metal plasma-facing components for fusion reactors, virtual engineering, and plasma-based space propulsion. The program is also the vehicle to recruit and train talented scientists and engineers with the new skills needed to perform the Laboratory’s mission. Many of the new hires through the program go on to become world-class scientists and engineers in their fields. This report provides descriptions and accomplishments of those LDRD projects that were completed during fiscal years 2018 through 2020.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Global optimization algorithms to compute thermodynamic equilibria in large complex systems with performance considerations

Several global optimization methods are reviewed that attempt to ensure that the integral Gibbs energy of a closed isothermal isobaric system is a global minimum to satisfy the necessary and sufficient conditions for thermodynamic equilibrium. In particular, the integral Gibbs energy function of a multicomponent system containing non-ideal phases may be highly non-linear and non-convex, which makes finding a global minimum a challenge. Consequently, a poor numerical approach may lead one to the false belief of equilibrium. Furthermore, confirming that one reaches a global minimum and that this is achieved with satisfactory computational performance becomes increasingly more challenging in systems containing many chemical elements and a correspondingly large number of species and phases. Several numerical methods that have been used for this specific purpose are reviewed with a benchmark study of three of the more promising methods using five case studies of varying complexity. A modification of the conventional Branch and Bound method is presented that is well suited to a wide array of thermodynamic applications, including complex phases with many constituents and sublattices, and ionic phases that must adhere to charge neutrality constraints. Also, a novel method is presented that efficiently solves the system of linear equations that exploits the unique structure of the Hessian matrix, which reduces the calculation from a O(N 3 ) operation to a O(N) operation. As a result, this combined approach demonstrates efficiency, reliability and capabilities that are favorable for integration of thermodynamic computations into multi-physics codes with inherent performance considerations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Data distribution service-based interoperability framework for smart grid testbed infrastructure

This study presents the design and implementation of a communication and control infrastructure for smart grid operation. The proposed infrastructure enhances the reliability of the measurements and control network. The advantages of utilizing the data-centric over message-centric communication approach are discussed in the context of smart grid applications. The data distribution service (DDS) is used to implement a data-centric common data bus for the smart grid. This common data bus improves the communication reliability, enabling distributed control and smart load management. These enhancements are achieved by avoiding a single point of failure while enabling peer-to-peer communication and an automatic discovery feature for dynamic participating nodes. The infrastructure and ideas presented in this paper were implemented and tested on the smart grid testbed. A toolbox and application programing interface for the testbed infrastructure are developed in order to facilitate interoperability and remote access to the testbed. This interface allows control, monitoring, and performing of experiments remotely. Furthermore, it could be used to integrate multidisciplinary testbeds to study complex cyber-physical systems (CPS).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

FPGAs as a Service to Accelerate Machine Learning Inference [PowerPoint]

Large-scale particle physics experiments face challenging demands for high-throughput computing resources both now and in the future. New heterogeneous computing paradigms on dedicated hardware with increased parallelization, such as Field Programmable Gate Arrays (FPGAs), offer exciting solutions with large potential gains. The growing applications of machine learning algorithms in particle physics for simulation, reconstruction, and analysis are naturally deployed on such platforms. We demonstrate that the acceleration of machine learning inference as a web service represents a heterogeneous computing solution for particle physics experiments that requires minimal modification to the current computing model. As examples, we retrain the ResNet-50 convolutional neural network to demonstrate state-of-the-art performance for top quark jet tagging at the LHC, and apply a ResNet-50 model with transfer learning for neutrino event classification. Using Project Brainwave b y Microsoft to accelerate the ResNet-50 image classification model, we achieve average inference times of 60 (10) milliseconds with our experimental physics software framework using Brainwave as a cloud (edge or on-premises) service, representing an improvement by a factor of approximately 30 (175) in model inference latency over traditional CPU inference in current experimental hardware. A single FPGA service accessed by many CPUs achieves a throughput of 600--700 inferences per second using an image batch of one, comparable to large batch-size GPU throughput and significantly better than small batch-size GPU throughput. Deployed as an edge or cloud service for the particle physics computing model, coprocessor accelerators can have a higher duty cycle and are potentially much more cost-effective.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Benchmarking the performance of uncertainty quantification methods for neural network-based interatomic potentials

Machine-learned interatomic potentials (ML-IAPs) continue to gain popularity as accurate, computationally efficient replacements for traditional, physics-based interatomic potentials and expensive ab initio methods. Uncertainty quantification (UQ) of ML-IAPs is a growing area of research as UQ is critical in many applications of IAPs, such as developing curated datasets, active learning-based data augmentation, self-improving models, and estimating the uncertainty of molecular dynamics simulations. In this paper, we construct and benchmark a series of different neural network potentials (NNPs) with varying network architectures to determine the performance of these models with respect to both the mean and uncertainty calibration error. Each NNP method is specifically designed to predict either epistemic or aleatoric uncertainty with particular focus on the differences in behavior between the epistemic and aleatoric uncertainty estimates. We benchmark these methods using multiple datasets common in the ML-IAP literature. The results show that the aleatoric uncertainty from single-shot model architectures is a competitive alternative to ensemble-based epistemic uncertainty predictions in regions of sufficient data-density. However, in regions where the representative data is sparse, aleatoric uncertainty models tend to overpredict and epistemic methods tend to underpredict the actual model error. We conclude that the type of UQ is crucial when discussing performance of probabilistic model results as different methods have different performance characteristics depending on the regime in which they are evaluated. Therefore, the type of UQ method should be carefully evaluated against both the data characteristics and requirements for the intended application.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Embedding hard physical constraints in neural network coarse-graining of three-dimensional turbulence

In recent years, deep learning approaches have shown much promise in modeling complex systems in the physical sciences. A major challenge in deep learning of partial differential equations is enforcing physical constraints and boundary conditions. In this work, we propose a general framework to directly embed the notion of an incompressible fluid into convolutional neural networks, and apply this to coarse-graining of turbulent flow. These physics-embedded neural networks leverage interpretable strategies from numerical methods and computational fluid dynamics to enforce physical laws and boundary conditions by taking advantage the mathematical properties of the underlying equations. Here, we demonstrate results on three-dimensional fully developed turbulence, showing that this technique drastically improves local conservation of mass, without sacrificing performance according to several other metrics characterizing the fluid flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Complexity and near extremal charged black branes

Highlights: • Trying to explore a remedy to overcome the puzzle in holographic complexity by assuming a cutoff behind the horizon and the corresponding boundary term. • Holographic complexity growth obtained by the WDW patch on shell action is insensitive to the UV cutoff at the late time. Therefore one may conclude that either on shell action does not compute complexity, or the late time is not given by the physical charges, or in our computations we have missed something! • We find that one needs to carefully compute on shell action taking into account the behind the horizon and the corresponding boundary term. Doing so, one arrives at a consistent picture. • We check this consistency by evaluating the late time complexity growth for near horizon limit of the near extremal black hole using two different methods. We would have obtained non-consistent picture, if we had not considered the behind the horizon cutoff as well as the corresponding boundary term. We compute holographic complexity of charged black brane solutions in arbitrary dimensions for the near horizon limit of near extremal case using two different methods. The corresponding complexity may be obtained either by taking the limit from the complexity of the charged black brane, or by computing the complexity for near horizon limit of near extremal solution. One observes that these results coincide if one assumes to have a cutoff behind horizon whose value is fixed by UV cutoff and also taking into account a proper counterterm evaluated on this cutoff. We also consider the situation for Vaidya charged black branes too.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Proton Isovector Helicity Distribution on the Lattice at Physical Pion Mass

We present a state-of-the-art calculation of the isovector quark helicity Bjorken-x distribution in the proton using lattice-QCD ensembles at the physical pion mass. We compute quasi-distributions at proton momenta P z ϵ {2.2; 2.6; 3.0} GeV on the lattice, and match them systematically to the physical parton distribution using large-momentum effective theory (LaMET).We reach an unprecedented precision through high statistics in simulations, large-momentum proton matrix elements, and control of excited-state contamination. The resulting distribution with combined statistical and systematic errors is in agreement with the latest phenomenological analysis of the spin-dependent experimental data; in particular, Δ$\bar{u}$(x) > Δ$\bar{d}$(x).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

How Data Becomes Physics: Inside the RACF

The RHIC & ATLAS Computing Facility (RACF) at the U.S. Department of Energy’s (DOE) Brookhaven National Laboratory sits at the center of a global computing network. It connects more than 2,500 researchers around the world with the data generated by millions of particle collisions taking place each second at Brookhaven Lab's Relativistic Heavy Ion Collider (RHIC, a DOE Office of Science User Facility for nuclear physics research), and the ATLAS experiment at the Large Hadron Collider in Europe. Watch this video to learn how the people and computing resources of the RACF serve these scientists to turn petabytes of raw data into physics discoveries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Intersystem crossing and exciton–defect coupling of spin defects in hexagonal boron nitride

Despite the recognition of two-dimensional (2D) systems as emerging and scalable host materials of single-photon emitters or spin qubits, the uncontrolled, and undetermined chemical nature of these quantum defects has been a roadblock to further development. Leveraging the design of extrinsic defects can circumvent these persistent issues and provide an ultimate solution. Here, we established a complete theoretical framework to accurately and systematically design quantum defects in wide-bandgap 2D systems. With this approach, essential static and dynamical properties are equally considered for spin qubit discovery. In particular, many-body interactions such as defect–exciton couplings are vital for describing excited state properties of defects in ultrathin 2D systems. Meanwhile, nonradiative processes such as phonon-assisted decay and intersystem crossing rates require careful evaluation, which competes together with radiative processes. From a thorough screening of defects based on first-principles calculations, we identify promising single-photon emitters such as Si VV and spin qubits such as Ti VV and Mo VV in hexagonal boron nitride. This work provided a complete first-principles theoretical framework for defect design in 2D materials.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Spin-neurons: A possible path to energy-efficient neuromorphic computers

Recent years have witnessed growing interest in the field of brain-inspired computing based on neural-network architectures. In order to translate the related algorithmic models into powerful, yet energy-efficient cognitive-computing hardware, computing-devices beyond CMOS may need to be explored. The suitability of such devices to this field of computing would strongly depend upon how closely their physical characteristics match with the essential computing primitives employed in such models. In this work, we discuss the rationale of applying emerging spin-torque devices for bio-inspired computing. Recent spin-torque experiments have shown the path to low-current, low-voltage, and high-speed magnetization switching in nano-scale magnetic devices. Such magneto-metallic, current-mode spin-torque switches can mimic the analog summing and “thresholding” operation of an artificial neuron with high energy-efficiency. Comparison with CMOS-based analog circuit-model of a neuron shows that “spin-neurons” (spin based circuit model of neurons) can achieve more than two orders of magnitude lower energy and beyond three orders of magnitude reduction in energy-delay product. The application of spin-neurons can therefore be an attractive option for neuromorphic computers of future.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Successive phase transitions and kink solutions in Φ⁸, Φ¹⁰, and Φ¹² field theories

We obtain exact solutions for kinks in Φ⁸, Φ¹⁰, and Φ¹² field theories with degenerate minima, which can describe a second-order phase transition followed by a first-order one, a succession of two first-order phase transitions and a second-order phase transition followed by two first-order phase transitions, respectively. Such phase transitions are known to occur in ferroelastic and ferroelectric crystals and in meson physics. In particular, we find that the higher-order field theories have kink solutions with algebraically-decaying tails and also asymmetric cases with mixed exponential-algebraic tail decay, unlike the lower-order Φ⁴ and Φ⁶ theories. Additionally, we construct distinct kinks with equal energies in all three field theories considered, and we show the co-existence of up to three distinct kinks (for a Φ¹² potential with six degenerate minima). We also summarize phonon dispersion relations for these systems, showing that the higher-order field theories have specific cases in which only nonlinear phonons are allowed. For the Φ¹⁰ field theory, which is a quasi-exactly solvable (QES) model akin to Φ⁶, we are also able to obtain three analytical solutions for the classical free energy as well as the probability distribution function in the thermodynamic limit.

36 MATERIALS SCIENCE↗
DOE OSTI2012

IER-163 Post-Experiment MCNP Calculations (U)

IER-163 has been modeled with high fidelity in MCNP6. The model k eff was high, as in other similar calculations. The fission ratio 238 U(n,f)/ 235 U(n,f) was 12.6% too small compared with measurements; the ratio 239 Pu(n,f)/ 235 U(n,f) was 11.5% too small compared with measurements; the iridium ratio 193 Ir(n,n´)/ 191 Ir(n,γ) was 16.4% too large; and the gold ratios 197 Au(n,2n)/ 197 Au(n,γ), 197 Au(n,2n)/ 235 U(n,f), and 197 Au(n,γ)/ 235 U(n,f) were within one standard deviation of the measured values. It is suggested that the calculated 235 U fission rate is too large and the calculated 238 U fission rate is too small.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2012

US/Russian Cooperative Program in Science and Technology

Developing a sophisticated theory to understand the electronic structure of 5f-metals is a great challenge to solid-state physics. Complicated electronic structures of 5f metals make their properties strongly sensitive to small energy changes produced by the addition of a small amount of alloy, impurities, or crystal structure defects caused by irradiation. A theoretical material science technique applicable to investigate these effects is atomistic simulation using Classical Molecular Dynamics (CMD). In contrast to ab initio techniques, CMD may include several million particles, so that there is a possibility of direct simulation of very low concentration impurities and defects (as well as phenomena such as plasticity and polymorphous transitions) under given conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Integral Criticality Estimators in MCATK(U)

The Monte Carlo Application Toolkit (MCATK) is a component-based software toolset for delivering customized particle transport solutions using the Monte Carlo method. Currently under development in the XCP-MCC group at Los Alamos National Laboratory, the toolkit has the ability to estimate the k eff and α eigenvalues for static geometries. This paper presents a description of the estimators and variance reduction techniques available in the toolkit and includes a preview of those slated for future releases. Along with the description of the underlying algorithms is a description of the available user inputs for controlling the iterations. The paper concludes with a comparison of the MCATK results with those provided by analytic solutions. The results match within expected statistical uncertainties and demonstrate MCATK's usefulness in estimating these important quantities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗