PBNC Presentation: Dynamic Modeling of a Latent Heat Thermal Battery
Presentation detailing dynamic modeling progress on a heat pipe integrated thermal battery system.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
SEARCH · Engineering Papers
Results for “thermal dynamics”
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.
Computational Fluid Dynamics Modeling of Solar Thermal Dry Reforming of Methane in a Parabolic Trough
Computational fluid dynamics simulations of solar-thermal dry reforming of methane using a parabolic trough configuration were performed. Parametric simulations of different combinations of gas flow rate, receiver tube emissivity, and geometric concentration ratio were conducted to determine configurations that could achieve the required catalyst temperatures of at least 700 °C to achieve high conversion of CH4 and CO2 to H2 and CO. Results showed that the concentration ratio of the parabolic trough collector had to be increased from ~70 to ~120 and the receiver-tube emissivity had to be reduced to ~0.2 to achieve bulk average catalyst temperatures of greater than 700 °C. Lower gas flow rates also reduced enthalpic heat losses and increased catalyst temperatures.
Ho, Clifford↗
Impact of Different Thermal Gradients on the Dynamics of Cylindrical Lithium-ion Cells Subject to Accelerated Aging and on Module Performance
This study investigates the impacts of applying different thermal gradient patterns to cylindrical lithium-ion cells in a module on cell dynamics (temperatures, current flows, state of charge), module performance (evolution of resistance, capacity, and energy versus cycle number), and module lifetime. The thermal gradients were generated using cooling plates (CPs) with three different flow-field designs, namely, straight, perpendicular, and U-turn. The study uses computational fluid dynamics (CFD), the pseudo-two-dimensional (P2D) battery model, capacity loss and increased impedance due to the growth of a solid-electrolyte-interphase, and the electric current distribution from module terminals to cells that depends on the series-parallel electrical connections among the cells. The impact of the thermal gradient (resulting from the CP designs) on the variability in resistance, current, state of charge, and voltage among the cells was analyzed and linked to differences in the module's performance. Applying a thermal gradient to parallel-connected strings of series-connected cells led to variation in the current through each parallel string and an imbalance in the voltage of series-connected cells. Module performance is poorer when the thermal gradient causes a voltage imbalance than when it causes a current imbalance. Module performance becomes the worst when both current variation and voltage imbalance happen together. For instance, the module's lifetime (estimated as reaching 80% of its initial capacity) varied by 5% to 17.5%, depending on the magnitude and pattern of the imposed thermal gradient. As the relative orientation between thermal gradients and cells' electrical connectivity influences the module's performance, appropriate consideration should be given to the choice of the CP, especially if large thermal gradients are allowed.
Battery thermal management↗
Prediction of Bi 2 Te 3 -Sb 2 Te 3 Interfacial Conductance and Superlattice Thermal Conductivity Using Molecular Dynamics Simulations
Bismuth telluride (Bi 2 Te 3 ) and its alloys with antimony telluride (Sb 2 Te 3 ) have long been considered to be the best room-temperature bulk thermoelectric (TE) materials. In recent decades, proof-of-concept demonstrations on Bi 2 Te 3 -Sb 2 Te 3 nanostructures have shown high TE performance due to reduction in lattice thermal conductivities. Particularly, ultra-low thermal conductivities have been observed in Bi 2 Te 3 -Sb 2 Te 3 1D superlattices, leading to thermoelectric figures of merit (ZT) as high as 2.4. In contrast, very few computational studies have been performed to provide insight into the phonon transport across these nanostructures. In this work, we use non-equilibrium molecular dynamics simulations with previously developed force fields to simulate thermal transport across Bi 2 Te 3 -Sb 2 Te 3 interfaces and superlattices. We first calculate the thermal conductance associated with a Bi 2 Te 3 -Sb 2 Te 3 interface across a temperature range of 200–400 K. Furthermore, the values are also compared with thermal conductances calculated by a modified Landauer transport formalism using phonon transmission coefficients obtained from the diffuse mismatch model. Our results show that inelastic scattering processes contribute to an increase in interfacial thermal conductance at higher temperatures. Finally, we calculate the thermal conductivities of Bi 2 Te 3 -Sb 2 Te 3 superlattices with varying period lengths from 2 to 18 nm. A minimum thermal conductivity of 0.27 W/mK is observed at a period length of 4 nm, which is attributed to the competition between incoherent and coherent phonon transport regimes. In comparison with previous experimental measurements in the literature, our results show good agreement with respect to the range of thermal conductivity values and the period length corresponding to the minimum superlattice thermal conductivity.
42 ENGINEERING↗
A molecular dynamics simulation study of thermal conductivity anisotropy in β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (β-HMX)
Molecular dynamics (MD) simulations were used to predict the thermal conductivity of β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (β-HMX) along directions normal to the (011), (110), and (010) crystal planes. These directions were selected based on the measured morphological importance of the corresponding crystal surfaces. A reverse non-equilibrium MD approach was used wherein a constant heat flux is imposed along a prescribed direction and the resulting steady-state temperature gradient determined. The coefficient of thermal conductivity λ is the quotient of heat flux and temperature gradient (i.e. Fourier's law). Finite-size effects and sensitivity to imposed heat flux were investigated. The results reveal a modest dependence of the conductivity on crystal orientation, significant finite-size effects, and low sensitivity to imposed flux so long as the Fourier's law analysis is limited to the spatial interval in the simulation cell for which the temperature gradient is constant. Infinite-length thermal conductivities were estimated by linear regression of λ -1 (L) versus reciprocal cell length L -1 for each direction. The predicted values are systematically larger, but within a factor of two, than most published experimental determinations, the latter of which were obtained for pressed-powder or composite samples rather than oriented single crystals.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A Dynamic Model of Refrigerator with Thermal Energy Storage for Demand Flexibility
A refrigerator is an essential appliance for both commercial and residential applications. Thermal (Cold) energy storage can effectively reduce the carbon emission from refrigerators. This preliminary study investigated the demand flexibility of a domestic refrigerator using thermal (cold) energy storage through dynamic modeling and controls. The models include both fresh and frozen compartments, a vapor compression cycle loop, and thermal storage panels. Results show that thermal energy storage can effectively achieve temperature control of the compartments within the 4-hour time window. The load reduction and cost savings (46.3%) were achieved through the TES using TOU rates.
Gao, Zhiming↗
Quantum turbulence, superfluidity, non-Markovian dynamics, and wave function thermalization
While quantum turbulence has been addressed both experimentally (predominantly for superfluid He 4 and He 3 ) and theoretically, the dynamics of various ensembles of quantized vortices has been followed in time only until the vortices have decayed into phonons. How this “thermalization” is achieved is still an unaddressed and thus an unelucidated question. The unitary Fermi gas (UFG) is a unique quantum system, which has no classical counterpart and is of relevance to neutron stars, cold atoms, condensed-matter and nuclear many-body systems. The non-Markovian evolution of an isolated UFG is put in evidence and its entire nonequilibrium evolution can be studied theoretically within a unified theoretical framework. The initial lattice of quantum vortices and antivortices evolves through a couple of vortex tangles and excitation of Kelvin waves, where vortices cross and reconnect, until very slowly thermalization sets in. Published by the American Physical Society 2024
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Gold–Thiolate Nanocluster Dynamics and Intercluster Reactions Enabled by a Machine Learned Interatomic Potential
Monolayer protected metal clusters comprise a rich class of molecular systems and are promising candidate materials for a variety of applications. While a growing number of protected nanoclusters have been synthesized and characterized in crystalline forms, their dynamical behavior in solution, including prenucleation cluster formation, is not well understood due to limitations both in characterization and first-principles modeling techniques. Recent advancements in machine-learned interatomic potentials are rapidly enabling the study of complex interactions such as dynamical behavior and reactivity on the nanoscale. Here, we develop an Au-S-C-H atomic cluster expansion (ACE) interatomic potential for efficient and accurate molecular dynamics simulations of thiolate-protected gold nanoclusters (Au n (SCH 3 ) m ). Trained on more than 30,000 density functional theory calculations of gold nanoclusters, the interatomic potential exhibits ab initio level accuracy in energies and forces and replicates nanocluster dynamics including thermal vibration and chiral inversion. Long dynamics simulations (up to 0.1 μs time scale) reveal a mechanism explaining the thermal instability of neutral Au 25 (SR) 18 clusters. Specifically, we observe multiple stages of isomerization of the Au 25 (SR) 18 cluster, including a chiral isomer. Additionally, we simulate coalescence of two Au 25 (SR) 18 clusters and observe series of clusters where the formation mechanisms are critically mediated by ligand exchange in the form of [Au-S] n rings.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Thermal and electric multidomain dynamic model for integration of power grid distribution with behind-the-meter devices
As renewable energy sources like solar and wind power become more integrated into the grid, coordinated control of behind-the-meter devices is crucial for enhancing grid flexibility and reliability and for meeting cost targets, with standardized models being developed to support this transition. The increasing flexibility and uncertainty of integrated renewable energy grids, along with interactions between various subsystems, make traditional steady-state modeling insufficient to capture transient and dynamic behaviors. Current models (e.g., composite load and battery equivalent models) focus on thermodynamic or electrical characteristics but overlook critical electromechanical interactions. This limits the ability to share performance information for grid services and hampers fast dynamic simulations. In addition, motor stalling is usually triggered by a fault event and attributed to the characteristics of the mechanical torque of the motor, resulting in absorption of a large amount of reactive power during the stalling period. Further, this significant withdrawal of reactive power will deteriorate the dynamic voltage stability of power grids and cause delayed voltage recovery. Therefore, an in-depth modeling of the thermodynamics or mechanical torque is essential to study the impacts of the realistic torque characteristics of those behind-the-meter devices on power system voltage stability. This study developed a dynamic multidomain model for building HVAC systems, such as air-source heat pumps, to simulate their thermal and electrical responses to grid transients. The model can accurately predict power metrics with a mean absolute percentage error of 10%, by validating against with power system computer-aided design performance data. Case studies demonstrate the model capability of capturing the transient response to sudden voltage changes, rapid load fluctuations, and system shutdowns respectively. During a sudden voltage drop (30% for 0.1s), a fully loaded heat pump’s motor speed dropped, continued declining, and shut down after 3.6s, with severe power oscillations and a torque spike. A partially loaded unit experienced temporary oscillations but stabilized. Under higher building loads, compressor speed increased from 64% to 100%, with power and torque rising before stabilizing. In safety-triggered shutdowns, power decreased after minor fluctuations, and torque briefly spiked before dropping to zero.
24 POWER TRANSMISSION AND DISTRIBUTION↗
Machine Learning-Accelerated First-Principles Molecular Dynamics Explains Anomalous Lattice Thermal Expansion in BaZr 0.78 Y 0.22 O 3-δ
Fuel cells are a vital clean energy technology that converts chemical energy directly into electricity with high efficiency, making them a cornerstone of a sustainable energy future. Herein we investigate the thermal and chemical lattice expansion behavior of hydrated BaZr 0.78 Y 0.22 O 3-δ using machine learning-accelerated ab initio molecular dynamics simulations. Here, our results reproduce the experimentally observed non-monotonic and anomalous temperature dependence of lattice expansion, which we attribute to the competing effects of thermal expansion and dehydration—two mechanisms that influence the lattice expansion in opposite directions. The importance of this work lies in its detailed demonstration of how advanced computational techniques can accurately capture complex environmental effects, providing a valuable framework for modeling similar phenomena in a variety of material systems and applications.
Proton conducting fuel cell↗
Ineffectiveness of formamidine in suppressing ultralow thermal conductivity in cubic hybrid perovskite FAPbI 3
Understanding lattice dynamics and thermal transport mechanisms in cubic hybrid organic–inorganic perovskites remain challenging due to strong anharmonicity and phase transitions. Here, we investigate the thermal transport behavior in benchmark cubic hybrid perovskite FAPbI 3 by coupling first principles-based anharmonic lattice dynamics with a linearized Wigner transport equation. Using the Temperature-Dependent Effective Potential (TDEP) method, we stabilize the negative soft modes, primarily dominated by organic FA + cations. Our calculations predict an ultra-low thermal conductivity of ~$0.63\,{\rm{W}}{{\rm{m}}}^{-1}{{\rm{K}}}^{-1}$ at 300 K, following a temperature dependence of T −0.740 . Contrary to common assumptions, we find that the [PbI 3 ] 1- units, rather than FA + cations, dominate thermal resistance. Furthermore, we demonstrate that anharmonic force constants are highly temperature-sensitive, relying on 0-K force constants significantly underestimates thermal conductivity. Our study not only elucidates the microscopic mechanisms governing thermal transport in FAPbI 3 but also provides a robust framework for modeling heat conduction in hybrid organic-inorganic compounds.
Zheng, Jiongzhi [Dartmouth College, Hanover, NH (U↗
Dynamic Nanocrystal Superlattices with Thermally Triggerable Lubricating Ligands
The size-dependent and collective physical properties of nanocrystals (NCs) and their self-assembled superlattices (SLs) enable the study of mesoscale phenomena and the design of metamaterials for a broad range of applications. However, the limited mobility of NC building blocks in dried NCSLs often hampers the potential for employing postdeposition methods to produce high-quality NCSLs. In this study, we present tailored promesogenic ligands that exhibit a lubricating property akin to thermotropic liquid crystals. The lubricating ability of ligands is thermally triggerable, allowing the dry solid NC aggregates deposited on the substrates with poor ordering to be transformed into NCSLs with high crystallinity and preferred orientations. The interplay between the dynamic behavior of NCSLs and the molecular structure of the ligands is elucidated through a comprehensive analysis of their lubricating efficacy using both experimental and simulation approaches. Coarse-grained molecular dynamic modeling suggests that a shielding layer from mesogens prevents the interdigitation of ligand tails, facilitating the sliding between outer shells and consequently enhancing the mobility of NC building blocks. Further, the dynamic organization of NCSLs can also be triggered with high spatial resolution by laser illumination. The principles, kinetics, and utility of lubricating ligands could be generalized to unlock stimuli-responsive metamaterials from NCSLs and contribute to the fabrication of NCSLs.
77 NANOSCIENCE AND NANOTECHNOLOGY↗
Contribution of cryogenic thermal cycling to the atomic dynamics in a La-based bulk metallic glass with different initial states
Structural rejuvenation in metallic glasses (MGs) induced by cryogenic thermal cycling has been intensively studied. However, the effect of thermal cycling on the atomic dynamics in MGs is still missing. In this work, we present a systematic study on the atomic dynamics in a La-based MG affected by sub-Tg annealing and thermal cycling. We find that the thermal cycling has little effect on the atomic dynamics in both as-cast and annealed MGs although it could significantly increase the free volume and change mechanical properties of the as-cast MG sample. In contrast, both the structure and mechanical properties of the well-relaxed sample are almost unchanged by thermal cycling. Further, the increased open volume, mainly in the form of atomic vacancies, does not greatly promote the atomic dynamics in the as-cast MG sample, indicating that the free volume induced by thermal cycling is not the key factor in stimulating the β-relaxation and the boson peak in the studied MG.
36 MATERIALS SCIENCE↗
Evolutionary multi-objective optimization and Pareto-frontal uncertainty quantification of interatomic forcefields for thermal conductivity simulations
Predictive Molecular Dynamics simulations of thermal transport require forcefields that can simultaneously reproduce several structural, thermodynamic and vibrational properties of materials like lattice constants, phonon density of states, and specific heat. This requires a multi-objective optimization approach for forcefield parameterization. Existing methodologies for forcefield parameterization use ad-hoc and empirical weighting schemes to convert this into a single-objective optimization problem. Here, we provide and describe software to perform multi-objective optimization of Stillinger–Weber forcefields (SWFF) for two-dimensional layered materials using the recently developed 3rd generation non-dominated sorting genetic algorithm (NSGA-III). NSGA-III converges to the set of optimal forcefields lying on the Pareto front in the multi-dimensional objective space. This set of forcefields is used for uncertainty quantification of computed thermal conductivity due to variability in the forcefield parameters. We demonstrate this new optimization scheme by constructing a SWFF for a representative two-dimensional material, 2H-MoSe 2 and quantifying the uncertainty in their computed thermal conductivity.
97 MATHEMATICS AND COMPUTING↗
An analysis of parameters affecting ampacity in aircraft bipolar MVDC power cables via coupled electrical, thermal, and computational fluid dynamic modelling
Abstract The next generation of aircraft, including more electric aircraft and all‐electric aircraft (AEA), requires electric power systems with high power density and low system mass specifications. Increasing the voltage of the system to the range of a few kV, medium voltage (MV), is a reasonable approach to achieving high‐power‐density and low‐system‐mass EPSs for aircraft applications. Higher voltages, however, pose many challenges for aviation MV power cables such as arcs and arc tracking, partial discharges (PDs), and thermal management. In this regard, thermal management is more challenging since heat transfer by convection is greatly reduced at wide‐body aircraft's cruising altitudes due to the reduced air pressure. In this paper, a finite element method (FEM) model is developed in COMSOL Multiphysics for an aircraft bipolar MVDC (±5 kV) power cable. Using the model, the maximum permissible cable current at a low pressure of 18.8 kPa (at an altitude of 12.2 km from sea level, the usual cruising altitude for wide‐body aircraft) is calculated. Also, an analytical model is developed based on analytical and proven empirical correlations governing conductive, radiative, and convective heat transfers at the steady state to estimate the ampacity of the bipolar cable system at reduced pressure. It was shown that the proposed analytical model can be used for atmospheric pressure and systems with a larger number of poles, expanding its range of applications. The results of the FEM and analytical models correlate at wide ranges of parameters such as ambient temperature, duct size, distance between the positive and negative pole cables, and the overall diameter of the cables. The influence of horizontal and vertical arrangement of poles is included in the analytical model. The results of this study can be used to design bipolar MVDC power cable systems for the envisaged wide‐body AEA.
Azizi, Arian↗
Thermal stability and coalescence dynamics of exsolved metal nanoparticles at charged perovskite surfaces
Exsolution reactions enable the synthesis of oxide-supported metal nanoparticles, which are desirable as catalysts in green energy conversion technologies. It is crucial to precisely tailor the nanoparticle characteristics to optimize the catalysts’ functionality, and to maintain the catalytic performance under operation conditions. We use chemical (co)-doping to modify the defect chemistry of exsolution-active perovskite oxides and examine its influence on the mass transfer kinetics of Ni dopants towards the oxide surface and on the subsequent coalescence behavior of the exsolved nanoparticles during a continuous thermal reduction treatment. Nanoparticles that exsolve at the surface of the acceptor-type fast-oxygen-ion-conductor SrTi 0.95 Ni 0.05 O 3–δ (STNi) show a high surface mobility leading to a very low thermal stability compared to nanoparticles that exsolve at the surface of donor-type SrTi 0.9 Nb 0.05 Ni 0.05 O 3–δ (STNNi). Our analysis indicates that the low thermal stability of exsolved nanoparticles at the acceptor-doped perovskite surface is linked to a high oxygen vacancy concentration at the nanoparticle-oxide interface. For catalysts that require fast oxygen exchange kinetics, exsolution synthesis routes in dry hydrogen conditions may hence lead to accelerated degradation, while humid reaction conditions may mitigate this failure mechanism.
25 ENERGY STORAGE↗
Grey-box modeling and application for building energy simulations - A critical review
Grey-box modeling, as one of the three fundamental modeling techniques for building energy models, has many advantages compared with black-box modeling and white-box modeling. Additionally, it has been widely applied to solve problems of building technologies, such as building load estimation, control and optimization, and building-grid integration. However, a thorough review of grey-box modeling is not available. This review study systematically investigated various aspects of grey-box modeling for buildings. First, the fundamental aspects of grey-box modeling are presented, including the theoretical background, modeling of building elements, modeling order, modeling diagram, and order reduction. Second, the detailed modeling approaches are discussed. Third, multiple applications of grey-box modeling are investigated for building energy domain, which are categorized into the following groups: heat dynamics analysis, thermal load estimation, building control and optimization, district/urban scale energy modeling, and building-grid integration. Finally, the available software packages for grey-box modeling are compared. Overall, the challenges of using grey-box modeling can be summarized as follows: (1) the theoretical limitations and assumptions of grey-box modeling are unclear; (2) grey-box model naming convention and structure are confusing; (3) grey-box model creation is vague; (4) suitable applications of grey-box models are unknown; and (5) grey-box models lack unified software solutions for wider adoption.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
Enhancing the dark-state bistability and thermal stability of zinc dynamic windows with nitrile electrolytes
Reversible metal electrodeposition (RME) for dynamic windows is a promising alternative to traditional electrochromics. We explore the reasons behind the bistability of practical RME dynamic windows under a wide temperature range, i.e., how well they maintain their dark state without external power. While Zn dynamic windows with dimethyl sulfoxide electrolytes exhibit adequate bistability at room temperature, their performance degrades at elevated temperatures due to the presence of small amounts of water. Here, to overcome this problem, we developed a set of Zn nitrile electrolytes. The integrated devices can maintain dark-state bistability over at least 1 month at 85°C with 0.1% transmission using a propionitrile (PrCN) electrolyte. Additionally, the PrCN-based Zn dynamic windows possess high optical contrast, color-neutral switching, and promising cyclability, making Zn RME in nitrile electrolytes a promising solution for achieving robust dynamic windows.
RME↗