SEARCH · Search NASA

Results for “nuclear level density”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Assessing Implicit and Explicit Polarizable Solvation Models for Nuclear–Electronic Orbital Systems: Quantum Proton Polarization and Solvation Energetics

Accurate simulations of many chemical processes require the inclusion of both nuclear quantum effects and a solvent environment. The nuclear–electronic orbital (NEO) approach, which treats electrons and select nuclei quantum mechanically on the same level, combined with a polarizable continuum model (PCM) for the solvent environment, addresses this challenge in a computationally practical manner. In this work, the NEO-PCM approach is extended beyond the IEF-PCM (integral equation formalism PCM) and C-PCM (conductor PCM) approaches to the SS(V)PE (surface and simulation of volume polarization for electrostatics) and ddCOSMO (domain decomposed conductor-like screening model) approaches. IEF-PCM, SS(V)PE, C-PCM, and ddCOSMO all exhibit similar solvation energies as well as comparable nuclear polarization within the NEO framework. The calculations show that the nuclear density does not leak out of the molecular cavity because it is much more localized than the electronic density. Finally, the polarization of quantized protons is analyzed in both continuum solvent and explicit solvent environments described by the polarizable MB-pol model, illustrating the impact of specific hydrogen-bonding interactions captured only by explicit solvation. Furthermore, these calculations highlight the relationship among solvation formalism, nuclear polarization, and energetics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Interplay of structural fluctuations and charge carrier dynamics is key for high performance of hybrid lead halide perovskites

The ability of methylammonium lead triiodide (MAPbI 3 ) to achieve photoelectric conversion efficiency that is on par with crystalline silicon has led to a surge of interest in perovskite photovoltaics. However, an in-depth understanding of how the ubiquitous coupling between the fast rovibrational movements of the organic cations and the phonon vibrations of the inorganic framework affects the relaxation and recombination of hot carriers remains largely elusive. Access to such knowledge is critical to guide design and increase efficiency of new classes of perovskite solar cells. In this report, we report a time-domain ab initio investigation of temperature dependent excited state dynamics in MAPbI 3 , with particular emphasis on nuclear anharmonic effects. The observed strong anharmonicity is attributed to softness of the material and unusual dynamical coupling between the organic and inorganic components. At an elevated temperature, the hydrogen bonding between MA and iodines is weakened, enhancing rotation of MA cations, which become more dynamically disordered. As a result, thermal vibrations of the inorganic Pb–I sublattice are suppressed, and the lattice anharmonicity is decreased. Thermal fluctuations of the electronic energy levels are found to follow the trend of anharmonic motions of Pb and I atoms, with holes relaxing faster to the band edges than electrons, due to higher density of the hole states. While elevated temperature accelerates intraband carrier cooling, it slows nonradiative carrier recombination. The latter result is important for performance, since solar cells and other devices heat up during operation. The reported signatures of coupled structural dynamics of the organic cations and inorganic framework unravel the interplay between anomalous structural fluctuations and charge carrier dynamics, which is of particular importance for fundamental understanding of the structure–property relationships in hybrid metal halide perovskites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Single-particle and collective excitations in Zn 66

Single-particle and collective excitations in 66 Zn have been investigated via the multinucleon transfer reaction, 26 Mg( 48 Ca, α4nγ) using the Gammasphere multidetector array and the Fragment Mass Analyzer. In addition to confirming and complementing the previously known low-spin structure, a new quasi-rotational band comprising several stretched E2 transitions has been established to high spins. However, due to fragmentary nature of its decay, it was not possible to link this sequence to the low-lying states and, thus, determine the absolute excitation energies, spins, and parities unambiguously. Large-scale shell-model calculations employing the JUN45 and jj44b effective interactions are able to successfully describe the low-spin structure and herewith confirm that it is dominated by single-particle excitations. Furthermore, the newly established rotational cascade is compared with known superdeformed bands in the A ≈ 60–70 mass region, and with results of calculations performed within the frameworks of the cranked shell model and the adiabatic and configuration-fixed constrained covariant density functional theory and the quantum particle-rotor model.

59 ≤ A ≤ 89↗

DOE OSTI2026

Helium co-implantation effects on dislocation loop evolution during in situ ion irradiations of a novel additive-manufactured Fe-9Cr steel alloy

This study examines the impact of helium injection rate on dislocation loop evolution during ion irradiation of a novel, additively manufactured (AM) 9Cr reduced-activation ferritic/martensitic (RAFM) steel. Through in situ transmission electron microscopy ion irradiations and post-irradiation characterization, the effects of simultaneous 1 MeV Kr²⁺ and 12 keV He⁺ ion irradiation at 400°C to 5 dpa with varying helium levels were examined. The findings reveal that helium co-injection significantly increases dislocation loop density, primarily through the formation of small defect features exhibiting contrast and projected morphologies consistent with a⟨100⟩ loops (∼1 nm), while suppressing larger a/2⟨111⟩ loops. The findings underscore the complex interplay among continuously increasing helium content, temperature, and alloy composition in the evolution of dislocation loops under irradiation in RAFM steels. The study highlights the need for further investigation into helium-mediated mechanisms, such as proposed interactions between self-interstitial atoms and helium atoms, to enhance the performance of RAFM steels in advanced nuclear reactors.

ANA2↗

DOE OSTI2020

Chiral effective field theory calculations of weak transitions in light nuclei

In this work, we report quantum Monte Carlo calculations of weak transitions in A≤10 nuclei, based on the Norfolk two- and three-nucleon chiral interactions, and associated one- and two-body axial currents. Furthermore, we find that the contribution from two-body currents is at the 2–3% level, with the exception of matrix elements entering the rates of 8Li, 8B, and 8He β decay. These matrix elements are suppressed in impulse approximation based on the (leading order) Gamow Teller transition operator alone; two-body currents provide a 20–30% correction, which is, however, insufficient to bring theory in agreement with experimental data. For the other transitions, the agreement with the data is satisfactory, and the results exhibit a negligible to mild model dependence when different combinations of Norfolk interactions are utilized to construct the nuclear wave functions. We report a complete study of two-body weak transition densities which reveals the expected universal behavior of two-body currents at short distances throughout the range of A=3 to A=10 systems considered here.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

High-resolution infrared spectroscopy of jet cooled CH 2 Br radicals: The symmetric CH stretch manifold and absence of nuclear spin cooling

Direct laser absorption of a slit supersonic discharge expansion provides the first high-resolution spectroscopic results on the symmetric CH stretch excitation (ν 1 ) of the bromomethyl (CH 2 Br) radical in the ground electronic state. Narrowband (<1 MHz) mid-infrared radiation is produced by difference-frequency generation of two visible laser beams, with the open shell halohydrocarbon radical generated by electron dissociative attachment of CH 2 Br 2 in a discharge and rapidly cooled to T rot = 18 ± 1 K in the subsequent slit-jet supersonic expansion. A rovibrational structure in the radical spectrum is fully resolved, as well as additional splittings due to spin-rotation effects and 79 Br/ 81 Br isotopologues in natural abundance. Spectroscopic constants and band origins are determined by fitting the transition frequencies to a non-rigid Watson Hamiltonian, yielding results consistent with a vibrationally averaged planar radical and an unpaired electron in the out-of-plane p π orbital. Additionally, extensive satellite band structure from a vibrational hot band is observed and analyzed. The hot band data is compared to CFOUR/VPT2 (CCSD(T)cc-pVQZ) ab initio anharmonic predictions of the vibration rotation alpha matrix, which permits unambiguous assignment to CH 2 symmetric-stretch excitation built on the singly excited CH 2 out-of-plane bending mode (ν 1 + ν 4 ← ν 4 ). Longitudinal cooling of the Doppler width in the slit-jet expansion geometry also reveals partially resolved hyperfine structure on transitions out of the lowest angular momentum states in excellent agreement with predictions based on microwave studies. High level ab initio MOLPRO calculations [CCSD(T)-f12b/VnZ-f12 (n = 3, 4, CBS)] are also performed with explicitly correlated f12 electron methods for the out-of-plane CH 2 bending mode over the halogen series CH 2 X (X = F, Cl, Br, I), which clearly reveals a non-planar geometry for X = F (with a ΔE ≈ 0.3 kcal/mol barrier) and yet planar equilibrium geometries for X = Cl, Br, and I. Finally, a detailed Boltzmann analysis of the transition intensities provides support for negligible collisional equilibration of the entangled H atom nuclear spin states on the few hundred microsecond time scale and high collision densities of a slit supersonic expansion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Bayesian inference of nuclear incompressibility from collective flow in mid-central Au+Au collisions at 400–1500 MeV/nucleon

The incompressibility K of symmetric nuclear matter (SNM) is determined through a Bayesian analysis of collective flow data from Au + Au collisions at beam energies $E = 400 -1500$ MeV/nucleon. This analysis utilizes a Gaussian process (GP) emulator applied to the isospin-dependent quantum molecular dynamics (IQMD) model for heavy-ion collisions, both with and without incorporating the momentum dependence of the single-nucleon potentials. Specifically, at the 68% confidence level, using rapidity and transverse velocity dependence of proton elliptic flow data with and without consideration of the momentum dependence, the inferred incompressibility values are $K=188.9^{+2.9}_{-4.5}$ MeV and $256.1^{+8.2}_{-8.7}$ MeV at $E = 400$ MeV/nucleon, respectively. When the transverse momentum dependence of proton-like directed flow data is included, the inferred incompressibility values become $K=222.3^{+9.0}_{-9.9}$ MeV and $K=285.5^{+6.7}_{-7.3}$ MeV, respectively. Furthermore, we found that the value of K derived from observables of proton elliptic flow increases with beam energy. Finally, this indicates that the equation of state (EoS) of nuclear matter hardens at higher densities and temperatures in reactions with higher beam energies.

Bayesian inference↗

DOE OSTI2026

Effects of Composition and Oxidation States on the Structures of Chromium-Containing Sodium Silicate Glasses: Molecular Dynamics Simulations using Machine Learning Interatomic Potentials

Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.

Puga, Christina L.↗

DOE OSTI2025

Atomic Orbital Energy Matching vs Overlap in Actinide-Ligand Dative Bonding

Prebonding conditions and dative bond formation in actinide(IV)hexachloride complexes with U, Np, and Pu are studied by density functional theory (DFT) calculations, focusing on the interplay of atomic orbital (AO) overlap and AO energy matching in the formation of molecular orbitals (MOs) and the subsequent identification of dative bonds. The extent of donation is tracked via population analysis, MO localization, and bond-order criteria. DFT Fock matrices are used to setup models in which varying numbers of valence AOs interact to form dative bonds. The results confirm, among other effects, a contribution to metal−ligand covalency from the actinide (An) 6p shell. Better matching of An(5f) and Cl(3p) energies occurs as the An(5f) level stabilizes with increasing An effective nuclear charge. A near-degeneracy occurs in the case of the α-spin orbitals of the Pu system, but it is inconsequential. Altogether, better AO energy matching for An(5f) and Cl(3p) along the series U, Np, Pu is counter-balanced by decreasing AO overlap and decreasing availability of 5f acceptor orbitals, leading to similar bond orders and extents of donation in the three systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Theoretical Examination of the Hydroxide Transport in Cobaltocenium-Containing Polyelectrolytes

Polymers incorporating cobaltocenium groups have received attention as promising components of anion-exchange membranes (AEMs), exhibiting a good balance of chemical stability and high ionic conductivity. In this work, we analyze the hydroxide diffusion in the presence of cobaltocenium cations in an aqueous environment based on the molecular dynamics of model systems confined in one dimension to mimic the AEM channels. In order to describe the proton hopping mechanism, the forces are obtained from the electronic structure computed at the density-functional tight-binding level. We find that the hydroxide diffusion depends on the channel size, modulation of the electrostatic interactions by the solvation shell, and its rearrangement ability. Hydroxide diffusion proceeds via both the vehicular and structural diffusion mechanisms with the latter playing a larger role at low diffusion coefficients. The highest diffusion coefficient is observed under moderate water densities (around half the density of liquid water) when there are enough water molecules to form the solvation shell, reducing the electrostatic interaction between ions, yet there is enough space for the water rearrangements during the proton hopping. Furthermore, the effects of cobaltocenium separation, orientation, chemical modifications, and the role of nuclear quantum effects are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Inferring relative dose-dependent color center populations in proton irradiated thoria single crystals using optical spectroscopy

In this work, we have utilized photoluminescence spectroscopy and optical ellipsometry to characterize the dose-dependence of the photoluminescence emission intensity and changes in optical absorption of thoria single crystals subject to irradiation with energetic protons at room- and elevated-temperatures. The photoluminescence peaks and the optical absorption bands are attributed to creation of new electronic states emerging from defects resulting from displacement damage. These bands are most likely associated with electrons trapped at the oxygen vacancy sites similar to color centers formed in other irradiated oxides and halides. Our experimental observations are supported by a standard density functional theory calculation of the electronic structure in pristine and oxygen vacancy-bearing thoria crystals. The dose-dependence of the intensity of the photoluminescence peaks is used to parameterize a rate theory model that estimates the concentration of color centers in the irradiated crystals. This parameterization provides optimized migration barrier parameters for oxygen interstitials and vacancies that simultaneously capture the optical response of the crystals irradiated at room- and elevated-temperature. These optical spectroscopy techniques offer a promising pathway to characterize the population of color centers formed at the sites of oxygen anion vacancies, particularly in irradiated nuclear fuels, where atomic-level defects cannot be readily imaged using electron microscopy. When combined with other direct and indirect characterization tools, our approach can provide new insight into defect formation and accumulation in energy materials over single atomic to extended length scales.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Using the Monte Carlo Method to Evaluate the Reliability of Screening Multifamily Housing for Radon

When screening for radon in a multifamily housing complex using a fixed sample density (e.g., testing 1 in 10 (10%) or 1 in 4 units (25%)), the statistical confidence is dependent upon the assumed elevated radon frequency. For example, testing 25% of units in a complex estimated to have eight units with elevated radon levels will provide 90% confidence. However, if it is assumed that, in the same complex, there are only three units with elevated radon levels, the confidence drops to around 58% for the same sample density. Furthermore, in the previous example, if elevated radon levels are not found during the screening, all that can be stated is that the screening provides 90% confidence that there are no more than seven units in the complex with elevated radon levels. To more fully illustrate this uncertainty, ten separate multifamily housing radon data sets with 1 to 10 units with radon levels ≥4 pCi/L (Table 1) were selected for analysis using the Monte Carlo statistical method. Unlike other mathematically based statistical approaches, the Monte Carlo statistical method relies on repeated analysis of randomly selected data. from a 100% sampled complex at various sample densities. Success for each simulation is defined as finding at least one unit with elevated radon levels. In this statistical method, confidence in the overarching conclusion can be greatly enhanced by repeating the simulated screening hundreds or even thousands of times. For this study, each of the ten data sets was simulated 1,000 times. For each simulation, the data set was randomized three times before the fixed percentage of data was selected.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2026

Measurement of the Neutron Elastic Electromagnetic Form Factor Ratio at Large Momentum Transfer

Exploring nucleon structure is vital both for understanding its origin and existence as well as for the advancement of the sciences. It helps us answer key questions such as how quark and gluon dynamics create 99% of the nucleon mass. Electron- nucleon scattering has been widely used for precision studies of the nucleon and nuclear structure since the Nobel Prize winning investigations by Robert Hofstadter and collaborators in the 1950s. These studies provide information about the spatial charge and current densities of the nucleon in terms of the electromagnetic form factors. The form factors are functions of four momentum transfer squared (Q2). Extending the electromagnetic form factor measurements to higher Q2 plays a critical role in furthering the understanding of nucleon structure. This motivated the Super BigBite Spectrometer (SBS) program at Jefferson Lab. The open nature of the spectrometers and the direct line of sight from the target to the tracking detector locations in experimental setups such as SBS creates high levels of background at the detectors. This necessitates the use of tracking detectors with high rate capability and good position resolution. Gas Electron Multiplier (GEM) detectors are an excellent choice for tracking detectors in such experiments. Understanding the performance of the GEM detectors is important not just for SBS experiments but also for future high-luminosity experiments. This thesis reports the exploratory results from the measurement of the neutron elastic electromagnetic form factor ratio (Gn E/Gn M) at high momentum transfer. A longitudinally polarized electron beam was scattered off a polarized 3He target, used as an effective polarized neutron target. In this experiment, the polarized 3He target achieved a world record polarization weighted luminosity at a beam current of 45 µA. Double spin asymmetry of the scattered neutron events is used to extract the neutron form factor ratio. Measurements were taken at Q2 = 3.0, 6.8, 9.8 (GeV/c)2. The lowest Q2 measurement is in good agreement with the existing world data, and the higher-Q2 measurements extend the Q2 reach well beyond the existing world data and are expected to remain unmatched for a long time.

Gamage, Vimukthi Haththotuwa [Univ. of Virginia, C↗

DOE OSTI2022

Design of a Scavenging Pyrrole Additive for High Voltage Lithium-Ion Batteries

We report 1-(dimethylamino) pyrrole (PyDMA) as an electrolyte additive for high voltage lithium-ion batteries based on LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC622)//Graphite with an upper cutoff voltage of 4.4 V. Density Functional Theory (DFT) modeling indicates that the unique structure of PyDMA could be effective in preventing the hydrolysis of LiPF 6 in a carbonate electrolyte, mitigating issues related to HF formation. The calculations also indicated that the additive would oxidize at lower potentials than typical electrolyte solvents, which could lead to protective films at the cathode surface. These expectations were tested using Nuclear Magnetic Resonance (NMR) and extensive electrochemical characterization. NMR studies confirmed the superb dehydrating capability of PyDMA, which successfully prevents HF formation even at high water content. Addition of 0.5 wt% PyDMA resulted in improved capacity retention in full-cells, and also in lower levels of transition metal dissolution from the cathode. Incremental capacity (dQ/dV) analysis indicates that benefits of PyDMA at low concentration (0.5–1 wt%) are associated with decreased rates of Li + -trapping reactions, and that higher concentrations of the additive can lead to isolation of cathode domains. Furthermore, our study indicates that PyDMA could be a promising electrolyte additive for high voltage lithium-ion batteries at a low concentration.

25 ENERGY STORAGE↗

DOE OSTI2024

First-principles investigation of the thermophysical properties of NaCl, PuCl 3 , and NaCl-PuCl 3 Molten salts

Molten salts have a variety of applications that span the nuclear and solar industries, and which involve thermal storage and heat transfer. There is a present knowledge gap in the thermophysical properties of molten salts, which limits the readiness level of molten salt applications. This is especially pertinent for molten salt reactors where the fissile material is dissolved within the molten salt. Ab initio molecular dynamics is a common method employed to investigate the structural and thermophysical properties at elevated temperatures. NaCl-PuCl 3 is a candidate fuel salt in molten salt reactors. The scope of this study is to investigate the seven unique compositions and a range of temperatures of NaCl-PuCl 3 and calculate the density, heat capacity, compressibility, enthalpy of mixing, and coefficient of thermal expansion. Within this work, the van der Waals (vdW) interactions are primarily handled with the vdW-DF2 functional but spot-checked with both the dDsC and DFT-D3 methods. The calculated densities are in good agreement with the NaCl and the eutectic compositions, and no experimental values exist for PuCl 3 at these temperatures. The densities are fit to a first-order Redlich-Kister expansion as a function of both composition and temperature. The heat capacity in the literature is scattered but the calculated values agree well with one of the experimental results. Here, the heat capacity increases as a linear function with respect to concentration at a rate of about 7.5 J/mol-K per 10 mol% PuCl 3 .

36 MATERIALS SCIENCE↗

DOE OSTI2022

Evolution of Urca Pairs in the Crusts of Highly Magnetized Neutron Stars

Abstract We report on the effects of strong magnetic fields on neutrino emission in the modified Urca process. We show that the effect of Landau levels on the various Urca pairs affects the neutrino emission spectrum and leads to an angular asymmetry in the neutrino emission. For low magnetic fields, the Landau levels have almost no effect on the cooling. However, as the field strength increases, the electron chemical potential increases resulting in a lower density at which Urca pairs can exist. For intermediate field strength, there is an interesting interference between the Landau level distribution and the Fermi distribution. For high enough field strength, the entire electron energy spectrum is eventually confined to a single Landau level producing dramatic spikes in the emission spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Heat transport in liquid water from first-principles and deep neural network simulations

In this work, we compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represented by a deep neural network (DNN) trained on DFT data, whereby the PES has an explicit local decomposition and the energy flux takes a particularly simple expression. By virtue of a gauge invariance principle, established by Marcolongo, Umari, and Baroni, the two approaches should be equivalent if the PES were reproduced accurately by the DNN model. We test this hypothesis by calculating the thermal conductivity, at the GGA (PBE) level of theory, using the direct formulation and its DNN proxy, finding that both approaches yield the same conductivity, in excess of the experimental value by approximately 60%. Besides being numerically much more efficient than its direct DFT counterpart, the DNN scheme has the advantage of being easily applicable to more sophisticated DFT approximations, such as meta-GGA and hybrid functionals, for which it would be hard to derive analytically the expression of the energy flux. We find in this way that a DNN model, trained on meta-GGA (SCAN) data, reduces the deviation from experiment of the predicted thermal conductivity by about 50%, leaving the question open as to whether the residual error is due to deficiencies of the functional, to a neglect of nuclear quantum effects in the atomic dynamics, or, likely, to a combination of the two.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Validation of the Porous Medium Approximation for Hydrodynamics Analysis in Compact Heat Exchangers

Compact heat exchangers (HXs) have gained attention in recent years in various fields such as solar and nuclear power generation, oil and gas, and refrigeration due to their low cost, high power density, and robustness in high-pressure and/or high-temperature environments. However, the large difference between a compact HX's overall dimensions (~m) and the much smaller scale of its channels (~mm) makes it challenging to model the entire HX at once, due to computational limitations. In this work, we treat the channeled region of a compact HX as a porous medium (PM) to circumvent the need to model/mesh each individual channel. This allows us to simulate the entire HX, including both the header and channeled regions while maintaining the computational cost at a practical level. Although the porous medium approach has been used to model heat exchangers, its validity is still questionable because (1) the resistance coefficients are heavily data-based and thus difficult to be applied to new heat exchangers and (2) the validation has been focused on matching the overall pressure drop in the channel region, which does not address whether such model can predict detailed pressure and velocity field. For the first time, this work addresses under what circumstances and with what uncertainty does the PM approach work for hydrodynamics modeling in compact HXs. By answering these questions, we introduce the PM approach as a powerful tool for HX hydrodynamics modeling that can predict not only the overall pressure drop but also the detailed pressure and velocity distributions.

Engineering↗