Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Thermodynamic limit”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Physical phenomena and the microgravity response

The living biological cell is not a sack of Newtonian fluid containing systems of chemical reactions at equilibrium. It is a kinetically driven system, not a thermodynamically driven system. While the cell as a whole might be considered isothermal, at the scale of individual macromolecular events there is heat generated, and presumably sharp thermal gradients exist at the submicron level. Basic physical phenomena to be considered when exploring the cell's response to inertial acceleration include particle sedimentation, solutal convection, motility electrokinetics, cytoskeletal work, and hydrostatic pressure. Protein crystal growth experiments, for example, illustrate the profound effects of convection currents on macromolecular assembly. Reaction kinetics in the cell vary all the way from diffusion-limited to life-time limited. Transport processes vary from free diffusion, to facilitated and active transmembrane transport, to contractile-protein-driven motility, to crystalline immobilization. At least four physical states of matter exist in the cell: aqueous, non-aqueous, immiscible-aqueous, and solid. Levels of order vary from crystalline to free solution. The relative volumes of these states profoundly influence the cell's response to inertial acceleration. Such subcellular phenomena as stretch-receptor activation, microtubule re-assembly, synaptic junction formation, chemotactic receptor activation, and statolith sedimentation were studied recently with respect to both their basic mechanisms and their responsiveness to inertial acceleration. From such studies a widespread role of cytoskeletal organization is becoming apparent.

Todd, Paul↗
NASA NTRS2016

Evaluating Secondary Inorganic Aerosols in Three Dimensions

The spatial distribution of aerosols and their chemical composition dictates whether aerosols have a cooling or a warming effect on the climate system. Hence, properly modeling the three-dimensional distribution of aerosols is a crucial step for coherent climate simulations. Since surface measurement networks only give 2-D data, and most satellites supply integrated column information, it is thus important to integrate aircraft measurements in climate model evaluations. In this study, the vertical distribution of secondary inorganic aerosol (i.e., sulfate, ammonium, and nitrate) is evaluated against a collection of 14 AMS flight campaigns and surface measurements from 2000 to 2010 in the USA and Europe. GISS ModelE2 is used with multiple aerosol microphysics (MATRIX, OMA) and thermodynamic (ISORROPIA II, EQSAM) configurations. Our results show that the MATRIX microphysical scheme improves the model performance for sulfate, but that there is a systematic underestimation of ammonium and nitrate over the USA and Europe in all model configurations. In terms of gaseous precursors, nitric acid concentrations are largely underestimated at the surface while overestimated in the higher levels of the model. Heterogeneous reactions on dust surfaces are an important sink for nitric acid, even high in the troposphere. At high altitudes, nitrate formation is calculated to be ammonia limited. The underestimation of ammonium and nitrate in polluted regions is most likely caused by a too simplified treatment of the NH3/NH4(+) partitioning which affects the HNO3/NO3(-) partitioning.

aerosols↗
DOE OSTI2023

High-Precision Calculation of Nanoparticle (Nanocrystal) Potentials of Mean Force and Internal Energies

Thermodynamic stability assessment of nanocrystal systems requires precise free energy calculations. This study highlights the importance of meticulous control over various factors, including the thermostat, time step, potential cutoff, initial configuration, sampling method, and overall simulation duration. Free energy computations in dry (solvent-free) systems are on the order of several hundred $k_BT$ but can be obtained with consistent accuracy. However, calculation of internal energies becomes challenging, as they are typically much larger in magnitude than free energies and exhibit significant noise and reduced reliability. To address this limitation, we propose a new internal energy estimate that drastically reduces the noise. Here, we also present formulas that enable the optimization of the parameters of the harmonic bias potential for optimal convergence. Finally, we discuss the implications of these findings for the computation of free energies in nanocrystal clusters and superlattices.

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

Structural Heterogeneity in Medium-Entropy AgMnSbPbTe 4 for Glassy Thermal Transport and High Thermoelectric Performance

Medium-entropy semiconductors represent a unique category of entropy-engineered materials. They possess a considerable level of randomness in atomic mixing, although this is not sufficient to conclusively achieve single-phase structure stabilization, in contrast to high-entropy materials. This introduces strong competition between the formation of different phases, which can potentially lead to structural heterogeneity. Here, in this work, we uncover endotaxial nanoprecipitates in the microscopically identified homogeneous medium-entropy semiconductor AgMnSbPbTe 4 . These nanoprecipitates initially crystallize in a cubic phase $(Fm\bar{3}m)$ within kinetically stabilized AgMnSbPbTe 4 , subsequently evolving into a thermodynamically stable monoclinic phase $(P2_1/c)$ during thermal annealing while maintaining an endotaxial relationship with the matrix lattice. This nanophase segregation and the resultant lattice mismatch at interfaces introduce strain fluctuations up to 5% at intervals of 20 nm across the entire microstructure. Within the matrix phase, atomic displacement of up to 23 pm was observed. This structural heterogeneity results in glass-like thermal transport behavior, achieving an ultralow lattice thermal conductivity κ L = 0.312 Wm –1 K –1 at 800 K, which is in accordance with the amorphous limit predicted by the Cahill model. The synergy of band convergence effect and well-maintained carrier mobility leads to a maximum ZT of 1.72 at 800 K and an average ZT avg of 1.02 over the temperature range of 300–825 K. This study highlights that the underexplored structural heterogeneity in medium-entropy semiconductors can potentially yield beneficial phenomena, such as the phonon-glass electron-crystal transport behavior in this case, which holds promise for advancing thermoelectric applications.

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

A tale of two polarization paradoxes: The diamagnetic polarization paradox

An accurate calculation of the total polarization charge density in a plasma is essential for a self-consistent determination of the electric field. Yet, the “diamagnetic polarization paradox” refers to the fact that there is a paradoxical factor of 1/2 difference between the pressure-driven “diamagnetic polarization” density calculated using real space drift theory vs action-angle space guiding center and gyrokinetic theory that has not been explained before. Here, we show that both results can be made consistent with one another. Half of the diamagnetic polarization is due to the transformation from the guiding center density to the real space density. The other half is due to the fact that, within the drift kinetic ordering assumptions, the guiding center density should be expressed as the gyroaverage of the density in the limit of vanishing Larmor radius. A comprehensive review of polarization is presented to complete the derivation, and then we derive results that are required for the polarization calculations to agree. Expressions for the diamagnetic polarization density are given that are accurate to first order in amplitude and all orders in gyroradius within the gyrokinetic theory for a constant magnetic field. Applications to Maxwell–Boltzmann particle distribution functions (PDFs), including anisotropic temperature, are presented. Local invariants, like total energy and toroidal momentum, do not generate net polarization effects; the electric and thermodynamic polarizations must precisely cancel. In contrast, anisotropic dependence on the magnetic moment generates a net polarization proportional to the temperature anisotropy. Finally, when sources are present, the equilibrium PDF is approximately the ratio of two orbit averages.

adiabatic theory↗
DOE OSTI2023

Thermodynamic stability of xenon-doped liquid argon detectors

Liquid argon detectors are employed in a wide variety of nuclear and particle physics experiments. The addition of small quantities of xenon to argon modifies its scintillation, ionization, and electroluminescence properties and can improve its performance as a detection medium. However, a liquid argon-xenon mixture can develop instabilities, especially in systems that require phase transitions or that utilize high xenon concentrations. In this work, we analyze the causes of these instabilities and describe a small (liter-scale) apparatus with a unique cryogenic circuit specifically designed to handle argon-xenon mixtures. The system is capable of condensing argon gas mixed with $\mathscr{O}$ (1%) xenon by volume and maintains a stable liquid mixture near the xenon saturation limit while actively circulating it in the gas phase. We also demonstrate control over instabilities that develop when the detector condition is allowed to deviate from optimized settings. This progress enables future liquid argon detectors to benefit from the effects of high concentrations of xenon doping, such as more efficient detection of low-energy ionization signals. In conclusion, this work also develops tools to study and mitigate instabilities in large argon detectors that use low concentration xenon doping.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Data-Driven Prediction of Formation Mechanisms of Lithium Ethylene Monocarbonate with an Automated Reaction Network

Interfacial reactions are notoriously difficult to characterize, and robust prediction of the chemical evolution and associated functionality of the resulting surface film is one of the grand challenges of materials chemistry. The solid-electrolyte interphase (SEI), critical to Li-ion batteries (LIBs), exemplifies such a surface film, and despite decades of work, considerable controversy remains regarding the major components of the SEI as well as their formation mechanisms. Here we use a reaction network to investigate whether lithium ethylene monocarbonate (LEMC) or lithium ethylene dicarbonate (LEDC) is the major organic component of the LIB SEI. Our data-driven, automated methodology is based on a systematic generation of relevant species using a general fragmentation/recombination procedure which provides the basis for a vast thermodynamic reaction landscape, calculated with density functional theory. The shortest pathfinding algorithms are employed to explore the reaction landscape and obtain previously proposed formation mechanisms of LEMC as well as several new reaction pathways and intermediates. For example, we identify two novel LEMC formation mechanisms: one which involves LiH generation and another that involves breaking the (CH 2 )O-C(=O)OLi bond in LEDC. Most importantly, we find that all identified paths, which are also kinetically favorable under the explored conditions, require water as a reactant. This condition severely limits the amount of LEMC that can form, as compared with LEDC, a conclusion that has direct impact on the SEI formation in Li-ion energy storage systems. Finally, the data-driven framework presented here is generally applicable to any electrochemical system and expected to improve our understanding of surface passivation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

State of the AGN: Progress Toward Understanding Black Hole Accretion Processes

Accretion of plasma onto black holes power some of the most powerful systems in the cosmos. Supermassive black holes at the centers of galaxies represent the high-mass limit of these objects and so account for the most luminous accretors. As a result, their influence spans vast spatial and temporal scales of cosmic phenomena: intracluster heating, intergalactic media, galactic feedback and star formation, kiloparsec-scale jets/outflows, variability over time scales of minutes to centuries, and luminous multi-wavelength electromagnetic emission extending all the way down to its event horizon. Their intrigue is heightened by the fact that they lie at the intersection of various physical laws---e.g., general relativistic gravity, magnetohydrodynamics, radiation, high-energy particle physics, thermodynamics, and photo-ionization---which all must be reconciled to arrive at a fundamental understanding and probe for new physics, like tests of general relativity. These physics ingredients must be incorporated into simulations performed somehow on dynamical scales ranging from that of the event horizon to parsec-scales. Fortunately, new computational and theoretical techniques---such as GPU computing, radiation transport, and and novel gridding techniques---are enabling progress to larger scales, more degrees of freedom, and even to binary systems. Some of the topics we will survey include recent progress on simulating the relationship between the disk-jet interaction, how tilted black holes behave, radiation-dominated flow, and how binary AGN affect the standard picture of black hole accretion. Along the way, we will highlight how new technologies have enabled these scientific rewards. Future directions and open questions will be provided to inspire discussion and interaction during the session.

LISA↗
DOE OSTI2023

Demonstrating non-isothermal vacuum membrane air dehumidification for efficient next-generation air conditioning

Vacuum membrane air dehumidification has gained significant interest in recent years as a highly efficient means of air dehumidification. Prior theoretical modeling work by the authors introduced the Active Membrane Energy Exchanger (AMX) concept, which combines active heat exchange and vacuum membrane dehumidification into one non-isothermal system, and found that it could outperform other air conditioning technologies under many conditions. However, no experimental literature exists on combining active heat exchangers and vacuum membrane dehumidification. Here, the goal of the present study is to evaluate the dehumidification (mass transfer) performance of the AMX concept relative to isothermal membrane dehumidification through three main methodologies: (1) membrane material-level characterization, (2) experimental prototype development and testing, and (3) computational fluid dynamics (CFD) simulations. The dense membranes used in this work showed up to a 40% increase in water vapor permeance at cooler temperatures, and the prototype system showed up to a 6% increase in humidity removal when the air was simultaneously cooled. Furthermore, the membrane module-integrated heat exchange tubes provided additional mixing and turbulence, leading to a 4–8% increase in humidity removal. The upper limit coefficient of performance was equal to approximately 2.2, necessitating advanced system designs to improve efficiency. This study is the first to highlight that combining the cooling and vacuum dehumidification processes can improve dehumidification performance.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Electrochemical Corrosion and Catalysis Dynamics of Tin Oxide during Water Oxidation

Metal oxide corrosion severely limits anodic electrocatalysis, particularly at high potentials in acidic environments, where degradation pathways remain poorly defined. This study establishes explicit connections between corrosion and electrocatalysis on tin oxide during water oxidation by examining the roles of lattice defects, reactive oxygen species, interfacial pH variations, and speciation of corroded tin in acid. We first demonstrate the presence of structural defects such as oxygen vacancies and substoichiometric Sn(II) species by integrating electron paramagnetic resonance spectroscopy, ultraviolet photoelectron spectroscopy, and Mott–Schottky analysis. Kohn–Sham density functional theory calculations reveal that explicit water structures thermodynamically stabilize reaction intermediates and lower reaction overpotentials. Moreover, we propose that water dissociation leads to hydrogen-bonding networks formed by H* and OH* intermediates, which may span the entire catalyst surface and decrease the interfacial pH to drive corrosion. In contrast, the electrochemical generation of reactive oxygen species is shown to play a minor role in catalyst corrosion during water oxidation using inductively coupled plasma mass spectrometry coupled with selective chemical scavengers. Square-wave voltammetry combined with rotating ring-disk electrodes is used to reveal that under open-circuit conditions, only Sn(IV) cations chemically dissolve from tin oxide, while both Sn(IV) and Sn(II) species electrochemically corrode during water oxidation. Our results unveil a dynamic and complicated interplay between corrosive and catalytic pathways on metal oxide electrocatalysts: a decrease in interfacial pH due to water oxidation exacerbates Sn(II)/Sn(IV) corrosion. Subsequently, the electrochemical corrosion of Sn(II)/Sn(IV) facilitates product formation from lattice oxygen, while the redeposition of corroded Sn(II) as Sn(IV) can enable oxygen exchange with water. By elucidating the roles of defects and interfacial chemistry, this work provides a roadmap for engineering improved electrocatalysts that balance activity and stability, a critical step toward scalable and durable energy technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Photon–photon chemical thermodynamics of frequency conversion processes in highly multimode systems

Abstract Frequency generation in highly multimode nonlinear optical systems is inherently a complex process, giving rise to an exceedingly convoluted landscape of evolution dynamics. While predicting and controlling the global conversion efficiencies in such nonlinear environments has long been considered impossible, here, we formally address this challenge even in scenarios involving a very large number of spatial modes. By utilizing fundamental notions from optical statistical mechanics, we develop a universal theoretical framework that effectively treats all frequency components as chemical reactants/products, capable of undergoing optical thermodynamic reactions facilitated by a variety of multi-wave mixing effects. These photon–photon reactions are governed by conservation laws that directly determine the optical temperatures and chemical potentials of the ensued chemical equilibria for each frequency species. In this context, we develop a comprehensive stoichiometric model and formally derive an expression that relates the chemical potentials to the optical stoichiometric coefficients, in a manner akin to atomic/molecular chemical reactions. This advancement unlocks new predictive capabilities that can facilitate the optimization of frequency generation in highly multimode photonic arrangements, surpassing the limitations of conventional schemes that rely exclusively on nonlinear optical dynamics. Notably, we identify a universal regime of Rayleigh–Jeans thermalization where an optical reaction at near-zero optical temperatures can promote the complete and entropically irreversible conversion of light to the fundamental mode at a target frequency. Our theoretical results are corroborated by numerical simulations in settings where second-harmonic generation, sum-frequency generation and four-wave mixing processes can manifest.

Optics↗
DOE OSTI2021

High-throughput design of high-performance lightweight high-entropy alloys

Developing affordable and light high-temperature materials alternative to Ni-base superalloys has significantly increased the efforts in designing advanced ferritic superalloys. However, currently developed ferritic superalloys still exhibit low high-temperature strengths, which limits their usage. Here we use a CALPHAD-based high-throughput computational method to design light, strong, and low-cost high-entropy alloys for elevated-temperature applications. Through the high-throughput screening, precipitation-strengthened lightweight high-entropy alloys are discovered from thousands of initial compositions, which exhibit enhanced strengths compared to other counterparts at room and elevated temperatures. The experimental and theoretical understanding of both successful and failed cases in their strengthening mechanisms and order-disorder transitions further improves the accuracy of the thermodynamic database of the discovered alloy system. This study shows that integrating high-throughput screening, multiscale modeling, and experimental validation proves to be efficient and useful in accelerating the discovery of advanced precipitation-strengthened structural materials tuned by the high-entropy alloy concept.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Improved reference system for the corrected rigid spheres equation of state model

The Corrected Rigid Spheres (CRIS) EOS model, developed from fluid perturbation theory using a hard sphere reference system, has been successfully used to calculate the EOS of many materials, including gases and metals. The radial distribution function (RDF) plays a pivotal role in choosing the sphere diameter, through a variational principle, as well as the thermodynamic response. Despite its success, the CRIS model has some shortcomings in that it predicts too large a temperature for liquid-vapor critical points, can break down at large compression, and is computationally expensive. In this paper, we first demonstrate that an improved analytic representation of the hard sphere RDF does not alleviate these issues. Relaxing the strict adherence of the RDF to hard spheres allows an accurate fit to the isotherms and vapor dome of the Lennard-Jones fluid using an arbitrary reference system. Second order correction are eliminated, limiting the breakdown at large compression and significantly reducing the computation cost. The transferrability of the new model to real systems is demonstrated on argon, with an improved vapor dome compared to the original CRIS model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1987

Polyphase alloys as rechargeable electrodes in advanced battery systems

The rechargeability of electrochemical cells is often limited by negative electrode problems. These may include loss of capacity, increased impedance, macroscopic shape change, dendrite growth, or a tendency for filamentary or whisker growth. In principle, these problems can be reduced or eliminated by the use of alloys that undergo either displacement or insertion reactions at reactant species activities less than unity, rather than pure elements. The fundamental reasons for some of these problems with elemental electrodes, as well as the basic principles involved in the different behavior of alloys, are briefly discussed. More information is now available concerning the thermodynamic and kinetic properties of a number of alloys of potential interest for use as electrodes in elevated temperature lithium battery systems. Recent results have extended these results down to ambient temperatures, indicating that some such materials may be of interest for use with new low temperature molten salt electrolytes, or with organic solvent electrolytes. The all solid mixed conductor matrix concept is also reviewed.

Huggins, Robert A.↗
NASA NTRS1990

A review of reaction rates and thermodynamic and transport properties for an 11-species air model for chemical and thermal nonequilibrium calculations to 30000 K

Reaction rate coefficients and thermodynamic and transport properties are reviewed and supplemented for the 11-species air model which can be used for analyzing flows in chemical and thermal nonequilibrium up to temperatures of 3000 K. Such flows will likely occur around currently planned and future hypersonic vehicles. Guidelines for determining the state of the surrounding environment are provided. Curve fits are given for the various species properties for their efficient computation in flowfield codes. Approximate and more exact formulas are provided for computing the properties of partially ionized air mixtures in a high energy environment. Limitations of the approximate mixing laws are discussed for a mixture of ionized species. An electron number-density correction for the transport properties of the charged species is obtained. This correction has been generally ignored in the literature.

Gupta, Roop N.↗
NASA NTRS1997

Kinetics and Thermochemistry of ClCO Formation from the Cl + CO Association Reaction

Laser flash photolysis of Cl2/CO/M mixtures (M = N2, CO, Ar, CO2) has been employed in conjunction with Cl((sup 2)P(sub J)) detection by time-resolved resonance fluorescence spectroscopy to investigate equilibration kinetics in the reactions Cl((sup 2)P(sub J)) + CO ClCO as a function of temperature (185-260 K) and pressure (14-200 Torr). The association and dissociation reactions are found to be in the low-pressure limit over the range of experimental conditions investigated. In N2 and/or CO buffer gases, the temperature dependences of the ClCO formation and dissociation reaction rate constants are described by the Arrhenius expressions k(sub 1) = (1.05 +/- 0.36) x 10(exp -34) exp[(810 +/- 70)/T] cm(exp 6)/molecules(exp 2).s and k(sub -1) = (4.1 +/- 3.1) x 10(exp -10) exp[(-2960 +/- 60)/T]cu cm/(molecule.s) (errors are 2 sigma). Second- and third-law analyses of the temperature dependence of the equilbrium constant (k/k-1) lead to the following thermodynamic parameters for the association reaction: Delta-H(sub 298) = -7.7 +/- 0.6 kcal/mol, Delta-H(sub 0) = -6.9 +/- 0.7 kcal/mol, Delta-S(sub 298) = -23.8 +/- 2.0 cal/mole.K, Delta-H(sub f,298)(ClCO) = 5.2 +/- 0.6 kcal/mol (errors are 2 sigma). The results repported in this study significantly reduce the uncertainties in all reported kinetic and thermodynamic parameters.

Nicovich, J. M.↗
NASA NTRS2015

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗
DOE OSTI2024

As‐Doped Polycrystalline CdSeTe: Localized Defects, Carrier Mobility and Lifetimes, and Impact on High‐Efficiency Solar Cells

Abstract The efficiency potential for single‐junction photovoltaics (PV) is described by the detailed balance model, which requires the elimination of nonradiative recombination and perfect minority carrier collection. Improvements in GaAs, Si, and perovskite PV follow this model. It might be more complex for CdTe, a leading thin‐film PV technology. While lifetime, passivation, and doping goals for 25% efficient CdTe solar cells are largely reached, voltage is ≈20% below the detailed balance limit. Why is that? In Se‐alloyed CdSe x Te 1‐x (Se is required for >20% efficiency) additional losses can occur due to electrostatic and bandgap fluctuations and due to electronic trap states. To understand mechanisms limiting CdSeTe solar cell performance and to suggest improvements, carrier dynamics, and transport in CdSe x Te 1‐x with variation in Se composition and as doping is analyzed. It is shown that trapping, likely due to anion‐site defects and their complexes, is correlated with low charge carrier mobility of 0.1–0.6 cm 2 (Vs) −1 . Even with 1000 ns charge carrier lifetimes, carrier diffusion length is less than the absorber thickness, reducing efficiency to ≈23%. Device simulations are used to analyze the performance of CdSe x Te 1‐x solar cells; thermodynamic models are not sufficient for absorbers with electronic disorder and trapping.

14 SOLAR ENERGY↗