Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “decoupling”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quantifying transport and electrocatalytic reaction processes in a gastight rotating cylinder electrode reactor via integration of Computational Fluid Dynamics modeling and experiments

Understanding the complexity of the multiple processes of mass, momentum, charge, and heat transport, and how these affect reaction kinetics at the electrode/electrolyte interface is one of the major challenges in the field of energy and catalysis. The rapid and rational scale-up of electrocatalytic systems to industrial scales require a detailed understanding of nonlinear transport-reaction processes, accessible only through the building of multi-physics models that capture with high fidelity the complexity of real-world devices. The gastight rotating cylinder electrode (RCE) reactor is a promising lab-scale tool that can decouple transport from intrinsic kinetics to generate data for first-principle models useful in the design of industrial, electrochemical reactors. Computational Fluid Dynamics (CFD) studies have previously been used to investigate the bulk flow in RCE reactors for simple corrosion and electroplating processes. However, the quantification of changes in local concentration within the viscous layer where catalysis takes place requires capturing the correct flow conditions inside the hydrodynamic boundary layer near the surface of the electrode. Further, this requires simulations with spatial resolution in the nm and μm scale and temporal resolutions between ms and s scales that are similar to the timescales for reactions on the electrode surface. In this study, experimental electrocatalysis is combined with CFD modeling to elucidate and parameterize the hydrodynamics in a gastight RCE reactor. CFD simulations of the electrochemical ferricyanide reduction reaction under mass transport limited conditions are used to evaluate the validity of the CFD model parameters by comparing calculated dimensionless mass transport descriptors to dimensionless correlations obtained experimentally. Justifications for assumptions and details of the simulation methods used in this study are presented to provide a detailed understanding of the effect that each model parameter has on the ability to accurately simulate electrocatalysis in RCE systems. The simulation methodology reported here is a first step towards the development of multi-scale models for the study of transport dependent electrocatalytic processes, such as the electrochemical transformation of CO 2 to fuels and chemicals.

42 ENGINEERING↗
DOE OSTI2025

Enhancing EnergyPlus capabilities to model dynamic building envelopes using python plugin

Nearly half of the energy consumption in the United States is related to buildings, resulting in an urgent need to develop innovative technologies to improve building energy efficiency. Dynamic building envelopes, comprising switchable insulation and thermal energy storage materials, have been proposed recently as a promising solution to reduce buildings' heating and cooling loads by thermally coupling the indoor environment with the ambient environment when beneficial while decoupling them when outdoor conditions are not favorable. Although various related technologies are still underway, the whole-building energy modeling tools, like EnergyPlus, do not have the capability to simulate the transient and dynamic nature of dynamic envelope materials and components to accurately predict their impact on building energy use. The objective of this study is to formulate a method in EnergyPlus simulation engine to model multilayer envelopes, comprising dynamic building materials with variable thermophysical properties, and discuss the changes made to the program using a Python plugin. Furthermore, the thermal performance of the dynamic envelopes using the proposed method is compared and verified with the results from a well-established commercial code, COMSOL Multiphysics. A parametric assessment is also conducted to evaluate the energy efficiency benefits of dynamic envelopes in a single-family residential building, demonstrating total annual energy savings up to 11.6 %, when a dynamic envelope operates alone, and up to 18.2 % when it is combined with a thin layer of phase change material as a thermal storage medium. Finally, a United States wide energy efficiency assessment is presented to showcase the geographical spread of the energy savings. The method designed and implemented in this study provides the researchers with the ability to implement their dynamic insulation methods in EnergyPlus and evaluate the whole building energy impact.

25 ENERGY STORAGE↗
DOE OSTI2026

Physics vs structure: A systematic benchmark of learning strategies for multi-zone building thermal dynamics

Recent advances in physics-informed and data-driven machine learning promise improved thermal models for advanced building control, yet there is limited quantitative evidence on when added physics structure and architectural complexity are beneficial. Here, this work presents a systematic benchmark of five representative system identification methods for modeling multi-zone building thermal dynamics: linear state-space models, multi-layer perceptrons, neural state-space models, neural ordinary differential equations, and physically-consistent neural networks. The methods are evaluated across multiple data regimes and zone coupling strategies. Using a high-fidelity multi-zone commercial building emulator, we examine short-term and long-term prediction accuracy, computational efficiency, and ease of development. Our results reveal critical trade-offs between prediction performance, model complexity, and physical consistency. We demonstrate that decoupled, nonlinear black-box models consistently outperform coupled physics-constrained architectures in both predictive accuracy and out-of-distribution robustness in majority of the test cases for the building type considered in the study. Our findings quantify the cost of complexity in building thermal modeling and provide concrete, actionable, scenario-based guidelines for selecting model classes for control-oriented applications.

Building thermal modeling↗
DOE OSTI2021

Energy-efficient sorption-based gas clothes dryer systems

Standard electric resistance and fuel-driven dehydration technologies exhibit a maximum coefficient of performance of well below 1 mainly due to enthalpy losses associated with the air leaving the dehydration system. To improve energy efficiency, condensing dryer systems condense the moisture captured from a product in a closed-loop air circulation cycle. Existing condensing dehydration systems including heat pump dryers, however, need to significantly cool the air to achieve dehumidification. The added cooling and subsequent heating to return the air to a desired drying temperature consume substantial energy and thus reduce drying performance. Here, an innovative sorption-based gas dehydration system is proposed to overcome barriers deteriorating energy efficiency in existing gas, electric, or heat pump dryer systems. Decoupling latent and sensible loads, the system employs a liquid-desiccant solution to directly capture air humidity, thereby allowing circulation of the air in a closed loop to achieve high drying energy efficiency. In other words, the system captures waste latent heat from the moisture produced during the dehydration process and reuses it to improve energy efficiency. This study focuses on a comprehensive quasi-steady-state thermodynamic modeling of the proposed sorption-based dehydration concept employed for a gas clothes dryer application to predict transient response and overall drying performance (i.e., time and energy metrics). The analysis indicates the proposed sorption-based gas clothes dryer system can deliver a specific moisture extraction rate of 1.71 kg of water per kWh (i.e., a combined energy factor of 3.167 kg (6.98 lbm) of dry cloth per kWh) with a drying time of 44 min. This is a 112% energy improvement compared with state-of-the-art gas clothes dryers exhibiting a combined energy factor of 1.50 kg (3.3 lbm) of dry cloth per kWh. The technology pursued here can potentially be employed as a platform for many fuel-driven equipment to take advantage of available waste thermal energy in the environment instead of simply burning a fuel.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

The reservoir network: A new network topology for district heating and cooling

Thermal district networks are effective solutions to substitute fossil fuels with renewable energy sources for heating and cooling. Moreover, thermal networking of buildings allows energy efficiency to be further increased. The waste heat from cooling can be reused for heating in thermal district systems. Because of bidirectional energy flows between prosumers, thermal networks require new hydraulic concepts. In this work, we present a novel network topology for simultaneous heating and cooling: the reservoir network. The reservoir network is robust in operation due to hydraulic decoupling of transfer stations, integrates heat sources and heat sinks at various temperature levels and is flexible in terms of network expansion. Furthermore, we used Modelica simulations to compare the new single-pipe reservoir network to a basecase double-pipe network, taking yearly demand profiles of different building types for heating and cooling. The electric energy consumed by the heat pumps and circulations pumps differs between the reservoir and base case networks by less than 1%. However, if the reservoir network is operated with constant instead of variable mass flow rate, the total electrical consumption can increase by 48% compared to the base case. As with any other network topology, the design and control of such networks is crucial to achieving energy efficient operation. Investment costs for piping and trenching depend on the district layout and dimensioning of the network. If a ring layout is applied in a district, the reservoir network with its single-pipe configuration is more economical than other topologies. For a linear layout, the piping costs are slightly higher for the reservoir network than for the base case because of larger pipe diameters.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Experimental study on a novel three-phase absorption thermal battery with high energy density applied to buildings

Buildings contribute to 75% of total electricity consumption and 80% of peak electricity demand in the United States. A significant portion of building electricity consumption goes to space conditioning. Therefore, using thermal energy storage (TES) to decouple the building electricity consumption and the fluctuating space-conditioning load has significant potential to shift peak electricity demand and improve the stability of the grid. To provide effective TES within the spatial constraints in buildings, a three-phase absorption thermal battery (TATB) system with very high energy storage density was designed and studied. The TATB stores low-temperature heat in hydrate crystals of salts and supports space conditioning through absorption systems. In this study, a benchtop TATB using lithium chloride (LiCl) hydrate crystals as the storage material was tested under practical operating conditions for its energy storage density and discharge rate. We report the TATB successfully generated LiCl hydrate crystals during charging mode and then dissolved the crystals during discharging mode. A material-based energy storage density of 300 kWh/m 3 (specific energy of 903 kJ/kg) was achieved with a discharge rate of up to 1.3 kW thermal energy. The test results verified the high performance of the TATB and encouraged further development of this technology.

25 ENERGY STORAGE↗
DOE OSTI2022

A fast and accurate modeling approach for water and steam thermodynamics with practical applications in district heating system simulation

In U.S. district heating (DH) systems, steam is the most common heat transport medium. Industry demand for new advanced modeling capabilities of complete steam DH systems is increasing; however, the existing models for water/steam thermodynamics are too slow for large system simulations because of computationally expensive algebraic loops that require the solution to nonlinear systems of equations. For practical applications, this work presents a novel split-medium approach that implements numerically efficient liquid water models alongside various water/steam models, breaking costly algebraic loops by decoupling mass and energy balance equations. New component models for steam DH systems are also presented. We implemented the models in the equation based Modelica language and evaluated accuracy and computing speed across multiple scales: from fundamental thermodynamic properties to complete districts featuring 10 to 200 buildings. Compared to district models with the IF97 water/steam model and equipment models from the Modelica Standard Library, the new implementation improves the scaling rate for large districts from cubic to quadratic with negligible compromise to accuracy. Additionally, for an annual simulation with 180 buildings, this translates to a computing time reduction from 33 to 1-1.5 h. These results are critically important for industry practitioners to simulate steam DH systems at large scales.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Anomalously high elastic modulus of a poly(ethylene oxide)-based composite electrolyte

The practical use of lithium metal anodes in solid-state batteries requires a polymer membrane with high lithium-ion conductivity, thermal/electrochemical stability, and mechanical strength. The primary challenge is to effectively decouple the ionic conductivity and mechanical strength of the polymer electrolytes. We report a remarkably facile single step synthetic strategy based on in-situ crosslinking of poly(ethylene oxide) (xPEO) in the presence of a woven glass fiber (GF). Such a simple method yields composite polymer electrolytes (CPE) of anomalously high elastic modulus up to 2.5 GPa over a broad temperature range (20 °C – 245 °C) that has never been previously documented. An unsupervised machine learning algorithm, K-mean clustering analysis, was implemented on the hyperspectral Raman mapping at the xPEO/GF interface. Using such a unique means, we show for the first time that the promoted mechanical strength originates from xPEO and GF interactions through dynamic hydrogen and ionic bonding. High ionic conductivity is achieved by the addition plasticizer (e.g. tetraglyme), where trifluoromethanesulfonate anions are tethered to the xPEO matrix and Li + cations are favorably transported through coordination with the plasticizer. Further, stringent galvanostatic cycling tests indicates the CPE can be stably cycled for >3000 h in a Li-metal symmetric cell at a moderate temperature (nearly 1500 Coulombs/cm 2 Li equivalents), outperforming most of the PEO-based electrolytes. The GF reinforced CPE reported here has multifunctional uses, such as solid electrolytes for all solid-state batteries and membranes for redox-flow batteries. Although the focus of this study is on lithium-based batteries, the results are equally promising for other alkali metal based batteries such as sodium and potassium.

25 ENERGY STORAGE↗
DOE OSTI2022

Mitigation of rapid capacity decay in silicon-LiNi 0.6 Mn 0.2 Co 0.2 O 2 full batteries

Silicon (Si)-based materials have been considered as the most promising anode materials for high-energy-density lithium-ion batteries because of their higher storage capacity and similar operating voltage, as compared to the commercial graphite (Gr) anode. But the use of Si anodes including silicon-graphite (Si-Gr) blended anodes often leads to rapid capacity decay in Si-Gr/LiNixMnyCo z O 2 (x+y+z=1) full cells, which has been attributed to surface instability of the Si component. In addition to stabilizing the surface, this work investigates the potential of the Si-Gr blended anodes in a full-cell configuration and its impact on the capacity contribution from active components. Using dQ/dV plots of the full cells, a powerful but simple-to-implement differential potential approach is developed to decouple the capacity contribution and degradation from the graphite and silicon components. Data collected from three-electrode cells confirm the results from the differential potential approach, which suggests a voltage slippage to a higher voltage at the blended anode side. Additionally, the voltage slippage causes a reduced utilization of the Gr component and exacerbates side reactions between the Si-Gr anode and carbonate electrolytes. Furthermore, based on these failure mechanisms, we adopted a mitigation strategy to tune the open circuit voltage of the prelithiated anode while stabilizing the surface. As a result, the full cells with the modified Si-Gr anodes (mass loading, 2.5 mAh/cm 2 ) offer a highly reversible full-cell energy density of 390 Wh/kg (based on the mass of both anode and cathode materials in a full cell) with a cycling CE of 99.9% over 200 cycles.

25 ENERGY STORAGE↗
DOE OSTI2020

Pervasive low-velocity layer atop the 410-km discontinuity beneath the northwest Pacific subduction zone: Implications for rheology and geodynamics

Regional triplication waveforms of five intermediate-depth events are modeled to simultaneously obtain the compressional (P) and shear (SH) wave velocity structure beneath northwestern Pacific subduction zone. Both the P- and SH-wave velocity models for three different sub-regions show a low-velocity layer (LVL) with a thickness of ~55-80 km lying above the 410-km discontinuity with a ~900 km lateral extent from the Japan Sea to the northeastern Asian continental margin. With the dihedral angle approaching to zero around 400 km, a minute amount of melt atop the 410-km discontinuity caused by the hydrous slab might completely wet olivine grain boundaries and result in a low seismic velocity layer in this specific subduction zone. This mechanism suggests that the 410-LVL is a low viscosity zone that would partially decouple the upper mantle from the transition zone. We infer that the widespread 410-LVL provides evidence for a water-bearing mantle transition zone beneath the western Pacific subduction zone.

58 GEOSCIENCES↗
DOE OSTI2026

Adaptation of virtual synchronous generators to dynamic conditions in power grids

Virtual synchronous generators (VSGs) are widely adopted as grid-forming controls for inverter-based resources. However, when grid conditions vary significantly as characterized by changes in short-circuit ratio (SCR) and the reactance-to-resistance (X/R) ratio, fixed-gain designs and the commonly used P–Q decoupling assumption can become inaccurate. Such conditions can degrade transient power performance, leading to oscillations, prolonged settling, and overshoot, particularly in stiff-grid operating points. This paper quantifies how grid strength and impedance-dependent coupling affect the active–reactive power dynamics of a conventional VSG over a broad range of SCR and X/R values. An adaptive VSG tuning framework is then developed by combining (i) a coupling-explicit, impedance-parameterized state-space model to enable systematic controller synthesis, (ii) a full-state-feedback law designed via pole placement to meet prescribed damping and settling-time specifications, and (iii) a physics-informed neural network (PINN)–based online grid-impedance estimator that updates controller gains in real time as grid conditions vary. Offline simulations in MATLAB/Simulink and real-time validation on an OPAL-RT platform show that the proposed method preserves consistent damping and settling behavior with reduced overshoot across wide SCR and X/R ranges, compared with fixed-gain VSG baselines.

Adaptive control↗
DOE OSTI2021

The quest for efficient oxygenated fuels: Examining interactions between lubricant components and oxygenates

This work empirically evaluated oxygenated blendstock-lubricant compatibility using a novel but simple method. The reactivity or compatibility was evaluated via gas chromatograph equipped with a flame ionization detector (GC-FID) before and after heating for 2h at 150 °C, which would roughly mimic conditions encountered in the combustion chamber of an Spark Ignition (SI) or Compression-Ignition (CI) engine. Five oxygenate blendstocks spanning several functional groups and chain lengths (prenol, isoprenol, butyl acetate, isohexanol, polyoxymethylene dimethyl ethers) were evaluated in a base fuel with nine variations of lubricant components of controlled concentrations, including fully formulated commercial lubricants. Baseline comparative examples were also generated, to decouple the effect of the oxygenate on the mixture, by heating only the base fuel with the lubricant component. Results demonstrate the relative compatibility between noted lubricants and oxygenates in base fuel upon heating. This was illustrated by minimal to no changes in the GC-FID traces. However, reactions of several of the chosen lubricants, especially Molybdenum di (2-ethylhexyl) phosphorodithioate (Molyvan L) and zinc dialkyldithiophosphate (ZDDP), with prenol and polyoxymethylene dimethyl ethers blendstocks showed substantial side reactions compared to corresponding baseline examples. Control experiments without oxygenates but comprising the same lubricant components showed no change in the mixture with heating. We concluded the incompatibility is due to the oxygenate and not the components of the base fuel. This preliminary work is valuable in establishing trends as well as provide useful information when a new component is introduced in the fuel.

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

Plasma shape and position control development for NSTX-U using the GSEvolve plasma simulator

NSTX-U discharges have been reproduced using the GSEvolve plasma simulator in order to improve their control and reduce oscillations during the plasma current (Ip) ramp-up phase, and with the ultimate goal of enabling reliable H-mode access with a double-null shape at high elongation, low internal inductance, and high plasma performance. The control objective is to eliminate undesired oscillations in the plasma shape and vertical position that start when the plasma diverts during the Ip ramp-up phase, and that happened recurrently during the first NSTX-U experimental campaign. Model-based control solutions have been developed, including alternative H-mode access recipes and a Multi-Input Multi-Output (MIMO) decoupling controller. These solutions have demonstrated their capability to reduce or eliminate the aforementioned oscillations in GSEvolve simulations, enabling reproducible H-mode access. The results of this work aim to inform NSTX-U experimental operations on potential control developments to maximize its scientific output, and highlight the powerful capabilities of high-fidelity control-simulation frameworks like GSEvolve.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Solid-phase Mn speciation in suspended particles along meltwater-influenced fjords of West Greenland

Manganese (Mn) is an essential micro-nutrient that can limit or, along with iron (Fe), co-limit phytoplankton growth in the ocean. Glacier meltwater is thought to be a key source of trace metals to high latitude coastal systems, but little is known about the nature of Mn delivered to glacially-influenced fjords and adjacent coastal waters. In this work, we combine in-situ dissolved Mn (dMn) measurements of surface waters with Mn K-edge X-ray absorption spectroscopy (XAS) data of suspended particles in four fjords of West Greenland. Data were collected from transects of up to 100 km in fjords with different underlying bedrock geology from 64 to 70°N. We found that dMn concentrations generally decreased conservatively with increasing salinity (from 80 to 120 nM at salinity < 8 to < 40 nM at salinities > 25). Dissolved Fe (dFe) trends in these fjords similarly declined with increasing distance from glacier outflows (declining from > 20 nM to < 8 nM). However, the dMn/dFe ratio increased rapidly likely due to the greater stability of dMn at intermediate salinities (i.e. 10–20) compared to rapid precipitation of dFe across the salinity gradient. The XAS data indicated a widespread presence of Mn(II)-rich suspended particles near fjord surfaces, with structures akin to Mn(II)-bearing phyllosilicates. However, a distinct increase in Mn oxidation state with depth and the predominance of birnessite-like Mn(IV) oxides was observed for suspended particles in a fjord with tertiary basalt geology. The similar dMn behaviour in fjords with different suspended particle Mn speciation (i.e., Mn(II)-bearing phyllosilicates and Mn(IV)-rich birnessite) is consistent with the decoupling of dissolved and particulate Mn and suggests that dMn concentrations on the scale of these fjords are controlled primarily by dilution of a freshwater dMn source rather than exchange between dissolved and particle phases. This work provides new insights into the Mn cycle in high latitude coastal waters, where small changes in the relative availabilities of dMn, dFe and macronutrients may affect the identity of the nutrient(s) proximally limiting primary production.

58 GEOSCIENCES↗
DOE OSTI2023

Rubidium and potassium isotopic variations in chondrites and Mars: Accretion signatures and planetary overprints

As moderately volatile elements, isotopes of Rb and K can trace volatilization processes in planetary bodies. Rubidium isotopic data are however very scarce, especially for non-carbonaceous meteorites. Here, in this study, we report combined Rb and K isotopic data (δ 87/85 Rb and δ 41/39 Κ) for 7 ordinary, 6 enstatite, and 4 Martian meteorite falls to understand the causes for the variations in volatile abundances and isotopic compositions. Bulk Rb and K isotopic compositions of planetary bodies are estimated to be (Table 1): Mars +0.10 ± 0.03 ‰ for Rb and -0.26 ± 0.05 ‰ for K, bulk OCs $-0.12^{+0.15}_{-0.24}$ ‰ for Rb and $-0.72^{+0.28}_{-0.41}$ ‰ for K, bulk ECs $-0.02^{+0.29}_{-0.26}$ ‰ for Rb and $-0.33^{+0.37}_{-0.23}$ ‰ for K. The bulk K isotopic compositions of subgroup OCs are estimated to be $-0.72^{+0.26}_{-0.55}$ ‰ for H chondrites, $-0.71^{+0.23}_{-0.39}$ ‰ for L chondrites, and $-0.77^{+0.63}_{-0.30}$ ‰ for LL chondrites. A broad correlation between the Rb and K isotopic compositions of planetary bodies is observed. The correlation follows a slope that is consistent with kinetic evaporation and condensation processes, suggesting volatility-controlled mass-dependent isotope fractionation (as opposed to nucleosynthetic anomalies). Individual ordinary and enstatite chondrites show large Rb and K isotopic variations (-1.02 to +0.29 ‰ for Rb and -0.91 to -0.15 ‰ for K). Samples of lower metamorphic grades display correlated elemental and isotopic fractionations between Rb and K, while samples of higher metamorphic grades show great scatter, suggesting that chondrite parent-body processes have decoupled the two elements and their isotopes at the sample scale. Several processes could have contributed to the observed isotopic variations of Rb and K, including (i) chondrule “nugget effect”, (ii) volatilization during parent-body thermal metamorphism (heat-induced vaporization and gas transport within parent bodies), (iii) thermal diffusion during parent-body metamorphism, and (iv) impact/shock heating. Quantitative modeling of the first two processes suggests that neither of them could produce isotopic variations large enough to explain the observed isotopic variations. Volatilization during parent-body thermal metamorphism [the scenario (ii)], which has been commonly invoked to explain the isotopic variations of volatile elements, is gas transport-limited and its effect on isotopic fractionations of moderately volatile elements should be negligible. Modeling of diffusion processes suggests that (iii) could produce K isotopic variation comparable to the observed variation. The large isotopic variations in non-carbonaceous meteorites are thus most likely due to diffusive redistribution of K and Rb during metamorphism and/or shock-induced heating and vaporization.

58 GEOSCIENCES↗
DOE OSTI2023

A hybrid data-driven and model-based approach for computationally efficient stochastic unit commitment and economic dispatch under wind and solar uncertainty

Stochastic unit commitment (UC) and economic dispatch (ED) are imperative in dealing with uncertainty in renewable forecast for power system operation and planning such that the overall expected production cost is minimized over the planning horizon. However, accurate calculation of the expected production cost requires assessment of a very large number of different scenarios of uncertain renewable resources, such as solar and wind, which is practically infeasible to simulate in real time. This article proposes a hybrid datadriven and physics-based model-predictive paradigm to efficiently solve for stochastic unit commitment and economic dispatch considering uncertainty in wind and solar power forecasts. Here, the novelty of the approach lies in decoupling the production cost estimation from the unit commitment and economic dispatch optimization problems under uncertainty without compromising on the fidelity of the solutions. A data-driven machine learning model is first developed to predict the mean optimal production cost. A physics-based inverse problem is then solved to get the stochastic UC and ED profiles from the expected cost. The presented approach considers, for the first time, solar uncertainty in UC/ED determination and enables efficient and accurate propagation of wind and solar uncertainty to estimate the statistics of the production cost. The effectiveness of the developed approach is demonstrated systematically on a stylized RTS-GMLC single-node system. The overall framework predicts the expected cost 62.5% more accurately than the existing state-of-the-art, on unforeseen days during the entire year, and yields, for the first time, the associated physically consistent UC and ED profiles. The solutions are also shown to be flexible in providing adequate daily reserves to address any statistical deviations from probabilistic power forecasts. The computational time associated with the presented method is only about 10 s compared to over 24 h needed for a conventional stochastic UC/ED determination under uncertainty on an Intel Core i9 processor with 32 GB of RAM.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Pulsating one-dimensional detonation in ammonia-hydrogen–air mixtures

In this study, one-dimensional detonations in ammonia/hydrogen-air mixtures are numerically investigated by solving the fully compressible Navier-Stokes equations with detailed chemistry. Pulsating instabilities with single-mode are observed during the detonation wave propagation, accompanied by periodic coupling and decoupling of the lead shock wave and the reactive front. The ratio between driver pressure and initial pressure determines the overdrive degree and thus the oscillatory mode of detonation for a premixture with certain composition. The effects of hydrogen dilution and mixture equivalence ratio on pulsating detonations are also examined under a constant driver pressure. The growing hydrogen fraction in fuel blends significantly increases the oscillation frequency. In addition, the pulsating detonation frequency rises with increasing equivalence ratio under fuel-lean conditions, peaks under stoichiometric conditions, and falls under fuel-rich conditions as the equivalence ratio increases further. Evolutions of reactants, main intermediate radicals, and products are analysed in both fuel-lean and fuel-rich conditions. A chemical explosive mode analysis further confirms the highly-autoignitive nature of the mixture in the induction zone between reaction front and shock front where thermal diffusion plays a negligible role.

08 HYDROGEN↗
DOE OSTI2023

Bayesian optimal experimental design for constitutive model calibration

Computational simulation is increasingly relied upon for high/consequence engineering decisions, which necessitates a high confidence in the calibration of and predictions from complex material models. However, the calibration and validation of material models is often a discrete, multi-stage process that is decoupled from material characterization activities, which means the data collected does not always align with the data that is needed. To address this issue, an integrated workflow for delivering an enhanced characterization and calibration procedure—Interlaced Characterization and Calibration (ICC)—is introduced and demonstrated. Further, this framework leverages Bayesian optimal experimental design (BOED), which creates a line of communication between model calibration needs and data collection capabilities in order to optimize the information content gathered from the experiments for model calibration. Eventually, the ICC framework will be used in quasi real-time to actively control experiments of complex specimens for the calibration of a high-fidelity material model. This work presents the critical first piece of algorithm development and a demonstration in determining the optimal load path of a cruciform specimen with simulated data. Calibration results, obtained via Bayesian inference, from the integrated ICC approach are compared to calibrations performed by choosing the load path a priori based on human intuition, as is traditionally done. The calibration results are communicated through parameter uncertainties which are propagated to the model output space (i.e. stress–strain). In these exemplar problems, data generated within the ICC framework resulted in calibrated model parameters with reduced measures of uncertainty compared to the traditional approaches.

42 ENGINEERING↗