Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reference 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 1,225 records · Page 68

DFT Calculation of Thermal Expansion within Debye-Grüneisen Framework Made Simple and Accurate

The density functional theory (DFT) based Debye-Grüneisen approach is a computational method used to calculate the thermodynamic properties of materials by relating the volume dependence of the Debye frequency to the Grüneisen parameter. In this study, we propose a new model to describe the volume dependence of the Grüneisen parameter. This model offers a simple yet accurate approach to predict the thermal expansion of materials within the framework of the Debye-Grüneisen theory. Specifically, for pure elements or species used as a reference system, we calibrated the parameters to reproduce the experimental coefficient of thermal expansion, heat capacity, and entropy. For compounds with arbitrary compositions, the parameters are determined by taking the arithmetic or geometric average of those obtained from the pure elements or species based on their composition. We then compared the model prediction with experiments. This allowed for a comprehensive characterization of the thermal expansion behavior of various materials.

Wang, Yi↗

A bi-level advanced control framework for large-scale control of buildings with system-level impact

Increased electricity consumption combined with new forms of generation is testing the reliability of our grid infrastructure. This work describes a method to improve the reliability of the grid through large-scale advanced building control. This paper develops a bi-level distributed control framework to shift the load of 153 buildings to achieve a system-level objective of tracking a power reference signal. This bi-level control is based on the previously-developed ANPV-MPC, a predictive controller that uses a Bayesian neural network to generate an accurate control model and adapt to changing conditions over time. By shifting the building electricity demand to better match the available power, the grid system supplying the buildings is more reliable as evidenced by the analysis of node voltages across an IEEE 13-bus distribution system. The proposed bi-level control framework tracks the system-level power reference with enough accuracy to regulate node voltages across the IEEE 13-bus distribution system within ANSI limits of ±5%. Additionally, the adaptive nature of ANPV-MPC allows each building across the system to adapt to changing conditions, further amplifying the system-level reliability.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Unified Model of Current-Limiting Grid-Forming Inverters for Large-Signal Analysis

Grid-forming (GFM) inverters can hardly withstand any overloading. As such, GFM inverters need a current limiter in their control system to avoid hardware damage during overloading or grid disturbances. When the current limiter is engaged for inverter protection, however, it dynamically changes the inverter output impedance, which can significantly affect the system behavior against large disturbances. For example, when the current limiter curtails the output power, the primary controller of the GFM inverter-which generates the voltage and angle reference-is prone to losing synchronization with the grid. This can cause large-signal instability in the network. This paper presents a unified GFM current-limiter model to gain a deeper understanding of the impact of the GFM inverter current limiting on large-signal instability and other system behaviors. The unified model captures several well-known current limiters for GFM inverters, and it reveals their large-signal equivalence characterized by the limiter angle. We use the unified model to present in-depth insights into the driving forces that cause large-signal instability in GFM inverter-driven grids. Further, this paper illustrates that the unified model can be used as a tool for assessing the system-level impacts of various GFM current-limiter designs during grid disturbances. Using numerical, analytical, electromagnetic transient models, and a hardware setup, the validity of the unified model is evaluated.

24 POWER TRANSMISSION AND DISTRIBUTION↗

CFD Simulations of Lower Plenum Mixing

Review of model development and validation performed in the Advanced Reactor Technologies (ART) program for thermal mixing at the outlet of High Temperature Gas Reactors (HTGRs). Understanding the mixing that occurs in the lower plenum in an HTGR is necessary to facilitate design improvements and to perform reactor safety analysis. Numerical models are one possible approach to gain a better understanding of mixing in the lower plenum. Given the complexity of the geometry and the intense mixing present, it is important to perform validation of numerical models. Three models have been developed during FY2025: a porous media with Pronghorn, a Reynolds Averaged Navier Stokes (RANS) with STAR-CCM+, and a Large Eddy Simulation (LES) with NekRS. The reference facility is a scaled-down version of the lower plenum of the High Temperature Gas-Cooled Reactor - Pebble-bed Module (HTR-PM) demonstration reactor. Preliminary results of the porous media and the RANS shows general good agreement against experimental benchmark data. Future work will leverage high-fidelity results obtained through LES to guide model selection and improvements to the lower-fidelity models, with particular attention to the Pronghorn porous media.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

A data-driven framework for permeability prediction of natural porous rocks via microstructural characterization and pore-scale simulation

Understanding the microstructure–property relationships of porous media is of great practical significance, based on which macroscopic physical properties can be directly derived from measurable microstructural informatics. However, establishing reliable microstructure–property mappings in an explicit manner is difficult, due to the intricacy, stochasticity, and heterogeneity of porous microstructures. In this paper, a data-driven computational framework is presented to investigate the inherent microstructure–permeability linkage for natural porous rocks, where multiple techniques are integrated together, including microscopy imaging, stochastic reconstruction, microstructural characterization, pore-scale simulation, feature selection, and data-driven modeling. A large number of 3D digital rocks with a wide porosity range are acquired from microscopy imaging and stochastic reconstruction techniques. A broad variety of morphological descriptors are used to quantitatively characterize pore microstructures from different perspectives, and they compose the raw feature pool for feature selection. Here high-fidelity lattice Boltzmann simulations are conducted to resolve fluid flow passing through porous media, from which reliable permeability references are obtained. The optimal feature set that best represents permeability is identified through a performance-oriented feature selection process, upon which a cost-effective surrogate model is rapidly fitted to approximate the microstructure-permeability mapping via data-driven modeling. This surrogate model exhibits great advantages over empirical/analytical formulas in terms of prediction accuracy and generalization capacity, which can predict reliable permeability values spanning four orders of magnitude. Besides, feature selection also greatly enhances the interpretability of the data-driven prediction model, from which new insights into the mechanism of how microstructural characteristics determine intrinsic permeability are obtained.

58 GEOSCIENCES↗

Temporal regulation of cold transcriptional response in switchgrass

Switchgrass low-land ecotypes have significantly higher biomass but lower cold tolerance compared to up-land ecotypes. Understanding the molecular mechanisms underlying cold response, including the ones at transcriptional level, can contribute to improving tolerance of high-yield switchgrass under chilling and freezing environmental conditions. Here, by analyzing an existing switchgrass transcriptome dataset, the temporal cis- regulatory basis of switchgrass transcriptional response to cold is dissected computationally. We found that the number of cold-responsive genes and enriched Gene Ontology terms increased as duration of cold treatment increased from 30 min to 24 hours, suggesting an amplified response/cascading effect in cold-responsive gene expression. To identify genomic sequences likely important for regulating cold response, machine learning models predictive of cold response were established using k -mer sequences enriched in the genic and flanking regions of cold-responsive genes but not non-responsive genes. These k -mers, referred to as putative cis -regulatory elements (pCREs) are likely regulatory sequences of cold response in switchgrass. There are in total 655 pCREs where 54 are important in all cold treatment time points. Consistent with this, eight of 35 known cold-responsive CREs were similar to top-ranked pCREs in the models and only these eight were important for predicting temporal cold response. More importantly, most of the top-ranked pCREs were novel sequences in cold regulation. Our findings suggest additional sequence elements important for cold-responsive regulation previously not known that warrant further studies.

60 APPLIED LIFE SCIENCES↗

Demonstrating SolarPILOT’s Python API Through Heliostat Optimal Aimpoint Strategy Use Case

SolarPILOT is a software package that generates solar field layouts and characterizes the optical performance of concentrating solar power (CSP) tower systems. SolarPILOT was developed by the National Renewable Energy Laboratory (NREL) as a stand-alone desktop application but has also been incorporated into NREL’s1 System Advisor Model (SAM) in a simplified format. Prior means for user interaction with SolarPILOT have included the application’s graphical interface, the SAM routines with limited configurability, and through a built-in scripting language called “LK.” This paper presents a new, full-featured, Python-based application programmable interface (API) for SolarPILOT, which we hereafter refer to as CoPylot. CoPylot provides access to all SolarPILOT’s capabilities to generate and characterize power tower CSP systems seamlessly through Python. Supported capabilities include (i) creating and destroying a model instance with message reporting tools; (ii) accessing and setting any SolarPILOT variable including custom land boundaries for field layouts; (iii) programmatically managing receiver and heliostat objects with varied attributes for systems with multiple receiver or heliostat types; (iv) generating, assigning, and modifying solar field layouts including the ability to set individual heliostat locations, aimpoints, soiling rates, and reflectivity levels; (v) simulating solar field performance; (vi) returning detailed results describing performance of individual heliostats, the aggregate field, and receiver flux distribution; and, (vii) exporting Python-based model instances to multiple file formats. CoPylot enables Python users to perform detailed CSP tower analysis utilizing either the Hermite expansion technique (analytical) or the SolTrace ray-tracing engine. In addition to CoPylot’s functionality, Python users have access to the over 100,000 open-source libraries to develop, analyze, optimize, and visualize power tower CSP research. This enables CSP researchers to perform analysis that was previously not possible through SolarPILOT’s existing interfaces. This paper discusses the capabilities of CoPylot and presents a use case wherein we demonstrate optimal solar field aiming strategies.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Demonstrating SolarPILOT's Python API Through Heliostat Optimal Aimpoint Strategy Use Case: Preprint

SolarPILOT is a software package that generates heliostat field layouts and characterizes the optical performance of concentrating solar power (CSP) tower systems. SolarPILOT was developed by the National Renewable Energy Laboratory (NREL) as a stand-alone desktop application but has also been incorporated into NREL's System Advisor Model (SAM) in a simplified format. Prior means for user interaction with SolarPILOT have included the application's graphical interface, the SAM routines with limited configurability, and through a built-in scripting language called "LK." This paper presents a new, full-featured Python-based application programmable interface (API) for SolarPILOT, which we hereafter refer to as CoPylot. CoPylot provides access to all SolarPILOT's capabilities to generate and characterize power tower CSP systems seamlessly through Python. Supported capabilities include (i) creating and destroying a model instances with message reporting tools; (ii) accessing and setting any SolarPILOT variable including custom land boundaries for field layout; (iii) programmatically managing receiver and heliostat objects with varied attributes for systems with multiple receiver or heliostat types; (iv) generating, assigning, and modifying heliostat field layouts including the ability to set individual heliostat locations, aimpoints, soiling rates, and reflectivity levels; (v) simulating heliostat field performance; (vi) returning detailed results describing performance of individual heliostats, the aggregate field, and receiver flux; and, (vii) exporting Python-based model instances to multiple file formats. CoPylot enables Python users to perform detailed tower CSP analysis utilizing either the Hermite expansion technique (analytical) or the SolTrace ray-tracing engine. In addition to CoPylot's functionality, Python users have access to the over 100,000 open-source libraries to develop, analyze, optimize, and visualize CSP tower research.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Eulerian thermo-mechanical simulations of heterogeneous solid propellants using an approximate projection method

We present a mathematical model that describes thermo-mechanical deformations and thermal gradients on the unsteady burning of a heterogeneous solid propellant. A scaling study shows that the deformations in solid at combustion timescales can be treated as quasi-static. The resulting thermo-mechanical formulation is formulated on a Cartesian grid and makes use of a weak form of Chorin-type projection method to deal with large difference in shear modulus of constituent materials. Further, a one-dimensional verification study is carried out by comparing numerical simulations with those of an analytical model. In addition, convergence studies for a two-dimensional propellant sandwich configuration are presented for the stress, velocity, and reference map components. Finally, simulations are carried out for a two-dimensional random propellant pack and the time-averaged burn rate and strains are reported.

42 ENGINEERING↗

North American energy system responses to natural gas price shocks

As of 2020, North American natural gas extraction and use in the electricity sector have both reached all-time highs. The combination of North America's increased reliance on natural gas with a potential disruption to the natural gas market has several energy security implications. Additionally, policymakers interested in economic resiliency will find this study's results useful for informing the implications of the energy sectors' long-term planning and investment decisions. This paper evaluates how both the electricity and natural gas sectors could respond to hypothetical gas price shocks under different system configurations. We impose unforeseen natural gas price shocks under reference and alternative configurations resulting from a renewable generation mandate or variations to renewable capacity costs. Results from several different models are presented for the electricity and natural gas sectors separately for Canada, Mexico, and the United States. Generally, the US becomes more (less) reliant on electricity imports from Canada given a high (low) gas price shock but increases (decreases) exports to Mexico. The renewable mandate is shown to buffer electricity price increases under high price shocks but price reductions under the low price shocks are dampened given less flexibility to take advantage of the low-priced natural gas. The United States is demonstrated to reduce natural gas production and net exports with high natural gas price shocks given a reduction in demand.

03 NATURAL GAS↗

Geothermal Play Fairway Analysis, Part 2: GIS methodology

Play Fairway Analysis (PFA) in geothermal exploration originates from a systematic methodology developed within the petroleum industry and is based on a geologic, geophysical, and hydrologic framework of identified geothermal systems. We tailored this methodology to study the geothermal resource potential of the Snake River Plain and surrounding region, but it can be adapted to other geothermal resource settings. We adapted the PFA approach to geothermal resource exploration by cataloging the critical elements controlling exploitable hydrothermal systems, establishing risk matrices that evaluate these elements in terms of both probability of success and level of knowledge, and building a code-based ‘processing model’ to process results. A geographic information system was used to compile a range of different data types, which we refer to as elements (e.g., faults, vents, heat flow, etc.), with distinct characteristics and measures of confidence. Discontinuous discrete data (points, lines, or polygons) for each element were transformed into continuous interpretive 2D grid surfaces called evidence layers. Because different data types have varying uncertainties, most evidence layers have an accompanying confidence layer which reflects spatial variations in these uncertainties. Risk layers, as defined here, are the product of evidence and confidence layers, and are the building blocks used to construct Common Risk Segment (CRS) maps for heat, permeability, and seal, using a weighted sum for permeability and heat, but a different approach with seal. CRS maps quantify the variable risk associated with each of these critical components. In a final step, the three CRS maps were combined into a Composite Common Risk Segment (CCRS) map, using a modified weighted sum, for results that reveal favorable areas for geothermal exploration. Additional maps are also presented that do not mix contributions from evidence and confidence (to allow an isolated view of evidence and confidence), as well as maps that calculate favorability using the product of components instead of a weighted sum (to highlight where all components are present). Our approach helped to identify areas of high geothermal favorability in the western and central Snake River Plain during the first phase of study and helped identify more precise local drilling targets during the second phase of work. By identifying favorable areas, this methodology can help to reduce uncertainty in geothermal energy exploration and development.

15 GEOTHERMAL ENERGY↗

Probability Density Function Control for Stochastic Nonlinear Systems using Monte Carlo Simulation

This paper presents a fast implementation framework of output probability density function (PDF) control for a class of stochastic nonlinear systems which are subjected to non-Gaussian noises. The statistical properties of the system outputs can be adjusted by shaping the dynamic output probability density function to track the reference stochastic distribution. However, the dynamic probability density function evolution is very difficult to obtain analytically even if the system model and the stochastic distributions of the noises are known. Motivated by Monte Carlo simulation, the dynamic probability density function can be estimated by sampling data which forms the contribution of this paper. In particular, the sampling points are generated following the stochastic distribution of the noise for each instant. These points go through the system would generate the histogram for system outputs, then the dynamic model can be established based on the dynamic histogram which reflects the randomness and the nonlinear dynamics of the investigated system. Based on the established model, the output probability density function tracking can be achieved and the simulation results and discussions show the effectiveness and benefits of the presented framework.

Zhang, Qichun↗

Informative and non-informative decomposition of turbulent flow fields

Not all the information in a turbulent field is relevant for understanding particular regions or variables in the flow. Here, we present a method for decomposing a source field into its informative Φ I (x, t) and residual Φ R (x, t) components relative to another target field. The method is referred to as informative and non-informative decomposition (IND). All the necessary information for physical understanding, reduced-order modelling and control of the target variable is contained in Φ I (x, t), whereas Φ R (x, t) offers no substantial utility in these contexts. The decomposition is formulated as an optimisation problem that seeks to maximise the time-lagged mutual information of the informative component with the target variable while minimising the mutual information with the residual component. The method is applied to extract the informative and residual components of the velocity field in a turbulent channel flow, using the wall shear stress as the target variable. We demonstrate the utility of IND in three scenarios: (i) physical insight into the effect of the velocity fluctuations on the wall shear stress; (ii) prediction of the wall shear stress using velocities far from the wall; and (iii) development of control strategies for drag reduction in a turbulent channel flow using opposition control. In case (i), IND reveals that the informative velocity related to wall shear stress consists of wall-attached high- and low-velocity streaks, collocated with regions of vertical motions and weak spanwise velocity. This informative structure is embedded within a larger-scale streak–roll structure of residual velocity, which bears no information about the wall shear stress. In case (ii), the best-performing model for predicting wall shear stress is a convolutional neural network that uses the informative component of the velocity as input, while the residual velocity component provides no predictive capabilities. Finally, in case (iii), we demonstrate that the informative component of the wall-normal velocity is closely linked to the observability of the target variable and holds the essential information needed to develop successful control strategies.

97 MATHEMATICS AND COMPUTING↗

Quantum Orbital Minimization Method for Excited States Calculation on a Quantum Computer

Herein we propose a quantum-classical hybrid variational algorithm, the quantum orbital minimization method (qOMM), for obtaining the ground state and low-lying excited states of a Hermitian operator. Given parametrized ansatz circuits representing eigenstates, qOMM implements quantum circuits to represent the objective function in the orbital minimization method and adopts a classical optimizer to minimize the objective function with respect to the parameters in ansatz circuits. The objective function has an orthogonality constraint implicitly embedded, which allows qOMM to apply a different ansatz circuit to each input reference state. We carry out numerical simulations that seek to find excited states of H 2 , LiH, and a toy model consisting of four hydrogen atoms arranged in a square lattice in the STO-3G basis with UCCSD ansatz circuits. Comparing the numerical results with existing excited states methods, qOMM is less prone to getting stuck in local minima and can achieve convergence with more shallow ansatz circuits.

97 MATHEMATICS AND COMPUTING↗

Development of a Transferable Density-Functional Tight-Binding Model for Organic Molecules at the Water/Platinum Interface

A computationally efficient and transferable approach for modeling reactions at metal/water interfaces could significantly accelerate our understanding and ultimately the development of new catalytic transformations, particularly in the context of the emerging field of biomass conversion. Here, we present a parametrization of Pt–X (X = H, O, C) density-functional tight-binding (DFTB) for addressing this need. We first constructed Pt–H, Pt–O, and Pt–C repulsive potential splines. These pairwise parameters were then augmented to include many-body interactions using the Chebyshev Interaction Model for Efficient Simulation (ChIMES). We compare the geometrical and energetic performances of both DFTB and DFTB/ChIMES methods with DFT reference data across a variety of organic molecules at the platinum surface from nanoparticles to single-crystal surfaces. DFTB shows limited transferability between extended crystal surfaces and small nanoparticles. This transferability is significantly improved through the introduction of three-body interactions with Pt in DFTB/ChIMES, which provides consistent results across various systems, with reductions in the RMSD from around 30 kcal/mol in DFTB to around 10 kcal/mol. We demonstrate the stability and reliability of the obtained parameters by performing metadynamic simulations for the adsorption of phenol on Pt(111). We observe that DFTB itself is undersolvating the surface, leading to only one or two chemisorbed water molecules in a c(4 × 6) unit cell. In contrast, DFTB/ChIMES leads to a coverage of about 0.5 ML and successfully captures the chemisorbed mode of phenol at both the solid/liquid and the solid/gas interfaces. Furthermore, in agreement with experimental measurements, the adsorption at the solid/liquid interface is significantly weaker than that at the solid/gas interface. As a result, we highlight that even with DFTB, where we can accumulate dynamics for more than 1 ns for a given system, the simulations are not fully converged.

Adsorption↗

A Large Ensemble Global Dataset for Climate Impact Assessments

We present a self-consistent, large ensemble, high-resolution global dataset of long-term future climate, which accounts for the uncertainty in climate system response to anthropogenic emissions of greenhouse gases and in geographical patterns of climate change. The dataset is developed by applying an integrated spatial disaggregation (SD) - bias-correction (BC) method to climate projections from the MIT Integrated Global System Model (IGSM). Four emission scenarios are considered that represent energy and environmental policies and commitments of potential future pathways, namely, Reference, Paris Forever, Paris 2 °C and Paris 1.5 °C. The dataset contains nine key meteorological variables on a monthly scale from 2021 to 2100 at a spatial resolution of 0.5°x 0.5°, including precipitation, air temperature (mean, minimum and maximum), near-surface wind speed, shortwave and longwave radiation, specific humidity, and relative humidity. We demonstrate the dataset’s ability to represent climate-change responses across various regions of the globe. This dataset can be used to support regional-scale climate-related impact assessments of risk across different applications that include hydropower, water resources, ecosystem, agriculture, and sustainable development.

54 ENVIRONMENTAL SCIENCES↗

Cost Analysis of Annual and Monthly Net Zero Energy Performance for Multifamily Buildings in the United States

As the United States continues its progress toward sustainable construction, net zero building design is becoming an increasingly important and popular topic. Two definitions of net zero energy performance, monthly and annual, provide different levels of energy autonomy and efficiency within a building. This article analyzes the viability and incremental cost for two-, three-, and four-story multifamily apartment buildings to reach both annual and monthly net zero energy performance throughout four climate zones in the United States using baseline reference buildings that represent current construction practices. Building size plays a large role in determining the capability for a building model to reach annual or monthly net zero. Two-story buildings are capable of reaching annual net zero with an increase in construction cost of about 4.4–5.6%. Three-story buildings in warmer climates can achieve annual net zero with an increase of about 5.1%, but models in colder climates cannot reach annual net zero performance. Four-story buildings cannot achieve net zero, owing to insufficient building area for solar arrays to produce electricity. Monthly net zero energy performance was significantly harder to achieve, with the only case reaching this goal being the two-story Houston based model, at an added cost of 8.2%. Generally speaking, only small buildings in warm climates will be able to achieve monthly net zero without vastly oversizing photovoltaic systems and increasing costs without adequate payback.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Helium metastable density determination in the COST reference source by absolutely calibrated optical emission spectroscopy

Helium metastable densities in the COST Reference Microplasma Jet are estimated for a variety of He/N2 admixtures and dissipated powers by applying a collisional-radiative model to absolutely calibrated optical emission spectroscopy measurements. This is accomplished by delineating the excitation mechanisms that result in the N2(C–B) and N2+(B–X) emission bands, the latter of which is strongly coupled to the presence of helium metastables. A number of other plasma parameters are established and discussed for each operating condition including the electron energy distribution function, reduced electric field, rate constants, and electron density. With these parameters, the reaction rates for the primary ionization pathways are also calculated, emphasizing the importance of helium metastables for discharge sustainment. Good agreement with the existing literature is found for most plasma parameters and for helium metastable densities, in particular. A clear [N2]−1 relationship between the nitrogen concentration and density of helium metastables is demonstrated, as has been identified in previous studies in analogous atmospheric pressure plasma jets. This validates the efficacy of this optical technique for determining helium metastable densities and establishes it as a viable, and in many cases, more accessible alternative to other means of quantifying helium metastables in low-temperature plasmas.

Physics↗