Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “optimal control theory”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Understanding Catalyst ‘Volcano’ Dependence Through Fermi-Level Controlled Kinetics Using Electronic Theory

The ubiquitous two-step Michaelis–Menten and Temkin–Boudart reaction mechanisms are extended to include the influence of the catalyst electronic subsystem in a 5-step mechanism. The resulting kinetic equation provides an alternative explanation for the well-known volcano-shaped dependence found in catalysis. The equilibrium constants of fast electronic steps are highlighted for their influence on adsorption and desorption through the relative concentration of charged versus neutral intermediates. This generalized concept can be widely applied to determine the optimal catalyst, based on the Fermi level of the material, for reactions proceeding via this universal reaction.

25 ENERGY STORAGE↗
DOE OSTI2023

Controlling Molecular Structure and Spin with Multiconfigurational Quantum Chemistry (Final Technical Report)

For many first-row transition metal complexes, structure-property relationships can be obtained from high-level molecular geometry optimizations and subsequent electronic structure studies. The project funded under this award utilized newly implemented fully internally contracted (FIC) nuclear gradients for extended multi-state (XMS) complete active space second-order multireference perturbation theory (CASPT2) to explore geometry changes in first-row transition metal coordination complexes. At the start of the project, only one full geometry optimization using FIC-CASPT2 analytical gradients had been reported (J. Chem. Theory Comput., 2016, 12 (8), 3781), and many open-questions regarding the performance and achievable accuracy in applying such computations to larger complexes persisted. Key deliverables in the report include the implementation of the numerical Hessian and subsequent vibrational analysis in the BAGEL program package. This includes both full Hessian vibrational analysis (FHVA) and a partial Hessian vibrational analysis (PHVA). The latter of which allows us to reduce the significant cost of computing the vibrational frequencies in a large molecule by focusing on the modes of highest interest. The gradient code also proved useful in the development of an approach to improve the Hubbard-U correction by removing bias in predicting spin-splitting in Fe(II) complexes via plane-wave density functional theory (DFT). Finally, we showed for three families of complexes (spin-crossover (SCO) complexes, metallocorroles complexes, and those with metal-metal bonds) that CASPT2 can result in good molecular geometries and established best practices in undertaking this work. We are now using this approach in collaborative projects with experimental groups, which would not have been possible at the start of this project. A common theme also arose through this work that static and dynamic correlation must be recovered in a balanced way to yield quantitative results. By systematically varying how both effects were included, we were able to make key chemical insights. This work supported the training of two postdoctoral scholars, one graduate student, and one undergraduate student in applying high-level wavefunction based methods to challenging systems.

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

Development of a Reactive Force Field for Simulating Photoinitiated Acrylate Polymerization

Light-driven and photo-curable polymer based additive manufacturing (AM) has enormous potential due to its excellent resolution and precision. Acrylated radical chain-growth polymerized resins are widely used in photopolymer AM due to their fast kinetics, and often serve as a departure point for developing other resin materials for photopolymer-based AM technologies. For successful control of the photopolymer resins, the molecular basis of the acrylate free-radical polymerization has to be understood in detail. We present an optimized reactive force field (ReaxFF) for molecular dynamics (MD) simulations of acrylate polymer resins that captures radical polymerization thermodynamics and kinetics. The force field is trained against an extensive training set including density functional theory (DFT) calculations of reaction pathways along the radical polymerization from methyl acrylate to methyl butyrate, bond dissociation energies, and structures and partial charges of several molecules and radicals. We also found that it was critical to train the force field against an incorrect, nonphysical reaction pathway observed in simulations that used parameters not optimized for acrylate polymerization. As a result, the parameterization process utilizes a parallelized search algorithm, and the resulting model can describe polymer resin formation, crosslinking density, conversion rate, and residual monomers of the complex acrylate mixtures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Understanding plasma turbulence through exact coherent structures

Plasma turbulence is a key challenge in understanding transport phenomena in magnetically confined plasmas. This work presents a generalized framework to analyze plasma turbulence that utilizes periodic orbit theory. In periodic orbit theory, doubly periodic solutions (coherent structures) of the governing equation(s) serve as building blocks of the considered turbulent dynamics. To illustrate the concept and method, the particularly simple Kuramoto–Sivashinsky (referred to here as LMRT for the original authors: LaQuey, Mahajan, Rutherford, and Tang) trapped-ion mode toy model is used. By applying numerical optimization techniques to the LMRT equation, we extract coherent spacetime patterns that represent the library of allowable fundamental structures of the equation. These structures provide a framework to systematically describe turbulence as a composition of recurrent solutions, revealing an underlying order within chaotic plasma motion. Although illustrated here using the simplified LMRT model for clarity, this framework provides a general strategy that can be extended to more complex and realistic models of plasma turbulence, including gyrokinetic systems. This offers a new method for predicting and potentially controlling transport processes in fusion plasmas by providing a bridge between nonlinear dynamical systems theory and plasma physics in the form of a generalized framework with which to analyze and understand spatially extended nonlinear partial differential equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Mitigating Green's function Monte Carlo signal-to-noise problems using contour deformations

The Green's function Monte Carlo (GFMC) method provides accurate solutions to the nuclear many-body problem and predicts properties of light nuclei starting from realistic two- and three-body interactions. Controlling the GFMC fermion sign problem is crucial, as the signal-to-noise ratio decreases exponentially with imaginary time, requiring significant computing resources. Inspired by similar scenarios in lattice quantum field theory and spin systems, in this work, we employ integration contour deformations to improve the GFMC signal-to-noise ratio. Machine learning techniques are used to select optimal contours with minimal variance from parametrized families of deformations. As a proof of principle, we consider the deuteron binding energies and Euclidean density response functions. Here, we only observe mild signal-to-noise improvement for the binding energy case. On the other hand, we achieve an order of magnitude reduction of the variance for Euclidean density response functions, paving the way for computing electron- and neutrino-nucleus cross sections of larger nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Factors influencing quantum evaporation of helium from polar semiconductors from first principles

While there is much indirect evidence for the existence of dark matter (DM), to date it has evaded detection. Current efforts focus on DM masses over ~GeV—to push the sensitivity of DM searches to lower masses, new DM targets and detection schemes are needed. In this work, we focus on the latter—a novel detection scheme recently proposed to detect 10–100 meV phonons in polar target materials. Previous work showed that well-motivated models of DM can interact with polar semiconductors to produce an athermal population of phonons. This new sensing scheme proposes that these phonons then facilitate quantum evaporation of 3 He from a van der Waals film deposited on the target material. However, a fundamental understanding of the underlying process is still unclear, with several uncertainties related to the precise rate of evaporation and how it can be controlled. In this work, we use ab initio density functional theory calculations to compare the adsorption energies of helium atoms on a polar target material, sodium iodide, to understand the underlying evaporation physics. We explore the role of surface termination, monolayer coverage, and elemental species on the rate of He evaporation from the target material. Using this, we discuss the optimal target features for He-evaporation experiments and their range of tunability through chemical and physical modifications such as applied field and surface termination.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The Compelling Need for a Mid-Scale Stellarator Facility

In the pursuit of the goal of commercial fusion as an abundant and safe source of energy, the stellarator is a leading concept with compelling attractiveness and demonstrated performance. A new mid-size stellarator is needed to retire risks and innovate towards a high performance, economically attractive, stellarator Fusion Pilot Plant. In this presentation we, as a community of US researchers from Universities, National Laboratories, and Private Industry, involved in studying the stellarator concept, lay out the programmatic and technical motivation for a new and modern mid-size stellarator research facility. A new mid-scale stellarator is needed to realize the potential predicted by a solid body of theory and simulation along with advances in computational tools for optimization and non-linear turbulence modeling. Notably, it is possible to combine the advantages of the stellarator (steady state, no current drive, no disruptions) with the good confinement regularly achieved in tokamaks. The top priorities for experimental work, and the motivation for a mid-scale stellarator experiment are: turbulence control, non-resonant divertor, MHD stability at large beta, confinement of fast particles, and coil simplification. A new mid-size quasi-symmetric stellarator, built as a user facility, would complement existing research at Wendelstein 7-X and LHD and strongly augment private industry. It would provide a program of innovative research, concept validation, theoretical advancement, and workforce development. Growing support and interest for stellarators by the fusion community and private industry affirms this rationale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Theoretical Optimization of Trapped-Bubble-Based Acoustic Metamaterial Performance

Acoustic metamaterials have proven to be a versatile tool for the precise control and manipulation of sound waves. One of the promising designs of acoustic metamaterials employ the arrays of bubbles and find applications for soundproofing, blast mitigation, and many others. An obvious advantage of bubble-based metamaterials is their ability to be relatively thin while absorbing low-frequency sound waves. The vast majority of theories developed to predict resonant behavior of bubble-based metamaterials capitalize on Minnaert frequency. Here, we propose a novel theoretical approach to characterize bubble-based metamaterials that are based on our previous findings for a single bubble trapped in circular cavity modeled as a thin clamped plate. We obtain analytical expressions for resonant frequencies of bubble metascreens using self-consistent approximation. Two geometry factors, distance between bubble centers and distance between bubble center and interface of acoustic impedance change, are taken into account. We demonstrate the existence of multiple bandgaps and possibility of switching between them via adjustment of geometry parameters and reflector properties.

Gritsenko, Dmitry (ORCID:0000000248181512)↗
DOE OSTI2023

Mitigating Green's function Monte Carlo signal-to-noise problems using contour deformations

The Green's function Monte Carlo (GFMC) method provides accurate solutions to the nuclear many-body problem and predicts properties of light nuclei starting from realistic two- and three-body interactions. Controlling the GFMC fermion-sign problem is crucial, as the signal-to-noise ratio decreases exponentially with Euclidean time, requiring significant computing resources. Inspired by similar scenarios in lattice quantum field theory and spin systems, in this work, we employ integration contour deformations to improve the GFMC signal-to-noise ratio. Machine learning techniques are used to select optimal contours with minimal variance from parameterized families of deformations. As a proof of principle, we consider the deuteron binding energies and Euclidean density response functions. We only observe mild signal-to-noise improvement for the binding energy case. On the other hand, we achieve an order of magnitude reduction of the variance for Euclidean density response functions, paving the way for computing electron- and neutrino-nucleus cross-sections of larger nuclei.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A General Framework for Bounding Approximate Dynamic Programming Schemes

For years, there has been interest in approximation methods for solving dynamic programming problems, because of the inherent complexity in computing optimal solutions characterized by Bellman’s principle of optimality. A wide range of approximate dynamic programming (ADP) methods now exists. It is of great interest to guarantee that the performance of an ADP scheme be at least some known fraction, say ß , of optimal. This letter introduces a general approach to bounding the performance of ADP methods, in this sense, in the stochastic setting. The approach is based on new results for bounding greedy solutions in string optimization problems, where one has to choose a string (ordered set) of actions to maximize an objective function. This bounding technique is inspired by submodularity theory, but submodularity is not required for establishing bounds. Instead, the bounding is based on quantifying certain notions of curvature of string functions; the smaller the curvatures the better the bound. The key insight is that any ADP scheme is a greedy scheme for some surrogate string objective function that coincides in its optimal solution and value with those of the original optimal control problem. The ADP scheme then yields to the bounding technique mentioned above, and the curvatures of the surrogate objective determine the value ß of the bound. The surrogate objective and its curvatures depend on the specific ADP.

discrete event systems↗
DOE OSTI2020

A software‐defined networks ‐based measurement method of network traffic for 6G technologies

Abstract In software‐defined networks (SDN) for 6G technologies, the controller should accurately and efficiently measure the flow traffic in switches for traffic engineering. Fine‐grained flow measurement can more accurately describe the traffic in the network, but it also consumes much more resources. To reduce the overhead incurred in the measurement process and obtain the approximate fine‐grained measurements, we propose a novel lightweight measurement scheme that runs in the controller. The novel lightweight measurement architecture consists of two parts: coarse‐grained measurement and interpolation‐optimization. In the first part, based on the SDN architecture, we use the pull‐based random sampling method to quickly obtain the coarse‐grained measurement of flow traffic through OpenFlow protocol. In the second part, we insert some discrete values into the coarse‐grained measurement with the interpolation theory, then we optimize interpolation results until finding the optimal fine‐grained flow traffic measurement by utilizing the multiconstraint method. We verified the feasibility of the proposed measurement method, and simulation results show that the measurement error of the proposed method is under 25%, but the flow measurement overhead of this method occupies only 3.3% compared with that of the fine‐grained method.

Huo, Liuwei↗
DOE OSTI2024

Critical role of solvation on CC13 porous organic cages for design of porous liquids

Efficient carbon capture requires the design of new materials with high CO 2 selectivity and gas adsorption capacity that can be incorporated into existing industrial processes. Porous liquids (PLs) are promising candidate materials that consist of a nanoporous host and a solvent forming a liquid with permanent porosity based on exclusion of the solvent from the interior of the nanoporous host. Stable PLs are based on solvent-nanoporous host interactions, which can be evaluated through molecular simulations. Here, time- and temperature-dependent density functional theory simulations were performed between four solvents, 2-bromophenol, 4-methylphenol, 2,4-dimethylphenol, and cyclohexanone and the CC13 porous organic cage (POC) as a prototypical PL composition. Overall, minimal reactions occurred in the PL including no changes in the POC structure. Additionally, POC-solvent coordination occurred through interactions of neighboring functional groups such as methyl/bromide and hydroxyl on the solvent molecules with the POC surface. Therefore, the location rather than the number of functional groups on the solvent molecule controls the POC-solvent interactions. Additionally, the POC pore window contracted or expanded up to 8% during solvation, which correlates with the experimental solubility and static solvent-POC binding, where solvents that caused less contraction of the POC pore window increased POC solubility. Finally, these results allow for the design of optimized POC-based PL compositions based on solvent-nanoporous host binding and variation in the pore window during solvation.

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

Suppressing CO formation in low-temperature methanol steam reforming via Ce-modified CuZnGa layered oxide catalysts

Cu-based layered double hydroxides (LDHs) are widely recognized as effective catalysts for low-temperature methanol steam reforming, yet achieving high hydrogen productivity together with near-complete suppression of CO formation remains challenging. Here, we report the synthesis and evaluation of a series of CuZnGa LDH-derived catalysts and Ce-modified analogues prepared via an aqueous miscible organic method, which enables high metal dispersion and precise structural control. The optimized CuZnGa catalyst exhibits a hydrogen production rate of 16.9 µmol H 2 ·g cat −1 ·s −1 at 180 °C with an H 2 /CO ratio exceeding 3500, outperforming many state-of-the-art low-temperature systems. Importantly, the incorporation of small amounts of Ce further suppresses CO formation while maintaining high hydrogen productivity. Combined spectroscopic characterization and density functional theory calculations reveal that Ce is incorporated into the LDH lattice by substituting Ga 3+ sites up to a critical threshold, beyond which highly dispersed CeO x species are formed. These species provide mobile lattice oxygen that participates in a Mars-van Krevelen-type pathway, selectively oxidizing CO and suppressing the reverse water-gas shift reaction. This study establishes a clear relationship between Ce speciation, oxygen mobility, and catalytic selectivity in LDH-derived systems. The resulting catalysts demonstrate the potential of interface-engineered Cu-based materials for efficient low-temperature hydrogen production with minimal CO contamination.

09 BIOMASS FUELS↗
DOE OSTI2023

Adaptive Time Step Control for Multirate Infinitesimal Methods

Multirate methods have been used for decades to temporally evolve initial-value problems in which different components evolve on distinct time scales, and thus use of different step sizes for these components can result in increased computational efficiency. Generally, such methods select these different step sizes based on experimentation or stability considerations. For problems that evolve on a single time scale, adaptivity approaches that strive to control local temporal error are widely used to achieve numerical results of a desired accuracy with minimal computational effort, while alleviating the need for manual experimentation with different time step sizes. However, there is a notable gap in the publication record on the development of adaptive time step controllers for multirate methods. In this paper, we extend the single-rate controller work of Gustafsson [ACM Trans. Math. Software, 20 (1994), pp. 496-517] to the multirate method setting. Specifically, we develop controllers based on polynomial approximations to the principal error functions for both the "fast" and "slow" time scales within multirate infinitesimal (MRI) methods. We additionally investigate a variety of approaches for estimating the errors arising from each time scale within MRI methods. We then numerically evaluate the proposed multirate controllers and error estimation strategies on a range of multirate test problems, comparing their performance against an estimated optimal performance. Through this work, we combine the most performant of these approaches to arrive at a set of multirate adaptive time step controllers that robustly achieve desired solution accuracy with minimal computational effort.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Chemical controls on the propagation and healing of subcritical fractures

Human activities involving subsurface reservoirs—resource extraction, carbon and nuclear waste storage—alter thermal, mechanical, and chemical steady-state conditions in these systems. Because these systems exist at lithostatic pressures, even minor chemical changes can cause chemically assisted deformation. Therefore, understanding how chemical effects control geomechanical properties is critical to optimizing engineering activities. The grand challenge in predicting the effect of chemical processes on mechanical properties lays in the fact that these phenomena take place at molecular scales, while they manifest all the way to reservoir scales. To address this fundamental challenge, we investigated chemical effects on deformation in model and real systems spanning molecular- to centimeter scales. We used theory, experiment, molecular dynamics simulation, and statistical analysis to (1) identify the effect of simple reactions, such as hydrolysis, on molecular structures in interfacial regions of stressed geomaterials; (2) quantify chemical effects on the bulk mechanical properties, fracture and displacement for granular rocks and single crystals; (3) develop initial understanding of universal scaling for individual displacement events in layered geomaterials; and (4) develop analytic approximations for the single-chain mechanical response utilizing asymptotically correct statistical thermodynamic theory. Taken together, these findings advance the challenging field of chemo-mechanics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Breaking Rayleigh’s Law with Spatially Correlated Disorder to Control Phonon Transport

Controlling thermal transport in insulators and semiconductors is crucial for many technological fields such as thermoelectrics and thermal insulation, for which a low thermal conductivity (κ) is desirable. A major obstacle for realizing low κ materials is Rayleigh’s law, which implies that acoustic phonons, which carry most of the heat, are insensitive to scattering by point defects at low energy. We demonstrate, with large scale simulations on tens of millions of atoms, that isotropic long-range spatial correlations in the defect distribution can dramatically reduce phonon lifetimes of important low-frequency heat-carrying modes, leading to a large reduction of κ—potentially an order of magnitude at room temperature. We propose a general and quantitative framework for controlling thermal transport in complex functional materials through structural spatial correlations, and we establish the optimal functional form of spatial correlations that minimize κ. Finally, we end by briefly discussing experimental realizations of various correlated structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Sensor and Actuator Attacks on Hierarchical Control Systems with Domain-Aware Operator Theory

Cyber-Physical Systems (CPSs) provide opportunities for cyber attacks to have physical impacts. Advanced Persistent Threats (APTs) are a subclass of cyber threats that act stealthily to avoid detection and enable long-term attacks. Here, we build on our past work in APT modelling to combine deception-based sensor bias attacks and direct actuator manipulations in attacks against a hierarchical control system. That past work used the Koopman operator to develop a data-driven, domain-aware, optimization-based attacker model. Using an expansion of this model, we compute several different attacks, including multiple simultaneous attacks, against a high-fidelity commercial building emulator and compare the impacts of those attacks to each other. One next step of interest is to construct a defender system, built on the same modelling approach, designed to detect and mitigate such attacks.

koopman operator, Cyber-Physical Security, machine↗
DOE OSTI2026

Methanol Partial Oxidation on Cu(111) and PtCu(111) Single-Atom Alloy Surfaces: Effect of Surface Oxygen Coverage on Selectivity

The selective oxidation of methanol to formaldehyde on Cu surfaces is an important and well-studied reaction. However, a systematic analysis of product selectivity as a function of oxygen coverage on Cu(111) and Cu-based single-atom alloys (SAAs) has not been previously reported. In this work, we present a comprehensive investigation of deuterated methanol (CD 3 OH) partial oxidation on Cu(111) and 1% PtCu(111) SAA surfaces as a function of preadsorbed oxygen coverage. Temperature-programmed desorption (TPD) and X-ray photoelectron spectroscopy (XPS) reveal that isolated Pt atoms inhibit the initial surface oxidation of Cu(111) under low oxygen exposures. Despite this difference in oxidation kinetics, the product selectivity remains largely unaffected: on both Cu(111) and PtCu(111), formaldehyde (CD 2 O) is the predominant partial oxidation product over a broad range of oxygen coverages. The selectivity toward formaldehyde peaks at intermediate oxygen coverages (∼0.3 monolayers, ML), indicating the existence of an optimal oxygen loading for partial oxidation. Notably, the similar product selectivities on Cu(111) and PtCu(111) over a range of surface oxygen coverage indicate that Pt single atoms do not significantly alter the reaction pathway or shift the optimal oxygen coverage for formaldehyde formation. Control experiments confirm that Cu(111) is unreactive toward methanol in the absence of oxygen, while PtCu(111) surfaces produce a small amount of formaldehyde even when oxygen is not preadsorbed, indicating that isolated Pt atoms facilitate O−H activation at below 150 K, leading to H 2 desorption, followed by C−D activation at higher temperatures (∼350 K). Density functional theory (DFT)-based calculations show that Pt atoms increase the O 2 dissociation barrier relative to Cu(111), consistent with the observed inhibition of oxidation at low exposures. Overall, this work provides the first detailed selectivity map for methanol oxidation on oxidized Cu(111) and PtCu(111) SAA surfaces. By linking classical mechanistic insights such as methoxy- and formate-mediated pathways with single-atom alloy catalyst design, this work demonstrates that while Pt substitution modulates the oxidation kinetics and oxygen binding, the overall selectivity toward formaldehyde is governed primarily by oxygen coverage. These findings underscore the potential of isolated dopants to tune surface oxidation behavior without compromising the intrinsic partial oxidation selectivity of copper-based catalysts.

36 MATERIALS SCIENCE↗