Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two phase flow”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Mode-multiplexed photonic integrated vector dot-product core from inverse design

Photonic computing has the potential to harness the full degrees of freedom (DOFs) of the light field, including the wavelength, spatial mode, spatial location, phase quadrature, and polarization, to achieve a higher level of computing parallelism and scalability than digital electronic processors. While multiplexing using the wavelength and other DOFs can be readily integrated on silicon photonics platforms with compact footprints, conventional mode-division multiplexed (MDM) photonic designs occupy areas exceeding tens to hundreds of microns for a few spatial modes, significantly limiting their scalability. Here, we utilize inverse design to demonstrate an ultracompact photonic computing core that calculates vector dot products based on MDM coherent mixing. Our dot-product core integrates the functionalities of two-mode multiplexers and one multimode coherent mixer within a nominal footprint of 5 μm x 3 μm . We have experimentally demonstrated computing examples on the fabricated dot-product core, including complex number multiplication and motion estimation using optical flow. The compact dot-product core design enables large-scale on-chip integration in a parallel photonic computing primitive cluster for high-throughput scientific computing and computer vision tasks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Rotary Inertia Friction Welding of Dissimilar High-Strength 422 Martensitic Stainless Steel and 4140 Low Alloy Steel for Heavy-Duty Engine Piston Fabrication

AISI 422 martensitic stainless steel with superior high temperature performance (oxidation resistance and strength) is under evaluation for replacing current heavy-duty piston crown materials, AISI 4140 martensitic steel and micro alloyed steel (MAS) 38MnSiVS5, to fabricate a multimaterial piston. This multimaterial piston concept further improved power density and fuel economy by allowing heavyduty diesel engines to operate at higher temperatures and pressures. Joining AISI 422 steel piston crowns with AISI 4140 steel piston skirts is a key manufacturing step for this multimaterial piston. However, the significant differences in strength, elevated temperature flow stress, alloy chemistry, and temper resistance between these two martensitic steels cause some weldability issues (cracking) and metallurgical challenges (alloying element migration/segregation) when using conventional fusion-based welding processes. Rotary inertia friction welding (RIFW), a solid-state welding process, has been the preferred method to join 4140 crowns to 4140 skirts (and MAS crowns to MAS skirts) in high-volume production of current heavy-duty diesel engine pistons. It has been used to join these two materials with relatively comparable alloy chemistry to fabricate pistons with MAS skirts and 4140 crowns. Meanwhile, RIFW has also been a preferred method of dissimilar metal welding. However, RIFW of dissimilar high-strength martensitic steels has yet to be widely pursued. The interfacial microstructure complexities created by the thermomechanical process and highly nonequilibrium phase transformations during RIFW are a significant challenge for understanding and predicting their joining behavior and have not been reported in detail. Here, in this work, defect-free AISI 422 steel-AISI 4140 multimaterial pistons were successfully fabricated using the RIFW process. The interfacial microstructure and mechanical properties of dissimilar 422/4140 steel RIFW in the as-welded condition were experimentally studied in detail. The results provide critical baseline information for understanding RIFW mechanisms and guiding subsequent postweld heat treatment (PWHT) practice.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Statistical distributions for transient transport

Here, this paper introduces the use of statistical distributions based on transport differential equations for clear distinction of transport modes within transient kinetic experiments. More specifically, novel techniques are developed for the transient data obtained through the Temporal Analysis of Products (TAP) reactor and are applicable to experiments where pulse response into a gas flow is used. The methodology allows distinguishing between two domains of diffusion transport in heterogeneous catalytic systems, i.e., Knudsen and non-Knudsen diffusion, using statistical fingerprints, and finding the transition domain. Two distribution parameters were obtained that directly result in coefficients that correspond to the concentration and the rate of transport. Using a linear relationship between the rate and concentration coefficients, Knudsen diffusion is revealed when the rate of transport is constant and non-Knudsen diffusion is confirmed when the rate of transport coefficient is a function of the concentration coefficient. As a result, accurate transport information can be extracted from experimental data even in the presence of comprising instrument drift or noise particularly when analyzing higher pressure pulse responses with complex transport. This enables more direct investigation of experiments influenced by gas-phase reactions.

TAP reactor↗
DOE OSTI2025

Glacial geology of the Hudson Mountains, Amundsen Sea sector, West Antarctica

The Hudson Mountains are situated in the eastern Amundsen Sea sector of the West Antarctic Ice Sheet, adjacent to Pine Island Glacier. They form a volcanic field of 17 stratovolcanoes and parasitic vents, preserved as nunataks. Two former tributaries of Pine Island Glacier (Larter and Lucchitta glaciers) flow through the mountains. Here we present a detailed study of the glacial geology of the area. We describe field observations and measurements of geomorphological features from 15 of the nunataks, meltwater ponds found on the surface of three nunataks and supraglacial features (ice dolines) from two sites near the present grounding line. Together these provide constraints on the past ice sheet extent, flow pathways and thermal regime, and enhance our understanding of the present hydrological regime – all of which are important as context for the observed modern ice sheet behaviour. We find evidence suggesting that all nunataks in the Hudson Mountains were covered by ice during the Last Glacial Maximum (defined here as 26.5-19 ka) and have since deglaciated. Faceted and polished erratic cobbles and boulders of exotic lithologies (syenites, alkali granites, granites, granodiorites, tonalites and gabbros) are numerous and perched on nunatak surfaces. A marked difference between the dominant erratic lithologies on nunataks adjacent to Pine Island Glacier (granite) and Lucchitta Glacier (granodiorite-tonalite) indicates that the ice sheet was transporting clasts from at least two distinct upstream source regions. The similarity in degree of weathering suggests, however, that all the erratics were transported by one phase of (warm-based) glaciation; their presence on or close to the summits of all except one nunatak indicates that the ice sheet during that time was at least 700 m thicker than present. These results are consistent with ice sheet model simulations which suggest that all nunataks in the Hudson Mountains were completely submerged by the Last Glacial Maximum ice sheet.

58 GEOSCIENCES↗
DOE OSTI2025

Geometric entropies and their Hamiltonian flows

In holographic theories, the Hubeny-Rangamani-Takayanagi (HRT) area operator plays a key role in our understanding of the emergence of semiclassical Einstein-Hilbert gravity. When higher derivative corrections are included, the role of the area is instead played by a more general functional known as the geometric entropy. It is thus of interest to understand the flow generated by the geometric entropy on the classical phase space. In particular, the fact that the associated flow in Einstein-Hilbert or Jackiw-Teitelboim (JT) gravity induces a relative boost between the left and right entanglement wedges is deeply related to the fact that gravitational dressing promotes the von Neumann algebra of local fields in each wedge to type II. This relative boost is known as a boundary-condition-preserving (BCP) kink-transformation. In a general theory of gravity (with arbitrary higher-derivative terms), it is straightforward to show that the flow continues to take the above geometric form when acting on a spacetime where the HRT surface is the bifurcation surface of a Killing horizon. However, the form of the flow on other spacetimes is less clear. In this paper, we use the manifestly-covariant Peierls bracket to explore such flows in two-dimensional theories of JT gravity coupled to matter fields with higher derivative interactions. The results no longer take a purely geometric form and, instead, demonstrate new features that should be expected of such flows in general higher derivative theories. We also show how to obtain the above flows using Poisson brackets.

AdS-CFT correspondence↗
DOE OSTI2025

Development of Hydrogen Burner for FT4000® Aeroderivative Engine - Final Report

This report details an effort to develop a retrofittable fuel/air injector for the FT4000® aeroderivative gas turbine that enables use of hydrogen as a carbon-free fuel for efficient power generation. The FT4000 engine’s low-NOx combustor was developed by Pratt & Whitney and RTX Technology Research Center with core technology from the Pratt & Whitney PW4000 turbofan aircraft engine. The current FT4000 production engine operates on either natural gas or No. 2 fuel oil with water injection to achieve high thermal efficiency and low emissions. This engine is fielded by Mitsubishi Power Aero and delivers 70 MW of power with a simple-cycle efficiency of over 41% when operating with wet compression. The work reported here advances the technology readiness level of the FT4000 combustor for operation with hydrogen, starting with an experimental assessment of the current production hardware with increasing hydrogen content mixed with natural gas and ending with improved nozzle concepts for fully robust operation with 100% hydrogen. High-pressure single-sector combustor rig tests have been completed, demonstrating the ability for the dual fuel nozzle to operate an FT4000 combustor on 100% hydrogen with low nitrogen oxide (NOx) emissions. Metal temperature measurements and video images of the flame structure from zero to 100% hydrogen highlight opportunities to improve the fuel nozzle robustness for high hydrogen conditions. The design of new fuel/air mixer concepts to improve durability and operability with high hydrogen levels was also completed. A total of eleven new concepts were developed and analyzed, ranging from modifications to the bill-of-materials nozzle to fully clean sheet designs. The concepts were evaluated with non-reacting and reacting flow evaluations to assess the performance of the new hardware designs. Non-reacting tests included Phase Doppler Particle Analysis (PDPA) for droplet size and velocity, mechanical patternation for liquid water flux, and planar laser induced fluorescence (PLIF) with acetone-seeding for gaseous fuel/air mixing characterization. Five scaled candidate fuel nozzle designs, in addition to a scaled bill-of-materials nozzle, were then successfully evaluated in an atmospheric pressure burner rig. The nozzles were evaluated for performance with natural gas, hydrogen/natural gas blends, and pure hydrogen. For 100% hydrogen, the nozzles were evaluated with and without water injection. Optical and infrared imaging of the flame and fuel nozzle was captured. NOx emissions were sampled from fixed emissions probes. The results show measurable differences between the various designs, and the data was used to down-select the two most promising designs to advance to future full pressure rig testing. Results from this study have cleared the current production FT4000 engines with dual fuel nozzles to operate at baseload power on blends of hydrogen mixed with natural gas and water. Two promising nozzle designs have been developed to enable fully robust operation with up to 100% hydrogen. These nozzles require validation at full baseload operating conditions before they can be introduced to the field.

03 NATURAL GAS↗
DOE OSTI2025

Multi-physics Topology OPtimization and Additive Manufacturing for High-temperature Heat Exchangers

This research significantly advances the understanding of high-temperature heat exchanger design through an integrated approach that combines topology optimization (TO), triply periodic minimal surface (TPMS) structures, additive manufacturing (AM) and thermohydraulic testing. Each of these components contributes uniquely to a unified, high-performance design, fabrication and testing workflow. Topology optimization serves as the foundation of the design methodology by providing a systematic way to determine the most effective material layout for separating hot and cold fluids while maximizing thermal performance. The researchers introduced a novel three-material optimization framework using two density fields to represent hot fluid, cold fluid, and solid domains. This approach enables automated discovery of optimal shapes and flow paths that cannot be intuitively designed, especially under constraints imposed by manufacturing technologies. Furthermore, constraints such as minimal wall thickness and overhang angles were embedded into the optimization process, ensuring that resulting designs are not only thermally efficient but also manufacturable using modern additive techniques. In parallel, the study delves into the use of Gyroid-based TPMS geometries for constructing the core of the heat exchanger. TPMS structures are known for their high surface area, excellent fluid mixing capabilities, and minimal pressure drop characteristics. The researchers applied a data-driven modeling framework using Heteroscedastic Sparse Gaussian Process Regression (HSGPR) combined with genetic algorithms. This allowed for the rapid evaluation and optimization of key geometric parameters such as frequency, iso-value, and phase shift. The result was a set of Gyroid structures tailored for high heat transfer and low flow resistance, demonstrating clear improvements over conventional straight-channel designs. After the designing process, additive manufacturing played a critical role by turning these highly complex, optimized geometries into physical components. Utilizing Laser Powder Bed Fusion (LPBF) with Haynes 282, the study demonstrated the feasibility of fabricating these heat exchangers at high precision. Post-processing methods, including dilation-erosion operations, were applied to ensure local features adhered to self-supporting constraints. The fabricated structures were then subjected to thermohydraulic testing under conditions representative of supercritical CO 2 Brayton cycles, validating the predicted performance and confirming the viability of the full design-to-fabrication pipeline. Finally, thermohydraulic testing across the above studies served as a crucial experimental validation of advanced heat exchanger. Under consistent high-temperature and high-pressure conditions using supercritical CO 2 , the testing demonstrated that both TO and Gyroid-based TPMS designs significantly outperformed conventional straight-channel HXs. The TO design achieved a 115% increase in UA and NTU and a 27.6% boost in gravimetric power density, while the data-driven optimized Gyroid design delivered a 166% increase in UA and NTU and improved effectiveness from 68.7% to 86.1%. These results validate the simulation models, confirm the manufacturability of complex geometries under AM constraints, and provide key insights into design-performance trade-offs, thereby advancing the development of high-efficiency, compact heat exchangers for extreme environments.

36 MATERIALS SCIENCE↗
DOE OSTI2026

MARS-Q modeling of low- n resistive kink-peeling modes and edge harmonic oscillations in DIII-D

Linear and quasilinear magnetohydrodynamic (MHD) modeling is carried out for two DIII-D discharges that both featured a transition from quiescent H-mode (QH) to wide-pedestal QH (WPQH). The MHD perturbations, associated with the edge harmonic oscillations (EHOs) observed during the QH-phase in both discharges, are identified as low-n (n is the toroidal mode number) resistive kink-peeling instabilities, with the ideal MHD counterpart remaining stable. The quasilinear model successfully simulates EHO-like perturbations during the QH phase in both discharges, confirming the experimental observations. A less intuitive finding is the EHO-like behavior involving the n = 2 perturbation, simulated for one of the discharges during the WPQH-phase. This result, while consistent with experimental observations, is obtained despite the fact that the initial perturbation is linearly stable. The eventual growth of the perturbation is solely due to nonlinear interaction between the MHD perturbation and the plasma toroidal flow. The occurrence of the EHO-like perturbation during the WPQH-phase is found to be sensitive to the initial profile of the plasma edge rotation. For DIII-D plasmas considered, the neoclassical toroidal viscosity, generated by three-dimensional low-n perturbations, is found to play a dominant role in modifying the edge flow during EHOs.

EHO↗
DOE OSTI2024

Solid State Quantum Refrigeration Superconducting, Absorption and Measurement Based (Final Technical Report)

During this DOE grant, DE-SC0017890, in place for the past six years, all proposed research was carried out and published in peer-reviewed papers, as well as other projects that emerged during the research. In that effort the research team accomplished all proposed research, as well as many closely related research projects discovered and conceived of during the grant. These works included “Efficient Quantum Measurement Engines”, a work published in Physical Review Letters, giving a theory of quantum measurement-based engines, which uses quantum measurement as a resource. These engines are designed to efficiently convert energy from the stochastic quantum measurement process into useful work. Further publications include “Experimental Realization of a Quantum Dot Energy Harvester”, a joint theory and experimental work in collaboration with the group of Charles Smith in Cambridge, UK, as well as long time theoretical collaborators, Rafael Sánchez and Björn Sothmann. This work, featured as an Editor’s Suggestion in Physical Review Letters, realized an earlier theoretical proposal of ours, whereby two resonant tunneling quantum dots are connected to a central electronic cavity that is heated by a hot energy source. We also published “Superconducting Quantum Refrigerator: Breaking and Rejoining Cooper Pairs with Magnetic Field Cycles” a work done in collaboration with experimentalist Francesco Giazotto from ENS Pisa, Italy, which also resulted in a patent. This paper, published in Phys. Rev. Applied, advanced the concept of a cyclic fridge based on the normal/superconducting phase transition together with layered materials separated by tunnel junctions. We also completed the proposed research on a heat transistor, publishing “Thermal transistor and thermometer based on Coulomb-coupled conductors”, carried out as a collaboration between my group and theorists Splettstoesser (Lund U., Sweden), Sothmann (U. Duisburg-Essen, Germany), and Sánchez (U. Autónoma de Madrid, Spain). We carried out an analysis of a quantum coupled to a quantum point contact as a sensitive thermometer and heat transistor. We found the optimal statistical estimator for the temperature and compared it with experiments on the same type of devices. We also investigated autonomous quantum absorption refrigerators using quantum dots to cool by using a very hot thermal reservoir to drive heat between two other reservoirs. In the article “Quantifying the quantum heat contribution from a driven superconducting circuit”, we demonstrated that for a driven superconducting circuit, we showed heat flow provided by a hot source to the qubit can be switched on and off by varying external parameters, the frequency and the intensity of the driving. In the work “Stochastic thermodynamic cycles of a mesoscopic thermoelectric engine”, we reconsidered the autonomous thermoelectric heat engine in terms of underlying cycles. Rather than periodic behavior, the cycles were stochastic in nature. Nevertheless, by undertaking a graph theoretical analysis of the elementary transport processed, great quantitative and qualitative insight could be found. We also considered the quantum measurement process and showed that a quantum version of Maxwell’s demon could be related to the work extraction of a quantum system, closely related to arrow-of-time measures for quantum measurement, as described in our article “Thermodynamics of quantum measurement and Maxwell's demon's arrow of time”. This work was selected in Phys. Rev. A as an Editor’s Suggestion. A recent preprint titled “Cyclic Superconducting Quantum Refrigerators Using Guided Fluxon Propagation” accomplished an important piece of this grant: to propose a new kind of quantum refrigerator using the dynamics of fluxons in a type II superconductor. This invention envisioned a race-track type geometry where fluxons are confined. By applying a gradient of magnetic field together with electric current in a Corbino geometry, the circulating fluxons can actively cool a cold reservoir, realizing a new type of cyclic superconducting refrigerator. We also investigated the possibility of thermal control from different points of view. The application of quantum measurement to the system gives a new kind of control on the system of interest – we have pioneered this approach and shown that measurement can boost the thermal power of quantum engines as described in “Continuous measurement boosted adiabatic quantum thermal machines”. The ability to have heat flows on demand is an outstanding challenge, and we have provided new solutions to this problem in Thermal control across a chain of electronic nanocavities” for a chain of electron cavities using gating voltage control. The control methods using qubit/qubit coupling to create absorption fridges at their most fundamental level have also been developed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

The Onset of Nucleate Boiling Prediction in a Mini Channel with the Eulerian Multiphase Flow

The prediction of boiling incipience is a critical issue for a reactor. Vapor lock causes operating instability and undesirable temperature rise, especially for small channels. This study uses the Eulerian multiphase flow boiling model to predict the onset of nucleate boiling in a millimeter scale rectangular channel. Two wall boiling models, namely the RPI (Rensselaer Polytechnic Institute) and the non-equilibrium sub cooled model are applied. Computational fluid dynamics is employed for this study to investigate the impact of operating conditions like liquid flow rate, operating pressure, heat flux, velocity profile, turbulence model, and inlet liquid temperature. Here, the effect of fundamental boiling parameters on the surface wall temperature and vapor volume fraction is also studied. The parameters under consideration are the bubble departure diameter, bubble departure frequency, nucleate site density, quenching time period, and interface heat transfer coefficients. A series of CFD calculations is conducted by varying the considered variables systematically in a wide range of flow conditions covering laminar, transition, and turbulent flows. The sub-component heat fluxes like the liquid convective, the vapor convective, the evaporative, and the quenching terms are monitored to infer the boiling dynamics in the transition region from the single-phase flow to the nucleate boiling zone. It is found that the CFD approach to detecting the boiling incipience point agrees reasonably with available experimental data. However, its limitations, like inaccuracy for lower flow rates and early transition to boiling flow, are also noticed due to the nature of the considered RPI model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Depinning, melting, and sliding of generalized Wigner crystals in Moiré systems

We numerically examine the depinning, sliding, and melting of commensurate and incommensurate Wigner crystals on two-dimensional hexagonal periodic substrates near fillings of 1 / 3 , 1 / 2 , and 2 / 3 to model the dynamics of generalized Wigner crystals in moiré heterostructures. At low temperatures where thermal fluctuations are irrelevant, for commensurate fillings we find a strongly pinned state that depins into a sliding crystal, while at incommensurate fillings the depinning threshold is strongly reduced. For fillings below 1 / 2 , the depinning occurs in a two-step process. Above the first depinning threshold, there is an extended range of drives where the conduction occurs via the sliding of antikinks along the charge stripes, followed by a second threshold where all the charges begin to slide. At finite temperatures, the external driving reduces the effective melting temperature at commensurate fillings and enhances creep at incommensurate fillings. We show that depinning into different sliding states, such as moving fluids or moving crystals, produces nonlinear features in the transport curves. We also show that transport is asymmetric on either side of a commensurate filling due to the different dynamics for interstitials versus holes in the commensurate structure. A variety of sliding states should be accessible for generalized Wigner crystals at finite temperatures even for low drives. If experiments can realize sliding dynamics in moiré systems, it would open another class of commensurate and incommensurate phases for study. Additionally, the sliding dynamics and antikink flow could provide unique functionalities for charge transport in moiré heterostructures. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

XRISM high-resolution X-ray spectroscopy of Cygnus X-1: Highly ionized iron absorption structures

We present the first high-resolution X-ray spectral analysis of Cygnus X-1 using XRISM. The observation wa3s carried out from 2024 April 7 to 10, covering the orbital phase range 0.65–0.17 during its low/hard state. Taking advantage of the exceptional energy resolution of the Resolve instrument, we examined highly ionized iron absorption lines and characterized the ionization states, column densities, and line-of-sight velocities of the absorbing plasma. Spectral analysis revealed an ionization parameter of $\xi \sim 3$, column densities of a few $\times 10^{21}$ cm$^{-2}$, and a blueshifted velocity of $\sim$100 km s$^{-1}$. The observation was divided into two phases: before and after orbital phase $\phi _{\rm {orb}} = 0.9$, corresponding to non-dipping and dipping intervals. While only weak absorption features were present before $\phi _{\rm {orb}} = 0.9$, strong absorption by He-like and H-like Fe appeared during the dipping phase. We measured equivalent widths of 2.3, 0.4, and 1.2 eV for He-like Fe K$\alpha$ and H-like Ly$\alpha _1$ and Ly$\alpha _2$, respectively—demonstrating the capability of XRISM Resolve to securely detect narrow absorption features of only a few eV. These measurements trace the motion of the absorbing material and offer insight into the kinematics and spatial distribution of the wind in the vicinity of the black hole. These findings enhance our understanding of wind-fed accretion in Cygnus X-1 and highlight the importance of continued high-resolution X-ray observations to further constrain the physical properties of winds and accretion flows in high-mass X-ray binaries.

Astronomy and AstroPhysics↗
DOE OSTI2026

Prototype Modeling for a Light-Trapping Planar-Cavity Enclosed Particle Solar Receiver

Concentrating solar thermal (CST) systems present a promising avenue for affordable and reliable energy production. Solar receivers are key components that determine the efficiency and longevity of these systems. Particle-based solar receivers have emerged as a compelling alternative to traditional technologies, offering several advantages that address limitations in current CST systems. This is especially true as next-generation CST technologies target applications including electricity generation, thermochemical processes, and industrial process heat, many of which necessitate higher operating temperatures than current commercial molten salt systems. Molten-salt thermal energy storage (TES) systems, commonly used in CSP, face challenges related to freezing and corrosion. Particle-based TES systems, in contrast, do not experience these issues, as particles are stable at high temperatures, exceeding 1000 degrees Celsius. This capability allows for a wider range of applications, including those requiring higher temperatures for industrial processes and efficient electricity generation. A novel innovation in particle-based solar receiver technology is the light-trapping planar cavity receiver (LTPCR) configuration developed by NREL. The LTPCR design consists of small cavity-like structures using opaque planar surfaces, enabling efficient capture and absorption of solar energy. A high incident flux concentration at the cavity aperture is absorbed on the receiver walls, and subsequently transferred to particles on the inside of cavities. The particles flow through the system, forming a fluidized bed inside of the receiver panels, effectively capturing the absorbed solar heat. Air is used as a fluidizing medium in this process to enhance particle heat transfer and mixing. The effectiveness of this design lies in its ability to manage solar flux conditions and ensure high solar-to-thermal receiver efficiency. A 100-kW prototype is currently being tested at the King Saud University in Saudi Arabia to assess the receiver performance. A range of modeling analyses for the optical, thermal, and mechanical effects were conducted to assess the performance of the receiver under on-sun conditions. The solar flux resulting from the KSU heliostat field was modeled using NREL SolTrace software and produced up to 600 kW/m2 at the receiver aperture. The solar flux absorbed on the receiver walls was then used within a computational fluid dynamics (CFD) model to predict wall temperature distributions along with radiation and convection loss. A two-phase CFD model was developed for the fluidized bed of silica sand inside the receiver panels to predict local wall-to-particle heat transfer coefficients, particle temperature distributions, and outlet temperature of the particles. We have also conducted analyses to understand the thermomechanical behavior of these innovative enclosed light-trapping solar receivers optimized for particle heating. We used finite element analysis (FEA) to predict the receiver's performance using temperature distributions obtained from CFD and based on the resulting stress profiles, evaluated creep-fatigue damage with a goal of achieving a 30-year service life. Analysis showed a significant impact of the particle-to-wall heat transfer coefficients (HTCs) on receiver performance, with higher HTCs resulting in reduced stress and increased lifespan. For instance, when using Inconel 740H, increasing the HTC from 800 W/m2 K to 1400 W/m2 K increased the creep life from 4,000 hours to over 100,000 hours. This highlights the importance of understanding and optimizing heat transfer in the design of high-efficiency receivers.

14 SOLAR ENERGY↗
DOE OSTI2024

Power and isotope effects in the ITER baseline scenario with tungsten and tungsten-equivalent radiators in DIII-D

Abstract Experiments in DIII-D document the ITER Baseline Scenario (IBS) at q 95 ∼ 3 and P IN / P LH ∼ 1–2, in both deuterium and hydrogen utilizing Kr and Xe as Tungsten-equivalent radiators. The power threshold for H-mode operation ( P LH ) was determined experimentally without added impurities and found to be about a factor of two higher than the scaling law. In recent IBS experiments in deuterium, intrinsic levels of metals such as Tungsten (W) or molybdenum and inconel are present that reduce the pedestal pressure by 20%–25%. A complete radiative collapse of deuterium IBS plasmas occurs at W core concentrations C W = 10 −5 . Simulations show that for core temperatures expected for ITER, the plasmas would not have a radiative collapse at C W = 1 × 10 −5 , moreover Q = 8–10 would still be achieved for C W up to 3 × 10 −5 . In contrast to deuterium, the IBS in hydrogen is not affected by intrinsic high-Z impurities, indicating that hydrogen H-modes in ITER may not inform the D-T phase with respect to W accumulation and discharge survival. Compared to deuterium, the pedestal pressure in hydrogen is ∼25% lower, with much higher ELM frequency of 150 Hz, decreasing with input power. Krypton was injected in a matrix scan of input power and impurity flow in IBS hydrogen discharges. Krypton impurity density profiles in hydrogen are similar to deuterium plasmas, but at Kr flows that are 2–3 times higher for the same input power. Krypton is transported into the core and affects the whole radius; at the highest injection rates a radiative collapse occurs at core radiation fractions of 0.3–0.35, consistent with the expected maximum W radiation fraction for ITER core plasmas. Comparing the results with previous International Tokamak Physics Activity database studies of the IBS confirms that at higher radiation fraction due to high-Z impurities, a drop in H 98 of >10% is observed. On the other hand, the results using Kr as a W-equivalent radiator indicate that metal (W) devices at lower core temperatures than ITER may provide overly pessimistic performance extrapolations to ITER for deuterium-tritium operation. The new DIII-D results support a more attractive option for the ITER Research Plan with a short hydrogen phase for system commissioning, transitioning to deuterium operations as soon as possible to provide relevant conditions for deuterium-tritium operations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Long-term socioeconomic trends and climate variability as drivers of virtual water scarcity in China

Water scarcity can have far-reaching sectoral impacts beyond its physical location through the propagation of virtual water flows. Socioeconomic and hydroclimatic changes affect local and virtual water scarcity by altering availability and demand. Yet most studies of this phenomenon focus on volumetric footprints, and the few on water scarcity risk have not examined hydroclimatic variability beyond long-term trends. In this study, we ask how gross domestic product (GDP) and population changes, long-term meteorological trends, sea surface temperature (SST) patterns, and interannual meteorological variability affect water scarcity in China, both locally (through the local water scarcity risk index, LWSR) and remotely (through the virtual water scarcity risk index, VWSR). Counterfactual scenarios were compared in a regression-and-simulation framework, with the socioeconomic and meteorological drivers varying over 1923–2019 and the multi-regional input–output structure staying fixed at 2017. Relative to a 5 year baseline centered on 2017, GDP and population changes have induced a cumulative 17%–50% increase in LWSR and 13%–21% increase in VWSR, outweighing the effect of long-term meteorological trends. phase change in one of two examined SST patterns induce 4%–13% differences in LWSR and 1%–4% differences in VWSR. Interannual meteorological variability induces 10%–20% standard deviations in LWSR and 3%–7% in VWSR. The findings highlight the importance of using longer time series to accurately assess local and virtual water scarcity situations. Water scarcity management should prioritize socioeconomic factors when planning at century-long timescales and prioritize hydroclimatic factors at multidecadal or shorter timescales. water managers should consider interannual variabilities in LWSR and VWSR and plan for potential occurrences of extreme conditions.

climate variability↗
DOE OSTI2024

A Review of Modeling Approaches for Predicting Frost Growth and Defrosting on Tube-Fin Heat Exchangers: Preprint

Frost formation and growth on the evaporator surface is a common process that deteriorates the air-refrigerant heat transfer and restricts airflow. This degrades the performance of the vapor compression system by increasing temperature lift and air-side pressure drop. To accurately predict these effects during coil frosting, as well as the energy use and duration of the defrost process, there is a need to estimate the heat and mass transfer, momentum transport, and solid-liquid and liquid-vapor phase change. Therefore, in the past few decades, continuous effort has been made to model frosting and defrosting processes using approaches ranging from empirical correlations to computational fluid dynamic models. To provide a clearer overview for researchers, engineers, and manufacturers in this field, this paper provides a comprehensive literature review for frosting and defrosting models. The paper begins with theoretical background of frost formation and defrost processes, and then reviews the common modeling approaches in literature and their underlying assumptions when trying to account for various physical phenomenon. Based on the literature review, the most critical modeling effort for frost formation is the determination of frost densification rate and frost growth rate. Various methods to predict these two parameters are reviewed. Empirical correlations commonly used for frost density and thermal conductivity are presented and compared. For the defrost process, various multi-stage models have been proposed with different assumptions. Some assume the presence of air gap between the tube wall and the frost, while others consider the melted frost flow due to gravity. We also review physics-based and empirical approaches to integrate defrost models into heat pump models. We conclude by identifying research gaps and providing recommendations.

defrost↗
DOE OSTI2024

Center for Tokamak Transients Simulations

The major goals for the Center for Tokamak Transients Simulations (CTTS) were to develop detailed numerical simulations of disruptive transients in tokamak experiments and to use them in theoretical studies of the thermal-quench (TQ) and current-quench (CQ) phases of disruption. The University of Wisconsin-Madison component of CTTS contributed relevant development for the NIMROD code (https://nimrodteam.org) and applied it to disruptions that involve global vertical displacement. Code-development contributions include boundary conditions that represent sheath effects which occur when tokamak plasma is in contact with the surrounding vessel wall during vertical displacement events (VDEs). We also implemented a reduced model of energetic runaway electrons (REs) to simulate the effect of these energetic electrons on the CQ phase of tokamak disruptions. Our computations of VDEs include 2D and 3D verification through benchmarks with two other widely used macroscale plasma simulation codes, M3D-C1 and JOREK. Over nonlinear evolution to the final termination of plasma current, the 2D simulation results on magnetic axis location, toroidal current, and halo current track each other well. With the 3D benchmark computation, all three codes found the same qualitative behavior of the asymmetric instabilities having increasing growth rate after contact with the wall, followed by nonlinear excitation of other toroidal harmonics. Our NIMROD computations for an idealized configuration show agreement with an analytical assessment by V. D. Pustovitov [Nucl. Fusion 55, 113032 (2015)] that the plasma and non-ideal wall remain in force-balance such that net forces only result from magnetic stresses over the outer surface of the wall. Our study of the influence of boundary conditions on VDEs shows that extended-MHD simulations of VDEs are sensitive to electron energy transport to the wall and that boundary conditions on flow are only significant to the extent that they influence the electron energy transport. Resistive-MHD NIMROD simulations of the tokamak discharges in the Madison Symmetric Torus show that magnetic perturbations of poloidal wavenumber m=3 lead to chaotic magnetic topology over the edge region, which helps deconfine energetic electrons. A second area of RE study shows that the coupled resistive MHD/reduced RE model reproduces the resistive-hose particle-beam instability. With respect to educational opportunities, this effort involved and supported four graduate students at the University of Wisconsin-Madison and one postdoctoral associate. The results of this study have been disseminated through journal publications, conference presentations, technical reports, and PhD dissertations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Testing- and Model- Based Optimization of Coal-fired Primary Heater Design for Indirect Supercritical CO 2 Power Cycles (Final Scientific and Technical Report)

The overall objective of this project was to perform the R&D necessary to mitigate the risk associated with the design of a primary heat exchanger for a solid-fired combustion system coupled with an indirect-fired closed-loop Brayton Cycle utilizing supercritical CO 2 . The key technological hurdle was the coupling of a solid-fuel firing system with the primary heater, which poses a singular challenge, which is the management of burner performance and operational conditions in a way to manage heat exchanger tube metal temperatures and temperature ramp rates in the absence of fluid phase change on the inside of the tubes. We designed and built the first ever pseudo power system employing a simple recuperated supercritical CO 2 closed-loop Brayton Cycle coupled to a solid-fuel fired system. Advanced coupled CFD and process modeling were used to design the primary heat exchanger (PHX), which consisted of both radiative and convective sections, to limit tube metal temperatures resulting from the heat release profile of the solid fuel flame near the radiative tubes. The heat exchanger was designed to produce finished CO 2 temperatures of 600 °C a pressure of 20.7 MPa and CO 2 flow of 5.5 kg/s. The constructed PHX was capable of 1.2 MWth heat uptake. During design of the PHX, the modeling showed that most variables influencing flame shape (burner stoichiometric ratio and register velocities and swirl) were not suitable to manage heat flux to the metal surfaces. This is because they substantially increased adiabatic flame temperature through the influence of localized stoichiometric ratio. Excess air and firing rate were the two most powerful variables that could be used to control tube surface temperatures. The coupled system was operated for a total of 407 hours, with the longest continuous run of 248 hours. For 62% of the operational time, the unit was unmanned and in automatic control. The fuels used for the testing included natural gas, two Utah Bituminous coals, woody biomass, and bagasse. During the testing we were able to verify the 1.2 MWth heat uptake and we operated at a finished CO 2 temperature of 607 °C and a pressure of 20.3 MPa simultaneously. The real-time corrosion rate of the Super 304H tube CO 2 surface in the region of the radiative section of the PHX were measured, at an approximate temperature of 550 °C. The two key variables related to corrosion rate are the pressure and flow rate of the CO 2 . A technoeconomic analysis was performed at a scale of 120 MWE. The updated analysis showed that the efficiency of an sCO 2 power producing plant will be related to the pressure drop of the PHX.

01 COAL, LIGNITE, AND PEAT↗