Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

ISO_Fortran_binding_m v0.1.0

The Fortran programming language standard defines a broad feature set supporting the interoperability of Fortran programs with program written according to the C programming language standard. Among Fortran's C-interoperability features is a a C header file "ISO_Fortran_binding.h" This header file defines the interface to various C data structures and functions that C programs may use to access Fortran data entities. The ISO_Fortran_binding_m software defines a native Fortran module that presents an interface to these same data structures and functions. ISO_Fortran_bind_m thus enables Fortran programs to access and manipulate Fortran entities in ways that precisely mirror what C programs can do using ISO_Fortran_binding.h. ISO_Fortran_binding_m facilitates writing portable standard-conforming Fortran programs that emulate non-interoperable features, e.g., dynamic polymorphism, in a standard-conforming interoperable way similar but broader than what is demonstrated in Berkeley Lab's Caffeine software [1]. ISO_Fortran_binding_m also enables a Fortran programmer to extend Fortran's capabilities to emulate certain C functionality such as memory address arithmetic or computing C's "sizeof" function. [1] https://github.com/BerkeleyLab/caffeine/blob/213e3df1c319f0663306354414f352acda42a24f/src/caffeine/collective_subroutines/co_reduce_s.f90#L88 [2] https://github.com/BerkeleyLab/ISO_Fortran_binding_m/blob/0c585362bb4f2c72cf9049c800a7115b529ec533/src/iso_fortran_binding_m.F90#L196

Rouson, Damian↗
DOE OSTI2025

PyDislocDyn: A Python code for calculating dislocation drag and other crystal properties

PyDislocDyn is a suite of python programs designed to perform various calculations for dislocation dynamics in the continuum limit. In particular, one of its main purposes is to calculate dislocation drag from phonon wind. Additional features include the averaging of elastic constants for polycrystals, the calculation of the dislocation field including its limiting velocities, and the calculation of dislocation self-energy and line tension.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Computation of the Seismic Stability of Earth Retaining Structures Final Report CRADA No. TSB-954-94

The purpose of this CRADA project was to evaluate the seismic stability of block retaining wall systems. Retaining wall systems are used extensively in private and commercial developments. This study was designed to develop and demonstrate a computer modeling technology to be used to predict the seismic stability of any block wall system design. The nonlinear finite element computer programs developed at LLNL and employed in the Computational Earthquake Initiative were utilized in this small business CRADA to analyze the seismic stability of the block retaining walls. The unique capability of the LLNL programs to rigorously model frictional contact in a dynamic analysis problem were used in a computer simulation of the dynamic interaction the block waiVsoil systems under seismic excitation.

58 GEOSCIENCES↗
DOE OSTI2020

Evaluation of a streaming readout solution for Jefferson Lab experiments

Jefferson Lab (JLab) operates the CEBAF accelerator that provides a high energy (up to 12 GeV) and high current (up to 150 µA) polarized electron beam for fixed target experiments. JLab's mission is to provide substantial progress in understanding the Quantum Chromo Dynamics (QCD), the current theory of the spectrum and the dynamics of hadronic matter. Given the broad scientific program a general-purpose detector, CLAS12, has been designed and constructed. Such a detector requires a sophisticated trigger and current experiments use an on-line FPGA-based system that relies upon custom firmware and electronics both of which are difficult reconfigure from one experiment to the next. To overcome these challenges transitioning to a streaming readout has been proposed since it would allow a more flexible, easier to debug, software trigger to be developed. To facilitate this JLab is gaining experience with streaming readout systems and is actively working to design and test a streaming data acquisition architecture. Significant work has already been done to make existing front-end electronics compatible with operation in streaming mode. Software-wise, much work remains to be done but an existing framework, TriDAS, which was developed for underwater neutrino telescopes can, in principle, be used with the existing JLab hardware. A test stand is already in place in the INDRA DAQ development lab at JLab. This paper will present the development of this TriDAS based test system and it?s operation in-beam using parts of the CLAS12 detector to evaluate streaming readout for a large general purpose detector.

Battaglieri, Marco↗
DOE OSTI2018

Decomposition of Nonlinear Dynamical Networks via Comparison Systems

In analysis and control of large-scale nonlinear dynamical systems, a distributed approach is often an attractive option due to its computational tractability and usually low communication requirements. Success of the distributed control design relies on the separability of the network into weakly interacting subsystems such that minimal information exchange between subsystems is sufficient to achieve satisfactory control performance. While distributed analysis and control design for dynamical network have been well studied, decomposition of nonlinear networks into weakly interacting subsystems has not received as much attention. In this article we propose a vector Lyapunov functions based approach to quantify the energy-flow in a dynamical network via a model of a comparison system. Introducing a notion of power and energy flow in a dynamical network, we use sum-of-squares programming tools to partition polynomial networks into weakly interacting subsystems. Examples are provided to illustrate the proposed method of decomposition.

Maruf, Abdullah Al↗
DOE OSTI2021

Short Range Correlations and Few Body Precision Measurements

This talk will discuss the exciting experimental program studying short range correlations and few body dynamics at Jefferson Lab in ?=3 nuclei. Quasi-elastic electron scattering on these light nuclei have readily calculable ground states but can already exhibit complex nuclear behavior including nucleon-nucleon (NN) short range correlations (SRCs). 3H and 3He, being mirror nuclei, exploit the maximum available isospin asymmetry and can be compared directly with ab initio calculations. This talk will summarize the recent Hall A 3H and 3He results. The first measurement of SRC nucleons in 3He using the CLAS detector from Hall B and identifying lead neutrons in ?? pairs at high missing momentum will be shown. These fascinating results will be extended in an approved experiment using the CLAS12 detector to measure high-statistics and absolute cross sections on 2H, 3He, and 3H with a wide range of kinematical coverage. This talk will discuss the recent and proposed measurements and their implications on nuclear models for few body nuclear dynamics and short range structure.

Szumila-Vance, Holly↗
DOE OSTI2024

Evaluation of Thermal Scattering Law for U 3 Si 2 Using Temperature-Dependent Phonon Density of States

The thermal scattering law for uranium silicide, a candidate material for high-density nuclear fuel, was evaluated from the temperature-dependent phonon density of states obtained using the quasi-harmonic approximation. For uranium silicide model systems of various volumes near the equilibrium volume, Hellmann-Feynman forces were calculated by performing ab initio lattice dynamics simulations based on density functional theory. Subsequently, the quasi-harmonic approximation was applied by interpolating the results of the harmonic approximation for each volume system. As a sensitivity test, the phonon density of states was calculated for the structures with temperature-dependent volumes at 298 K and 1000 K, obtained by the quasi-harmonic approximation, and the corresponding thermal neutron scattering cross section was evaluated. The temperature-dependent volume and heat capacity of uranium silicide obtained by the quasi-harmonic approximation showed reasonable agreement with experimental data while lattice thermal expansion showed deviation from experimental observations. Thermal scattering cross sections evaluated from 0 K and temperature-dependent phonon density of states exhibited differences of up to 4.1% at 1000 K. This study examines the impact of temperature on the thermal scattering law, employing ab initio lattice dynamics simulations through the quasi-harmonic approximation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Interactive Visualization of Complex Seismic Data and Models Using Bokeh

Visualizing multidimensional data and models becomes more challenging as the volume and resolution of seismic data and models increase. But thanks to the development of powerful and accessible computer systems, a model web browser can be used to visualize complex scientific data and models dynamically. In this paper, we present four examples of seismic model visualization using an open-source Python package Bokeh. One example is a visualization of a surface-wave dispersion data set, another presents a view of three-component seismograms, and two illustrate methods to explore a 3D seismic-velocity model. Unlike other 3D visualization packages, our visualization approach has a minimum requirement on users and is relatively easy to develop, provided you have reasonable programming skills. Finally, utilizing familiar web browsing interfaces, the dynamic tools provide us an effective and efficient approach to explore large data sets and models.

58 GEOSCIENCES↗
DOE OSTI2022

Gaussian processes meet NeuralODEs: a Bayesian framework for learning the dynamics of partially observed systems from scarce and noisy data

We present a machine learning framework (GP-NODE) for Bayesian model discovery from partial, noisy and irregular observations of nonlinear dynamical systems. The proposed method takes advantage of differentiable programming to propagate gradient information through ordinary differential equation solvers and perform Bayesian inference with respect to unknown model parameters using Hamiltonian Monte Carlo sampling and Gaussian Process priors over the observed system states. This allows us to exploit temporal correlations in the observed data, and efficiently infer posterior distributions over plausible models with quantified uncertainty. The use of the Finnish Horseshoe as a sparsity-promoting prior for free model parameters also enables the discovery of parsimonious representations for the latent dynamics. A series of numerical studies is presented to demonstrate the effectiveness of the proposed GP-NODE method including predator–prey systems, systems biology and a 50-dimensional human motion dynamical system. This article is part of the theme issue ‘Data-driven prediction in dynamical systems’.

Science & Technology - Other Topics↗
DOE OSTI2020

Dynamic Modeling of Full Converter Adjustable-speed Pumped Storage Hydropower (FC AS-PSH)

Full converter adjustable-speed pumped storage hydropower (FC AS-PSH) technology, as one of advanced-PSH technology, is developed from wind turbine technology. By making the synchronous machine connect to the grid through a full-size converter, FC AS-PSH has a wider adjustment range of speed and a better reactive power control capability compared with a doubly-fed asynchronous generator AS-PSH technology. When it plays as an energy backup in the power system, FC AS-PSH can provide a much faster response than conventional-PSH (C-PSH) which makes this technology provide better ancillary service for a high renewable penetrated system. In this paper, the dynamic modeling of FC AS-PSH is fully studied. We develop a detailed model of this technology in the IEEE 14-bus system based on GE Positive Sequence Load Flow (PSLF) platform. Especially, the first governor model is developed based on the Engineer’s Program Control Language (EPCL) user-defined model in this platform. All operation modes are validated and studied under a system contingency. Besides, comparison cases between FC AS-PSH and C-PSH are studied to show advantages providing from FC AS-PSH when it works with renewable energy.

50 EE - Wind and Water Power Program - Water (EE-4↗
DOE OSTI2025

Dynamic Capacitive Wireless Power Transfer System for Electric Vehicles (CRADA Final Report)

Road transportation accounts for 23% of our nation's total energy consumption, 59% of our petroleum consumption and 22% of our total emissions. Electric vehicles (EVs) have much higher well-to-wheel efficiency compared to gasoline vehicles and can reduce our dependence on oil, minimize emissions, improve local air quality and provide a platform for technological innovation and economic growth. While major progress has been made in the development of electric vehicles, their penetration remains under 1 The main hurdles in their widespread adoption are high cost, limited range and long charging times, due to limitations in battery technology. An approach to overcome these hurdles is to substantially reduce the on-board energy storage and instead deliver power wirelessly to the vehicle while it is in motion dynamic wireless power transfer (WPT). For EV charging applications, researchers have traditionally focused on inductive WPT. An issue with inductive WPT systems is that for magnetic flux guidance and shielding, they require ferrite cores, making them expensive, bulky and difficult to embed in the roadway. Under an ARPA-E IDEAS grant, the CU Boulder team has recently developed an alternative approach to WPT based on capacitive coupling that is much smaller, lighter, less expensive, and easier to embed in the roadway than inductive WPT; and have demonstrated a 1.2-kW 6.78-MHz version of this capacitive WPT system on a stationary platform. Under this seed grant program, preliminary work will be carried out on a dynamic capacitive WPT system to prepare a strong proposal for the ARPA-E OPEN 2018 FOA.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Optimal Equilibrium Selection of Price-Maker Agents in Performance-Based Regulation Market

This paper analyzes the oligopolistic equilibrium of multiple price-maker agents in performance-based regulation (PBR) markets. In these markets, there are price-maker agents representing some frequency regulation (FR) providers and a number of independent price-taker FR providers. An equilibrium problem with equilibrium constraints (EPEC) model is employed in this paper to study the equilibria of a PBR market in the presence of price-maker agents and pricetaker FR providers. Due to the incorporation of the FR providers' dynamics, the proposed model is reformulated as a mixed-integer linear programming (MILP) problem over innovative mathematical techniques. An optimal equilibrium point is also selected for the market, where none of the agents is unique deviator and the dynamic performance of power system is improved simultaneously. The effectiveness of proposed optimal equilibrium point is evaluated in the numerical results section by comparing the outputs with conventional optimal dispatches of the FR providers.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

XPlacer/Tracer

XPlacer/Tracer is a dynamic analysis tool that finds bad memory access patterns in heterogeneous codes. XPlacer/Tracer consists of two components. (1) A configurable ROSE plugin for source code instrumentation of heterogeneous code written in the C/C++ and CUDA programming languages. (2) A sample configuration that instruments memory accesses to dynamic memory in CPU and GPU codes, and a runtime library that tracks these memory accesses at runtime. The collected information is reported in form of a textual summary or as memory map that can be converted to images.

Pirkelbauer, PeterM.↗
DOE OSTI2023

Up up down down left right left right B A Start for the catalytic hackers of programmable materials

Catalysts have advanced over the last century to accelerate and control reactions based on static active sites. More effective catalysis can be achieved using catalysts that change with time over the course of a reaction, providing a dynamic free energy landscape that is tailored to each step in the reaction sequence. Here, the catalyst is modulated via an input program that directs the surface to change physically or electronically with time, providing information regarding the extent and duration of change optimized for specific combinations of chemistry and programmable catalyst surfaces. While in its infancy, programmable catalysis is advancing with parallel efforts to establish fundamental principles of dynamic catalysts, design of programmable materials, and strategies to design input programs that will control catalysis for faster and more selective reactions.

catalyst↗
DOE OSTI2020

Pilgrim : A thermal rate constant calculator and a chemical kinetics simulator

Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

METHODOLOGY AND APPLICATION OF PHYSICAL SECURITY EFFECTIVENESS BASED ON DYNAMIC FORCE-ON-FORCE MODELING

This paper describes ongoing work within the Light Water Reactor Sustainability (LWRS) Program at Idaho National Laboratory (INL) to optimize security and cost of nuclear power plants (NPPs). It reviews the conservatisms in conventional physical security posture and regulations. It introduces the dynamic risk assessment tool developed at INL, Event Modeling Risk Assessment using Linked Diagrams (EMRALD). The dynamic assessment methodology leverages EMRALD to process results of force-on-force (FOF) simulations and crediting safety mitigation actions from probabilistic risk assessment (PRA) models as well as diverse and flexible coping strategies (FLEX) mitigation strategies. Timing information from these simulations are compared against the available time to perform mitigations obtained from Reactor Excursion and Leak Analysis Program (RELAP5) simulations. To illustrate the methodology, a station blackout (SBO) attack scenario was modeled in commercially available FOF simulation tools. The simulation results provide valuable insights into possible attack outcomes and as the probabilistic risk of a core damage event given these outcomes. Safety mitigation procedures were modeled in EMRALD, and were dependent on the attack outcomes by considering human operator uncertainties. RELAP5 simulations incorporating human and hardware uncertainties were performed to estimate the distribution of time-to-core damage. The results demonstrate that, even in the extreme case of a successful adversarial attack, plant mitigation strategies provide significantly high-likelihood of preventing radiological release. The proposed modeling and simulation framework of integrating FLEX equipment with FOF models enables the NPPs to credit FLEX portable equipment in the plant security posture, resulting in an efficient and optimized physical security.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Defect-Driven Redox Interplay on Anatase TiO 2 : Surface-Structure Dependent Activation for CO 2 Hydrogenation Catalysis

Titanium dioxide (TiO 2 ) is one of the most extensively studied oxides as an active catalyst or catalyst support, particularly in energy and environmental applications, but the atomistic mechanisms governing its dynamic response to reactive environments and their correlation to reactivity remain largely elusive. Using in situ environmental transmission electron microscopy (ETEM), synchrotron X-ray diffraction (XRD), ambient-pressure X-ray photoelectron spectroscopy (AP-XPS), temperature-programmed reduction (TPR), reactivity measurements, and theoretical modeling, we reveal the dynamic interplay between oxygen loss and replenishment of anatase TiO 2 under varying reactive conditions. Under H 2 exposure, anatase TiO 2 undergoes surface reduction via lattice oxygen loss, forming Ti 3 O 5 . In contrast, CO 2 exposure induces oxygen replenishment, reversing stoichiometry. In mixed H2/CO 2 environments, the reverse water–gas shift (RWGS) reaction proceeds selectively on stepped and high-indexed TiO 2 surfaces, whereas the thermodynamically stable TiO 2 (101) surface remains inactive and intact. Critically, H 2 pretreatment generates oxygen vacancies on TiO 2 (101), transforming it into an active Ti 3 O 5 or defect-rich surface that catalyzes RWGS. By correlating surface structure, defect dynamics, and gas-phase interactions, this work deciphers the competition between H 2 -driven reduction and CO 2 -driven oxidation pathways at the atomic scale. Furthermore, these insights establish defect engineering as a strategic lever to activate inert TiO 2 facets, advancing the design of adaptive catalysts for sustainable fuel synthesis technologies.

36 MATERIALS SCIENCE↗