Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “optimal start”

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 361 records · Page 20

Iterative subspace algorithms for finite-temperature solution of Dyson equation

One-particle Green’s functions obtained from the self-consistent solution of the Dyson equation can be employed in the evaluation of spectroscopic and thermodynamic properties for both molecules and solids. However, typical acceleration techniques used in the traditional quantum chemistry self-consistent algorithms cannot be easily deployed for the Green’s function methods because of a non-convex grand potential functional and a non-idempotent density matrix. Moreover, the optimization problem can become more challenging due to the inclusion of correlation effects, changing chemical potential, and fluctuations of the number of particles. In this paper, we study acceleration techniques to target the self-consistent solution of the Dyson equation directly. We use the direct inversion in the iterative subspace (DIIS), the least-squared commutator in the iterative subspace (LCIIS), and the Krylov space accelerated inexact Newton method (KAIN). We observe that the definition of the residual has a significant impact on the convergence of the iterative procedure. Based on the Dyson equation, we generalize the concept of the commutator residual used in DIIS and LCIIS and compare it with the difference residual used in DIIS and KAIN. The commutator residuals outperform the difference residuals for all considered molecular and solid systems within both GW and GF2. For a number of bond-breaking problems, we found that an easily obtained high-temperature solution with effectively suppressed correlations is a very effective starting point for reaching convergence of the problematic low-temperature solutions through a sequential reduction of temperature during calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Supply chain optimization of sustainable aviation fuel from carinata in the Southeastern United States

We report carinata is a purpose-grown oilseed feedstock for renewable fuels, including sustainable aviation fuel (SAF) that can replace conventional aviation fuel (CAF). Given carinata is a new crop in the Southeastern United States, it is crucial to analyze its sustainability from a supply chain perspective. This study developed a mixed-integer linear programming (MILP) model and simulated it for 20 years, starting from a farm (county-level data) and ending at the airport. About 2.06 million ha in Alabama, Florida, and Georgia combined were found suitable for carinata production. Given the three-year rotation period, about 0.69 million hectares can be cultivated annually, approximately 65% of which was in Georgia. About 2.4% of the combined SAF annual demand of four major airports (about 210 thousand t) in the study area is satisfied at that level of carinata cultivation. However, all available SAF was supplied to the Atlanta airport as this decision minimizes the supply chain cost. A total of 1343 storage units, one oil extraction mill, and one biorefinery were needed to meet the overall demand. We found that SW Georgia is the top supplier of carinata seeds. The unit cost of production and carbon intensity were estimated to be $\$$0.89 L -1 (or $\$$26.79 GJ -1 ) and 0.91 kg CO 2 e L -1 (or 27.28 kg CO 2 e GJ -1 ), respectively. This carbon intensity of carinata-based SAF was 67.8% lower than that of CAF. With variations included in SAF demand, yield, and soil carbon sequestration, carbon savings remained between 66.5% and 67.8%. Given the GHG advantage of SAF over CAF, there is justification for subsidies required to make SAF competitive.

10 SYNTHETIC FUELS↗

Identification of Small-Molecule Inhibitors of Fibroblast Growth Factor 23 Signaling via In Silico Hot Spot Prediction and Molecular Docking to α-Klotho

Fibroblast growth factor 23 (FGF23) is a therapeutic target for treating hereditary and acquired hypophosphatemic disorders, such as X-linked hypophosphatemic (XLH) rickets and tumor-induced osteomalacia (TIO), respectively. FGF23-induced hypophosphatemia is mediated by signaling through a ternary complex formed by FGF23, the FGF receptor (FGFR), and α-Klotho. Currently, disorders of excess FGF23 are treated with an FGF23-blocking antibody, burosumab. Small-molecule drugs that disrupt protein/protein interactions necessary for the ternary complex formation offer an alternative to disrupting FGF23 signaling. In this study, the FGF23:α-Klotho interface was targeted to identify small-molecule protein/protein interaction inhibitors since it was computationally predicted to have a large fraction of hot spots and two druggable residues on α-Klotho. We further identified Tyr433 on the KL1 domain of α-Klotho as a promising hot spot and α-Klotho as an appropriate drug-binding target at this interface. Subsequently, we performed in silico docking of ~5.5 million compounds from the ZINC database to the interface region of α-Klotho from the ternary crystal structure. Following docking, 24 and 20 compounds were in the final list based on the lowest binding free energies to α-Klotho and the largest number of contacts with Tyr433, respectively. Herein, five compounds were assessed experimentally by their FGF23-mediated extracellular signal-regulated kinase (ERK) activities in vitro, and two of these reduced activities significantly. Both these compounds were predicted to have favorable binding affinities to α-Klotho but not have a large number of contacts with the hot spot Tyr433. ZINC12409120 was found experimentally to disrupt FGF23:α-Klotho interaction to reduce FGF23-mediated ERK activities by 70% and have a half maximal inhibitory concentration (IC50) of 5.0 ± 0.23 μM. Molecular dynamics (MD) simulations of the ZINC12409120:α-Klotho complex starting from in silico docking poses reveal that the ligand exhibits contacts with residues on the KL1 domain, the KL1–KL2 linker, and the KL2 domain of α-Klotho simultaneously, thereby possibly disrupting the regular function of α-Klotho and impeding FGF23:α-Klotho interaction. ZINC12409120 is a candidate for lead optimization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Effect of Geometry, Spin, and Orbital Optimization in Achieving Accurate, Correlated Results for Iron–Sulfur Cubanes

Iron–sulfur clusters comprise an important functional motif in the catalytic centers of biological systems, capable of enabling important chemical transformations at ambient conditions. This remarkable capability derives from a notoriously complex electronic structure that is characterized by a high density of states that is sensitive to geometric changes. The spectral sensitivity to subtle geometric changes has received little attention from correlated, large active space calculations, owing partly to the exceptional computational complexity for treating these large and correlated systems accurately. Here, to provide insight into this aspect, we report the first Complete Active Space Self Consistent Field (CASSCF) calculations for different geometries of the [Fe(II/III) 4 S 4 (SMe) 4 ] -2 clusters using two complementary, correlated solvers: spin-pure Adaptive Sampling Configuration Interaction (ASCI) and Density Matrix Renormalization Group (DMRG). We find that the previously established picture of a double-exchange driven magnetic structure, with minute energy gaps (<1 mHa) between consecutive spin states, has a weak dependence on the underlying geometry. However, the spin gap between the singlet and the spin state 2S + 1 = 19, corresponding to a maximal number of Fe-d electrons being unpaired and of parallel spin, is strongly geometry dependent, changing by a factor of 3 upon slight deformations that are still within biologically relevant parameters. The CASSCF orbital optimization procedure, using active spaces as large as 86 electrons in 52 orbitals, was found to reduce this gap compared to typical mean-field orbital approaches. Our results show the need for performing large active space calculations to unveil the challenging electronic structure of these complex catalytic centers and should serve as accurate starting points for fully correlated treatments upon inclusion of dynamical correlation outside the active space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Particle-based fast jet simulation at the LHC with variational autoencoders

Abstract We study how to use deep variational autoencoders (VAEs) for a fast simulation of jets of particles at the Large Hadron Collider. We represent jets as a list of constituents, characterized by their momenta. Starting from a simulation of the jet before detector effects, we train a deep VAE to return the corresponding list of constituents after detection. Doing so, we bypass both the time-consuming detector simulation and the collision reconstruction steps of a traditional processing chain, speeding up significantly the events generation workflow. Through model optimization and hyperparameter tuning, we achieve state-of-the-art precision on the jet four-momentum, while providing an accurate description of the constituents momenta, and an inference time comparable to that of a rule-based fast simulation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Plateau Phenomenon in Gradient Descent Training of RELU Networks: Explanation, Quantification, and Avoidance

The ability of neural networks to provide ‘best in class’ approximation across a wide range of applications is well-documented. Nevertheless, the powerful expressivity of neural networks comes to naught if one is unable to effectively train (choose) the parameters defining the network. In general, neural networks are trained by gradient descent type optimization methods,a stochastic variant thereof. In practice, such methods result in the loss function decreases rapidly at the beginning of training but then, after a relatively small number of steps, significantly slow down. The loss may even appear to stagnate over the period of a large number of epochs, only to then suddenly start to decrease fast again for no apparent reason. This so-called plateau phenomenon manifests itself in many learning tasks. The present work aims to identify and quantify the root causes of plateau phenomenon.analysis is carried out in the setting of univariate ReLU networks. No assumptions are made on the number of neurons relative to the number of training data, and our results hold for both the lazy and adaptive regimes. Here, the main findings are: plateaux correspond to periods during which activation patterns remain constant, where activation pattern refers to the number of data points that activate a given neuron; quantification of convergence of the gradient flow dynamics; and, characterization stationary points in terms solutions of local least squares regression lines over subsets of the training data. Based on these conclusions, we propose a new iterative training method, the Active Neuron Least Squares (ANLS), characterised by the explicit adjustment of the activation pattern at each step, which is designed to enable a quick exit from a plateau. Illustrative numerical examples are included throughout.

97 MATHEMATICS AND COMPUTING↗

Uncovering the mechanism of selective stabilization of high-energy diastereoisomers via inclusion

Abstract Supramolecular systems may be used to stabilize otherwise unstable isomers to find alternative synthetic pathways. It has been reported that cucurbit[8]uril can stabilize trans -I and trans -II Cu $$^{{\textrm{II}}}$$ II cyclam, whereas trans -III is the only non-substituted trans Cu $$^{{\textrm{II}}}$$ II cyclam diastereoisomer found outside of the host molecule experimentally. Quantum chemistry methods can provide valuable insight into the intermolecular interactions involved in these inclusion complexes. All five possible trans diastereoisomers of Cu $$^{{\textrm{II}}}$$ II cyclam were studied within the host molecule to calculate the interaction energy and free energy of association for each complex. The relative free energies of the five free cyclams confirm that trans -I and trans -II are the most energetically accessible diastereoisomers from the initial trans -III starting point. Energy decomposition analysis was used to identify the attractive and repulsive interactions between cyclam and cucurbit[8]uril and showed that trans -II encounters repulsive forces almost three times greater than trans -I, which may explain the 7:3 ratio of trans -I to trans -II within cucurbit[8]uril that occurs experimentally. Optimized complex geometries with trans -III, IV, and V show that the cyclams protrude out of cucurbit[8]uril, whereas trans -I and trans -II become more encapsulated and elongate the host, suggesting that the position of the cyclam is extremely important when forming non-covalent interactions. Our results agree with the experimental findings and provide greater insight into why the most stable isolated cyclam diastereoisomer, trans -III, does not form a complex.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nucleon-pair truncation of the shell model for medium-heavy nuclei

Background: Establishing computationally tractable models of atomic nuclei is a long-time goal of nuclear structure physics. A flexible framework which easily includes excited states and many-body correlations is the configuration-interaction shell model (SM), but the exponential growth of the basis means one needs an efficient truncation scheme, ideally one that includes both deformation and pairing correlations. Purpose: Here, we propose an efficient truncation scheme of the SM: starting from a pair condensate variationally defined by Hartree-Fock single-particle states and the particle-number conserved Bardeen-Cooper-Schrieffer (NBCS) approximation, we carry out projection of states with good angular momentum. Methods: After generating Hartree-Fock single-particle states with Kramers degeneracy in a SM space, we optimize the pair amplitudes in the NBCS by minimizing the energy, and then use linear algebra projection (LAP) of states with good angular momentum. Both NBCS and LAP are computationally fast. Results: Our calculations yield good agreement with full configuration-interaction SM calculations for low-lying states of transitional and rotational nuclei with axially symmetric and triaxial deformation in medium- and heavy-mass regions: 44,46,48 Ti, 48,50 Cr, 52 Fe, 60,62,64 Zn, 66,68 Ge, 68 Se, and 108,110 Xe. We predict low-lying states of 112–114 Ba and 116–120 Ce, nuclei difficult to reach by large-scale SM calculations. Conclusions: Both pair correlation and the configuration mixing between different intrinsic states play a key role in reproducing collectivity and shape coexistence, demonstrating the utility of this truncation scheme of the SM to study transitional and deformed nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

UNRAVEL: A Robust Modular Velocity Dealiasing Technique for Doppler Radar

Unfold Radar Velocity (UNRAVEL) is an open-source modular Doppler velocity dealiasing algorithm for weather radars. UNRAVEL is an algorithm that does not need external reference velocity data, making it easily applicable. The proposed algorithm includes 11 core modules and 2 dealiasing strategies. UNRAVEL is an iterative algorithm. The goal is to build the dealiasing results starting with the strictest possible continuity tests in azimuth and range and, after each step, relaxing the parameters to include more results from a progressively growing number of reference points. UNRAVEL also has modules that perform 3D continuity checks. Thanks to this modular design, the number of dealiasing strategies can be expanded in order to optimize the dealiasing results. While the first driver dealiases Doppler velocity from each tilt independently from one another, the second driver also performs a three-dimensional continuity check of the velocity using successive elevations. The proposed dealiasing algorithm is tested using severe weather data from an S-band Doppler radar that have been aliased to mimic aliased radial velocity patterns that would be observed by a C-band Doppler radar. Artificially aliasing S-band data permits creation of a reference to which the performance of various dealiasing techniques can be compared. Comparisons show that UNRAVEL consistently outperforms other established dealiasing algorithms for the test period selected in this work.

47 OTHER INSTRUMENTATION↗

Control of the Plasma-Material Interface for Long Pulse Optimization in EAST and KSTAR

This is the Final Technical Report for the Grant DE-SC0016322. The major goals of the project are to use HIDRA to evaluate and optimize the FLiLi and LiMIT concepts, in preparation for the final designs of liquid Li PFCs for testing in EAST long pulse discharges. In the first year the goal is to fully commission HIDRA and begin liquid Li tests. In the second year it is to test the FLiLi and LiMIT designs in HIDRA and start the production of the first LiMIT test plate for EAST.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Building a Computational and Experimental Rapid Response Pipeline to Counter the Coronavirus Disease 2019 Outbreak and Emerging Biothreats

The COVID-19 pandemic underscored the promise of monoclonal antibody-based prophylactic and therapeutic drugs, especially where protective candidates can be rapidly identified or developed for emerging biothreats and escape variants. Current cutting-edge technologies for this purpose still rely on pathogen-exposed convalescent volunteers and a large screening effort to find a proverbial needle in a haystack. Computational design of protective antibodies based on pre-existing templates skips those requirements and allows for greater control over the breadth and target epitope, while also co-optimizing for potency and developability or other biophysical characteristics. We approached this problem by building and expanding an in vitro experimental rapid antibody production and characterization pipeline to support development of an autonomous, closed loop, active learning software system based on structural simulation and ground truth experimental data to design and evaluate antibody antigen interactions. Starting from early in the pandemic, we targeted SARS-CoV-2, the causative agent of COVID-19, by re-purposing neutralizing antibodies against SARS-CoV-1 that had been identified in the wake of that outbreak in the early 2000’s. We successfully re-targeted three different anti-SARS-CoV-1 antibodies to neutralize SARS-CoV-2 in vitro, where the antibodies were generated externally and tested through conventional binding and neutralization assays internally or with collaborators. As antibodies were identified from the blood of humans infected with SARS-CoV-2, we shifted to collaborate with academic partners to develop improved versions of their human-derived antibodies. This work reached its most important stage in rapid response to the emergence of the Omicron variant of concern (VOC) in late 2021. In a matter of weeks, enabled by on demand innovation to our screening pipeline, we computationally designed and experimentally characterized derivative antibodies of COV2-2130, one of two antibodies from Vanderbilt that form the basis of the AstraZeneca Evusheld prophylactic drug product. This drug product suffers a serious loss of efficacy against Omicron BA.1 and BA.1.1, the first Omicron strains. Due to tight integration of computational design and experimental evaluation, we were able to identify a pair of designs with potent neutralization of the main targets Omicron BA.1 and BA.1.1; but also the earlier Delta variant, and subsequent Omicron strains including BA.2, BA.4, BA.5, and BA.2.75, demonstrating that our multi-target design process can, by its nature, produce robust antibody designs that strictly improve over the parental antibody. These results, recognized by a 2022 Director’s Science and Technology award, have enabled the follow-on GUIDE program, to commence in FY23. While earlier design campaigns were substantially outsourced, we have engineered better and faster processes internally to better compliment, calibrate, and speed computational designs. As part of the follow-on GUIDE program, we will stand up a rapid and high-throughput antibody production and characterization facility staffed with the expertise and capabilities to foster our current collaboration across PLS and ENG as well as other partnerships toward computational design of biologics.

59 BASIC BIOLOGICAL SCIENCES↗

A Two-Level Model Predictive Control-Based Approach for Building Energy Management including Photovoltaics, Energy Storage, Solar Forecasting and Building Loads

This paper uses a two-level model predictive control-based approach for the coordinated control and energy management of an integrated system that includes photovoltaic (PV) generation, energy storage, and building loads. Novel features of the proposed local controller include (1) the ability to simultaneously manage building loads and energy storage to achieve different operational objectives such as energy efficiency, economic cost efficiency, demand response and grid optimization through the design of specific power trajectory tracking performance functionals, (2) an energy trim function that minimizes the impact of solar forecasting errors on system performance, and (3) the design of a state of charge controller that uses day-ahead forecast of solar power and building loads to intialize energy storage at the start of each day. The local controller is tested in simulation using an exemplary system with PV generation, energy storage and dispatchable building loads. Two sample days with different PV forecasts and multiple case scenarios are considered, and the performance of the algorithm in managing the real and reactive net building load trajectories and the ramp rate of PV injections into the utility network are evaluated. The simulations are based on actual forecasted and measured PV data, and the results show that the local controller meets the tracking requirements for real and reactive power within the operating constraints of the building.

14 SOLAR ENERGY↗

STREAM: A technology planning and capacity expansion model for the industrial sector

The Strategic Technology Roadmapping and Energy, Environmental, and Economic Analysis Model—STREAM—is an optimization-based modeling tool and analysis framework to assist with strategic planning and technology investments of the industrial sector. This open-source framework is written in Julia using the JuMP package, which enables users to model future “pathways” for incumbent and future production technologies, costs, fuels and energy carriers, and energy and non-energy environmental impacts from industries as they transform in pursuit of a robust and competitive manufacturing sector. The model starts with an initial stock of industrial production technologies and assets at a facility level and then determines pathways that minimize cost, subject to an array of possible constraints on demand, market shares, environmental flows, and other exogenously specified operational considerations such as capacity utilization rates or regional energy costs. Key features of the framework include flexibility to model a wide range of industries and industrial technologies/processes at varying levels of granularity, ability to perform parametric sensitivity analyses, and ability to visualize model results using visualization objects.

capacity expansion↗

csphothybv3: A hybrid sky radiance scan scenario version 3 data

This hybrid scan starts at the position of the Sun and first follows the great circle that is orthogonal to the Principal Plane until it reaches a viewing zenith angle of 75 degrees and then follows the almucantar keep that viewing zenith angle constant. The scan is then repeated in the opposite azimuth direction. This mode was designed to optimize the range of scattering angles and significantly improves the retrieval of inversion parameters especially at low solar zenith angles.

54 ENVIRONMENTAL SCIENCES↗

Global kinetic modeling of NH 3 -SCR with two sites of NH 3 storage on Cu-SSZ-13

A comprehensive selective catalytic reduction (SCR) model is developed to detail the global kinetic reactions in the catalyst which is involved in two type of sites for high- and low-temperature NH 3 adsorption/desorption, as well as a low temperature H 2 O storage. The model accounts for the formation of ammonium nitrate in fast SCR and the decomposition formation of ammonium nitrate in the absence of NO, and addresses N 2 O relative to NO and NO 2 reactions with NH 3 . The model has been validated against commercial Cu-SSZ-13 measurements of NH 3 inventory, NH 3 oxidation, NO oxidation, standard SCR, fast SCR and NO 2 SCR at both comprehensive steady-state and transient conditions. In addition, the application of the SCR model for emissions control over a transient drive cycle at cold start was demonstrated based on diesel oxidation catalyst (DOC) and SCR aftertreatment system for a passenger car. In conclusion, the results show that the SCR model can provide the detailed emissions estimations for steady-state and transient conditions, can assist with understanding the impact of the SCR reaction mechanism, and can assist with the design and optimization of SCR catalysts to develop innovative technologies for co-optimization of engine fuel economy and emissions control technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Engineering intricacies of implementing single-atom alloy catalysts for low-temperature electrocatalytic CO 2 reduction

Catalysts research for electrocatalytic CO 2 reduction reactions (CO 2 R) has undergone rapid growth in the last decade. Single-atom alloy catalysts (SAAs) featuring atomically dispersed metal dopants on host metal surfaces have shown promises in boosting CO 2 R yield by optimizing the structure and synergy of the catalytic metals at the atomic scale. Despite the exciting development of SAAs for CO 2 R in fundamental science, dedicated studies for its engineering implementation have been absent. We use this perspective to discuss our non-exhaustive engineering considerations for implementing SAAs for CO 2 R. Here, the perspective starts with a brief overview of the current research status for SAAs in CO 2 R, followed by focal points on structure uncertainties associated with catalyst manufacturing, catalyst layer degradation during reaction, and possibilities for SAAs to mitigate the salt precipitation issue at the device level. We hope our opinions will engage increasing attention toward the engineering catalysis research for applying SAAs to CO 2 R at scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Robustness of energy performance of Zero-Net-Energy (ZNE) homes

Zero-net-energy (ZNE) homes produce an adequate amount of energy on-site to meet their energy demand based on source energy for an entire year. California building energy efficiency standards require new residential buildings started in 2020 to be ZNE. For various reasons, a home designed as ZNE may not achieve ZNE performance in real operation. Here we aimed to quantify the robustness of the energy performance of ZNE homes due to weather variability, climate change, and the uncertainty of occupant behavior. A single-family ZNE house, based on the optimal cost-effective design in three California climate zones, was used to develop the EnergyPlus simulation models. Weather variations were considered from a combination of the historical 30 years’ actual meteorological year (AMY) weather data, typical year weather data in TMY3, and future weather data based on Intergovernmental Panel on Climate Change scenarios. Three scenarios of occupant behavior from the energy perspective were defined to represent the uncertainty about occupants’ activities, comfort requirements, and their adaptive interactions with buildings and systems. In terms of annual source energy, the simulation results of the ZNE homes showed: (1) a decrease of 23–38 percent for occupants with energy austerity behavior and an increase of 120–130 percent for occupants with energy wasteful behavior, compared with the baseline assumption of normal occupants; (2) a variation range of –15 percent to +14 percent for the results using 30-year AMY weather data compared with the baseline results using TMY3 weather data; (3) an increase of 10–13 percent with future weather in Fresno and Riverside and a decrease of 15 percent with San Francisco; and (4) climate change can reduce the gap between the austerity and wasteful consumption. These findings provide insights into how ZNE homes may perform in reality and inform architects, engineers, occupants, and policymakers to pay more attention to occupant behavior on design, operation, and regulations of ZNE homes to ensure energy performance robustness.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

EIC Detector Overview

The Electron Ion Collider will have two interaction regions that can be instrumented with detectors. Here, the first region will be instrumented as part of the project and needs to be capable of delivering the physics that has been outlined by the National Academy of Sciences and ready at the start of beam commissioning near the end of this decade. Plans for a second complementary detector to be located at a second interaction region are already in progress and will hopefully come to fruition just few years after the first detector comes online. While the basic parameters of these detectors are being selected using conventional approaches, the optimization of the detectors is already being enhanced by making use of advanced optimization techniques.

47 OTHER INSTRUMENTATION↗