Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Critical Fluids”

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 451 records · Page 25

A Study of Nucleate Boiling with Forced Convection in Microgravity

The ultimate objective of basic studies of flow boiling in microgravity is to improve the understanding of the processes involved, as manifested by the ability to predict its behavior. This is not yet the case for boiling heat transfer even in earth gravity, despite the considerable research activity over the past 30 years. The elements that constitute the nucleate boiling process - nucleation, growth, motion, and collapse of the vapor bubbles (if the bulk liquid is subcooled) - are common to both pool and flow boiling. It is well known that the imposition of bulk liquid motion affects the vapor bubble behavior relative to pool boiling, but does not appear to significantly influence the heat transfer. Indeed, it has been recommended in the past that empirical correlations or experimental data of pool boiling be used for design purposes with forced convection nucleate boiling. It is anticipated that such will most certainly not be possible for boiling in microgravity, based on observations made with pool boiling in microgravity. In earth gravity buoyancy will act to remove the vapor bubbles from the vicinity of the heater surface regardless of how much the imposed bulk velocity is reduced, depending, of course, on the geometry of the system. Vapor bubbles have been observed to dramatically increase in size in pool boiling in microgravity, and the heat flux at which dryout took place was reduced considerably below what is generally termed the critical heat flux (CHF) in earth gravity, depending on the bulk liquid subcooling. However, at heat flux levels below dryout, the nucleate pool boiling process was enhanced considerably over that in earth gravity, in spite of the large vapor bubbles formed in microgravity and perhaps as a consequence. These large vapor bubbles tended to remain in the vicinity of the heater surface, and the enhanced heat transfer appeared to be associated with the presence of what variously has been referred to as a liquid microlayer between the bubble and the heater surface. The enhancement of the boiling process with low velocities in earth gravity for those orientations producing the formation of a liquid macrolayer described above, accompanied by "sliding" vapor bubbles, has been demonstrated. The enhancement was presented as a function of orientation, subcooling, and heated length, while a criterion for the heat transfer for mixed natural/forced convection nucleate boiling was given previously. A major unknown in the prediction and application of flow boiling heat transfer in microgravity is the upper limit of the heat flux for the onset of dryout (or critical heat flux - CHF), for given conditions of fluid-heater surfaces, including geometry, system pressure and bulk liquid subcooling. It is clearly understood that the behavior in microgravity will be no different than on earth with sufficiently high flow velocities, and would require no space experimentation. However, the boundary at which this takes place is still an unknown. Previous results of CHF measurements were presented for low velocity flow boiling at various orientations in earth gravity as a function of flow velocity and bulk liquid subcooling, along with preliminary measurements of bubble residence times on a flat heater surface. This showed promise as a parameter to be used in modeling the CHF, both in earth gravity and in microgravity. The objective of the work here is to draw attention to and show results of current modeling efforts for the CHF, with low velocities in earth gravity at different orientations and subcoolings. Many geometrical possibilities for a heater surface exist in flowing boiling, with boiling on the inner and outer surfaces of tubes perhaps being the most common. If the vapor bubble residence time on and departure size from the heater surface bear a relationship to the CHF, as results to be given indicate, it is important that visualization of and access to vapor bubble growth be conveniently available for research purposes. In addition, it is desirable to reduce the number of variables as much as possible in a fundamental study. These considerations dictated the use of a flat heater surface, which is rectangular in shape, 1.91 cm by 3.81 cm (0.75 x 1.5 inches), consisting either of a 400 Angstrom thick semi-transparent gold film sputtered on a quartz substrate which serves simultaneously as a heater and a resistance thermometer, or a copper substrate of the same size. The heater substrate is a disc which can be rotated so that the heated length in the flow direction can be changed from 1.91 to 3.81 cm (0.75 to 1.5 inches). The fluid is R-113, and the velocities can be varied between 0.5 cm/s and 60 cm/s. For a sufficiently low velocity the CHF can be modeled reasonably well at various orientations by the correlation for pool boiling corrected for the influence of bulk liquid subcooling, multiplied by the square root of q, the angle relative to horizontal. This arises from equating buoyancy and drag forces in the inverted positions where the vapor bubbles are held against the heater surface as they slide. A distortion of the measurements relative to pool boiling occurs as the flow velocity increases. In modeling this effect at different levels of subcooling it appeared appropriate to estimate the volumetric rate of vapor generation, using measurements of bubble frequency (or residence time), void fraction and average bubble boundary layer thickness. These were determined with the use of a platinum hot wire probe 0.025 mm in diameter by 1.3 mm long, applying a constant current to distinguish between contact with liquid or vapor. Two-dimensional spatial variations are obtained with a special mechanism to resolve displacements in increments of 0.025 mm. From such measurements it was determined that the fraction of the surface heat transfer resulting in evaporation varies inversely with the subcooling correction factor for the CHF. The measured inverse bubble residence time is normalized relative to that predicted for an infinite horizontal flat plate at the CHF, and is correlated well with the CHF normalized relative to that for pool boiling, for various orientation angles and subcooling levels. This correspondence is then combined with a normalizing factor for the energy flux leaving the heater surface at the CHF and the computed bubble radius at departure, determined from the balance between the outward velocity of the interface due to evaporation and the buoyance induced velocity of the center of mass of the bubble. The product of the CHF and the corresponding residence time was determined to be a constant for all orientations at a given bulk flow velocity and liquid subcooling, and must be determined empirically for each velocity and subcooling at present. It then becomes possible to predict the CHF for the different orientations, velocities, and subcoolings. These are compared with normalized measurements of the CHF for velocities ranging from 4 cm/s to 55 cm/s, subcoolings from 2.8 to 22.2 K, over orientations angles of 360 degrees.

Merte, Herman, Jr.↗

Evaluation of an Interferometric Sensor for In-Space Detection of Gas Leaks

Space mission planning often involves long-term storage of volatile liquids or high-pressure gases. These may include cryogenic fuels and oxidizers, high-pressure gases, and life-support-critical consumables. The risk associated with the storage of fluids and gases in space systems has long been an issue and the ability to retain these fluids is often tied to mission success. A leak in the storage or distribution system can cause many different problems, including a simple, but mission endangering, loss of inventory or, in severe cases, unbalanced thrust loads on a flight vehicle. Cryogenic propellants are especially difficult to store, especially over a long duration. The propellant can boil off and be lost through the insulating walls of the tank or simple thermal cycling of the fittings, valves, and propellant feed lines may unseat seals allowing the fluid to escape. Current NASA missions call for long-duration in-space storage of propellants, oxidizers, and life support supplies. Leaks of a scale detectable through a pressure drop in the storage tank are often catastrophic and have long been the focus of ground-based mitigation efforts where redundant systems are often employed. However, there is presently no technology available for detecting and monitoring low-level, but still mission-endangering, gas leaks in space. Standard in-space gas detection methods either have a very limited pressure range over which they operate effectively or are limited to certain gases. Mass spectrometer systems are able to perform the detection tasks, but their size, mass and use of high voltage, which could potentially lead to an arc that ignites a combustible propellent, severely limit their usefulness in a space system. In this paper, we present results from testing of the light-based interferometric gas monitoring and leak detection sensor shown in Fig. 1. The output of the sensor is an interference fringe pattern that is a function of the gas density, and commensurate index of refraction, in the sample region. Changes in the density of gas cause the interference fringes to move across a photodiode detector, providing a temporal history of the leak. The sensor is fiber coupled and constructed from solid optics, allowing for placement almost anywhere on the spacecraft. It is also advantageous in that it consumes very little power and does not introduce an ignition source. Data are presented demonstrating the capability of the sensor to measure density variations in different gas species. In addition, the transient response of the sensor in vacuum is demonstrated. These data extend and improve upon the results previously presented by the authors in Ref. [1].

Polzin, Kurt A.↗

Green Liquid Monopropellant Thruster

Physical Sciences, Inc. (PSI), and Orbital Technologies Corporation (ORBITEC) are developing a unique chemical propulsion system for next-generation NASA science spacecraft and missions. The system is compact, lightweight, and can operate with high reliability over extended periods of time and under a wide range of thermal environments. The system uses a new storable, low-toxicity liquid monopropellant as its working fluid. In Phase I, the team demonstrated experimentally the critical ignition and combustion processes for the propellant and used the data to develop thruster design concepts. In Phase II, the team developed and demonstrated in the laboratory a proof-of-concept prototype thruster. A Phase III project is envisioned to develop a full-scale protoflight propulsion system applicable to a class of NASA missions.

Joshi, Prakash B.↗

Optimized thermodynamic properties of REE aqueous species (REE 3+ and REEOH 2+ ) and experimental database for modeling the solubility of REE phosphate minerals (monazite, xenotime, and rhabdophane) from 25 to 300 °C

Rare earth elements (REE) are critical elements found in monazite, xenotime, and hydrated REE phosphates which typically form in hydrothermal mineral deposits. Accurate predictions of the solubility of these REE phosphates and the speciation of REE in aqueous fluids are both key to understanding the controls on the transport, fractionation, and deposition of REE in natural systems. Previous monazite and xenotime solubility experiments indicate the presence of large discrepancies between experimentally derived solubility constants versus calculated solubilities by combining different data sources for the thermodynamic properties of minerals and aqueous species at hydrothermal conditions. In this study, these discrepancies were resolved by using the program GEMSFITS to optimize the standard partial molal Gibbs energy of formation (Δ f G° 298 ) of REE aqueous species (REE 3+ and REE hydroxyl complexes) at 298.15 K and 1 bar while keeping the thermodynamic properties fixed for the REE phosphates. A comprehensive experimental database was compiled using solubility data available between 25 and 300 °C. The latter permits conducting thermodynamic parameter optimization of Δ f G° 298 for REE aqueous species. Optimal matching of the rhabdophane solubility data between 25 and 100 °C requires modifying the Δ f G° 298 values of REE 3+ by 1–6 kJ/mol, whereas matching of the monazite solubility data between 100 and 300 °C requires modifying the Δ f G° 298 values of both REE 3+ and REEOH 2+ by ~15–31 kJ/mol and ~2–10 kJ/mol, respectively. For xenotime, adjustments of Δ f G° 298 values by 1–26 kJ/mol are only necessary for the REE 3+ species. The optimizations indicate that the solubility of monazite in acidic solutions is controlled by the light (L)REE 3+ species at <150 °C and the LREEOH 2+ species at >150 °C, whereas the solubility of xenotime is controlled by the heavy (H)REE 3+ species between 25 and 300 °C. Based on the optimization results, we conclude that the revised Helgeson-Kirkham-Flowers equation of state does not reliably predict the thermodynamic properties of REE 3+ , REEOH 2+ , and likely other REE hydroxyl species at hydrothermal conditions. We therefore provide an experimental database (ThermoExp_REE) as a basic framework for future updates, extensions with other ligands, and optimizations as new experimental REE data become available. As a result, the optimized thermodynamic properties of aqueous species and minerals are available open access to accurately predict the solubility of REE phosphates in fluid-rock systems.

58 GEOSCIENCES↗

Prediction of Self-Diffusion in Binary Fluid Mixtures Using Artificial Neural Networks

Artificial neural networks (ANNs) were developed to accurately predict the self-diffusion constants for individual components in binary fluid mixtures. The ANNs were tested on an experimental database of 4328 self-diffusion constants from 131 mixtures containing 75 unique compounds. The presence of strong hydrogen bonding molecules may lead to clustering or dimerization resulting in non-linear diffusive behavior. To address this, self- and binary association energies were calculated for each molecule and mixture to provide information on intermolecular interaction strength and were used as input features to the ANN. An accurate, generalized ANN model was developed with an overall average absolute deviation of 4.1%. Forward input feature selection reveals the importance of critical properties and self-association energies along with other fluid properties. Additional ANNs were developed with subsets of the full input feature set to further investigate the impact of various properties on model performance. The results from two specific mixtures are discussed in additional detail: one providing an example of strong hydrogen bonding and the other an example of extreme pressure changes, with the ANN models predicting self-diffusion well in both cases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Critical exponent for the viscosity of carbon dioxide and xenon

The viscosities of carbon dioxide and xenon have been measured near their critical points and the critical exponent y characterizing the asymptotic divergence has been determined. Both fluids yielded exponents in the range y = 0.041 + or - 0.001 and thus also fell in the range y = 0.042 + or - 0.002 from an earlier study of four binary liquids. This agreement between experiments is the first evidence that pure fluids and binary liquids are in the same dynamic universality class. A recent theoretical value for y is 0.032. The 30 percent discrepancy is much greater than the combined errors from experiment and theory. The torsion oscillator viscometer operated at low frequency and low shear rate to avoid systematic errors caused by critical slowing down. Far from T(c) the analysis accounted for the crossover from critical to noncritical temperature dependence, where the latter was obtained from previously published correlations. Corrections for gravitational stratification were included close to T(c).

Berg, R. F.↗

Liquid Hydrogen Sensor Considerations for Space Exploration

The on-orbit management of liquid hydrogen planned for the return to the moon will introduce new considerations not encountered in previous missions. This paper identifies critical liquid hydrogen sensing needs from the perspective of reliable on-orbit cryogenic fluid management, and contrasts the fundamental differences in fluid and thermodynamic behavior for ground-based versus on-orbit conditions. Opportunities for advanced sensor development and implementation are explored in the context of critical Exploration Architecture operations such as on-orbit storage, docking, and trans-lunar injection burn. Key sensing needs relative to these operations are also examined, including: liquid/vapor detection, thermodynamic condition monitoring, mass gauging, and leak detection. Finally, operational aspects of an integrated system health management approach are discussed to highlight the potential impact on mission success.

Moran, Matthew E.↗

Fluid Expulsion, Habitability, and the Search for Life on Mars

Habitability assessments are critical for identifying settings in which potential biosignatures could exist in quantities large enough to be detected by rovers. Habitability depends on 1) the potential for long-lived liquid water, 2) conditions affording protection from surface processes destructive to organic biomolecules, and 3) a source of renewing nutrients and energy. Of these criteria, the latter is often overlooked. Here we present an analysis of a large "ghost" crater in northern Chryse Planitia [1] that appears to have satisfied each of these requirements, with several processes providing potential sources of nutrient/energy renewal [1-2]. This analysis can serve as a model for identifying other localities that could provide similarly favorable settings in which to seek evidence of life on Mars.

Oehler, Dorothy Z.↗

Reliability of Copper Inverse Opal Surfaces for Extreme-Heat-Flux Micro-Coolers in Low-Global-Warming-Potential Refrigerant R-1233zd Pool Boiling Experiments

This presentation provides a brief snapshot of the InterPACK paper InterPACK2023-113781. The paper explores copper inverse opal (CIO) surface reliability in pool boiling experiments in water and a new, low-global-warming-potential (GWP = 1) hydrofluoroolefin (HFO) refrigerant R-1233zd. The CIO-based structure is intended to develop enhanced two-phase heat transfer surfaces for extreme-heat-flux (approximately 1 kW/cm2) micro-coolers. In this study, a limited number of pool boiling experiments were performed using water and HFO-1233zd fluid, and the reliability of the CIO-based surfaces was evaluated. Critical heat flux (CHF) values in HFO-1233zd at 40 Degrees Celsius to 45 Degrees Celsius saturation temperatures and the corresponding saturation pressures were also measured. The CHF values with the refrigerant are significantly lower compared to those with water, but the refrigerant allows for a wider usable temperature range in the end application of the micro-coolers and is not limited to data centers with controlled ambient conditions. Reliability experiments with CIO surface samples - involving pool boiling with water on the CIO surfaces for approximately 48 hours and with HFO-1233zd for 144 hours - showed no structural degradation of the enhanced surface or any significant performance drop in heat transfer coefficients. The CIO surface samples in water were oxidized, most likely due to the presence of air in water and in the experimental vessel.

critical heat flux↗

Thermoconvective Instability in a Rotating Magnetic Field

The effect of a rotating magnetic field (RMF) on the stability of a fluid contained in a cylindrical column and heated from below is investigated. The RMF increases the critical Rayleigh number for asymmetric flow modes but does not affect the onset of instability for axisymmetric modes. The critical Rayleigh number is dependent upon the relative penetration of the magnetic field into the cylinder and the Prandtl number of the fluid. Instability first develops in the form of a single asymmetric meridional roll rotating around the axis of the cylinder, driven by the azimuthal component of the magnetic field.

Volz, M. P.↗

Preparing an Incompressible-Flow Fluid Dynamics Code for Exascale-Class Wind Energy Simulations

The U.S. Department of Energy has identified exascale-class wind farm simulation as critical to wind energy scientific discovery. A primary objective of the ExaWind project is to build high-performance, predictive computational fluid dynamics (CFD) tools that satisfy these modeling needs. GPU accelerators will serve as the computational thoroughbreds of next-generation, exascale-class supercomputers. Here, we report on our efforts in preparing the ExaWind unstructured mesh solver, Nalu-Wind, for exascale-class machines. For computing at this scale, a simple port of the incompressible-flow algorithms to GPUs is insufficient. To achieve high performance, one needs novel algorithms that are application aware, memory efficient, and optimized for the latest-generation GPU devices. The result of our efforts are unstructured-mesh simulations of wind turbines that can effectively leverage thousands of GPUs. In particular, we demonstrate a first-of-its-kind, incompressible-flow simulation using Algebraic Multigrid solvers that strong scales to more than 4000 GPUs on the Summit supercomputer.

algebraic multigrid↗

Two-phase flows within systems with ambient pressure

In systems where the design inlet and outlet pressures are maintained above the thermodynamic critical pressure, it is often assumed that two phase flows within the system cannot occur. Designers rely on this simple rule of thumb to circumvent problems associated with a highly compressible two phase flow occurring within the supercritical pressure system along with the uncertainties in rotordynamics, load capacity, heat transfer, fluid mechanics, and thermophysical property variations. The simple rule of thumb is adequate in many low power designs but is inadequate for high performance turbomachines and linear systems, where two phase regions can exist even though outlet pressure is greater than critical pressure. Rotordynamic-fluid-mechanic restoring forces depend on momentum differences, and those for a two phase zone can differ significantly from those for a single-phase zone. Using the Reynolds equation the angular velocity, eccentricity, geometry, and ambient conditions are varied to determine the point of two phase flow incipience.

Hendricks, R. C.↗

Carbon Dioxide Absorption Heat Pump

A carbon dioxide absorption heat pump cycle is disclosed using a high pressure stage and a super-critical cooling stage to provide a non-toxic system. Using carbon dioxide gas as the working fluid in the system, the present invention desorbs the CO2 from an absorbent and cools the gas in the super-critical state to deliver heat thereby. The cooled CO2 gas is then expanded thereby providing cooling and is returned to an absorber for further cycling. Strategic use of heat exchangers can increase the efficiency and performance of the system.

Jones, Jack A.↗

A qualitative view of cryogenic fluid injection into high speed flows

The injection of supercritical pressure, subcritical temperature fluids, into a 2-D, ambient, static temperature and static pressure supersonic tunnel and free jet supersonic nitrogen flow field was observed. Observed patterns with fluid air were the same as those observed for fluid nitrogen injected into the tunnel at 90 deg to the supersonic flow. The nominal injection pressure was of 6.9 MPa and tunnel Mach number was 2.7. When injected directly into and opposing the tunnel exhaust flow, the observed patterns with fluid air were similar to those observed for fluid nitrogen but appeared more diffusive. Cryogenic injection creates a high density region within the bow shock wake but the standoff distance remains unchanged from the gaseous value. However, as the temperature reaches a critical value, the shock faded and advanced into the supersonic stream. For both fluids, nitrogen and air, the phenomena was completely reversible.

Hendricks, R. C.↗

Thermal Management Challenges and Technology Options for Lunar Regenerative Fuel Cells

Regenerative Fuel Cells (RFCs) enable lunar surface hardware, from pressurized rovers to surface habitats, to survive the lunar night by providing the required energy storage with a lower mass than state of the art battery systems. RFCs consist of fuel cell and electrolyzer electrochemical stack(s), a fluidic balance of plant, and electrical balance of plant. Fuel cells operate by converting the stored chemical energy of reactants into DC electrical power, heat, and water. When supplied DC power from an available external source, electrolyzers regenerate the reactants through electrolysis. The balance of plant is a critical part of the integrated system consisting of many individual components and provides fluid pressure management and conditioning, fluid storage, and thermal management. Component technology development is required to create an RFC system that will be capable of surviving the harsh lunar environment and maximize specific energy (W-hr/kg). One critical component requiring development includes a kW-scale high flux thermal switch. The primary requirement for the thermal switch is the ability to maintain low thermal resistance with the RFC heat rejection interface during the lunar day to maximize heat rejection while switching to a high thermal resistance during the lunar night such that the RFC can conserve as much thermal energy as possible. Different technology including freeze tolerant pumped loop, passive louvers, and variable conductance heat pipes were assessed as potential solutions. NASA is exploring development options through internal and external partnerships to develop this technology critical for a spaceflight RFC system.

Energy Storage↗

Exploring numerical blow-up phenomena for the Keller–Segel–Navier–Stokes equations

The Keller–Segel–Navier–Stokes system governs chemotaxis in liquid environments. This system is to be solved for the organism and chemoattractant densities and for the fluid velocity and pressure. It is known that if the total initial organism density mass is below 2π there exist globally defined generalised solutions, but what is less understood is whether there are blow-up solutions beyond such a threshold and its optimality. Motivated by this issue, a numerical blow-up scenario is investigated. Approximate solutions computed via a stabilised finite element method founded on a shock capturing technique are such that they satisfy a priori bounds as well as lower and L 1 (Ω) bounds for the organism and chemoattractant densities. In particular, these latter properties are essential in detecting numerical blow-up configurations, since the non-satisfaction of these two requirements might trigger numerical oscillations leading to non-realistic finite-time collapses into persistent Dirac-type measures. Our findings show that the existence threshold value 2π encountered for the organism density mass may not be optimal and hence it is conjectured that the critical threshold value 4π may be inherited from the fluid-free Keller–Segel equations. Additionally it is observed that the formation of singular points can be neglected if the fluid flow is intensified.

97 MATHEMATICS AND COMPUTING↗

Dynamic scaling of order parameter fluctuations in model B

Brookhaven National Laboratory has performed a dedicated set of experiments, called the Beam Energy Scan (BES), with the goal of location a possible critical point in the phase diagram of nuclear matter, This point is analogous to the endpoint of the liquid-gas phase transition in water, where the liquid phase corresponds to a phase of hadrons, and the gas phase to a plasma of quarks and gluons. In order to describe the data generated by these experiments we have studied numerical simulations of the stochastic diffusion equation with a conserved charge. Here, we focus on the dynamics in the vicinity of a critical point in the Ising universality class. The model we consider is expected to describe the critical dynamics near a possible QCD critical point if the coupling of the order parameter to the momentum density of the fluid can be neglected. The simulations are performed on a spatial lattice, and the time evolution is performed using a Metropolis algorithm. We determine the dynamical critical exponent z ≃ 3.972(2), which agrees with predictions of the epsilon expansion. We also study non-equilibrium sweeps of the reduced temperature and observe approximate Kibble-Zurek scaling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Performance Portability Evaluation of Fluid-Structure Interaction Simulations on Heterogeneous Platforms

The rapid proliferation of heterogeneous programming languages and multi-vendor hardware has underscored the critical need to evaluate the performance portability of scientific applications. In this work, we present the systematic porting and optimization of a massively parallel fluid-structure interaction code across multiple heterogeneous programming frameworks for deployment on leadership-class supercomputers from major vendors. Our analysis focuses on at-scale performance for simulations involving hundreds of millions of deformable cells, executed on a combination of CPUs and GPUs spanning thousands of nodes on exascale machines. We benchmark the performance of each implementation, highlighting the trade-offs inherent in adopting diverse programming models. Key insights regarding the portability of CUDA on multi-vendor platforms, the superior multi-core CPU performance from SYCL, and architectural considerations on performance optimization are distilled from our experience, offering guidance to other users of high performance computing based on our findings.

Martin, Aristotle [Duke University]↗