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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Approximate Dynamic Programming With Enhanced Off-Policy Learning for Coordinating Distributed Energy Resources

Herein this paper proposes an innovative approximate dynamic programming (ADP) method for distributed energy resource coordination with the loss of life of battery energy storage system (BESS) explicitly modeled. The dispatch policy is designed to account for both calendrical and cyclical aging effects on BESS, explicitly modeling the impacts of ambient temperature on BESS lifespan. The proposed ADP employs an adaptive critic method and enhanced off-policy deterministic policy gradient (DPG) strategy, addressing the limitations of the on-policy gradient-based ADP approaches, including inadequate exploration, low data usage, and computational complexity. In particular, a customized policy is proposed to guide the algorithm to explore some promising decisions and thereby improve exploration capability and learning efficiency compared to conventional DPG-based learning approaches, which may struggle to find a global optimum due to random noisy action-based exploration or require expert demonstration with extra effort. The proposed method is illustrated using the IEEE 123-node system and compared with the existing ADP methods to prove solution accuracy and demonstrate the effects of incorporating degradation models into control design. Case studies showed that the proposed ADP effectively coordinates DERs with a 10 times smaller optimization gap compared to existing methods, and the incorporation of the BESS life loss model ensures the expected lifespan and results in significant cost savings.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Physics-Informed Graph Neural Networks for Collaborative Dynamic Reconfiguration and Voltage Regulation in Unbalanced Distribution Systems

Network reconfiguration has long been employed as a strategic approach to minimize power distribution system losses and effectively regulate voltage levels. Tap-changing voltage regulators are also critical for controlling bus voltages, especially in accommodating the increasing integration of distributed energy resources (DERs) with intermittent outputs. This paper introduces novel methodologies to address the challenges of dynamic reconfiguration and optimal tap setting in unbalanced three-phase distribution systems. We propose an approximated mixed-integer quadratically constrained program (MIQCP) to model dynamic reconfiguration, along with a pioneering formulation for voltage regulator (VR) tap-setting based on Special Ordered Set type 1 (SOS1). To mitigate computational complexity, we propose a physics-informed spatial-temporal graph convolutional network (STGCN) with an integrated link classifier. The proposed approach enables efficient solution generation by fixing specific variables in the MIQCP instance and solving the simplified sub-MIP using an MIP solver. Numerical studies demonstrate the superior prediction accuracy of our STGCN model compared to baseline neural network models, resulting in reduced DER curtailment and voltage deviation with shorter computation time.

dynamic reconfiguration↗
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 OSTI2025

GPR_calculator: An on-the-fly surrogate model to accelerate massive nudged elastic band calculations

We present GPR_calculator, a package based on Python and C++ programming languages to build an on-the-fly surrogate model using Gaussian Process Regression (GPR) to approximate computationally expensive electronic structure calculations. The key idea is to dynamically train a GPR model during the simulation that can accurately predict energies and forces with uncertainty quantification. When the uncertainty is high, the costly electronic structure calculation is performed to obtain the ground truth data, which is then used to update the GPR model. To illustrate the effectiveness of GPR_calculator, we demonstrate its application in Nudged Elastic Band (NEB) simulations of surface diffusion and reactions, achieving 3-10 times acceleration compared to pure ab initio calculations. The source code is available at https://github.com/MaterSim/GPR_calculator.

Gaussian process regression↗
DOE OSTI2023

Plant-specific Model and Data Analysis using Dynamic Security Modeling and Simulation

The requirements for U.S. nuclear power plants to maintain a large on-site physical security force contribute to their high operational costs. The cost of maintaining the current physical security posture is approximately 10% of the overall operation and maintenance budget for commercial nuclear power plants. The goal of the Light Water Reactor Sustainability (LWRS) program’s physical security pathway is to develop tools, methods, and technologies and provide the technical basis for an optimized physical security posture. The conservatisms built into current security postures may be analyzed and minimized in order to reduce security costs while still ensuring adequate security and operational safety. The research performed at Idaho National Laboratory within LWRS program’s physical security pathway has successfully developed a dynamic force-on-force modeling framework using various computer simulation tools and integrating them with the dynamic assessment Event Modeling Risk Assessment using Linked Diagrams (EMRALD) tool. This document provides an update on the progress in applying a dynamic computational framework that links results from a commercially available force-on-force simulation tool, a commercially available thermal-hydraulic tool, and EMRALD to an operating commercial nuclear power plant. This report is only a summary of the progress and does not contain specific modeling results as those contain sensitive security information. This process of including plant procedures and multiple analysis results is being called Modeling and Analysis for Safety Security using Dynamic EMRALD Framework or MASS-DEF. Previous reports described how a user could integrate their plant-specific force-on-force models with the dynamic simulation tool EMRALD, model operator actions, integrate with probabilistic risk assessment tools, such as CAFTA (Computer Aided Fault Tree Analysis System) or SAPHIRE (Systems Analysis Programs for Hands-on Integrated Reliability Evaluations), and with thermal-hydraulic tools, such as RELAP-5. Previous reports applied various combinations of available simulations codes with EMRALD using generic plant models to demonstrate how to perform the analysis. This report documents the results of applying the dynamic computational framework to an actual nuclear facility using their security scenarios and timelines. This report does not contain any plant's sensitive information and/or Safeguards Information. The purpose of this study was to verify that results achieved using generic models are similar to actual plant results and to refine our guidance on the use of the framework. This assessment enables further analysis, such as what-if scenarios and staff-reduction evaluation, thereby optimizing physical security at plants.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

An Evaluation of The Dynamic Physical Security Risk Assessment Methodology for Fleet-Wide Applications

The requirements for U.S. nuclear power plants to maintain a large onsite physical security force contribute to their high operational costs. The cost of maintaining the current physical security posture is approximately 10% of the overall operation and maintenance budget for commercial nuclear power plants. The goal of the Light Water Reactor Sustainability (LWRS) program’s physical security pathway is to develop tools, methods, and technologies and provide the technical basis for an optimized physical security posture. The conservatisms built into current security postures may be analyzed and minimized to reduce security costs while still ensuring adequate security and operational safety. The research performed at Idaho National Laboratory within LWRS program’s physical security pathway has successfully developed a dynamic force-on-force modeling framework using various computer simulation tools and integrating them with the dynamic assessment Event Modeling Risk Assessment using Linked Diagrams (EMRALD) tool. This integrated process for physical security analysis is named Modeling and Analysis for Safety Security using Dynamic EMRALD Framework (MASS-DEF). This document provides an update on the progress in applying the MASS-DEF process to an operating commercial nuclear power plant as well as additional industry feedback regarding use of the tool for other physical security risk-informed topics. This report is only a summary of the progress and does not contain specific modeling results as those contain sensitive security information. Previous reports described how a user could integrate their plant-specific force-on-force models with the dynamic simulation tool EMRALD, model operator actions, and integrate with probabilistic risk assessment tools, such as CAFTA (Computer Aided Fault Tree Analysis System) or SAPHIRE (Systems Analysis Programs for Hands-on Integrated Reliability Evaluations), and with thermal-hydraulic tools, such as RELAP-5 or MAAP. Previous reports applied various combinations of available simulations codes with EMRALD using generic plant models to demonstrate how to perform the analysis. This report is an update the progress of applying the dynamic computational framework to an actual nuclear facility using their security scenarios and timelines. This report also provides an update to the procedural guidance for the MASS-DEF process and an overview of the generic models available for use by utilities. This report does not contain any plant’s sensitive information and/or safeguards information. This study’s purpose was to verify that the results achieved using generic models are similar to actual plant results and refine our guidance on the use of the framework. This assessment enables further analysis, such as what-if scenarios and staff-reduction evaluation, thereby optimizing physical security at plants.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

The TDHF code Sky3D version 1.2

The Sky3D code has been widely used to describe nuclear ground states, collective vibrational excitations, and heavy-ion collisions. The approach is based on Skyrme forces or related energy density functionals. The static and dynamic equations are solved on a three-dimensional grid, and pairing is been implemented in the BCS approximation. This updated version of the code aims to facilitate the calculation of nuclear strength functions in the regime of linear response theory, while retaining all existing functionality and use cases. The strength functions are benchmarked against available RPA codes, and the user has the freedom of choice when selecting the nature of external excitation (from monopole to hexadecapole and more). Some utility programs are also provided that calculate the strength function from the time-dependent output of the dynamic calculations of the Sky3D code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗
DOE OSTI2025

A tri-level distribution locational marginal price-based demand response framework

Here, in this paper, we propose a tri-level, nested, two-stage price-based demand response (PBDR) framework that considers distribution locational marginal price (DLMP) as DR enabler between load-serving entities (LSE), demand response providers (DRPs), and customers in the day-ahead distribution market. It enables LSE and customer interactions by using multiple DRPs, positioned in-between, and independently optimizes their objectives. The problem is formulated using linear power flow with approximated power losses and its application in DLMP as DR pricing. The tri-level problem is solved using a nested reformulation & decomposition (R&D) method and tested on the real Indian-108 bus distribution system under various dynamic pricings. Further, the temporal–spatial variations in DLMPs are assessed using fairness criteria. Numerical analyses demonstrate that DLMP applications can effectively improve economic efficiency, and transparency in DR programs valuation with a favorable fairness margin. The results show that DLMP as DR pricing signal induces (0-2) % variation in DLMP for DR participation up to 10 %. Further, it gives over 90 % fairness over temporal–spatial variation for all the customers.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

A development of wide dynamic-range halo monitor for 8 GeV proton beams at FNAL

The FNAL accelerator complex has been upgrading in increasing beam intensity and beam quality. A new beam halo diagnostic device is required in the beam transport line between booster and Recycler. For this purpose, it was decided to introduce the wide dynamic range monitor technique that was developed in 2012 and has been in operation at the J-PARC beam transport line. The device is a two-dimensional beam profile monitor, and it has a dynamic range of approximately six digits of magnitude by using of Optical Transition Radiation and fluorescence screens. Eliminating harmful beam halos is the most important technique for high-intensity proton accelerators. Therefore, beam halo diagnosis is indispensable and becomes more and more important. New FNAL device has been manufactured in a collaboration between J-PARC and FNAL as a part of U.S.-Japan Science and Technology Cooperation Program in High Energy Physics. The equipment will be manufactured at J-PARC and will be shipped to FNAL in 2025.We designed the device to satisfy FNAL specifications: the beam energy, intensity, and size. Currently, most of the equipments are under construction. The large-aperture optical system has been completed and its optical characteristics are being evaluated at J-PARC. We have been also investigating measurement methods corresponding to FNAL bunch trains. This paper reports on the current status of these developments.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Population assignment of migratory Westslope Cutthroat Trout (WCT) in the Clark Fork–Pend Oreille River basin

Abstract Objective The Clark Fork–Pend Oreille River basin of northeastern Washington and the Idaho Panhandle historically supported a robust metapopulation of the Westslope Cutthroat Trout (WCT) Oncorhynchus lewisi, a western native salmonid of high cultural and economic value. The construction of impassible hydroelectric dams and smaller instream barriers has prevented the return of migratory WCT to spawning tributaries, leading to the fragmentation of this metapopulation over the past 100 years. One such impassible barrier is Albeni Falls Dam (AFD) near Newport, Washington, which was completed without fish passage capabilities in 1955. We sought to examine large-scale genetic patterns in the study area and determine the most likely spawning tributary of origin for migratory WCT captured below AFD. Methods We created a genetic baseline representative of populations within the Clark Fork–Pend Oreille River basin from upstream and downstream of the dam using 191 biallelic single-nucleotide polymorphism genetic markers. Our data set included 124 collections, which allowed for an examination of population structure and hatchery influence across the study area and provided a robust tool for population assignment. Population assignment tests were conducted using the program RUBIAS. Result Population assignment tests were successful for all pure WCT of unknown origin despite potential influence from hatchery lineages across the study area. Of 83 migratory WCT captured below AFD, approximately 80% were assigned to tributaries upstream of AFD with a posterior assignment probability of at least 90%. Only one fish was assigned to a tributary downstream of AFD. Conclusion Our results indicate that AFD disrupts the natural metapopulation dynamics of WCT populations in the basin. Passage for WCT at this barrier would reestablish metapopulation connectivity within the basin by allowing migratory individuals to make genetic contributions to populations upstream of the dam.

Wells, Craig D. (ORCID:0000000258513250)↗
DOE OSTI2026

Margin to onset of nucleate boiling and flow instability studies for preliminary MITR design-demonstration element thermal-hydraulics

This study covers the onset of flow instability (OFI) preliminary results obtained from leveraging correlations, in addition to the preliminary thermal hydraulics results such as pressure, flow velocity, temperature, and oxide layer over the design demonstration experiment (DDE) for the Massachusetts Institute of Technology Reactor (MITR). Current computational fluid dynamics (CFD) models in fluid structure interaction (FSI) have added the capability of assessing margins to onset of nucleate boiling (ONB). This study initiates the capability to model the margin to OFI and ONB presented for the MITR. Such study is supportive of the United States High Performance Research Reactor (USHPRR) program. Previous studies provided preliminary thermal-hydraulic and mechanical analyses of the hydrodynamic effects in the MITR DDE under conservative approximations for plate power distribution. In conclusion, this study focuses on providing insights into the OFI future research direction optimizing the transport of thermal energy, mass-flow rates, flow-channel geometries, and boundary conditions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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 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 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 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 OSTI2026

Near-Optimal Performance of Stochastic Model Predictive Control

Here, this article presents a regret analysis for stochastic model predictive control (SMPC) in linear systems with quadratic performance index and additive and multiplicative uncertainties. Under a finite support assumption, the problem can be cast as a finite-dimensional quadratic program, but the problem becomes quickly intractable as the problem size grows exponentially in the horizon length. SMPC aims to compute approximate solutions by solving a sequence of problems with truncated prediction horizons and committing the solution in a receding-horizon fashion. Although this approach is widely used in practice, its performance relative to the optimal solution is not well understood. This article reports for the first time a rigorous near-optimal performance guarantee of SMPC: under stabilizability and detectability conditions, the regret of SMPC is exponentially small in the prediction horizon length, allowing SMPC to achieve near-optimal performance at a substantially reduced computational expense.

93E20, 93B45↗
DOE OSTI2023

Atomistic Simulations of Polydisperse Lignin Melts Using Simple Polydisperse Residue Input Generator

Understanding the physics of lignin will help rationalize its function in plant cell walls as well as aiding practical applications such as deriving biofuels and bioproducts. Here, in this work, we present SPRIG (Simple Polydisperse Residue Input Generator), a program for generating atomic-detail models of random polydisperse lignin copolymer melts i.e., the state most commonly found in nature. Using these models, we use all-atom molecular dynamics (MD) simulations to investigate the conformational and dynamic properties of polydisperse melts representative of switchgrass (Panicum virgatum L.) lignin. Polydispersity, branching and monolignol sequence are found to not affect the calculated glass transition temperature, T g . The Flory–Huggins scaling parameter for the segmental radius of gyration is 0.42 ± 0.02, indicating that the chains exhibit statistics that lie between a globular chain and an ideal Gaussian chain. Below T g the atomic mean squared displacements are independent of molecular weight. In contrast, above T g , they decrease with increasing molecular weight. Therefore, a monodisperse lignin melt is a good approximation to this polydisperse lignin when only static properties are probed, whereas the molecular weight distribution needs to be considered while analyzing lignin dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗