Search NASA⌕ Search

SEARCH · Search NASA

Results for “approximate dynamic programming”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

The potential of positron emission tomography for intratreatment dynamic lung tumor tracking: A phantom study

Purpose: This study aims to evaluate the potential and feasibility of positron emission tomography for dynamic lung tumor tracking during radiation treatment. The authors propose a center of mass (CoM) tumor tracking algorithm using gated-PET images combined with a respiratory monitor and investigate the geometric accuracy of the proposed algorithm. Methods: The proposed PET dynamic lung tumor tracking algorithm estimated the target position information through the CoM of the segmented target volume on gated PET images reconstructed from accumulated coincidence events. The information was continuously updated throughout a scan based on the assumption that real-time processing was supported (actual processing time at each frame ≈10 s). External respiratory motion and list-mode PET data were acquired from a phantom programmed to move with measured respiratory traces (external respiratory motion and internal target motion) from human subjects, for which the ground truth target position was known as a function of time. The phantom was cylindrical with six hollow sphere targets (10, 13, 17, 22, 28, and 37 mm in diameter). The measured respiratory traces consisted of two sets: (1) 1D-measured motion from ten healthy volunteers and (2) 3D-measured motion from four lung cancer patients. The authors evaluated the geometric accuracy of the proposed algorithm by quantifying estimation errors (Euclidean distance) between the actual motion of targets (1D-motion and 3D-motion traces) and CoM trajectories estimated by the proposed algorithm as a function of time. Results: The time-averaged error of 1D-motion traces over all trajectories of all targets was 1.6 mm. The error trajectories decreased with time as coincidence events were accumulated. The overall error trajectory of 1D-motion traces converged to within 2 mm in approximately 90 s. As expected, more accurate results were obtained for larger targets. For example, for the 37 mm target, the average error over all 1D-motion traces was 1.1 mm; and for the 10 mm target, the average error over all 1D-motion traces was 2.8 mm. The overall time-averaged error of 3D-motion traces was 1.6 mm, which was comparable to that of the 1D-motion traces. There were small variations in the errors between the 3D-motion traces, although the motion trajectories were very different. The accuracy of the estimates was consistent for all targets except for the smallest. Conclusions: The authors developed an algorithm for dynamic lung tumor tracking using list-mode PET data and a respiratory motion signal, and demonstrated proof-of-principle for PET-guided lung tumor tracking. The overall tracking error in phantom studies is less than 2 mm.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Traversable wormhole dynamics on a quantum processor

The holographic principle, theorized to be a property of quantum gravity, postulates that the description of a volume of space can be encoded on a lower-dimensional boundary. The anti-de Sitter (AdS)/conformal field theory correspondence or duality is the principal example of holography. The Sachdev–Ye–Kitaev (SYK) model of N >> 1 Majorana fermions has features suggesting the existence of a gravitational dual in AdS 2 , and is a new realization of holography. Here, we invoke the holographic correspondence of the SYK many-body system and gravity to probe the conjectured ER=EPR relation between entanglement and spacetime geometry through the traversable wormhole mechanism as implemented in the SYK model. A qubit can be used to probe the SYK traversable wormhole dynamics through the corresponding teleportation protocol. This can be realized as a quantum circuit, equivalent to the gravitational picture in the semiclassical limit of an infinite number of qubits. Here we use learning techniques to construct a sparsified SYK model that we experimentally realize with 164 two-qubit gates on a nine-qubit circuit and observe the corresponding traversable wormhole dynamics. Despite its approximate nature, the sparsified SYK model preserves key properties of the traversable wormhole physics: perfect size winding, coupling on either side of the wormhole that is consistent with a negative energy shockwave, a Shapiro time delay, causal time-order of signals emerging from the wormhole, and scrambling and thermalization dynamics. Our experiment was run on the Google Sycamore processor. By interrogating a two-dimensional gravity dual system, our work represents a step towards a program for studying quantum gravity in the laboratory. Future developments will require improved hardware scalability and performance as well as theoretical developments including higher-dimensional quantum gravity duals and other SYK-like models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2010

Simulated combined abnormal environment fire calculations for aviation impacts.

Aircraft impacts at flight speeds are relevant environments for aircraft safety studies. This type of environment pertains to normal environments such as wildlife impacts and rough landings, but also the abnormal environment that has more recently been evidenced in cases such as the Pentagon and World Trade Center events of September 11, 2001, and the FBI building impact in Austin. For more severe impacts, the environment is combined because it involves not just the structural mechanics, but also the release of the fuel and the subsequent fire. Impacts normally last on the order of milliseconds to seconds, whereas the fire dynamics may last for minutes to hours, or longer. This presents a serious challenge for physical models that employ discrete time stepping to model the dynamics with accuracy. Another challenge is that the capabilities to model the fire and structural impact are seldom found in a common simulation tool. Sandia National Labs maintains two codes under a common architecture that have been used to model the dynamics of aircraft impact and fire scenarios. Only recently have these codes been coupled directly to provide a fire prediction that is better informed on the basis of a detailed structural calculation. To enable this technology, several facilitating models are necessary, as is a methodology for determining and executing the transfer of information from the structural code to the fire code. A methodology has been developed and implemented. Previous test programs at the Sandia National Labs sled track provide unique data for the dynamic response of an aluminum tank of liquid water impacting a barricade at flight speeds. These data are used to validate the modeling effort, and suggest reasonable accuracy for the dispersion of a non-combustible fluid in an impact environment. The capability is also demonstrated with a notional impact of a fuel-filled container at flight speed. Both of these scenarios are used to evaluate numeric approximations, and help provide an understanding of the quantitative accuracy of the modeling methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Impact of Heaters on Molten Salt Reactors Dynamics

Recently, there has been renewed interest in the Molten Salt Reactor (MSR) concept. This interest is mainly due to its advantages over Light Water Reactors (LWRs). Among these advantages are the flexibility in the fuel choice and the possibility to burn actinides [1]. Currently, two molten research reactors are under development. A 2 MW molten salt experimental reactor has been constructed and is planned for operation in China. A 1 MW Molten Salt Research Reactor is under construction in the United States [2]. Thus, Research and Development (R&D) programs are needed for the MSR technology. Reactor dynamics studies are crucial safety evaluation studies. The development of reactor dynamics tools for MSRs started during the Molten Salt Reactor Program (MSRP) at ORNL, focusing on the Molten Salt Reactor Experiment (MSRE) [3]. Although various dynamic simulation tools are available, R&D efforts to develop specialized tools for MSRs are still ongoing [4]. The primary motivation for this interest is the goal of commercializing MSRs. Developing a dynamics tool for MSRs is crucial for analyzing their safety and supporting their demonstration efforts. A crucial safety issue in MSRs is ensuring that the fuel salt in the primary loop remains molten. This concern becomes particularly important when the reactor operates at lower power levels. Most chloride and fluoride salts used as fuel have a melting point of approximately 450 °C [5]. Therefore, maintaining fuel salt temperatures above this melting point is essential to prevent it from freezing. Freezing of the fuel salt can lead to volume expansion, potentially causing damage to reactor components [6]. Moreover, it can create local blockages within fuel salt channels, reducing cooling efficiency and resulting in hot spots [7]. Various methods can prevent fuel salt from freezing at low power levels. These include reducing the fuel mass flow rate in the primary loop, decreasing heat exchange between the primary and secondary loops in the reactor, and using electric fuel salt heaters. Reducing the mass flow rate of fuel salt in the primary loop can increase the fuel salt transit time in the core, enhancing heat production and increasing temperature. However, it can potentially lead to challenges. These include increasing the potential for fuel salt deposition on the reactor channel walls [8], affecting the overall heat transfer processes. Thus, alternative methods should be used [9]. This research investigates the effects of fuel salt heaters on the dynamics and stability during MSR operation. The investigation includes studying the reactor stability with the fuel salt heaters on and off. In addition, the reactor response to transients is compared with and without the operation of fuel salt heaters. These transients include reactivity insertion, primary and secondary pump failure, and heat sink temperature increase. The results of the reactivity insertion transient are presented in the summary.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Analysis of Dynamic Stress in The Reactor Core Barrel Due to the Pressure Pulsations

This article presents vibrations analysis of the reactor core barrel caused by pressure pulsations induced by the main coolant pump. For this purpose, the calculations of the pressure distribution in the annulus between the core barrel and the reactor pressure vessel, bounded above by a separating ring were performed. Using transfer matrix method is obtained the solution of two-dimensional problem of pressure pulsations in the annulus between reactor core barrel and reactor vessel. The calculation results are compared with the pulsation pressure measurements made at commissioning unit 2 of the South Ukraine Nuclear Power Station. The distribution of pressure over the height of core barrel was obtained, which makes possible to estimate its strength for variant deformation of the core barrel as a beam, and in the case of deformation of the core barrel as a shell. The calculation results are used to assess the reliability of core barrel pre-load, which clamps the core barrel flange in place at the top, at full power operating. (author)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Model-based and Model-free Designs for an Extended Continuous-time LQR with Exogenous Inputs

We present an extended linear quadratic regulator (LQR) design for continuous-time linear time-invariant (LTI) systems in the presence of exogenous inputs. We first propose a model-based solution with cost minimization guarantees for states and inputs using dynamic programming (DP). The control law consists of a combination of the optimal state feedback and an additional optimal term dependent on the exogenous inputs. The control gains for the two components are obtained by solving a set of matrix differential equations. We provide these solutions for both finite horizons and steady-state cases. In the second part of the paper, we formulate a reinforcement learning (RL) based algorithm which does not need any model information except the input matrix, and can compute an approximate steady-state LQR gain using measurements of the states, the control inputs, and the exogenous inputs. Both model-based and data-driven optimal control algorithms are tested with a numerical example under different exogenous inputs showcasing the effectiveness of the designs.

Mukherjee, Sayak↗
DOE OSTI2016

A performance model for the communication in fast multipole methods on high-performance computing platforms

Exascale systems are predicted to have approximately 1 billion cores, assuming gigahertz cores. Limitations on affordable network topologies for distributed memory systems of such massive scale bring new challenges to the currently dominant parallel programing model. Currently, there are many efforts to evaluate the hardware and software bottlenecks of exascale designs. It is therefore of interest to model application performance and to understand what changes need to be made to ensure extrapolated scalability. The fast multipole method (FMM) was originally developed for accelerating N-body problems in astrophysics and molecular dynamics but has recently been extended to a wider range of problems. Its high arithmetic intensity combined with its linear complexity and asynchronous communication patterns make it a promising algorithm for exascale systems. In this paper, we discuss the challenges for FMM on current parallel computers and future exascale architectures, with a focus on internode communication. We focus on the communication part only; the efficiency of the computational kernels are beyond the scope of the present study. We develop a performance model that considers the communication patterns of the FMM and observe a good match between our model and the actual communication time on four high-performance computing (HPC) systems, when latency, bandwidth, network topology, and multicore penalties are all taken into account. To our knowledge, this is the first formal characterization of internode communication in FMM that validates the model against actual measurements of communication time. The ultimate communication model is predictive in an absolute sense; however, on complex systems, this objective is often out of reach or of a difficulty out of proportion to its benefit when there exists a simpler model that is inexpensive and sufficient to guide coding decisions leading to improved scaling. The current model provides such guidance.

Computer Science↗
DOE OSTI2024

LoVoCCS. II. Weak Lensing Mass Distributions, Red-sequence Galaxy Distributions, and Their Alignment with the Brightest Cluster Galaxy in 58 Nearby X-Ray-luminous Galaxy Clusters

The Local Volume Complete Cluster Survey is an ongoing program to observe nearly a hundred low-redshift X-ray-luminous galaxy clusters (redshifts 0.03 < z < 0.12 and X-ray luminosities in the 0.1–2.4 keV band L X500c > 10 44 erg s −1 ) with the Dark Energy Camera, capturing data in the u, g, r, i, z bands with a 5σ point source depth of approximately 25th–26th AB magnitudes. Here, we map the aperture masses in 58 galaxy cluster fields using weak gravitational lensing. These clusters span a variety of dynamical states, from nearly relaxed to merging systems, and approximately half of them have not been subject to detailed weak lensing analysis before. In each cluster field, we analyze the alignment between the 2D mass distribution described by the aperture mass map, the 2D red-sequence (RS) galaxy distribution, and the brightest cluster galaxy (BCG). We find that the orientations of the BCG and the RS distribution are strongly aligned throughout the interiors of the clusters: the median misalignment angle is 19° within 2 Mpc. We also observe the alignment between the orientations of the RS distribution and the overall cluster mass distribution (by a median difference of 32° within 1 Mpc), although this is constrained by galaxy shape noise and the limitations of our cluster sample size. These types of alignment suggest long-term dynamical evolution within the clusters over cosmic timescales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

FY21 Progress Report: SRNL Analysis of ICCWR LCM and WAMS data for Corrosion and Cracking

The development of algorithms for machine learning and data analysis for the 3013 Surveillance Program is a collaborative effort by the Savannah River National Laboratory (SRNL) and the University of South Carolina (USC). For corrosion detection, Laser Confocal Microscope (LCM) or Wide Area 3D Measurement System (WAMS) data is extracted from large binary files, with software written to convert the data to physical attributes (e.g., height, color and grayscale values; all as functions of a location in a plane projection). A user-friendly Matlab Graphical User Interface (GUI) that reads data from either LCM or WAMS files was developed to integrate input data with software developed for processing and evaluation. The GUI can selectively download binary data, interrogate data attributes, label data, flag significant features, execute Machine Learning (ML) algorithms, output parameters for trained ML algorithms, report ML model accuracy with respect to labeled data, and generate graphical representations for various analyses. Features can be called out by user-specified thresholds, manual labeling or machine learning algorithms when they have been completed. The ability to rapidly label data is important because of the volume of data required for training machine learning algorithms. The GUI has the flexibility to allow addition of improved ML algorithms, methods for data visualization, and statistical computations. Statistical analyses via the GUI include areas of pits within a defined range of pit depths, correlations between Red-Green-Blue (RGB) or grayscale intensity and relative surface height, covariances between values associated with features, and feature histograms. The development of supervised machine learning algorithms, however, has been hindered by a lack of training data. The machine learning algorithms for crack identification are being refined but require improvements to the true positive rate for crack detection. This shortcoming is an artifact of the limited training data currently available, perhaps more so than the structure of the neural networks. At present, the best results are had from a consensus over an ensemble of randomly generated Deep Neural Network (DNN) or Convolutional Neural Network (CNN) algorithms. Although the consensus accuracy method has yielded optimum true positive and true negative rates in excess of 80%, additional validation testing is necessary. In addition to the suite of LCM data that was initially used, and which represents the majority of the work presented in this report, WAMS image data was also reviewed at a preliminary level. The review included a comparison between image resolution and dynamic range for each method. WAMS (ZON file) image data was found to have a pixel pitch of 3.69μm compared to 1 μm for the LCM (vk4 file) data, which implies a lower resolution for the WAMS images. Conversely, the ratio of dynamic range of the WAMS data to the LCM data was approximately 41:20 for height data, suggesting that information from WAMS should more accurately determine the depth of pits. At present, the significance of the greater dynamic range of the WAMS data relative to the LCM data has not yet been evaluated.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

ELECTRONIC STRUCTURE METHODS AND PROTOCOLS WITH APPLICATION TO DYNAMICS, KINETICS AND THERMOCHEMISTRY

Hydrocarbon combustion involves the reaction dynamics of a tremendous number of species beginning with many-component fuel mixtures and proceeding via a complex system of intermediates to form primary and secondary products. Combustion conditions corresponding to new advanced engines and/or alternative fuels rely increasingly on autoignition and low-temperature-combustion chemistry. In these regimes various transient radical species such as HO2, ROO·, ·QOOH, HCO, NO2, HOCO, and Criegee intermediates play important roles in determining the detailed as well as more general dynamics. A clear understanding and accurate representation of these processes is needed for effective modeling. Given the difficulties associated with making reliable experimental measurements of these systems, computation can play an important role in developing these energy technologies. Accurate calculations have their own challenges since even within the simplest dynamical approximations such as transition state theory, the rates depend exponentially on critical barrier heights and these may be sensitive to the level of quantum chemistry. Moreover, it is well-known that in many cases it is necessary to go beyond statistical theories and consider the dynamics. Quantum tunneling, resonances, radiative transitions, and non-adiabatic effects governed by spin-orbit or derivative coupling can be determining factors in those dynamics. Building upon progress made during a period of prior support through the DOE Early Career Program, this project combines developments in the areas of potential energy surface (PES) fitting and multistate multireference quantum chemistry to allow spectroscopically and dynamically/kinetically accurate investigations of key molecular systems (such as those mentioned above), many of which are radicals with strong multireference character and have the possibility of multiple electronic states contributing to the observed dynamics. An ongoing area of investigation is to develop general strategies for robustly convergent electronic structure theory for global multichannel reactive surfaces including diabatization of energy and other relevant surfaces such as dipole transition. Combining advances in ab initio methods with automated interpolative PES fitting allows the construction of high-quality PESs (incorporating thousands of high-level data) to be done rapidly through parallel processing on high-performance computing (HPC) clusters. In addition, new methods and approaches to electronic structure theory will be developed and tested through applications. This project will explore limitations in traditional multireference calculations (e.g., MRCI) such as those imposed by internal contraction, lack of high-order correlation treatment and poor scaling. Methods such as DMRG-based extended active-space CASSCF and various Quantum Monte Carlo (QMC) methods will be applied (including VMC/DMC and FCIQMC). Insight into the relative significance of different orbital spaces and the robustness of application of these approaches on leadership class computing architectures will be gained. Synergy with other components of this research program such as automated PES fitting and multireference quantum chemistry will be used to address challenges encountered by the standard approaches to computational thermochemistry (those being single-reference quantum chemistry and perturbative treatments of the anharmonic vibrational energy, which break down for some cases of electronic structure or floppy strongly coupled vibrational modes).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Simulating non-native cubic interactions on noisy quantum machines

As a milestone for general-purpose computing machines, we demonstrate that quantum processors can be programed to efficiently simulate dynamics that are not native to the hardware. Moreover, on noisy devices without error correction, we show that simulation results are significantly improved when the quantum program is compiled using modular gates instead of a restricted set of standard gates. We demonstrate the general methodology by solving a cubic interaction problem, which appears in nonlinear optics, gauge theories, as well as plasma and fluid dynamics. To encode the non-native Hamiltonian evolution, we decompose the Hilbert space into a direct sum of invariant subspaces in which the nonlinear problem is mapped to a finite-dimensional Hamiltonian simulation problem. Furthermore, in a three-states example, the resultant unitary evolution is realized by a product of approximately 20 standard gates, using which approximately ten simulation steps can be carried out on state-of-the-art quantum hardware before results are corrupted by decoherence. In comparison, the simulation depth is improved by more than an order of magnitude when the unitary evolution is realized as a single cubic gate, which is compiled directly using optimal control. Alternatively, parametric gates may also be compiled by interpolating control pulses. Modular gates thus obtained provide high-fidelity building blocks for quantum Hamiltonian simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Notes from the Unit Commitment Linear Programming Project

The Unit Commitment Problem (UCP) is a computational bottleneck in electrical grid optimization. The problem bears a close resemblance to a similar problem in High Performance Computing: the Dynamic Volt age/Frequency Scheduling (DVFS) problem. In my earlier work I successfully tranformed the DVFS into a closely related problem amenable to a linear programming solution and was thus able to schedule problem instances several orders of magnitude larger than any that had been at tempted before. This technical report details my unsuccessful efforts to derive a similar transformation for the UCP. The problem proved impossible not because generators were either on or off, but because the state of an earlier timeslice constrains the states of latter timeslices in ways that cannot be sufficiently well approximated for useful solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Cryogenic Design and Thermal Analysis of EIC Central Detector (ePIC) Solenoid Magnet (MARCO)

The Electron Ion Collider (EIC) physics program utilizes a 2.0 T superconducting magnet at the heart of its ePIC detector system. This approximately 3.5 m long and 2.84 m diameter warm bore magnet has a 20 tons cold mass which is conduction-cooled using liquid helium at 4.5 K. A closed loop active thermosiphon system is chosen to facilitate the cooling and to maintain a minimum of 2 K temperature margin on the peak operating temperature (4.7 K) of the superconductor. Here, this paper presents the cryogenic design of the cooling system and the thermal analysis of the solenoid. A Computational fluid dynamics (CFD) model was developed to quantify the performance of the two–phase closed thermosiphon system and predict the temperature gradient on the cold mass.

Gopinath, Sandesh [Thomas Jefferson National Accel↗
DOE OSTI2019

Identification of Ebola Virus Inhibitors Targeting GP2 Using Principles of Molecular Mimicry

A key step in the Ebola virus (EBOV) replication cycle involves conformational changes in viral glycoprotein 2 (GP2) which facilitate host-viral membrane fusion and subsequent release of the viral genome. Ebola GP2 plays a critical role in virus entry and has similarities in mechanism and structure to the HIV gp41 protein for which inhibitors have been successfully developed. In this work, a putative binding pocket for the C-terminal heptad repeat in the N-terminal heptad repeat trimer was targeted for identification of small molecules that arrest EBOV-host membrane fusion. Two computational structure-based virtual screens of ~1.7 M compounds were performed (DOCK program) against a GP2 five-helix bundle, resulting in 165 commercially available compounds purchased for experimental testing. Based on assessment of inhibitory activity, cytotoxicity, and target specificity, four promising candidates emerged with 50% inhibitory concentration values in the 3 to 26 μM range. Molecular dynamics simulations of the two most potent candidates in their DOCK-predicted binding poses indicate that the majority of favorable interactions involve seven highly conserved residues that can be used to guide further inhibitor development and refinement targeting EBOV. The most recent Ebola virus disease outbreak, from 2014 to 2016, resulted in approximately 28,000 individuals becoming infected, which led to over 12,000 causalities worldwide. The particularly high pathogenicity of the virus makes paramount the identification and development of promising lead compounds to serve as inhibitors of Ebola infection. To limit viral load, the virus-host membrane fusion event can be targeted through the inhibition of the class I fusion glycoprotein of. In the current work, several promising small-molecule inhibitors that target the glycoprotein GP2 were identified through systematic application of structure-based computational and experimental drug design procedures.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2011

Conservative self-force correction to the innermost stable circular orbit: Comparison with multiple post-Newtonian-based methods

Barack and Sago [Phys. Rev. Lett. 102, 191101 (2009)] have recently computed the shift of the innermost stable circular orbit (ISCO) of the Schwarzschild spacetime due to the conservative self-force that arises from the finite-mass of an orbiting test-particle. This calculation of the ISCO shift is one of the first concrete results of the self-force program, and provides an exact (fully relativistic) point of comparison with approximate post-Newtonian (PN) computations of the ISCO. Here this exact ISCO shift is compared with nearly all known PN-based methods. These include both 'nonresummed' and 'resummed' approaches (the latter reproduce the test-particle limit by construction). The best agreement with the exact (Barack-Sago) result is found when the pseudo-4PN coefficient of the effective-one-body (EOB) metric is fit to numerical relativity simulations. However, if one considers uncalibrated methods based only on the currently known 3PN-order conservative dynamics, the best agreement is found from the gauge-invariant ISCO condition of Blanchet and Iyer [Classical Quantum Gravity 20, 755 (2003)], which relies only on the (nonresummed) 3PN equations of motion. This method reproduces the exact test-particle limit without any resummation. A comparison of PN methods with the ISCO in the equal-mass case (computed via sequences of numerical relativity initial-data sets) is also performed. Here a (different) nonresummed method also performs very well (as was previously shown). These results suggest that the EOB approach - while exactly incorporating the conservative test-particle dynamics and having several other important advantages - does not (in the absence of calibration) incorporate conservative self-force effects more accurately than standard PN methods. I also consider how the conservative self-force ISCO shift, combined in some cases with numerical relativity computations of the ISCO, can be used to constrain our knowledge of (1) the EOB effective metric, (2) phenomenological inspiral-merger-ringdown templates, and (3) 4PN- and 5PN-order terms in the PN orbital energy. These constraints could help in constructing better gravitational-wave templates. Lastly, I suggest a new method to calibrate unknown PN terms in inspiral templates using numerical-relativity calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Automated Dynamic Demand Response Implementation on a Micro-grid

In this paper, we describe a system for real-time automated Dynamic and Sustainable Demand Response with sparse data consumption prediction implemented on the University of Southern California campus microgrid. Supply side approaches to resolving energy supply-load imbalance do not work at high levels of renewable energy penetration. Dynamic Demand Response (D 2 R) is a widely used demand-side technique to dynamically adjust electricity consumption during peak load periods. Our D 2 R system consists of accurate machine learning based energy consumption forecasting models that work with sparse data coupled with fast and sustainable load curtailment optimization algorithms that provide the ability to dynamically adapt to changing supply-load imbalances in near real-time. Our Sustainable DR (SDR) algorithms attempt to distribute customer curtailment evenly across sub-intervals during a DR event and avoid expensive demand peaks during a few sub-intervals. It also ensures that each customer is penalized fairly in order to achieve the targeted curtailment. We develop near linear-time constant-factor approximation algorithms along with Polynomial Time Approximation Schemes (PTAS) for SDR curtailment that minimizes the curtailment error defined as the difference between the target and achieved curtailment values. Our SDR curtailment problem is formulated as an Integer Linear Program that optimally matches customers to curtailment strategies during a DR event while also explicitly accounting for customer strategy switching overhead as a constraint. We demonstrate the results of our D 2 R system using real data from experiments performed on the USC smartgrid and show that 1) our prediction algorithms can very accurately predict energy consumption even with noisy or missing data and 2) our curtailment algorithms deliver DR with extremely low curtailment errors in the 0.01-0.05 kWh range.

Kuppannagari, Sanmukh R.↗
DOE OSTI2015

Quantum dynamics study on the CHIPR potential energy surface for the hydroperoxyl radical: The reactions O + OH⇋O{sub 2} + H

Quantum scattering calculations of the O({sup 3}P)+OH({sup 2}Π)⇌O{sub 2}({sup 3}Σ{sub g}{sup −})+H({sup 2}S) reactions are presented using the combined-hyperbolic-inverse-power-representation potential energy surface [A. J. C. Varandas, J. Chem. Phys. 138, 134117 (2013)], which employs a realistic, ab initio-based, description of both the valence and long-range interactions. The calculations have been performed with the ABC time-independent quantum reactive scattering computer program based on hyperspherical coordinates. The reactivity of both arrangements has been investigated, with particular attention paid to the effects of vibrational excitation. By using the J-shifting approximation, rate constants are also reported for both the title reactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Computer simulation on spatial resolution of X-ray bright-field imaging by dynamical diffraction theory for a Laue-case crystal analyzer

Recently, dark-field imaging (DFI) and bright-field imaging (BFI) have been proposed and applied to visualize X-ray refraction effects yielded in biomedical objects. In order to clarify the spatial resolution due to a crystal analyzer in Laue geometry, a program based on the Takagi-Taupin equation was modified to be used for carrying out simulations to evaluate the spatial resolution of images coming into a Laue angular analyzer (LAA). The calculation was done with a perfect plane wave for diffraction wave-fields, which corresponded to BFI, under the conditions of 35 keV and a diffraction index 440 for a 2100 {mu}m thick LAA. As a result, the spatial resolution along the g-vector direction showed approximately 37.5 {mu}m. 126 {mu}m-thick LAA showed a spatial resolution better than 3.1 {mu}m under the conditions of 13.7 keV and a diffraction index 220.

36 MATERIALS SCIENCE↗