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 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Finite Element Solver for Fission Dynamics

FELIX is a physics computer code used to model fission fragment mass distributions in a fully quantum-mechanical, microscopic framework that only relies on our current knowledge of nuclear forces [1]. It is an implementation of the time-dependent generator coordinate method (TDGCM), which simulates the dynamics of a collective quantum wave-packet assuming the motion is adiabatic. In typical applications of the TDFCM, the nuclear collective wavepacket is obtained as a superposition of wavefunctions obtained by solving the Hartree-Fock-Bogoliubov equations of nuclear density functional theory (DFT). The program calculates at each time step the coefficients of that superposition.

Regnier, D.↗
DOE OSTI2014

TAS: 89 0227: TAS Recovery Act - Optimization and Control of Electric Power Systems: ARRA

The name SuperOPF is used to refer several projects, problem formulations and soft-ware tools intended to extend, improve and re-define some of the standard methods of optimizing electric power systems. Our work included applying primal-dual interior point methods to standard AC optimal power flow problems of large size, as well as extensions of this problem to include co-optimization of multiple scenarios. The original SuperOPF problem formulation was based on co-optimizing a base scenario along with multiple post-contingency scenarios, where all AC power flow models and constraints are enforced for each, to find optimal energy contracts, endogenously determined locational reserves and appropriate nodal energy prices for a single period optimal power flow problem with uncertainty. This led to example non-linear programming problems on the order of 1 million constraints and half a million variables. The second generation SuperOPF formulation extends this by adding multiple periods and multiple base scenarios per period. It also incorporates additional variables and constraints to model load following reserves, ramping costs, and storage resources. A third generation of the multi-period SuperOPF, adds both integer variables and a receding horizon framework in which the problem type is more challenging (mixed integer), the size is even larger, and it must be solved more frequently, pushing the limits of currently available algorithms and solvers. The consideration of transient stability constraints in optimal power flow (OPF) problems has become increasingly important in modern power systems. Transient stability constrained OPF (TSCOPF) is a nonlinear optimization problem subject to a set of algebraic and differential equations. Solving a TSCOPF problem can be challenging due to (i) the differential-equation constraints in an optimization problem, (ii) the lack of a true analytical expression for transient stability in OPF. To handle the dynamics in TSCOPF, the set of differential equations can be approximated or converted into equivalent algebraic equations before they are included in an OPF formulation. In Chapter 4, a rigorous evaluation of using a predefined and fixed threshold for rotor angles as a mean to determine transient stability of the system is developed. TSCOPF can be modeled as a large-scale nonlinear programming problem including the constraints of differential-algebraic equations (DAE). Solving a TSCOPF problem can be challenging due to (i) the differential-equation constraints in an optimization problem, (ii) the lack of a true analytical expression for transient stability constraint in OPF. Unfortunately, even the current best TSCOPF solvers still suffer from the curse of dimensionality and unacceptable computational time, especially for large-scale power systems with multiple contingencies. In chapter 5, thse issues will be addressed and a new method to incorporate the transient stability constraints will be presented.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Bi-Level Linear Programming Model for Automatic Load Shedding: A Distributed Wide-Area Measurement System-based Solution

Load shedding is currently implemented as a two-step based approach. In the first step, manual load shedding is taken place, were system operators, using estimates, inform distribution utilities of predicted stressful conditions. Information provided include the potential use of energy reserves, as well as load shedding amount. In a second step, automatic load shedding is done. The latter is realized using protection relays. While considering frequency variation, pre-defined values of load to be shed and correspondent number of stages for such to be realized are transformed into relay settings. Under-frequency protection relays use only local measurements towards decision making, thus operate in a decentralized architecture. Decision making is done in milliseconds plus breaker time. While this approach has provided much system reliability, considering the new smart grid paradigm, where system dynamics are much faster due to increasing renewable resources penetration, in some operating conditions it will generate sub-optimal solutions, such as islanding. Phasor measurement units provide a source of information which can be useful for this problem. Centralized architecture-based solutions for automatic load shedding, as present in the state-of-the-art, require though total processing times which are not acceptable for real-life implementation. In this work, considering the above, a bi-level linear programming model is presented. The model is implemented considering a distributed architecture while leveraging phasor measurement units data. The upper-level model estimates the current system state. Results of this model are embedded in a lower-level model, which decision variables are the location and load value to be shed. Easy-to-implement model, built-on the classic weighted least squares solution, highlight potential aspects towards real-life applications.

Bretas, Arturo Suman↗
DOE OSTI2023

TEAMER: Experimental performance characterization of a shrouded axial-flow turbine

Sitkana has developed a shrouded hydrokinetic turbine with a modular, low-cost design that can be scaled to meet the needs of remote communities. With technical support from the University of Washington, Sitkana sought to experimentally characterize the mechanical power and structural loads of various 1:3.3 scale rotor geometries. In all, 11 different rotor geometries were characterized with variations in height-to-diameter ratio, blade number, and blade type (foiled versus flat). All tests were conducted in Reynolds-independent flow conditions in the Alice C. Tyler Flume at the University of Washington. Results allow Sitkana to (1) refine the optimal rotor geometry, (2) validate numerical models, and (3) predict power output for a full-scale system. This project is part of the TEAMER RFTS 8 (request for technical support) program.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Modeling Distributed Energy Resources for Analyzing Distribution Systems with High Renewable Penetration: Preprint

Increasing levels of inverter-based distributed energy resources (IBDERs) impact the legacy overcurrent distribution protection systems. The fault current injections of IBDERs are limited to between 1-2 p.u. of the rated current. The combination of the limited fault current and varying system load makes it increasingly difficult to reliably set the overcurrent protection. To address this challenge, researchers are proposing adaptive overcurrent protection mechanisms that can adapt to changing system conditions. Accurate modeling of IBDERs is required to develop reliable protection schemes. This paper presents a photovoltaics-based IBDER model in OpenDSS for adaptive overcurrent protection. The paper presents the validation of the proposed IBDER model against a detailed model in electromagnetic transient programs. Protection analysis of the Electric Power Research Institute J1 feeder with the proposed IBDER model is also presented.

adaptive↗
DOE OSTI2015

Correlation between inter-spin interaction and molecular dynamics of organic radicals in organic 1D nanochannels

One-dimensional (1D) molecular chains of 4-substituted-2,2,6,6-tetramethyl-1-piperidinyloxyl (4-X-TEMPO) radicals were constructed in the crystalline 1D nanochannels of 2,4,6-tris(4-chlorophenoxy)-1,3,5-triazine (CLPOT) used as a template. The ESR spectra of CLPOT inclusion compounds (ICs) using 4-X-TEMPO were examined on the basis of spectral simulation using EasySpin program package for simulating and fitting ESR spectra. The ESR spectra of [(CLPOT){sub 2}-(TEMPO){sub 1.0}] IC were isotropic in the total range of temperatures. The peak-to-peak line width (ΔB{sub pp}) became monotonically narrower from 2.8 to 1.3 mT with increase in temperature in the range of 4.2–298 K. The effect of the rotational diffusion motion of TEMPO radicals in the CLPOT nanochannels for the inter-spin interaction of the [(CLPOT){sub 2}-(TEMPO){sub 1.0}] IC was found to be smaller than the case of [(TPP){sub 2}−(TEMPO){sub 1.0}] IC (TPP = tris(o-phenylenedioxy)cyclotriphosphazene) reported in our previous study. The ΔB{sub pp} of the [(CLPOT){sub 2}-(TEMPO){sub 1.0}] IC in the whole range of temperatures was much narrower than the estimation to be based on the Van Vleck’s formula for the second moment of the rigid lattice model where the electron spin can be considered as fixed; 11 mT of Gaussian line-width component. This suggests the possibility of exchange narrowing in the 1D organic-radical chains of the [(CLPOT){sub 2}-(TEMPO){sub 1.0}] IC. On the other hand, the ESR spectra of [(CLPOT){sub 2}-(MeO-TEMPO){sub 0.41}] IC (MeO-TEMPO = 4-methoxy-TEMPO) were reproduced by a superposition of major broad isotropic adsorption line and minor temperature-dependent modulated triplet component. This suggests that the IC has the part of 1D organic-radical chains and MeO-TEMPO molecules isolated in the CLPOT nanochannels.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Enhanced Dual Confocal Measurement System

A measurement instrument utilizing dual, chromatic, confocal, distance sensors has been jointly developed by General Atomics (GA) and Sandia National Laboratories (SNL) for thickness and flatness measurement of target components used in dynamic materials properties (DMP) experiments on the SNL Z-Machine (Z). Compared to previous methods used in production of these types of targets, the tool saves time and yields a 4x reduction in thickness uncertainty which is one of the largest sources of error in equation of state (EOS) measurements critical to supporting the NNSA’s Stockpile Stewardship program and computer modeling of high energy density experiments. It has numerous differences from earlier instruments operating on the dual confocal sensor principle to accommodate DMP components including larger lateral travel, longer working distance, ability to measure flatness in addition to thickness, built-in thickness calibration standards for quickly checking calibration before and after each measurement, and streamlined operation. Thickness and flatness of 0.2mm-3.3mm thick sections of diamond machined copper and aluminum can be measured to “sub-micron” accuracy. Sections up to 6mm thick can be measured with as-yet undermined accuracy. Furthermore, samples must have one surface which is flat to within 300µm, lateral dimensions of no more than 50mm x 50mm, and height less than 40mm.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Large-scale parametric modeling of spent nuclear fuel dynamics in the 30 cm package drop scenario

Packages used to transport spent nuclear fuel (SNF) are required by the U.S. Code of Federal Regulations 10 CFR 71.71 to demonstrate satisfactory performance during a drop scenario. While the CFR is meant to ensure safe package function, it does not evaluate survival of the SNF within. The U.S. Department of Energy Spent Fuel and Waste Science and Technology program is working on closing the knowledge gap related to the response of SNF to external mechanical loads, including the hypothetical 30 cm package drop scenario in the CFR. In support of this effort, LS-DYNA finite element simulations were developed by Pacific Northwest National Laboratory (PNNL) to model generic drop scenarios at both the package and fuel assembly level. The models were validated against one-third scale package and full scale fuel assembly drop test data and were exercised to predict fuel cladding strains in a narrow range of model configurations. This work describes a large-scale parametric study conducted by PNNL using the previously developed and validated PWR finite element model, with the addition of a new generic BWR assembly model. The motivation for the parametric study was to characterize the broad range of SNF responses in the 30 cm package drop scenario. This was accomplished by varying the drop orientation, fuel assembly type (17x17 PWR and 10x10 BWR), burnup, cladding temperature, spacer grid buckling load, package mass, impact limiter stiffness, and mechanical gap conditions within the basket. A MATLAB framework was developed to automate LS-DYNA model generation and execution on PNNL institutional computing resources. In total, over 2000 simulations were performed. For each simulation, the SNF response was quantified in terms of permanent grid deformation, fuel rod contact pressure, and strains within the fuel rods, guide tubes, and water rods. The results provide valuable insight into the range of responses that could be reasonably expected from SNF in the hypothetical drop scenario, as well as the sensitivity to each input parameter. The results of this parametric study are a key component of the testing and modeling strategy the Spent Fuel and Waste Science and Technology program is using to close the external loads knowledge gap.

Kadooka, Kevin↗
DOE OSTI2022

Improving Milling Performance on New Computer Numerical Control Machining Centers with Machining Dynamics

This project is directed at improving the efficiency and productivity of the existing installed US machine tool base by correcting a fundamental system-wide flaw in the existing programmable machine tool infrastructure. The project will catalyze a shift from a geometry-only mindset to a geometry -and-machine -performance mindset. The project will lead to a dramatic overhaul of the existing ecosystem for programming of Computer Numerically Controlled (CNC) machine tools, and will result in energy savings, cost savings, reduced scrap, and improved competitiveness of US manufacturers.

42 ENGINEERING↗
DOE OSTI2023

Pseudospectral convex optimization for on-ramp merging control of connected vehicles

It can be a daunting task for human drivers to merge into highways because of the intricate vehicle negotiations and potential risk within limited time and space. Connected vehicle (CV) technologies could be a solution to this problem and offer many benefits to the road safety, traffic mobility, and energy efficiency. However, real-time optimal control of CVs is still an open challenge, due to the nonlinear vehicle dynamics, non-convex fuel consumption model, and highly dynamic uncertain inter-vehicle interactions. To tackle these issues, a novel real-time optimal control approach that balances the computational efficiency and solution optimality is proposed for the purpose of onboard application. To this end, the pseudospectral collocation method is integrated with a sequential convex programming approach to develop two new optimization algorithms, which are implemented within a model predictive control (MPC) framework to allow for real-time generation of optimal merging speed profiles. One algorithm leverages the line search technique to improve convergence, and the other benefits from the trust region method for better computational efficiency. The optimality and convergence process of both proposed algorithms are investigated by comparing their solutions with a popular non-linear solver. Furthermore, simulation results show that the proposed methods outperform the benchmark in terms of computational cost, fuel consumption, and traffic efficiency. In particular, the proposed fuel-economy merging rule can save 57.1% fuel consumption on average on four different traffic volumes. Meanwhile, the proposed optimal control algorithms can reduce 2.2% travel time on average comparing to the “first-in-first-out” merging rule.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2019

Elucidation of the Reaction Mechanism for High-Temperature Water Gas Shift over an Industrial-Type Copper–Chromium–Iron Oxide Catalyst

The water gas shift (WGS) reaction is of critical importance for the chemical industry, as it constitutes, coupled with methane reforming, the main industrial route to produce hydrogen. Copper–chromium–iron oxide-based catalysts have been widely used for the high-temperature WGS reaction industrially. The WGS reaction mechanism by the CuCrFeO x catalyst has been debated for years, mainly between a “redox” mechanism involving the participation of atomic oxygen from the catalyst and an “associative” mechanism proceeding via a surface formate-like intermediate. Here, advanced in situ characterization techniques (infrared spectroscopy, temperature-programmed surface reaction (TPSR), near-ambient pressure XPS (NAP-XPS), and inelastic neutron scattering (INS)) were applied to determine the nature of the catalyst surface and identify surface intermediate species under WGS reaction conditions. The surface of the CuCrFeO x catalyst is found to be dynamic and becomes partially reduced under WGS reaction conditions, forming metallic Cu nanoparticles on Fe 3 O 4 . Neither in situ IR not INS spectroscopy detect the presence of surface formate species during WGS. TPSR experiments demonstrate that the evolution of CO 2 and H 2 from the CO/H 2 O reactants follows different kinetics than the evolution of CO 2 and H 2 from HCOOH decomposition (molecule mimicking the associative mechanism). Steady-state isotopic transient kinetic analysis (SSITKA) (CO + H 2 16 O → CO + H 2 18 O) exhibited significant 16 O/ 18 O scrambling, characteristic of a redox mechanism. Computed activation energies for elementary steps for the redox and associative mechanism by density functional theory (DFT) simulations indicate that the redox mechanism is favored over the associative mechanism. The combined spectroscopic, computational, and kinetic evidence in the present study finally resolves the WGS reaction mechanism on the industrial-type high-temperature CuCrFeO x catalyst that is shown to proceed via the redox mechanism.

SSITKA↗
DOE OSTI2022

Coolant Channel Design for Additively Manufactured Reactor Cores

Additive manufacturing (AM) methods are currently being explored for applications in nuclear reactors to make advanced reactors more efficient, safe, and reliable. The Transformational Challenge Reactor (TCR) program has explored AM for nuclear by designing a high-temperature gas reactor (HTGR) using an AM silicon carbide fuel form with uranium nitride–tristructural isotropic fuel. In this study, we detail the design process for the TCR fuel form’s coolant channels using computational fluid dynamics models with conjugate heat transfer. Additionally, this work discusses how these models were interfaced with other design teams, project milestones, and the agile design method used to mature the reactor design. The methodology deployed was able to create a channel design with lower maximum fuel temperatures and thermal stresses in the fuel form over traditional channel designs that can be manufactured subtractively. These results were achieved with only small manufacturing penalties. Results are discussed and presented on lessons learned for designing AM components for nuclear reactors. Finally, areas of opportunity are discussed for advanced design tools to further automate design activities and optimize reactors with fewer built-in assumptions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Experimental targets for dark photon dark matter

Ultralight dark photon dark matter features distinctive cosmological and astrophysical signatures and is also supported by a burgeoning direct-detection program searching for its kinetic mixing with the ordinary photon over a wide mass range. Dark photons, however, cannot necessarily constitute the dark matter in all of this parameter space. In minimal models where the dark photon mass arises from a dark Higgs mechanism, early-Universe dynamics can easily breach the regime of validity of the low-energy effective theory for a massive vector field. In the process, the dark sector can collapse into a cosmic string network, precluding dark photons as viable dark matter. We establish the general conditions under which dark photon production avoids significant backreaction on the dark Higgs and identify regions of parameter space that naturally circumvent these constraints. After surveying implications for known dark photon production mechanisms, we propose novel models that set well-motivated experimental targets across much of the accessible parameter space. We also discuss complementary cosmological and astrophysical signatures that can probe the dark sector physics responsible for dark photon production.

Dark matter↗
DOE OSTI2012

Real-time Pricing Demand Response in Operations

Abstract—Dynamic pricing schemes have been implemented in commercial and industrial application settings, and recently they are getting attention for application to residential customers. Time-of-use and critical-peak-pricing rates are in place in various regions and are being piloted in many more. These programs are proving themselves useful for balancing energy during peak periods; however, real-time (5 minute) pricing signals combined with automation in end-use systems have the potential to deliver even more benefits to operators and consumers. Besides system peak shaving, a real-time pricing system can contribute demand response based on the locational marginal price of electricity, reduce load in response to a generator outage, and respond to local distribution system capacity limiting situations. The US Department of Energy (DOE) is teaming with a mid-west electricity service provider to run a distribution feeder-based retail electricity market that negotiates with residential automation equipment and clears every 5 minutes, thus providing a signal for lowering or raising electric consumption based on operational objectives of economic efficiency and reliability. This paper outlines the capability of the real-time pricing system and the operational scenarios being tested as the system is rolled-out starting in the first half of 2012.

Demand Response↗
DOE OSTI2016

TH-EF-BRB-04: 4π Dynamic Conformal Arc Therapy Dynamic Conformal Arc Therapy (DCAT) for SBRT

Purpose: To develop an efficient and effective trajectory optimization methodology for 4π dynamic conformal arc treatment (4π DCAT) with synchronized gantry and couch motion; and to investigate potential clinical benefits for stereotactic body radiation therapy (SBRT) to breast, lung, liver and spine tumors. Methods: The entire optimization framework for 4π DCAT inverse planning consists of two parts: 1) integer programming algorithm and 2) particle swarm optimization (PSO) algorithm. The integer programming is designed to find an optimal solution for arc delivery trajectory with both couch and gantry rotation, while PSO minimize a non-convex objective function based on the selected trajectory and dose-volume constraints. In this study, control point interaction is explicitly taken into account. Beam trajectory was modeled as a series of control points connected together to form a deliverable path. With linear treatment planning objectives, a mixed-integer program (MIP) was formulated. Under mild assumptions, the MIP is tractable. Assigning monitor units to control points along the path can be integrated into the model and done by PSO. The developed 4π DCAT inverse planning strategy is evaluated on SBRT cases and compared to clinically treated plans. Results: The resultant dose distribution of this technique was evaluated between 3D conformal treatment plan generated by Pinnacle treatment planning system and 4π DCAT on a lung SBRT patient case. Both plans share the same scale of MU, 3038 and 2822 correspondingly to 3D conformal plan and 4π DCAT. The mean doses for most of OARs were greatly reduced at 32% (cord), 70% (esophagus), 2.8% (lung) and 42.4% (stomach). Conclusion: Initial results in this study show the proposed 4π DCAT treatment technique can achieve better OAR sparing and lower MUs, which indicates that the developed technique is promising for high dose SBRT to reduce the risk of secondary cancer.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Emulating Multiple Inheritance in Fortran 2003/2008

Although the high-performance computing (HPC) community increasingly embraces object-oriented programming (OOP), most HPC OOP projects employ the C++ programming language. Until recently, Fortran programmers interested in mining the benefits of OOP had to emulate OOP in Fortran 90/95. The advent of widespread compiler support for Fortran 2003 now facilitates explicitly constructing object-oriented class hierarchies via inheritance and leveraging related class behaviors such as dynamic polymorphism. Although C++ allows a class to inherit from multiple parent classes, Fortran and several other OOP languages restrict or prohibit explicit multiple inheritance relationships in order to circumvent several pitfalls associated with them. Nonetheless, what appears as an intrinsic feature in one language can be modeled as a user-constructed design pattern in another language. The present paper demonstrates how to apply the facade structural design pattern to support a multiple inheritance class relationship in Fortran 2003. The design unleashes the power of the associated class relationships for modeling complicated data structures yet avoids the ambiguities that plague some multiple inheritance scenarios.

97 MATHEMATICS AND COMPUTING↗