SEARCH · Engineering Papers
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.
We present our SciDAC FASTMath-HEP partnership results for tuning generative adversarial models (GANs) for high energy physics applications. The GANs are used in hybrid simulations to accelerate otherwise time-consuming computations. We optimize for both, prediction accuracy and variability with the goal to find GAN architectures that are reliable and robust.
Searches for new phenomena in complex scientific data are predominantly model-dependent, optimized for specific hypotheses, and therefore limited in their coverage of the space of possible signals. Recently, new AI-based model-agnostic search strategies, many of which have been pioneered in high-energy physics, have been proposed which provide a complementary paradigm, prioritizing broad exploration over tailored analyses. These techniques offer an opportunity to enhance the overall discovery potential of modern experiments, especially in regimes where theoretical guidance is scarce. In this document, we review the conceptual framework behind the main classes of AI-based model-agnostic strategies. We discuss the potential pitfalls of these methods, and strategies for their validation and interpretation. We aim for this document to serve as a useful reference both for practitioners and for researchers interested in learning more about these model-agnostic search strategies.
DOE and NETL have created an advanced, open-source computational platform through the Institute for the Design of Advanced Energy Systems (IDAES). The IDAES platform is a very extensive modeling environment that can be used for a broad range of power plant and process applications. MapEx Software is developing and commercializing a software application that makes it easier to set-up and run IDAES-based analyses. The MapEx-developed software application will replace the need for custom Python language coding with a user-friendly graphical user interface where the user can construct a flowsheet diagram of the IDAES model by inserting icons representing the plant equipment onto the screen. This application will make the implementation of modeling and optimization of existing fossil-fired power plants more straight-forward and less time-consuming. The effort focuses on performance monitoring and optimization of plant operations for the existing coal-fired power plant fleet but is built on a structure that allows expansion into the broad range of applications where IDAES methods may be applied.
A fiber-optic visible–near-infrared absorption spectroscopy system in a glove box was demonstrated for remote quantification of Am(III) (0–500 µM) and HNO 3 (0.1–9 M) using partial least squares regression (PLSR) models. The sensor platform, featuring a simple plug-and-play spectrophotometer, can enable noninvasive, real-time monitoring of actinide process solutions. To establish a flexible PLSR model calibration strategy, a D-optimal design developed using Nd(III) in previous studies was successfully extended to an actinide system with Am(III) to effectively minimize sample set size while maintaining robust prediction performance. The results suggest strong spectral similarities between Nd(III) and Am(III) and validate Nd(III) as an effective optical surrogate for trivalent actinide species. This work also supports the generalizability of a D-optimal training set selection approach for two-factor systems. The PLS1 models for Am(III) and HNO 3 outperformed a PLS2 model and maintained reasonable performance in the presence of interfering U(VI). The resulting sensor system and multivariate approach provides a flexible and scalable solution for process monitoring, control, and safety in diverse nuclear applications.
The following article details a model predictive control (MPC) to improve grid resilience when faced with variable generation resources. This topic is of significant interest to utility power systems where distributed intermittent energy sources will increase significantly and be relied on for electric grid ancillary services. Previous work on MPCs has focused on narrowly targeted control applications such as improving electric vehicle (EV) charging infrastructure or reducing the cost of integrating Energy Storage Systems (ESSs) into the grid. In contrast, this article develops a comprehensive treatment of the construction of an MPC tailored to electric grids and then applies it integration of intermittent energy resources. To accomplish this, the following article includes a description of a reduced order model (ROM) of an electric power grid based on a circuit model, an optimization formulation that describes the MPC, a collocation method for solving linear time-dependent differential algebraic equations (DAEs) that result from the ROM, and an overall strategy for iteratively refining the behavior of the MPC. Next, the algorithm is validated using two separate numerical experiments. First, the algorithm is compared to an existing MPC code and the results are verified by a numerically precise simulation. It is shown that this algorithm produces a control comparable to existing algorithms and the behavior of the control carefully respects the bounds specified. Second, the MPC is applied to a small nine bus system that contains a mix of turbine-spinning-machine-based and intermittent generation in order to demonstrate the algorithm’s utility for resource planning and control of intermittent resources. This study demonstrates how the MPC can be tuned to change the behavior of the control, which can then assist with the integration of intermittent resources into the grid. The emphasis throughout the paper is to provide systematic treatment of the topic and produce a novel nonlinear control compatible design framework applicable to electric grids and the control of variable resources. This differs from the more targeted application-based focus in most presentations.
Efficient analytical and computational tools for simultaneous optimal design of the structural and control components of aeroservoelastic systems are presented. The optimization objective is to achieve aircraft performance requirements and sufficient flutter and control stability margins with a minimal weight penalty and without violating the design constraints. Analytical sensitivity derivatives facilitate an efficient optimization process which allows a relatively large number of design variables. Standard finite element and unsteady aerodynamic routines are used to construct a modal data base. Minimum State aerodynamic approximations and dynamic residualization methods are used to construct a high accuracy, low order aeroservoelastic model. Sensitivity derivatives of flutter dynamic pressure, control stability margins and control effectiveness with respect to structural and control design variables are presented. The performance requirements are utilized by equality constraints which affect the sensitivity derivatives. A gradient-based optimization algorithm is used to minimize an overall cost function. A realistic numerical example of a composite wing with four controls is used to demonstrate the modeling technique, the optimization process, and their accuracy and efficiency.
Explore the source record for details and available documents.
ABSTRACT Introduction Extensive trauma, commonly seen in wounded military Service Members, often leads to a severe sterile inflammation termed systemic inflammatory response syndrome (SIRS), which can progress to multiple organ dysfunction syndrome (MODS) and death. MODS is a serious threat to wounded Service Members, historically causing 10% of all deaths in trauma admissions at a forward deployed combat hospital. The importance of this problem will be exacerbated in large-scale combat operations, in which evacuation will be delayed and care of complex injuries at lower echelons of care may be prolonged. The main goal of this study was to optimize an existing mouse model of lethal SIRS/MODS as a therapeutic screening platform for the evaluation of immunomodulatory drugs. Materials and Methods Male C57BL/6 mice were euthanized, and the bones and muscles were collected and blended into a paste termed tissue–bone matrix (TBX). The TBX at 12.5%–20% relative to body weight of each recipient mouse was implanted into subcutaneous pouches created on the dorsum of anesthetized animals. Mice were observed for clinical scores for up to 48 hours postimplantation and euthanized at the preset point of moribundity. To test effects of anesthetics on TBX-induced mortality, animals received isoflurane or ketamine/xylazine (K/X). In a separate set of studies, mice received TBX followed by intraperitoneal injection with 20 mg/kg or 40 mg/kg Eritoran or a placebo carrier. All Eritoran studies were performed in a blinded fashion. Results We observed that K/X anesthesia significantly increased the lethality of the implanted TBX in comparison to inhaled anesthetics. Although all the mice anesthetized with isoflurane and implanted with 12.5% TBX survived for 24 hours, 60% of mice anesthetized with K/X were moribund by 24 hours postimplantation. To mimic more closely the timing of lethal SIRS/MODS following polytrauma in human patients, we extended observation to 48 hours. We performed TBX dose–response studies and found that as low as 15%, 17.5%, and 20% TBX caused moribundity/mortality in 50%, 80%, and 100% mice, respectively, over a 48-hour time period. With 17.5% TBX, we tested if moribundity/mortality could be rescued by anti-inflammatory drug Eritoran, a toll-like receptor 4 antagonist. Neither 20 mg/kg nor 40 mg/kg doses of Eritoran were found to be effective in this model. Conclusions We optimized a TBX mouse model of SIRS/MODS for the purpose of evaluating novel therapeutic interventions to prevent trauma-related pathophysiologies in wounded Service Members. Negative effects of K/X on lethality of TBX should be further evaluated, particularly in the light of widespread use of ketamine in treatment of pain. By mimicking muscle crush, bone fracture, and necrosis, the TBX model has pleiotropic effects on physiology and immunology that make it uniquely valuable as a screening tool for the evaluation of novel therapeutics against trauma-induced SIRS/MODS.
Motion primitives allow for application of discrete search algorithms to rapidly produce trajectories in complex continuous space. The maneuver automaton (MA) provides an elegant formulation for creating a primitive library based on trims and maneuvers. However, performance is fundamentally limited by the contents of the primitive library. If the library is too sparse, performance can be poor in terms of path cost, whereas a library that is too large can increase run time. This work outlines new methods for using genetic algorithms to prune a primitive library. The proposed methods balance the path cost and planning time while maintaining the reachability of the MA. The genetic algorithm in this paper evaluates and mutates populations of motion primitive libraries to optimize both objectives. Here, we illustrate the performance of these methods with a simulated study using a nonlinear medium-fidelity F-16 model. We optimize a library with the presented algorithm for obstacle-free navigation and a nap-of-the-Earth navigation task. In the obstacle-free navigation task, we show a tradeoff of a 10.16% higher planning cost for a 96.63% improvement in run time. In the nap-of-the-Earth task, we show a tradeoff of a 9.712% higher planning cost for a 92.06% improvement in run time.
Although the random forest (RF) model is a powerful machine learning tool that has been utilized in many wind speed/power forecasting studies, there has been no consensus on optimal RF modeling strategies. This study investigates three basic questions which aim to assist in the discernment and quantification of the effects of individual model properties, namely: (1) using a standalone RF model versus using RF as a correction mechanism for the persistence approach, (2) utilizing a recursive versus direct multi-step forecasting strategy, and (3) training data availability on model forecasting accuracy from one to six hours ahead. These questions are investigated utilizing data from the FINO1 offshore platform and Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) C1 site, and testing results are compared to the persistence method. At FINO1, due to the presence of multiple wind farms and high inter-annual variability, RF is more effective as an error-correction mechanism for the persistence approach. The direct forecasting strategy is seen to slightly outperform the recursive strategy, specifically for forecasts three or more steps ahead. Finally, increased data availability (up to ~8 equivalent years of hourly training data) appears to continually improve forecasting accuracy, although changing environmental flow patterns have the potential to negate such improvement. We hope that the findings of this study will assist future researchers and industry professionals to construct accurate, reliable RF models for wind speed forecasting.
Metallic fuels have seen increased interest for future sodium fast reactors due to their material properties: high thermal conductivities and advantageous neutronic properties allow for greater fission densities. One drawback to typical metallic fuels is zirconium redistribution, which impacts this advantageous material and its neutronic properties. Unfortunately, the processes behind zirconium migration behavior are understood using first principles, so before these fuels are implemented in future fast reactors, characterization and fuel qualification regimes must be completed. These activities can be supported through the use of robust modeling using the most accurate empirical models currently available to fuel researchers around the world. The tool that allows researchers to model this complex coupled thermo-mechanical behavior and nuclear properties is BISON. Additionally, BISON model parameters need to be compared against PIE measurements. The current work utilizes two fuel pins from EBR-II experiment X441 to optimize various model parameters, including porosity correction factor, thermal conductivity, phase transition temperature, and diffusion coefficient multipliers, before implementing the final model for seven fuel pins with differing characteristics. To properly evaluate the BISON simulations, the results are compared to PIE metallography data for each fuel pin, to ensure the zirconium redistribution is properly reflected in the simulation results. Six out of seven analyzed fuel pins demonstrate good agreement between the metallography images and BISON results, showing alignment of the Zr-rich, Zr-depleted, and moderately Zr-enriched zones at various axial heights along the fuel pins. Further work is needed to refine the model parameters for general pin use.
One of the main charges of beam physics is to improve the overall beam quality delivered to experiments and applications, which is generally quantified by the beam brightness. High brightness beams have numerous applications, including x-ray free electron lasers and ultrafast electron diffraction (UED). To these ends, this thesis details efforts to use high-performance models --- high-fidelity models that execute quickly --- for controls and diagnostics. To demonstrate the generality of these techniques, this thesis focuses on three of the most successful photoinjector designs currently used worldwide: the UCLA/SLAC/BNL-type high-gradient S-band gun, the continuous-wave high-repetition rate VHF APEX gun and the L-band DESY-PITZ-type gun. Work will be shown from Pegasus (UCLA), HiRES (LBNL) and FAST (FNAL).Specifically, data-driven models for online virtual diagnostics are presented, in this case, in the context of UED at HiRES, leading to a temporal resolution improvement. Methods for improving the fidelity of physics-based models are discussed, with examples at HiRES, Pegasus and FAST. Markov-chain Monte Carlo analysis is applied to match simulations, in the context of photocathode studies at HiRES and Pegasus. Lastly, the augmentation of online model-based predictions with model-independent optimization is explored in a fluctuating environment at Pegasus and HiRES. A central theme of this dissertation is working with parameter fluctuations when modeling an accelerator beamline. Long-term drifts and shot-to-shot jitter exist in every accelerator to a varying degree and therefore play an important role in every chapter of this dissertation. This thesis attempts to address the issues associated with these fluctuations when trying to develop faithful model representations of the system.
The problem of satellite synthesis can be described as optimally allotting locations and sometimes frequencies and polarizations, to communication satellites so that interference from unwanted satellite signals does not exceed a specified threshold. In this report, mathematical programming models and optimization methods are used to solve satellite synthesis problems. A nonlinear programming formulation which is solved using Zoutendijk's method and a gradient search method is described. Nine mixed integer programming models are considered. Results of computer runs with these nine models and five geographically compatible scenarios are presented and evaluated. A heuristic solution procedure is also used to solve two of the models studied. Heuristic solutions to three large synthesis problems are presented. The results of our analysis show that the heuristic performs very well, both in terms of solution quality and solution time, on the two models to which it was applied. It is concluded that the heuristic procedure is the best of the methods considered for solving satellite synthesis problems.
Most crystallographic data processing methods use pixel integration. In serial femtosecond crystallography (SFX), the intricate interaction between the reciprocal lattice point and the Ewald sphere is integrated out by averaging symmetrically equivalent observations recorded across a large number (10 4 -10 6 ) of exposures. Although sufficient for generating biological insights, this approach converges slowly, and using it to accurately measure anomalous differences has proved difficult. This report presents a novel approach for increasing the accuracy of structure factors obtained from SFX data. A physical model describing all observed pixels is defined to a degree of complexity such that it can decouple the various contributions to the pixel intensities. Model dependencies include lattice orientation, unit-cell dimensions, mosaic structure, incident photon spectra and structure factor amplitudes. Maximum likelihood estimation is used to optimize all model parameters. The application of prior knowledge that structure factor amplitudes are positive quantities is included in the form of a reparameterization. The method is tested using a synthesized SFX dataset of ytterbium(III) lysozyme, where each X-ray laser pulse energy is centered at 9034 eV. This energy is 100 eV above the Yb 3+ L -III absorption edge, so the anomalous difference signal is stable at 10 electrons despite the inherent energy jitter of each femtosecond X-ray laser pulse. This work demonstrates that this approach allows the determination of anomalous structure factors with very high accuracy while requiring an order-of-magnitude fewer shots than conventional integration-based methods would require to achieve similar results.
The main objective of the developed software is to reduce the cost per foot during drilling, in other words, optimize the drilling operational parameters in achieving optimum ROP while avoiding critical operational parameters due to either low ROP, drillstring vibration, accelerated cutter wear, or low MSE. The developed software can also be used for post-well analysis to provide insight and lessons learned for future drilling operations. Several functions are available in the software to help the user perform drilling analysis, optimization, and simulation.
Seismic tomography harnesses earthquake data to explore the inaccessible structure of the Earth. Adjoint waveform tomography (AWT), a method of seismic tomography, updates the tomographic model by optimizing the fit between observed earthquake data and synthetic waveforms. The synthetic data are calculated by solving the wave equation through a given 3D model. An important requirement to calculating synthetics is the source information (location, centroid time, depth, and moment tensor). Errors in source information affect the quality of the synthetics produced, which in turn can limit how structure can be inferred in the AWT workflow. Here, to test the effect of updating source information, we used MTTime (Chiang, 2020), a time-domain full-waveform moment tensor inversion code, to calculate the moment tensors and depths of 118 earthquakes that occurred in California and Nevada over a 20-yr period. We calculated 3D Green’s functions using a 3D seismic wavespeed model of California and Nevada (Doody et al., 2023b). We show that the inverted solutions provide better waveform fits than the Global Centroid Moment Tensor catalog and increase usable, well-correlated data by up to 7%. Therefore, we argue that recalculating source parameters should be considered in AWT workflows, particularly for smaller magnitude events (M w > 5.0).
Two-step solar thermochemical cycling (STC) based on nonstoichiometric oxides is an ideal means of solar fuel (e.g., H 2 , CO) production. Screening of nonstoichiometric oxides with excellent thermodynamic performance is key to achieving high solar-to-fuel efficiency. However, application-driven materials assessment intended for reactor-level solar fuel production performance requires mimicking realistic operating conditions of on-sun tests in a laboratory setting, which makes oxides assessment and screening an onerous task to accomplish experimentally. In this work, a rapid assessment and screening model of nonstoichiometric oxides for two-step solar thermochemical cycling assuming fixed-bed flow pattern and quasi-equilibration of the solid with the flowing gas phase is developed, with solar-to-fuel efficiency being the target function of optimization. The model accounts for the thermodynamic parameters of oxide materials and typical operating conditions of experimental thermochemical cycling. This study employed the model to explore and compare two groups of typical nonstoichiometric oxides (CeO 2 - and (LaSr)MnO 3 -based) for their maximum efficiency under their uniquely optimized cycling conditions. The results show that CeO 2 can reach a maximum efficiency of 12.9% at reduction temperature of 1500 °C, which is superior to other candidate materials, including 20 mol% Zr-doped CeO 2 (10.1%), La 0.6 Sr 0.4 MnO 3 (2.5%) and La 0.8 Sr 0.2 MnO 3 (3.3%). Even when the reduction temperature is lowered to 1350 °C, ceria yields the highest efficiency amongst the candidate STC materials considered. The optimal cycling strategy depends on the inherent thermodynamic properties of the oxides. This approach serves as a framework for assessing the maximum efficiency and optimal conditions of candidate thermochemical materials within a range of constraints rather than comparing materials under arbitrary cycling conditions, which may inherently favor one material over another. For oxidation temperatures below 800–1000 °C, the model could be further improved by considering reaction kinetics.
Scientific machine learning (SciML) often operates in ill-conditioned, weakly identifiable regimes due to limited data or indirect observations. In such settings, optimization and inference are highly sensitive to the starting point, making initialization--often under-reported--a consequential degree of freedom. Random initialization is not a neutral default as it induces an implicit prior over candidate solutions and can systematically bias the result, producing large run-to-run variability. Here, we formalize this view by treating initialization as a hidden confounder in SciML and develop a unifying theory for structure-aware initialization via numerical continuation, constructing warm starts from related problem instances. Across representative tasks, including physics-informed neural networks, maximum likelihood estimation, and variational inference, warm starts have been shown to consistently reduce optimization effort and improve reliability.