Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “critical cooling rate”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Computational Fluid Dynamics (CFD) applications in rocket propulsion analysis and design

Computational Fluid Dynamics (CFD) has been used in recent applications to affect subcomponent designs in liquid propulsion rocket engines. This paper elucidates three such applications for turbine stage, pump stage, and combustor chamber geometries. Details of these applications include the development of a high turning airfoil for a gas generator (GG) powered, liquid oxygen (LOX) turbopump, single-stage turbine using CFD as an integral part of the design process. CFD application to pump stage design has emphasized analysis of inducers, impellers, and diffuser/volute sections. Improvements in pump stage impeller discharge flow uniformity have been seen through CFD optimization on coarse grid models. In the area of combustor design, recent CFD analysis of a film cooled ablating combustion chamber has been used to quantify the interaction between film cooling rate, chamber wall contraction angle, and geometry and their effects of these quantities on local wall temperature. The results are currently guiding combustion chamber design and coolant flow rate for an upcoming subcomponent test. Critical aspects of successful integration of CFD into the design cycle includes a close-coupling of CFD and design organizations, quick turnaround of parametric analyses once a baseline CFD benchmark has been established, and the use of CFD methodology and approaches that address pertinent design issues. In this latter area, some problem details can be simplified while retaining key physical aspects to maintain analytical integrity.

Mcconnaughey, P. K.↗
NASA NTRS1995

Nonlinear stability of non-stationary cross-flow vortices in compressible boundary layers

The nonlinear evolution of long wavelength non-stationary cross-flow vortices in a compressible boundary layer is investigated and the work extends that of Gajjar (1994) to flows involving multiple critical layers. The basic flow profile considered in this paper is that appropriate for a fully three-dimensional boundary layer with O(1) Mach number and with wall heating or cooling. The governing equations for the evolution of the cross-flow vortex are obtained and some special cases are discussed. One special case includes linear theory where exact analytic expressions for the growth rate of the vortices are obtained. Another special case is a generalization of the Bassom & Gajjar (1988) results for neutral waves to compressible flows. The viscous correction to the growth rate is derived and it is shown how the unsteady nonlinear critical layer structure merges with that for a Haberman type of viscous critical layer.

Gajjar, J. S. B.↗
DOE OSTI2024

Weathering Incongruence in Mountainous Mediterranean Climates Recorded by Stream Lithium Isotope Ratios

Lithium isotope ratios (δ 7 Li) of rivers are increasingly serving as a diagnostic of the balance between chemical and physical weathering contributions to overall landscape denudation rates. Here, we show that intermediate weathering intensities and highly enriched stream δ 7 Li values typically associated with lowland floodplains can also describe small upland watersheds subject to cool, wet climates. This behavior is revealed by stream δ 7 Li between +22.4 and +23.5‰ within a Critical Zone observatory located in the Cévennes region of southern France, where dilute stream solute concentrations and significant atmospheric deposition otherwise mask evidence of incongruence. The water-rock reaction pathways underlying this behavior are quantified through a multicomponent, isotope-enabled reactive transport model. Using geochemical characterization of soil profiles, bedrock, and long-term stream samples as constraints, we evolve the simulation from an initially unweathered granite to a steady state weathering profile which reflects the balance between (a) fluid infiltration and drainage and (b) bedrock uplift and soil erosion. Enriched stream δ 7 Li occurs because Li is strongly incorporated into actively precipitating secondary clay phases beyond what prior laboratory experiments have suggested. Chemical weathering incongruence is maintained despite relatively slow reaction rates and moderate clay accumulation due to a combination of two factors. First, reactive primary mineral phases persist across the weathering profile and effectively “shield” the secondary clays from resolubilization due to their greater solubility. Second, the clays accumulating in the near-surface profile are relatively mature weathering byproducts. These factors promote characteristically low total dissolved solute export from the catchment despite significant input of exogenous dust.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Capillary-Enhanced Two-Phase Micro-Cooler Using Copper-Inverse-Opal Wick with Silicon Microchannel Manifold for High-Heat-Flux Cooling Application

In this work, we demonstrate a two-phase capillary-fed boiling micro-cooler that consists of a ~ 25-..mu..m-thick copper inverse opal (CIO) porous wicking structure for high-heat-flux boiling and a silicon 3D-manifold for distributed liquid delivery and vapor extraction across a 0.5 cm x 0.5 cm heated area. At low inlet water mass flow rates of 1.5 to 1.9 g(min)-1, the micro-cooler displays nearly two-phase boiling with exit vapor quality ~ 1 and a high critical heat flux (CHF) of 253 to 320 W cm-2 with low superheat of ~ 10 degrees C resulting in a thermal resistance of boiling ~ 0.025 cm2 degrees C W-1 or heat transfer coefficient of 0.4 MW m-2 degrees C-1. For higher flow rates of 5, 10, and 15 g(min)-1, the micro-cooler exhibits a hybrid single-phase and two-phase cooling regime where the contribution of the sensible heat (single-phase) cooling is linearly added to that of the two-phase cooling. For the highest flow rate of 15 g(min)-1, the CHF is increased to ~ 500 W cm-2 resulting in an overall thermal resistance of ~ 0.18 cm2 degrees C W-1. However, the two-phase heat transfer effectiveness, which estimates the utilization level of the inlet mass flow rate for two-phase boiling, is reduced to ~ 0.11. To achieve the best cooling system performances, the micro-cooler must operate entirely within the two-phase boiling regime (exit vapor quality or two-phase heat transfer effectiveness ~ 1). Ideally, the "coolant" should be delivered near its saturation temperature (~ 100 degrees C for water), which provides significant advantages for the energy-efficient operation of data centers and power electronics. We present detailed analysis with Infrared and high speed camera images at various inlet flow rates and heat fluxes to understand complex heat transfer in the micro-cooler. Furthermore, a conjugate thermofluidic simulation model, which incorporates the physics of capillary-fed boiling in a porous copper layer, agrees well with the experimental data.

capillary-enhanced boiling↗
DOE OSTI2022

Analysis of precooling optimization for residential buildings

To reduce peak demand and/or energy cost for residential buildings, optimal precooling strategies are becoming important as an alternative to rule-based precooling strategies that are intuitive but may not be optimal. Since precooling optimization is heavily dependent on a variety of factors such as the home thermal properties, HVAC system, weather, thermal comfort criteria, and utility rate structure, the individual and collective impact of these factors on precooling performance needs to be analyzed. In addition, since the indoor air temperature is affected by heat transfer to and from the interior wall surface, performance analysis in view of the interior wall surface temperature is also essential. Therefore, in this paper, an optimal precooling strategy that accounts for the aforementioned factors and utilizes a second-order thermal network model, is proposed. With this strategy, the HVAC on/off control signal that minimizes 24-hour energy cost while maintaining thermal comfort, is determined. Through extensive simulations, it is found that the proposed optimal precooling strategy is able to adapt to changing conditions and that having a sufficiently low interior wall surface temperature during precooling is critical for avoiding expensive on-peak operation. Here, the reason for the latter is that such a temperature indicates that enough “cooling energy” has been stored. It is also found that weather has the most dominant impact on the precooling performance, followed by home thermal condition, with the rated cooling capacity and utility rate structure having the least impact.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Analysis of Conduction Cooling Strategies for Wire Arc Additive Manufacturing

Metal additive manufacturing (AM) processing consists of numerous parameters which take time to optimize for various geometries. One aspect of the metal AM process that continues to be explored is the control of thermal energy accumulation during component manufacturing due to the melting and solidification of the feedstock. Excessive energy accumulation causes thermal failure of the component while minimal energy accumulation causes lack of fusion with the build plate or previous layer. The ability to simulate the thermal response of an AM component can increase research efficiency by reducing the time to optimize thermal energy accumulation. This paper presents an effective implementation of finite element analysis to determine the thermal response of a wire arc additive manufactured component with various build plate sizes and cooling methods including, integral build plate cooling, oversized build plates with passive cooling, and non-integral build plate cooling. The use of integral build plate cooling channels was shown to decrease the interpass temperature at the conclusion of the build process by 55% and build plate temperature by 96% compared to the conventionally deposited sample with 20 second dwell time. The use of a tall build plate with passive cooling was shown to reduce the interpass temperature by 32% as compared to the conventionally deposited sample with 20 second dwell time. Each cooling strategy evaluated decreased the interpass temperature within a range of 20–55% which enables higher deposition rates and decreased dwell times during depositions. The cooling strategies are designed to be implemented in a hybrid or retrofit AM platform to mitigate concerns of the thermal input from the additive process having detrimental effects on the precision of the machining process. This paper shows that accurate simulations of all strategies can be used to accurately predict the thermal response of the various strategies discussed. These cooling strategies will allow for increased deposition rates with comparable interpass temperature and decreased dwell time, increasing deposition efficiency. This model and these simulations are verified by experimental results. It is concluded that passive strategies, such as the over-sized tall build plate, can be used when liquid coolant in the AM environment could negatively affect the deposition process. Active cooling strategies, such as the integral build plate cooling could be used if low thermal conductivity materials are deposited or higher material deposition rates are desired. This paper discusses the use of active and passive cooling used during AM and shows how a simulation model can be used to make design choices for cooling strategies. The model also enables verification of select critical process parameters such as dwell times for a desired interpass temperature.

Heinrich, Lauren↗
DOE OSTI2021

The Fundamental Creep Behavior Model of Gr.91 Alloy by ICME Approach (Final Report)

The current project mainly focused on the fundamental creep cracking mechanism of the Grade 91 system during the operation conditions in the advanced technologies FE power plants and build the link among Composition-Processing Parameters-Phase Stability-Microstructure-Creep Resistance. A model based on computational thermodynamics and diffusion kinetics will be developed to provide guidance on how to improve the creep resistance of the alloy system. The PI mainly investigate the Gr.91 base alloy and weldment with the Integrated Computational Materials Engineering (ICME) approach. The long-term goal of the proposed program is to develop a model for creep and fatigue resistant alloys with different elemental systems, compositions, and processing parameters of weldment and heat-treatment. In addition, new creep and fatigue resistant alloys are to be designed based on the predictions from the model. The computational as well as experimental results include three sections: 1) The compositional optimization and secondary phases evaluation regarding the creep resistance in Grade 91 steel through the CALPHAD approach. In this section, the formation of the critical secondary phases, M23C6, MX and Z-phase were well predicted based on many conditions like the temperature as well as compositions. Meanwhile, the role of different alloying elements was also considered as they will directly affect the formation of the secondary phases mentioned above based on the compositional evolution of the alloying elements. 2) An investigation of Creep Resistance in Grade 91 Steel through Computational Thermodynamics. We systematically considered the evolution of the four critical temperatures, Ac3, Ac1, the threshold of M23C6 and Z-phase, on basis of the computational approach. Meanwhile, the equilibrium cooling as well as Scheil simulations were also considered to perform further predictions based on the different cooling rates. 3) Creep lifetime and microstructural evolution of the Grade 91 steel. In this part, mainly experimental approach was performed on the real creep test from the beginning of the commercial Gr.91 alloys. The alloys were normalized, tempered, welded (Gleeble) and PWHT and finally subjected to the creep lifetime analyses. Finally, all the samples after each process were further characterized based on the OM, XRD and SEM technique to draw final conclusions.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Compatibility of LaFe13−x−yMnxSiyH1.6 and Eutectic Liquid GaInSn Alloy

The heat transfer rate of magnetocaloric regenerators is a topic of extensive research and the cyclability of these regenerators is critical to the operation of systems with a high coefficient of performance (e.g., potentially >22, significantly higher than typical vapor compression cooling technologies). To enable a high operating frequency that will result in a high specific cooling power, the heat transfer fluid should have high thermal conductivity and lower specific heat, i.e., higher thermal diffusivity. Eutectic metal alloys possess these qualities, such as gallium–indium–tin (Galinstan), whose thermal diffusivity has been found to be approximately an order of magnitude higher than water. For this study, the effects of eutectic liquid Galinstan exposure on the phase stability of LaFe13−x−yMnxSiyH1.6 magnetocaloric powders in an active magnetic regenerator device were investigated. The powders were characterized before and after exposure to Galinstan using X-ray diffraction, in which the phases were determined using the Rietveld refinement technique and X-ray fluorescence. It was found that after Galinstan exposure, hydrogen containing phases were present in the powder, suggesting that the hydrogen was lost from the magnetocaloric phase. The magnetocaloric phase degradation indicates that the powder was incompatible with the Galinstan metal in an environment with moisture.

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

Indicators of Global Climate Change 2023: annual update of key indicators of the state of the climate system and human influence

Intergovernmental Panel on Climate Change (IPCC) assessments are the trusted source of scientific evidence for climate negotiations taking place under the United Nations Framework Convention on Climate Change (UNFCCC). Evidence-based decision-making needs to be informed by up-to-date and timely information on key indicators of the state of the climate system and of the human influence on the global climate system. However, successive IPCC reports are published at intervals of 5–10 years, creating potential for an information gap between report cycles. We follow methods as close as possible to those used in the IPCC Sixth Assessment Report (AR6) Working Group One (WGI) report. We compile monitoring datasets to produce estimates for key climate indicators related to forcing of the climate system: emissions of greenhouse gases and short-lived climate forcers, greenhouse gas concentrations, radiative forcing, the Earth's energy imbalance, surface temperature changes, warming attributed to human activities, the remaining carbon budget, and estimates of global temperature extremes. The purpose of this effort, grounded in an open-data, open-science approach, is to make annually updated reliable global climate indicators available in the public domain. As they are traceable to IPCC report methods, they can be trusted by all parties involved in UNFCCC negotiations and help convey wider understanding of the latest knowledge of the climate system and its direction of travel. The indicators show that, for the 2014–2023 decade average, observed warming was 1.19 [1.06 to 1.30] °C, of which 1.19 [1.0 to 1.4] °C was human-induced. For the single-year average, human-induced warming reached 1.31 [1.1 to 1.7] °C in 2023 relative to 1850–1900. The best estimate is below the 2023-observed warming record of 1.43 [1.32 to 1.53] °C, indicating a substantial contribution of internal variability in the 2023 record. Human-induced warming has been increasing at a rate that is unprecedented in the instrumental record, reaching 0.26 [0.2–0.4] °C per decade over 2014–2023. This high rate of warming is caused by a combination of net greenhouse gas emissions being at a persistent high of 53±5.4 Gt CO 2 e yr -1 over the last decade, as well as reductions in the strength of aerosol cooling. Despite this, there is evidence that the rate of increase in CO 2 emissions over the last decade has slowed compared to the 2000s, and depending on societal choices, a continued series of these annual updates over the critical 2020s decade could track a change of direction for some of the indicators presented here.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2017

Updates on Modeling the Water Cycle with the NASA Ames Mars Global Climate Model

Global Circulation Models (GCMs) have made steady progress in simulating the current Mars water cycle. It is now widely recognized that clouds are a critical component that can significantly affect the nature of the simulated water cycle. Two processes in particular are key to implementing clouds in a GCM: the microphysical processes of formation and dissipation, and their radiative effects on heating/ cooling rates. Together, these processes alter the thermal structure, change the dynamics, and regulate inter-hemispheric transport. We have made considerable progress representing these processes in the NASA Ames GCM, particularly in the presence of radiatively active water ice clouds. We present the current state of our group's water cycle modeling efforts, show results from selected simulations, highlight some of the issues, and discuss avenues for further investigation.­

Global Circulation Models (GCMs)↗
NASA NTRS1994

The hard X-ray emission spectra from accretion columns in intermediate polars

We consider the hard (greater than 2 keV) X-ray emission from accretion columns in an intermediate polar system, GK Per, using a simple settling solution. The rate of photon emission per logarithmic energy interval can be fitted with a power law, E(exp -gamma), with gamma approximately 2.0, in agreement with observations. This index is only weakly dependent on the mass accretion rate, dot-M, for dot-M in the range of a few times 10(exp 16-18) g/s. The peak energy of the photon spectra (after photoelectric absorption) is expected to be E(sub p) approximately (5 keV) gamma(exp -1/3) (N(sub H)/10(exp 23)/sq cm)(exp 1/3) where N(sub H) is the hydrogen column density along the line of sight. The observed spectra of GK Per and possibly of V1223 Sgr suggest N(sub H) approximately 10(exp 23)/sq cm. This large N(sub H) may be due to partially ionized preshock column material. Alternatively, we also consider absorption by the cool outer parts of an accretion disk. In this case the photoelectric absorption depth in the disk is a sensitive function of inclination. For GK Per the required inclination is approximately 83 deg. For mass accretion rates larger than a critical rate of approximately 10(exp 18) g/s, X-ray emission from the column accretion is significantly affected by radiation drag. Although the mass accretion rate increases dramatically during outbursts, the observed hard (greater than 2 keV) X-ray luminosity will not rise proportionately. The slope and peak energy of the outburst spectra are only weakly affected. We conclude that the observed X-ray spectra can be explained by this simple analytic solution and that the production of hard X-rays from the accretion shock at the magnetic poles in the intermediate polars is in general agreement with the observations. However, since the X-ray emission and absorption depend on the mass accretion rate in a complicated manner, observed hard X-ray luminosities (greater than 2 keV) are not a good indicator of the mass accretion rate in the X-ray emitting regions.

Yi, Insu↗
DOE OSTI2023

County-level assessment of behind-the-meter solar and storage to mitigate long duration power interruptions for residential customers

Customer concerns over electric system resilience could drive early adoption of behind-the-meter solar-plus-storage (BTM PVESS), especially as wildfire, hurricane, and other climate-driven risks to electric grids become more pronounced. However, the resilience benefits of BTM PVESS are poorly understood, especially for residential customers, owing to lack of data and methodological challenges, making it difficult to forecast adoption trends. In this paper, we develop a methodology to model the performance of BTM PVESS in providing backup power across a wide range of customer types, geography / climate conditions, and long duration power interruption scenarios, considering both whole-building backup and backup of specific critical loads. We combine novel, disaggregated end-use load profiles across the continental United States with temporally and geospatially aligned solar generation estimates. We then implement a PVESS dispatch algorithm to calculate the amount of load served during interruptions. We find that PVESS with 10 kWh of storage can meet a limited set of critical loads in most United States counties during any month of the year, though this capability drops to meeting only 86% of critical load, averaged across all counties and months, when heating and cooling are considered critical. Backup performance is lowest in winter months where electric heat is common (southeast and northwest U.S.) and in summer months in places with large cooling loads (southwest and southeast U.S.). Winter backup performance varies by roughly 20% depending on infiltration rates, while summer performance varies by close to 15% depending on the efficiency of the central air-conditioning system. Differences in temperature set-points in Harris County correspond to a 40% range in winter backup performance and a 20% range in summer performance. Economic calculations show that a customer’s resilience value of PVESS must be high to motivate adoption of these systems.

14 SOLAR ENERGY↗
DOE OSTI2025

Evaluation of PBR Spent Fuel Criticality and Dose Rate Compliance for Storage and Transportation

Spent tri-structural isotropic (TRISO)–based fuels have a strong track record in storage and transportation without documented incidents. This work seeks to reduce uncertainty to aid in more informed spent fuel management of TRISO-based fuels by modeling both fresh and spent pebble bed reactor (PBR) fuel and comparing the results to the regulatory standards from 10 CFR 71. SCALE was used for all modeling due to it having fast and accurate methods for handling PBR fuel modeling, as well as having an efficient method for shielding calculations in monaco with automated variance reduction using importance calculations (MAVRIC), which utilizes the consistent adjoint-driven importance sampling (CADIS) and the forward-weighted consistent adjoint-driven importance sampling (FW-CADIS) methods. KENO-VI was used for all criticality calculations, TSUNAMI was used for uncertainty quantification on k-effective, TRITON and the Oak Ridge isotope generation code (ORIGEN) were both used for depletion of the fuel, and MAVRIC was used for shielding calculations. For criticality assessments, this study focused on the requirement that the value of the neutron multiplication factor, k-effective (k-eff), would not exceed a peak value of 0.95, including uncertainty, with 95% confidence. Criticality was initially examined by modeling fresh fuel from three different designs—HTR-10 fuel, PBMR-400 fuel, and demonstration fuel representative of a TRISO-fueled modern high-temperature gas reactor (HTGR) design, henceforth referred to as Demo HTGR—and placing them into various sized containers with conditions described in 10 CFR 71 to quantify the peak k-eff state. When the peak value of 0.95 k-eff was exceeded, mitigation methods were examined in those scenarios. Burnup credit, pebble displacement in areas of strong neutron multiplication, and random pebble replacement using pebbles of various compositions and replacement fractions were examined. In summary, the criticality of PBR fuels can be well accounted for by restricting container size, taking credit for burnup, or by displacing/replacing pebbles. Uncertainty of the k-eff due to nuclear data uncertainties was recorded at ~0.6644%Δk/k, or roughly 664% mil (pcm). The nuclear data–induced uncertainty was relatively small and should not require significant modification in the design to be accounted for. Revisions to the evaluated nuclear data file values have been shown to have a larger impact than nuclear data–induced uncertainty. For dose rate aspects, U.S. Nuclear Regulatory Commission regulations require a maximum dose rate of 10 millirem per hour (mrem/h) at 2 meters. In examining the dose rate behavior of spent PBR fuel, the representative Demo HTGR fuel was modeled exclusively due to it possessing the highest target burnup of the examined fuels. Equilibrium cycle modeling methods were used to produce a higher-fidelity discharge isotopic composition than simple assumptions, such as reflected pebbles. The discharge composition was used as a source term in the fixed-source transport shielding calculations, and dose rates were calculated at 2 m for the shortest possible cooling time. The low concentration of fuel material led to dose rates that were in line with regulatory limits, despite the high burnup when compared to traditional light water reactor fuels. In conclusion, the methods employed in this study would require more work to further verify and validate and are limited to the criticality and dose rate analyses performed.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Film condensation with high heat fluxes and scaled experiments using pure steam for reactor containment cooling

Condensation tests were performed using a newly developed test facility for scaling the passive containment cooling system (PCCS) to a small modular reactor (SMR). The PCCS of the SMR plays a pivotal role in ensuring greater safety, reliability, and compactness than what is afforded by traditional reactors. Therefore, a well-designed PCCS is essential to SMRs. However, previous studies and test data were unsuitable for scaling, due to high variation in the test geometry and operating conditions. This study intends to close this research gap by using a novel designed scaled test facility consisting of vertical condensing test sections featuring 1-, 2-, and 4-inch-diameter condensing tubes with annular water cooling, and by applying superheated and saturated steam with different steam mass flow ranges of 5–25 g/s. Further, the primary test data, including axial temperatures, mass flow rates, and pressures, were used in conjunction with a standard data reduction method to estimate critical parameters such as heat fluxes, heat transfer coefficients, and condensation rates. These scaled test data would support improving empirical correlations and validating condensation models to identify scaling distortion for SMR PCCSs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2000

Thermocapillary Flow and Aggregation of Bubbles on a Solid Wall

During the electrolytic evolution of oxygen bubbles forming on a vertically oriented transparent tin oxide electrode, bubbles were found to be mutually attractive. The mechanism of the aggregation had never been explained satisfactorily until Guelcher et al. attributed it to thermocapillary flow. The gradient of surface tension of the liquid at the bubble's surface, which was established because of reaction heat and ohmic heat loss at the electrode wall, drives flow of the liquid adjacent to each bubble; the bubble "pumps" fluid along its surface away from the wall. Fluid flows toward the bubble to conserve mass and entrains nearby bubbles in the flow pattern. The same logic would apply when two bubbles of equal size are adjacent to each other on a warm wall. Each bubble drives thermocapillary flow and hence entrains the other in its flow pattern, which drives the aggregation. Our objective here is to perform experiments where the temperature gradient at the wall is well known and controlled. The theory can be quantitatively tested by studying aggregation of bubble pairs of equal size, and by varying system parameters such as temperature gradient, bubble size and fluid viscosity. The results are then compared with the theory in a quantitatively rigorous manner. We demonstrate that the theory without adjustable parameters is capable of quantitatively modeling the rate of aggregation of two bubbles. The equations governing the thermocapillary flow around a single stationary bubble on a heated or cooled wall in a semi-infinite domain were solved. Both Reynolds number and Marangoni number were much less than unity. The critical result is that liquid in the vicinity of a warm wall flows toward a stationary collector bubble. Consequently the thermocapillary flow around the stationary bubble entrains another bubble toward itself. The bubbles undergo hindered translation parallel to the wall with velocity U while the fluid flow field is described with u. Two velocities were equated by using a wall hindrance parameter q: U = qu [1] which shows the velocity of bubble is proportional to the entraining velocity. The hindrance parameter q can experimentally be measured independently. q can also be calculated by solving the equations of motion for a bubble translating parallel to a solid wall. The experimental cell is cylindrical with an ID of 10 cm and consists of a 1 cm deep main cell filled with silicone oil and flanked by two thermal reservoirs. The upper thermal reservoir was heated and the lower thermal reservoir was cooled so that the bubbles aggregate. Two types of silicone oil (eta = 0.02 and 0.50 Pa s) were used. Two equal sized air bubbles were injected into the cell with a syringe. The center-to-center distance of bubbles was observed through a microscope. Bubble radius ranged from 0.40 mm to 0.65 mm and the temperature gradients along with the cell ranged from 1400 to 5000 K/m. The bubbles aggregated when heat flows from the wall to the fluid. The velocities of bubbles were in the range of 1 - 10 microns/s. The separation r decreased more quickly when the temperature gradient was higher, bubble size was larger, and the oil viscosity was lower. r decreased more rapidly as the bubbles approached each other. Dimensionless time was arbitrarily set to be zero when the dimensionless center-to-center distance between the bubbles was 4. All the bubble trajectories fall onto one line, especially in the range of dimensionless distance from 4 to 3. This means the relative movement of the bubble pair is proportional to the temperature gradient and bubble size and it is inversely proportional to the viscosity of the oil. This result strongly suggests that the thermocapillary flow-based aggregation mechanism is correct. A value of q can be estimated by fitting the scaled data to Eq. [1]. A best fit value of q was obtained as q = 0.26 with a standard deviation of 0.03. Independent experimental results for q for a 0.5 mm radius bubble, give values of q in the range 0.11 to 0.23. The value of q obtained from solving the equations of motion reveals q has values in the range 0.23 - 0.30. Since the full scale of possible values of q is zero to one, the maximum deviation of independently determined values of q from the best fit value was 15% of this full scale. Thus reasonable quantitative agreement between theory and experiment has been obtained.

Kasumi, Hiroki↗
DOE OSTI2022

Thermo-economic analysis of a solid oxide fuel cell-gas turbine hybrid with commercial off-the-shelf gas turbine

The thermodynamic and economic performance of a natural gas-fueled solid oxide fuel cell (SOFC)-gas turbine (GT) hybrid system with a commercial off-the-shelf GT is investigated. Until today, commercial GTs are primarily engineered for the direct use of energy dense fuels, such as natural gas, and the integration of commercial GTs into a SOFC-GT hybrid remains challenging. In this work technically feasible and economically viable GT operating modes are identified to gauge the technology’s competitiveness on the free market. Steady state, full load operating conditions were established for various GT operating modes: I) constant spool speed operation, II) variable spool speed operation and III) partially closed compressor inlet guide vanes. For each GT operating mode, the performance was investigated over a range of SOFC fuel utilization factors, while considering physical constraints inside the SOFC, such as local temperature gradients in flow direction as well as overall cell temperature differences. The results of the thermodynamic evaluation served as inputs for the economic analysis of the SOFC-GT hybrid power plant. The results show that the integration of an off-the-shelf GT not only results in a significant derating of the GT, but also substantially impacts the SOFC operation, the main power producer in this SOFC-GT hybrid. To maximize the SOFC power output it is desirable operate the GT in a region of high air mass flow rates and low pressure ratios, which increases the number of stacks that can be accommodated and reduces the SOFC cooling requirement. The lowest costs of electricity are obtained at constant spool speed operation, while the highest efficiencies are reached at variable spool speed operation. Critical for the integration of off-the-shelf GTs, which historically have been designed for natural gas, is the surge margin. The largest surge margins are obtained by closing the compressor inlet guide vanes. Furthermore, higher SOFC fuel utilization factors are shown to increase the surge margin as the turbine firing temperature is decreased.

25 ENERGY STORAGE↗
NASA NTRS2004

A New Test Method for Material Flammability Assessment in Microgravity and Extraterrestrial Environments

The standard oxygen consumption (cone) calorimeter (described in ASTM E 1354 and NASA STD 6001 Test 2) is modified to provide a bench-scale test environment that simulates the low velocity buoyant or ventilation flow generated by or around a burning surface in a spacecraft or extraterrestrial gravity level. The Equivalent Low Stretch Apparatus (ELSA) uses an inverted cone geometry with the sample burning in a ceiling fire (stagnation flow) configuration. For a fixed radiant flux, ignition delay times for characterization material PMMA are shown to decrease by a factor of three at low stretch, demonstrating that ignition delay times determined from normal cone tests significantly underestimate the risk in microgravity. The critical heat flux for ignition is found to be lowered at low stretch as the convective cooling is reduced. At the limit of no stretch, any heat flux that exceeds the surface radiative loss at the surface ignition temperature is sufficient for ignition. Regression rates for PMMA increase with heat flux and stretch rate, but regression rates are much more sensitive to heat flux at the low stretch rates, where a modest increase in heat flux of 25 kW/m2 increases the burning rates by an order of magnitude. The global equivalence ratio of these flames is very fuel rich, and the quantity of CO produced in this configuration is significantly higher than standard cone tests. These results [2] demonstrate the ELSA apparatus allows us to conduct normal gravity experiments that accurately and quantifiably evaluate a material s flammability characteristics in the real-use environment of spacecraft or extra-terrestrial gravitational acceleration. These results also demonstrate that current NASA STD 6001 Test 2 (standard cone) is not conservative since it evaluates materials flammability with a much higher inherent buoyant convective flow.

Olson, S. L.↗
DOE OSTI2022

Evaluation of mitigation strategies for radioactive fission gases in fluoride-salt-cooled high-temperature reactors

Here, the Fluoride-salt-cooled High-temperature Reactor (FHR) is a promising Generation IV nuclear reactor due to its improved safety design, near atmospheric working pressure, and high electric power conversion efficiency. However, control of radioactive fission gases, such as tritium, is a critical issue in FHRs due to the significantly larger production rate and potentially larger leakage rate of tritium compared to those in Light Water Reactors (LWRs). Therefore, this study proposes two mitigation strategies for tritium in FHRs: Double-Wall Heat eXchanger with a Tritium Carrier (DWHX-TC) option and Single-Wall Heat eXchanger with a Tritium Barrier (SWHX-TB) option. For the DWHX-TC option, two tube surface configurations, i.e., plain and fluted tubes, and four tritium carriers, i.e., helium, FLiBe, FLiNaK, and KF-ZrF 4 , are initially considered, among which the best tube configuration and promising tritium carrier are then selected and used for design optimization of the DWHX. For the SWHX-TB option, two cases, tritium barriers with and without cracks, are investigated in this study. In addition, the two mitigation strategies are evaluated by an in-house, one-dimensional, coupled HEat and MAss Transfer (HEMAT) code developed to predict the leakage rate of tritium in nuclear systems. Our analysis shows that (1) a double-wall configuration using fluted tubes as both the inner and outer tubes is better than other double-wall configurations investigated; (2) helium, as a tritium carrier, is superior to other candidates, including FLiBe, FLiNaK, and KF-ZrF 4 ; (3) the leakage rate of tritium (into the atmosphere) is reduced from 3400 Ci/day, taking the Advanced High-Temperature Reactor (AHTR) as a reference, to 4.0 Ci/day for an optimum design of the DWHX-TC option using helium as the tritium carrier and 1.7 Ci/day for an optimum design of the SWHX-TB option using 30-µm thick silicon carbide (SiC) as the tritium barrier (no cracks); and (4) potential cracks in tritium barriers may significantly deteriorate mass transfer performance of the SWHX-TB option. For example, a cracking rate of 0.1% leads to an increase of 694.1% in the leakage rate under the conditions investigated.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗