Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal system modeling”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

System efficiency of packed bed TES with radial flow vs. axial flow – Influence of aspect ratio

Here, this paper compares the net system efficiency, including thermal efficiency and pressure drop effects, of radial versus axial flow packed beds for thermal energy storage. The traditional packed bed system is a cylindrical geometry where fluid flows axially from one end to another. However, issues of thermal stratification and high-pressure drop have led to recent studies on radial flow systems. One potential benefit is the reduced pressure drop in a radial flow system. This paper compares the performance of radial flow and axial flow systems at a range of aspect ratios (AR = H/D bed ) from 0.21 to 1.92 using a numerical model where the storage volume is held constant in all cases. When the radial flow bed is at a low aspect ratio (short/wide), the thermal front is improved but the pressure drop is high. At a high aspect ratio, the velocity is reduced in radial flow, leading to decreased pressure drop but an increased spreads in the thermal front that lowers thermal efficiency. The opposite trends are noted in axial flow. Thermal efficiencies of 83–91 % were noted for radial flow, while they ranged from 85 to 94 % in axial flow. Net efficiencies including pressure drop ranged from 74 to 82 % for radial flow and 80–87 % for axial flow. In both systems, a peak net efficiency was noted between the highest and lowest aspect ratio. While some aspect ratios with radial flow outperform axial flow from a net efficiency perspective, the results show that the highest net efficiency from axial flow is higher than that from radial flow. Overall, this paper highlights the importance of innovative TES designs and their potential to improve energy efficiency.

25 ENERGY STORAGE↗
DOE OSTI2026

Achieving 20.1% Efficiency in Organic Solar Cells Through Interconnected Fibrillar Networks via Local Molecular Stacking

ABSTRACT The performance of organic solar cells (OSCs) is governed by how molecular packing evolves into interconnected networks that facilitate exciton dissociation and charge transport. Using an all‐small‐molecule blend DR3TSBDT:Y6 as a model system, we study how local molecular stacking evolves into performance‐relevant morphology during solvent vapor annealing (SVA) and subsequent thermal annealing (TA). SVA promotes end‐to‐end stacking of amorphous acceptors to form interconnected fibrils, while TA compacts inter‐fibril spacing without disrupting favorable local order. Such molecular‐to‐morphological refinements broaden light absorption, enhance charge transport, and markedly improve device efficiency. Extending this approach to additional blend systems (D18:Y6, D18:L8‐BO, and DR3TSBDT:L8‐BO) yields similar structural evolution and performance gains, with the D18:L8‐BO system achieving up to 20.10% PCE. Our study establishes control over local stacking in amorphous acceptors into fibrillar networks as a general and effective route to realize high‐performance OSCs.

Wu, Junying↗
DOE OSTI2020

On the prediction of critical heat flux using a physics-informed machine learning-aided framework

The critical heat flux (CHF) corresponding to the departure from nucleate boiling (DNB) crisis is essential to the design and safety of a two-phase flow boiling system. Despite the abundance of predictive tools available to the thermal engineering community, the path for an accurate, robust CHF model remains elusive due to lack of consensus on the DNB triggering mechanism. This work aims to apply a physics-informed machine learning (ML)-aided hybrid framework to achieve superior predictive capabilities. Such a hybrid approach takes advantage of existing understanding in the field of interest (i.e., domain knowledge) and uses ML to capture undiscovered information from the mismatch between the actual and domain knowledge-predicted target. A detailed case study is carried out with an extensive DNB-specific CHF database to demonstrate (1) the improved performance of the hybrid approach as compared to traditional domain knowledge-based models, and (2) the hybrid model’s superior generalization capabilities over standalone ML methods across a wide range of flow conditions. The hybrid framework could also readily extend its applicability domain and complexity on the fly, showing an elevated level of flexibility and robustness. Based on the case study conclusions, the window-type extrapolation mapping methodology is further proposed to better inform high-cost experimental work.

42 ENGINEERING↗
DOE OSTI2022

Technoeconomic Cost Analysis of NREL Concentrating Solar Power Gen3 Liquid Pathway

The National Renewable Energy Laboratory is leading the liquid (molten salt) power tower pathway for the U.S. Department of Energy's concentrating solar power Gen3 initiative. The Gen3 liquid pathway required updated designs to three major components: the tower and receiver, the thermal energy storage tanks, and the power cycle. We assume a 100 MWe net system output and used the System Advisor Model (SAM) to complete a technoeconomic cost analysis of the Gen3 liquid pathway design and estimate its levelized cost of electricity. This paper summarizes the methodology and results of that analysis.

chloride salt↗
DOE OSTI2025

Multi-device analysis of energy loss duration and pellet penetration with implications for shattered pellet injection in ITER

A robust disruption mitigation system (DMS) requires accurate characterization of key disruption timescales, one of the most notable being the thermal quench (TQ). Recent modeling of shattered pellet injection (SPI) into ITER plasmas, using JOREK and INDEX, suggests long TQ durations (6–10 ms) and slow cold front propagation due to the large plasma size. If validated, these predictions would have an impact on the desired pellet parameters and mitigation strategies for the ITER DMS. To resolve these questions, a database of SPI experiments from several small-to-large sized devices (J-TEXT, KSTAR, AUG, DIII-D, and JET) has been compiled under the auspices of the International Tokamak Physics Activity MHD, disruptions, and control topical group. Analysis of the energy loss duration (proxy for the TQ duration) with machine size is presented for both mixed neon/deuterium (Ne/D) SPI and pure deuterium (D) SPI. Several metrics for the energy loss onset (e.g. soft x-ray signal drop, I p dip, and radiation flash) were considered as the conventional metric, electron cyclotron emission, is often cut-off during SPI. Several scalings with different onset metrics showed an increase in energy loss duration with machine size. The energy loss duration was additionally shown to be a function of the ratio between the number of SPI neon atoms injected and the stored energy. Analysis of the pellet shard position relative to the cold front found that in larger devices, pellets are typically found inboard of the q = 2 surface at the energy loss onset. Lastly, the delay between the pellet shards hitting the q = 2 surface and the energy loss onset was additionally found to increase with machine size. This suggests that the pellet shards in large devices will penetrate faster and further than the cooling front.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Applying Infrared Thermography as a Method for Online Monitoring of Turbine Blade Coolant Flow

As gas turbine engine manufacturers strive to implement condition-based operation and maintenance, there is a need for blade monitoring strategies capable of early fault detection and root-cause determination. Given the importance of blade cooling flows to turbine blade health and longevity, there is a distinct lack of methodologies for coolant flowrate monitoring. The present study addresses this identified opportunity by applying an infrared thermography system on an engine-representative research turbine to generate data-driven models for prediction of blade coolant flowrate. Thermal images were used as inputs to a linear regression and regularization algorithm to relate blade surface temperature distribution with blade coolant flowrate. Additionally, this study investigates how coolant flowrate prediction accuracy is influenced by the number and breadth of diagnostic measurements. Here, the results of this study indicate that a source of high-fidelity training data can be used to predict blade coolant flowrate within about six percent error. Furthermore, identification of prioritized sensor placement supports application of this technique across multiple sensor technologies capable of measuring blade surface temperature in operating gas turbine engines, including spatially resolved and point-based measurement techniques.

42 ENGINEERING↗
DOE OSTI2022

Reconfiguring band-edge states and charge distribution of organic semiconductor–incorporated 2D perovskites via pressure gating

Two-dimensional (2D) semiconductor heterostructures are key building blocks for many electronic and optoelectronic devices. Reconfiguring the band-edge states and modulating their interplay with charge carriers at the interface in a continuous manner have long been sought yet are challenging. Here, using organic semiconductor–incorporated 2D halide perovskites as the model system, we realize the manipulation of band-edge states and charge distribution via mechanical—rather than chemical or thermal—regulation. Compression induces band-alignment switching and charge redistribution due to the different pressure responses of organic and inorganic building blocks, giving controllable emission properties of 2D perovskites. We propose and demonstrate a “pressure gating” strategy that enables the control of multiple emission states within a single material. We also reveal that band-alignment transition at the organic-inorganic interface is intrinsically not well resolved at room temperature owing to the thermally activated transfer and shuffling of band-edge carriers. This work provides important fundamental insights into the energetics and carrier dynamics of hybrid semiconductor heterostructures.

42 ENGINEERING↗
DOE OSTI2022

Deterministic Transport: An Introduction [Slides]

In this lecture, we cover deterministic particle transport. We focus on neutral particles (neutrons, gamma-rays, thermal photons). Our goal is to predictively and accurately model how free particles interact with their environments and properties of systems involving those particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Multi-fidelity information fusion and resource allocation

This project created and demonstrated a framework for the efficient and accurate prediction of complex systems with only a limited amount of highly trusted data. These next generation computational multi-fidelity tools fuse multiple information sources of varying cost and accuracy to reduce the computational and experimental resources needed for designing and assessing complex multi-physics/scale/component systems. These tools have already been used to substantially improve the computational efficiency of simulation aided modeling activities from assessing thermal battery performance to predicting material deformation. This report summarizes the work carried out during a two year LDRD project. Specifically we present our technical accomplishments; project outputs such as publications, presentations and professional leadership activities; and the project’s legacy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Deterministic Transport: An Introduction [Slides]

In this lecture, we cover deterministic particle transport. We focus on neutral particles (neutrons, gamma-rays, and thermal photons). Our goal is to predictively and accurately model how free particles interact with their environments and properties of systems involving those particles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

High Flux Isotope Reactor Low-Enriched Uranium High Density Silicide Fuel Preliminary Design Update: System Transient Analysis

As a part of conversion efforts from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel under direction of the National Nuclear Security Administration of the U.S. Department of Energy, multiple proposed designs of the High Flux Isotope Reactor (HFIR) have been created and assessed regarding reactor physics performance metrics, including designs utilizing uranium silicide dispersion fuel (U3Si2-Al). This report updates the previous analyses that evaluated the nuclear safety performance of LEU fuel designs with respect to selected accident events from the HFIR Safety Analysis Report (SAR). Both the Low Density (LD) and High Density (HD) Optimized designs’ reactivity initiated accident fuel performance improved relative to the HEU fuel, attributed to greater 238 U negative Doppler feedback. However, the thermal margins for primary coolant system accidents were reduced with some acceptance criteria unable to be met. The need to resolve reduced thermal margin, open modeling items, and unresolved assumptions was identified.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

LIM1TR: Lithium-Ion Modeling with 1-D Thermal Runaway (V.1.0)

LIM1TR (Lithium-Ion Modeling with 1-D Thermal Runaway) is an open-source code that uses the finite volume method to simulate heat transfer and chemical kinetics on a quasi 1-D domain. The target application of this software is to simulate thermal runaway in systems of lithium-ion batteries. The source code for LIM1TR can be found at https://github.com/ajkur/lim1tr. This user guide details the steps required to create and run simulations with LIM1TR starting with setting up the Python environment, generating an input file, and running a simulation. Additional details are provided on the output of LIM1TR as well as extending the code with custom reaction models. This user guide concludes with simple example analyses of common battery thermal runaway scenarios. The corresponding input files and processing scripts can be found in the “Examples” folder in the on-line repository, with select input files included in the appendix of this document.

25 ENERGY STORAGE↗
DOE OSTI2025

Oklo Sponsored Testing using PELICAN: System, Subchannel, and High-Fidelity Software Validation (Final CRADA Report)

Argonne National Laboratory (the Contractor), located in Lemont IL, and Oklo, Inc., (the Participant), headquartered in Santa Clara, CA, propose to enter into a Cooperative Research and Development Agreement (CRADA) to perform a gap analysis of thermal hydraulic data, perform prototypical fuel assembly pressure drop and cavitation model validation, generate the experimental data as well as the corresponding uncertainties for this matrix and, finally, develop the validation models with the Argonne system level code SAS4A/SASSYS-1, the Argonne subchannel analysis code DASSH, the Argonne high fidelity code Nek5000, and/or the Idaho National Laboratory code Pronghorn Subchannel. The work outlined below will significantly improve the experimental and validation database currently available for liquid metal fast reactors, thus making it a viable part of a comprehensive reactor design and licensing suite to be used by the participant.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Material Plasma Exposure eXperiment (MPEX) High Heat Flux Bellows for Component Alignment in High Microwave Environment

The Material Plasma Exposure eXperiment (MPEX) device is a linear plasma device developed to perform plasma material interaction experiments under the conditions prototypic of a fusion reactor divertor. MPEX has multiple systems that must be precisely aligned to the plasma axis, including an electron cyclotron heating system that emits up to 400 kW of microwave power into the vacuum vessel. Five distinct systems require precise alignment on the MPEX device, thus requiring four bellows, all of which are adjacent to the plasma at a relatively high heat flux of approximately 47 kW/m2 and microwave power regions. The MPEX high heat flux bellows (HHFB) is designed to deliver 6 degrees of freedom positioning. The HHFB includes titanium–zirconium–molybdenum (TZM) inserts that are brazed into a Glidcop AL-15 body using a high-temperature braze alloy, thus blocking direct line of sight to an edge-welded bellows from the plasma and microwave screen to block microwaves from the bellows. Further, a custom ConFlat knife edge is machined into the Glidcop AL-15, so the vacuum flanges do not need a braze or weld joint on the vacuum interface. Fingerstock or copper mesh is used to restrict microwave power from entering the interstitial space between the water-cooled Glidcop AL-15 body and edge-welded bellows. Glidcop AL-15 was selected as the material of choice for the water-cooled body because it can maintain mechanical integrity at elevated temperatures, and it also allows for a high-temperature braze. TZM was selected for its machinability and compatibility with vacuum and plasma requirements. A test article is also being considered for a similarly shaped component, the MPEX limiter, that will demonstrate the integrity of the braze joint under high thermal load. Results from this testing will be extrapolated to deduce the lifetime and integrity of the HHFB design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Thermal Overloading Risk Mitigation With a Semi-Analytical Probabilistic Model on Branch Current

A semi-analytical formulation is presented in this paper for the probability computation of branch current in multiphase systems. The developed formula is derived based on the linear power flow model in rectangular coordinates. The system uncertainty injections can be renewable energy resources or loads and are modeled using a Gaussian mixture model (GMM). The developed formula can be used to compute the line current violation probability as well as integrate into optimal power flow problem as chance-constraint relaxation. Here, the proposed formula is first compared with the Matlab embedded numerical integration function to show its performance. Besides, the semi-analytical formula is validated and compared with the Monte Carlo simulation method using the IEEE 123-bus system, EPRI Ckt5, and Ckt7 systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Spectral form factor of a quantum spin glass

It is widely expected that systems which fully thermalize are chaotic in the sense of exhibiting random-matrix statistics of their energy level spacings, whereas integrable systems exhibit Poissonian statistics. In this paper, we investigate a third class: spin glasses. These systems are partially chaotic but do not achieve full thermalization due to large free energy barriers. We examine the level spacing statistics of a canonical infinite-range quantum spin glass, the quantum p-spherical model, using an analytic path integral approach. We find statistics consistent with a direct sum of independent random matrices, and show that the number of such matrices is equal to the number of distinct metastable configurations — the exponential of the spin glass “complexity” as obtained from the quantum Thouless-Anderson-Palmer equations. We also consider the statistical properties of the complexity itself and identify a set of contributions to the path integral which suggest a Poissonian distribution for the number of metastable configurations. Our results show that level spacing statistics can probe the ergodicity-breaking in quantum spin glasses and provide a way to generalize the notion of spin glass complexity beyond models with a semi-classical limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Thermomechanical Modeling of Copper

The plastic deformation of metals is a dissipative process. Some fraction of the plastic work is converted to heat which, given the temperature dependent response of metals, produces a thermal-mechanical coupling. In various cases, for instance when the loading is dynamic, this interaction can impact the resulting response of a material and/or system. Thus, appropriately capturing the heat generation from plastic work is necessary for various solid mechanics analysis. Determination of the fraction of work converted to heat has been long studied. Recent developments have demonstrated that the fraction is not constant but depends on various state variables. Resolving these features requires combined modeling and experimental studies. To this end, 304L stainless steel – a poor thermal conductor – was recently subjected to such an investigation. Advanced modeling capabilities were deployed to assess novel thermomechanically coupled experiments. As a complement to that study, in the current work a similar investigation is performed on copper – a good thermal conductor – to assess performance on the opposite end of the spectrum. The current document discusses these modeling efforts.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Arbitrage and Capacity Firming in Coordination with Day-Ahead Bidding of a Hybrid PV Plant

A hybrid PV plant (HPP) combines a photovoltaic (PV) plant with a battery energy storage system (BESS), which is considered a promising step towards the future of renewable power plants by the U.S. Department of Energy. When the renewable penetration reaches a significant level, a hybrid PV plant can bid in as a controllable thermal plant in the future electricity market. In this study, a bidding and BESS scheduling model is proposed for the HPP. The robust optimization (RO) technique has been utilized to identify the worst-case scenario of uncertainties during the bidding process. To address the overly conservative issue of the single-stage RO, we have decoupled the BESS schedule for arbitrage and PV capacity firming by a two-stage RO formulation. By comparing the output of single-stage RO and two-stage RO, the two-stage RO bids and schedules in a more aggressive manner, which increases the income of HPP. Also, the penalty of under-generation is considered in our model so that the day-ahead bidding decision and arbitrage schedules can be adjusted based on the potential UNDER-GENERATION penalty. Because the proposed model is non-convex and contains multi-stages, the Column-and-Constraint Generation (C&CG) algorithm is applied to the model as the solution. The proposed model has shown better economic performance compared to a state-of-art single-stage bidding method in case studies.

BESS scheduling↗