Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flexible models”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Multiscale Electricity Modeling for Evaluating Carbon Capture and Sequestration Technologies (MEME-CCS)

This effort employs and adapts an existing, rigorous multiscale electricity modeling platform at the National Renewable Energy Laboratory (NREL) to evaluate carbon capture and sequestration (CCS) and negative emissions technologies (NET) from the ARPA-E FLECCS program. NREL's modeling platform includes the Regional Energy Deployment System (ReEDS) electric sector capacity expansion model, which projects future electricity generation mixes at sub-state-level resolution that are downscaled to the unit-level to enable hourly, zonal or nodal electricity production cost modeling in the PLEXOS model. The resulting hourly price data from PLEXOS is provided to technology developers under the ARPA-E FLECCS program to enable technology-specific economic analysis. The ReEDS-PLEXOS modeling suite is well-established for examining electric sector futures with high renewable energy penetrations, energy storage, electrification, and distributed generation. The key advancement proposed herein utilizes collaboration with CCS experts at the University of Wyoming and the FLECCS teams to create innovative methods for representing CCS and NET in the ReEDS and PLEXOS models. Expanded technology options, new operational parameterizations, and detailed data defining CO2 capture, transportation, and storage systems are being integrated into these models to allow an unprecedented combination of scope and resolution for exploring the future of CCS and NET. The resulting capabilities will take advantage of high-performance computing resources to permit wide-ranging scenario analysis to assess CCS and NET deployment under alternative CO2 prices, fossil fuel prices, electricity demand growth, and other electric sector characteristics. Final outcomes will include publicly available hourly grid operation and price data for any U.S. region of interest along with open-access capacity expansion tools for evaluating CCS/NET systems. These products will help an emerging CCS/NET industry in the United States by providing economics-driven guidance to technology developers while informing policy and investment decisions in the public and private sectors.

air capture↗
DOE OSTI2025

Diagnostic-free onboard battery health assessment

Diverse usage patterns induce complex and variable aging behaviors in lithiumion batteries, complicating accurate health diagnosis and prognosis. Separate diagnostic cycles are often used to untangle the battery’s current state of health from prior complex aging patterns. However, these same diagnostic cycles alter the battery’s degradation trajectory, are time-intensive, and cannot be practically performed in onboard applications. Here, in this work, we leverage portions of operational measurements in combination with an interpretable machine learning model to enable rapid, onboard battery health diagnostics and prognostics without offline diagnostic testing and the requirement of historical data. We integrate mechanistic constraints within an encoder-decoder architecture to extract electrode states in a physically interpretable latent space and enable improved reconstruction of the degradation path. The health diagnosis model framework can be flexibly applied across diverse application interests with slight fine-tuning.

battery aging reconstruction↗
DOE OSTI2025

Neural chaos: A spectral stochastic neural operator

Building surrogate models for operators with uncertainty quantification capabilities is essential for many engineering applications where randomness–such as variability in material properties, boundary conditions, and initial conditions–is unavoidable. Polynomial Chaos Expansion (PCE) is widely recognized as a go-to method for constructing stochastic surrogates in both intrusive and non-intrusive ways, and it has recently been used in the context of operator learning. However, its application becomes challenging for complex or high-dimensional processes, as achieving accuracy requires higher-order polynomials, which can increase computational demand and/or the risk of overfitting. Furthermore, PCE requires specialized treatments to manage random variables that are not independent, and these treatments may be problem-dependent or may fail with increasing complexity. Here, in this work, we adopt the same formalism as the spectral expansion used in PCE; however, we replace the classical polynomial basis functions with neural network (NN) basis functions to leverage their expressivity. To achieve this, we propose an algorithm that identifies NN-parameterized basis functions in a purely data-driven manner, without any prior assumptions about the joint distribution of the random variables involved, whether independent or dependent, or about their marginal distributions. The proposed algorithm identifies each NN-parameterized basis function sequentially, ensuring they are orthogonal with respect to the data distribution. The basis functions are constructed directly on the joint stochastic variables without requiring a tensor product structure or assuming independence of the random variables. This approach may offer greater flexibility for complex stochastic models, while simplifying implementation compared to the tensor product structures typically used in PCE to handle random vectors. This is particularly advantageous given the current state of open-source packages, where building and training neural networks can be done with just a few lines of code and extensive community support. We demonstrate the effectiveness of the proposed scheme through several numerical examples of varying complexity and provide comparisons with classical PCE.

Polynomial chaos expansion↗
DOE OSTI2025

The influence of wind veer and drivetrain flexibility on fatigue loading for large floating wind turbines

To reduce costs, offshore wind turbines are expected to be designed with significantly increased rotor diameters. Larger turbines become more flexible and span a larger portion of the atmospheric boundary layer. With these changes, the validity of traditional modeling assumptions should be investigated. This work challenges two common assumptions: (1) that the drivetrain can be considered rigid (except in torsion) and does not couple with the rotor and tower and 2) that wind directional change with height (veer) does not greatly influence the fatigue damage in the tower, blades and drivetrain. Two large semi-submersible floating wind turbines are considered: a 15 and a 22 MW reference turbine. Both use direct-drive generators. Aero-hydro-servo-elastic simulations are performed using OpenFAST, with drivetrain bending flexibility and main bearing response implemented in the coupled analysis. The turbines are subjected to a set of load cases at below-, near- and above-rated mean wind speeds, assembled based on the 3 km Norwegian reanalysis (NORA3) hourly wind and wave hindcast data for Utsira Nord, off the coast of Norway. In each load case, conditions with and without veer are simulated to evaluate the influence of veer on damage equivalent loads (DELs) of the turbine tower, blades and main bearings. Further, these load cases are applied to evaluate the influence of drivetrain flexibility on aero-elastic turbine response. The results indicate that, depending on the veer gradient, mean wind speed, operating regime and turbine size, veer can be very important for tower-top DELs and the fluctuations of main bearing radial loads, while main bearing and blade-root flapwise DELs are less affected. Considering these specific load cases and turbine models, drivetrain flexibility is found to significantly influence tower-top DELs of the largest turbine: the tower-top fore-aft and torsional damage equivalent moments of the 22 MW turbine are reduced by more than 20 % at near-rated wind speeds when the drivetrain is modeled as flexible.

17 WIND ENERGY↗
DOE OSTI2022

Toward Addressing the Challenge to Predict the Heat Capacities of RDX and HMX Energetic Materials

Availability of heat capacity as function of pressure and temperature is an essential prerequisite for development of a computational multiscale strategy capable to address the evolution of microstructure and energy release in advanced high energy density materials. In the case of 1,3,5-trinitro-1,3,5-triazinane (RDX) and 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) systems as two of the most studied energetic materials, there are substantial gaps in experimental data, with available heat capacities values distributed only in a region close to standard ambient conditions. In this study we demonstrate how these major experimental limitations can be addressed in the case of the RDX and HMX systems based on the combined use of classical and quantum mechanical calculations. We show that by considering ideal gas properties evaluated using quantum mechanical methods, and residual properties obtained from molecular simulations using fully flexible atomistic force field models, excellent agreement can be obtained for the predicted heat capacities to the most recent experimental values. An important advantage of the current computational methodology is that it allows evaluation of both constant-volume and constant-pressure heat capacities for a broad interval of temperatures and pressures, which encompasses solid and liquid phases conditions. In the case of the solid α and γ phases of RDX and the β phase of HMX, the predicted results follow closely both the available experimental data at standard ambient conditions and the results obtained using density functional theory calculations at high pressures, a regime where experimental data are not available. A perspective to expand the current methodology is also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Long‐Range Electron Transport Rates Depend on Wire Dimensions in Cytochrome Nanowires

Abstract The ability to redirect electron transport to new reactions in living systems opens possibilities to store energy, generate new products, or probe physiological processes. Recent work by Huang et al. showed that 3D crystals of small tetraheme cytochromes (STC) can transport electrons over nanoscopic to mesoscopic distances by an electron hopping mechanism, making them promising materials for nanowires. However, fluctuations at room temperature may distort the nanostructure, hindering efficient electron transport. Classical molecular dynamics simulations of these fluctuations at the nano‐ and mesoscopic scales allowed us to develop a graph network representation to estimate maximum electron flow that can be driven through STC wires. In longer nanowires, transient structural fluctuations at protein‐protein interfaces tended to obstruct efficient electron transfer, but these blockages are ameliorated in thicker crystals where alternative electron transfer pathways become more efficient. The model implies that more flexible proteinprotein interfaces limit the required minimum diameter to carry currents commensurate with conventional electronics.

Chemistry↗
DOE OSTI2021

Pyomo.GDP: an ecosystem for logic based modeling and optimization development

We present three core principles for engineering-oriented integrated modeling and optimization tool sets—intuitive modeling contexts, systematic computer-aided reformulations, and flexible solution strategies—and describe how new developments in Pyomo.GDP for Generalized Disjunctive Programming (GDP) advance this vision. We describe a new logical expression system implementation for Pyomo.GDP allowing for a more intuitive description of logical propositions. The logical expression system supports automated reformulation of these logical constraints to linear constraints. We also describe two new logic-based global optimization solver implementations built on Pyomo.GDP that exploit logical structure to avoid “zero-flow” numerical difficulties that arise in nonlinear network design problems when nodes or streams disappear. These new solvers also demonstrate the capability to link to external libraries for expanded functionality within an integrated implementation. We present these new solvers in the context of a flexible array of solution paths available to GDP models. Finally, we present results on a new library of GDP models demonstrating the value of multiple solution approaches.

42 ENGINEERING↗
DOE OSTI2022

Qualitative trend analysis based on a mixed-integer representation

Shape constrained spline fitting is a useful method to impose prior knowledge onto flexible semi-parametric models during parameter estimation. Most typically, the function shape is imposed through order restrictions on the regression coefficients. The intended shape is considered known or selected based on heuristic rules. In this study, we present a method to estimate the optimal set of order restrictions to segment a univariate data series into episodes with distinct shapes. This is also known as the qualitative trend analysis (QTA) problem. The obtained solution uses a trade-off between lack-of-fit and model complexity. Further, our practical implementation takes inspiration from the generalized order restricted information criterion (GORIC) for inequality-constrained model selection. From this, one learns (a) that QTA can be formulated as a mixed-integer quadratic program (MIQP) and (b) that the newly proposed mixed order restricted information criterion (MORIC) enables optimal segmentation. This is illustrated through didactic case studies.

42 ENGINEERING↗
DOE OSTI2024

Polyconvex neural network models of thermoelasticity

Machine-learning function representations such as neural networks have proven to be excellent constructs for constitutive modeling due to their flexibility to represent highly nonlinear data and their ability to incorporate constitutive constraints, which also allows them to generalize well to unseen data. Here, in this work, we extend a polyconvex hyperelastic neural network framework to (isotropic) thermo-hyperelasticity by specifying the thermodynamic and material theoretic requirements for an expansion of the Helmholtz free energy expressed in terms of deformation invariants and temperature. Different formulations which a priori ensure polyconvexity with respect to deformation and concavity with respect to temperature are proposed and discussed. The physics-augmented neural networks are furthermore calibrated with a recently proposed sparsification algorithm that not only aims to fit the training data but also penalizes the number of active parameters, which prevents overfitting in the low data regime and promotes generalization. The performance of the proposed framework is demonstrated on synthetic data, which illustrate the expected thermomechanical phenomena, and existing temperature-dependent uniaxial tension and tension-torsion experimental datasets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

An Ab Initio -Derived Force Field for Amorphous Silica Interfaces for Use in Molecular Dynamics Simulations

Here, we present a classical interatomic force field, silica-DDEC, to describe the interactions of amorphous and crystalline silica surfaces, parametrized using density functional theory-based charges. Charge schemes for silica surfaces were developed using the density-derived electrostatic and chemical (DDEC) method, which reproduces atomic charges of the periodic models as well as the electrostatic potential away from the atom sites. Lennard–Jones parameters were determined by requiring the correct description of (i) the amorphous silica density, coordination defects, and local coordination geometry, relative to experimental measurements, and (ii) water-silica interatomic distances compared with ab initio results. Deprotonated surface silanol sites are also described within the model based on DDEC charges. The result is a general electronic structure-derived model for describing fully flexible amorphous and crystalline silica surfaces and interactions of liquids with silica surfaces of varying structure and protonation state.

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

Energy, power, and infrastructure demands from electrifying airport ground support equipment at United States airports

As the airline industry seeks to reduce costs and transition to clean energy, electric ground support equipment is emerging as a favorable option. The integration of electric ground support equipment into airport operations requires careful planning for vehicle deployment, charging infrastructure, and grid impacts. We develop a flexible, bottom-up modeling framework to assess energy and infrastructure needs across more than 300 U.S. airports. Our analysis estimates site-specific power and energy demand, equipment counts by type, charger requirements, and costs. Here we quantify the magnitude of new electrical loads created by the electrification of airport ground support equipment, finding that peak power demand at the largest airports can reach up to 20 megawatts, with annual electricity consumption approaching 51,000 megawatt-hours. We further show that behind-the-meter battery energy storage systems and solar photovoltaic systems can reduce peak load and lower total system costs by as much as 10 million dollars.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

TDCOSMO - IV. Hierarchical time-delay cosmography – joint inference of the Hubble constant and galaxy density profiles

The H0LiCOW collaboration inferred via strong gravitational lensing time delays a Hubble constant value of H0 = 73.3−1.8+1.7 km s−1 Mpc−1, describing deflector mass density profiles by either a power-law or stars (constant mass-to-light ratio) plus standard dark matter halos. The mass-sheet transform (MST) that leaves the lensing observables unchanged is considered the dominant source of residual uncertainty in H0. We quantify any potential effect of the MST with a flexible family of mass models, which directly encodes it, and they are hence maximally degenerate with H0. Our calculation is based on a new hierarchical Bayesian approach in which the MST is only constrained by stellar kinematics. The approach is validated on mock lenses, which are generated from hydrodynamic simulations. We first applied the inference to the TDCOSMO sample of seven lenses, six of which are from H0LiCOW, and measured H0 = 74.5−6.1+5.6 km s−1 Mpc−1. Secondly, in order to further constrain the deflector mass density profiles, we added imaging and spectroscopy for a set of 33 strong gravitational lenses from the Sloan Lens ACS (SLACS) sample. For nine of the 33 SLAC lenses, we used resolved kinematics to constrain the stellar anisotropy. From the joint hierarchical analysis of the TDCOSMO+SLACS sample, we measured H0 = 67.4−3.2+4.1 km s−1 Mpc−1. This measurement assumes that the TDCOSMO and SLACS galaxies are drawn from the same parent population. The blind H0LiCOW, TDCOSMO-only and TDCOSMO+SLACS analyses are in mutual statistical agreement. The TDCOSMO+SLACS analysis prefers marginally shallower mass profiles than H0LiCOW or TDCOSMO-only. Without relying on the form of the mass density profile used by H0LiCOW, we achieve a ∼5% measurement of H0. While our new hierarchical analysis does not statistically invalidate the mass profile assumptions by H0LiCOW – and thus the H0 measurement relying on them – it demonstrates the importance of understanding the mass density profile of elliptical galaxies. The uncertainties on H0 derived in this paper can be reduced by physical or observational priors on the form of the mass profile, or by additional data.Key words: gravitational lensing: strong / galaxies: general / galaxies: kinematics and dynamics / distance scale / cosmological parameters / cosmology: observations⋆ The full analysis is available at https://github.com/TDCOSMO/hierarchy_analysis_2020_public.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A New 1D Model for Thermal Mixing and Stratification in Advanced Reactor Transients

Thermal mixing and stratification in large pools and enclosures play a critical role in the safety and performance of pool-type nuclear reactors, particularly during transient scenarios involving significant temperature differences between incoming and bulk coolant. Accurate modeling of these phenomena is essential for predicting system behavior and supporting passive safety features such as natural circulation. Here, this paper presents a new 1D model for thermal mixing and stratification, developed and implemented in the SAM code. The model represents a large pool as 1D coolant jet channels and zero-dimensional bulk pool volumes, enabling the simulation of a wide range of flow configurations, including hot and cold jet interactions, stratified layers, and the influence of complex geometries such as ceilings, free surfaces, and internal obstacles. Heat exchange between jet and pool regions is governed by closure relations calibrated against 3D computational fluid dynamics (CFD) simulations. The model improves upon earlier approaches by incorporating time-dependent jet characteristics and capturing the associated delay effects more accurately. Code-to-code comparisons and validation against experimental data from the Thermal Stratification Test Facility demonstrate the model’s accuracy and flexibility. This work offers two key contributions: (1) an efficient and robust method for simulating thermal mixing and stratification at the system level, eliminating the need for external coupling between system analysis codes and CFD, and (2) a significant enhancement of SAM’s capabilities to analyze thermal stratification phenomena in advanced reactor systems.

SAM↗
DOE OSTI2026

Microscopic constraints for the equation of state and structure of neutron stars: A Bayesian model mixing framework

Bayesian model mixing (BMM) is a statistical technique that can combine constraints from different regions of an input space in a principled way. Here we extend our BMM framework for the equation of state (EOS) of strongly interacting matter from symmetric nuclear matter to asymmetric matter, specifically focusing on zero-temperature, charge-neutral, 𝛽-equilibrated matter. We use Gaussian processes (GPs) to infer constraints on the neutron-star matter EOS at intermediate densities from two different microscopic theories: chiral effective-field theory (𝜒⁢EFT) at baryon densities around nuclear saturation, 𝑛 𝐵 ∼ 𝑛 0 , and perturbative QCD at asymptotically high baryon densities, 𝑛 𝐵 ⩾ 20⁢𝑛 0 . The uncertainties of the 𝜒⁢EFT and pQCD EOSs are obtained using the BUQEYE truncation error model. We demonstrate the flexibility of our framework through the use of two categories of GP kernels: conventional stationary kernels and a nonstationary changepoint kernel. We use the latter to explore potential constraints on the dense matter EOS by including exogenous data representing theory predictions and heavy-ion collision measurements at densities ⩾ 2⁢𝑛 0 . We also use our EOSs to obtain neutron-star mass-radius relations and their uncertainties. Finally, our framework, whose implementation will be available through a GitHub repository, provides a prior distribution for the EOS that can be used in large-scale neutron-star inference frameworks.

Bayesian methods↗
DOE OSTI2024

Urban-Scale Control of School Bus Fleet Charging and Discharging Strategies Using Single and Multi-Stage Optimization

This paper presents a dual-strategy approach to optimizing charging and discharging schedules for school bus fleets, using the limited charging infrastructure effectively. We aim to ensure that each bus is fully charged for daily operations and aids in grid stability during peak demand. The first strategy utilizes linear programming to schedule overnight charging at available station sockets and strategic discharging during peak periods, efficiently coordinating limited resources. The second strategy employs metaheuristic techniques for continuous optimization, focusing on precise power requirements and offering greater flexibility than the linear model.

Selim, Alaa↗
DOE OSTI2021

Robust Deep Gaussian Process-Based Probabilistic Electrical Load Forecasting Against Anomalous Events

The abnormal events, such as the unprecedented COVID-19 pandemic, can significantly change the load behaviors, leading to huge challenges for traditional short-term forecasting methods. This article proposes a robust deep Gaussian processes (DGP)-based probabilistic load forecasting method using a limited number of data. Since the proposed method only requires a limited number of training samples for load forecasting, it allows us to deal with extreme scenarios that cause short-term load behavior changes. In particular, the load forecasting at the beginning of abnormal event is cast as a regression problem with limited training samples and solved by double stochastic variational inference DGP. The mobility data are also utilized to deal with the uncertainties and pattern changes and enhance the flexibility of the forecasting model. The proposed method can quantify the uncertainties of load forecasting outcomes, which would be essential under uncertain inputs. Extensive comparison results with other state-of-the-art point and probabilistic forecasting methods show that our proposed approach can achieve high forecasting accuracies with only a limited number of data while maintaining the excellent performance of capturing the forecasting uncertainties.

anomalous events↗
DOE OSTI2020

Performance simulation of the soft gamma-ray concentrator

The soft gamma-ray concentrator is a telescope mission concept utilizing a suitable arrangement of bent multilayer structures of alternating low- and high-density materials. This lens is able to channel gamma-ray photons via total external reflection and concentrate the incident radiation to a point. The channeling technique offers the potential for concentrating gamma rays with focal lengths <10 m and energies >100 keV, beyond the reach of current grazing-incidence hard x-ray mirrors. For the performance estimation of such an instrument, we have developed a flexible set of computer modeling tools to compute the optical properties of multilayer structures, predict the channeling efficiency for a given multilayer configuration, and aid in the optimization of potential gamma-ray concentrator-based telescope designs. This modeling includes the multilayer optical properties calculated by the IMD software, the ray tracing using an IDL code, and the focal plane detector simulation by MEGAlib. We illustrate the potential of this approach by presenting simulated astronomical observations from a balloon-borne platform. The final result, including simulated effective area, instrument sensitivity, and polarization performance, shows that the gamma-ray concentrator will provide greatly increased sensitivity for next-generation soft gamma-ray missions with modest cost and complexity.

47 OTHER INSTRUMENTATION↗