Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat Flux”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Isotope effects under the influence of global radial electric fields in a helical configuration

Isotope effects under the influence of a radial electric field are examined in a helical magnetic field configuration. We perform global gyrokinetic simulations with additional poloidal rotations to estimate quasi-linear heat flux due to ion temperature gradient mode under the mixing length model. In single-ion-species plasmas, the mass number dependency of heat flux agrees with gyro-Bohm scaling in the absence of a radial electric field. Favorable mass number dependencies violating gyro-Bohm scaling are observed in the presence of a global radial electric field or a heavy hydrogen component in multi-ion-species plasmas. The radial electric field and the heavy hydrogen component affect the heat flux through an increase of wavelength as well as mode stabilization. Poloidal Mach number characterizes the transition from unfavorable to favorable mass number dependency under radial electric fields. While the heat flux is independent of mass number for a given poloidal Mach number, the heat flux decreases for higher mass numbers in a given radial electric field. The heat flux is also independent of average mass number in multi-ion-species plasmas because the heavy hydrogen component effectively enhances the light hydrogen heat flux. The present results are potentially relevant to the violation of gyro-Bohm scaling observed in the recent deuterium experiments in the Large Helical Device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2008

Mobile Lid Convection Beneath Enceladus' South Polar Terrain

Enceladus' south polar region has a large heat flux, 55-110 milliwatts per square meter (or higher), that is spatially associated with cryovolcanic and tectonic activity. Tidal dissipation and vigorous convection in the underlying ice shell are possible sources of heat; however, prior predictions of the heat flux carried by stagnant lid convection range from F(sub conv) ~15 to 30 milliwatts per square meter, too low to explain the observed heat flux. The high heat flux and increased cryovolcanic and tectonic activity suggest that near-surface ice in the region has become rheologically and mechanically weakened enough to permit convective plumes to reach close to the surface. If the yield strength of Enceladus' lithosphere is less than ~1-10 kPa, convection may instead occur in the mobile lid" regime, which is characterized by large heat fluxes and large horizontal velocities in the near-surface ice. I show that model ice shells with effective surface viscosities between 10(exp 16) and 10(exp 17) Pa s and basal viscosities between 10(exp 13) and 10(exp 15) Pa s have convective heat fluxes comparable to that observed by the Cassini Composite Infrared Spectrometer. If this style of convection is occurring, the south polar terrain should be spreading horizontally with v~1-10 millimeter per year and should be resurfaced in ~0.1-10 Ma. On the basis of Cassini imaging data, the south polar terrain is ~0.5 Ma old, consistent with the mobile lid hypothesis. Maxwell viscoelastic tidal dissipation in such ice shells is not capable of generating enough heat to balance convective heat transport. However, tidal heat may also be generated in the near-surface along faults as suggested by Nimmo et al. and/or viscous dissipation within the ice shell may occur by other processes not accounted for by the canonical Maxwell dissipation model.

Barr, Amy C.↗
NASA NTRS1995

Electron energy transport in the solar wind: Ulysses observations

The electron heat flux in the solar wind has been measured by the Ulysses solar wind plasma experiment in the ecliptic from 1 to 5 AU and out of the ecliptic during the recently completed pass over the solar south pole and the ongoing pass over the solar north pole. Although the electron heat flux contains only a fraction of the kinetic energy of the solar wind. the available energy is sufficient to account for the non-adiabatic expansion of the solar wind electrons. The Ulysses measurements indicate that the electron heat flux is actively dissipated in the solar wind. The exact mechanism or mechanisms is unknown. but a model based on the whistler heat flux instability predicts radial gradients for the electron heat flux in good agreement with the data. We will present measurements of the correlation between wave activity measured by the unified radio and plasma experiment (URAP) and the electron heat flux throughout the Ulysses mission. The goal is to determine if whistler waves are a good candidate for the observed electron heat flux dissipation. The latitudinal gradients of the electron heat flux. wave activity. and electron pressure will be discussed in light of the changes in the magnetic field geometry from equator to poles.

Scime, Earl↗
NASA NTRS1998

Design, Construction, and Qualification of a Microscale Heater Array for Use in Boiling Heat Transfer

Boiling heat transfer is an efficient means of heat transfer because a large amount of heat can be removed from a surface using a relatively small temperature difference between the surface and the bulk liquid. However, the mechanisms that govern boiling heat transfer are not well understood. Measurements of wall temperature and heat flux near the wall would add to the database of knowledge which is necessary to understand the mechanisms of nucleate boiling. A heater array has been developed which contains 96 heater elements within a 2.5 mm square area. The temperature of each heater element is held constant by an electronic control system similar to a hot-wire anemometer. The voltage that is being applied to each heater element can be measured and digitized using a high-speed A/D converter, and this digital information can be compiled into a series of heat-flux maps. Information for up to 10,000 heat flux maps can be obtained each second. The heater control system, the A/D system and the heater array construction are described in detail. Results are presented which show that this is an effective method of measuring the local heat flux during nucleate and transition boiling. Heat flux maps are obtained for pool boiling in FC-72 on a horizontal surface. Local heat flux variations are shown to be three to six times larger than variations in the spatially averaged heat flux.

Rule, T. D.↗
NASA NTRS1996

A Computational Fluid Dynamic and Heat Transfer Model for Gaseous Core and Gas Cooled Space Power and Propulsion Reactors

A computational model based on the axisymmetric, thin-layer Navier-Stokes equations is developed to predict the convective, radiation and conductive heat transfer in high temperature space nuclear reactors. An implicit-explicit, finite volume, MacCormack method in conjunction with the Gauss-Seidel line iteration procedure is utilized to solve the thermal and fluid governing equations. Simulation of coolant and propellant flows in these reactors involves the subsonic and supersonic flows of hydrogen, helium and uranium tetrafluoride under variable boundary conditions. An enthalpy-rebalancing scheme is developed and implemented to enhance and accelerate the rate of convergence when a wall heat flux boundary condition is used. The model also incorporated the Baldwin and Lomax two-layer algebraic turbulence scheme for the calculation of the turbulent kinetic energy and eddy diffusivity of energy. The Rosseland diffusion approximation is used to simulate the radiative energy transfer in the optically thick environment of gas core reactors. The computational model is benchmarked with experimental data on flow separation angle and drag force acting on a suspended sphere in a cylindrical tube. The heat transfer is validated by comparing the computed results with the standard heat transfer correlations predictions. The model is used to simulate flow and heat transfer under a variety of design conditions. The effect of internal heat generation on the heat transfer in the gas core reactors is examined for a variety of power densities, 100 W/cc, 500 W/cc and 1000 W/cc. The maximum temperature, corresponding with the heat generation rates, are 2150 K, 2750 K and 3550 K, respectively. This analysis shows that the maximum temperature is strongly dependent on the value of heat generation rate. It also indicates that a heat generation rate higher than 1000 W/cc is necessary to maintain the gas temperature at about 3500 K, which is typical design temperature required to achieve high efficiency in the gas core reactors. The model is also used to predict the convective and radiation heat fluxes for the gas core reactors. The maximum value of heat flux occurs at the exit of the reactor core. Radiation heat flux increases with higher wall temperature. This behavior is due to the fact that the radiative heat flux is strongly dependent on wall temperature. This study also found that at temperature close to 3500 K the radiative heat flux is comparable with the convective heat flux in a uranium fluoride failed gas core reactor.

Anghaie, S.↗
DOE OSTI2024

Optimization of nonlinear turbulence in stellarators

We present new stellarator equilibria that have been optimized for reduced turbulent transport using nonlinear gyrokinetic simulations within the optimization loop. The optimization routine involves coupling the pseudo-spectral GPU-native gyrokinetic code GX with the stellarator equilibrium and optimization code DESC . Since using GX allows for fast nonlinear simulations, we directly optimize for reduced nonlinear heat fluxes. To handle the noisy heat flux traces returned by these simulations, we employ the simultaneous perturbation stochastic approximation (SPSA) method that only uses two objective function evaluations for a simple estimate of the gradient. We show several examples that optimize for both reduced heat fluxes and good quasi-symmetry as a proxy for low neoclassical transport. Finally, we run full transport simulations using the T3D stellarator transport code to evaluate the changes in the macroscopic profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2024

Defining the Operational Envelope for Air Flows in the Miniature Arc-Jet Research Chamber (mARC II)

The second-generation 30 kW miniature Arc-jet Research Chamber (mARC II) at NASA Ames Research Center produces high enthalpy flows relevant for entry systems ground testing. The mARC II facility has recently undergone upgrades, including the installation of a new vacuum system to address the issues preventing it from maintaining underexpanded flow under test conditions. In this work, we present data obtained from an Integrated Systems Testing campaign and provide an initial assessment of arc-jet performance following the upgrades. Air is used as the working gas for the standard mARC II arc-heater configuration with two constrictor disks. Seven runs were investigated for five test conditions to assess the lowest achievable stagnation point heat fluxes for air flow rates of 0.15 or 0.25 g s −1 . The heat flux was measured using a water-cooled Gardon gauge (Ø4.76 mm, 3/16" hemispherical) at 70 mm from the nozzle exit plane. The new vacuum system produced test box pressures in the medium (fine) vacuum range (~0.03 torr, 4 Pa) prior to gas addition and successfully maintained underexpanded flow after gas addition. The upgrade yielded a ~4X reduction in heat flux relative to the previous system for the same set test conditions. We report the lowest heat fluxes measured in mARC II to date, ranging from 26 to 81 W cm −2 , for sonic flow enthalpies of 4–14 MJ kg −1 . Bulk enthalpies estimated using an energy balance method (EB 2 ) are reported for the first time using mARC II.Initial data suggests EB 2 generally estimates lower enthalpies than sonic flow methods for mARC II. Lastly, laminar axisymmetric Navier–Stokes simulations were performed using the NASA DPLR code. Numerical heat flux results show good agreement with experiments at low arc powers (3–15% difference at the minimum set arc current), but discrepancy increases with arc power (49% difference at the maximum set arc current).

Pitot tube↗
NASA NTRS2010

Multiphysics Modeling for Dimensional Analysis of a Self-Heated Molten Regolith Electrolysis Reactor for Oxygen and Metals Production on the Moon and Mars

The technology of direct electrolysis of molten lunar regolith to produce oxygen and molten metal alloys has progressed greatly in the last few years. The development of long-lasting inert anodes and cathode designs as well as techniques for the removal of molten products from the reactor has been demonstrated. The containment of chemically aggressive oxide and metal melts is very difficult at the operating temperatures ca 1600 C. Containing the molten oxides in a regolith shell can solve this technical issue and can be achieved by designing a self-heating reactor in which the electrolytic currents generate enough Joule heat to create a molten bath. In a first phase, a thermal analysis model was built to study the formation of a melt of lunar basaltic regolith irradiated by a focused solar beam This mode of heating was selected because it relies on radiative heat transfer, which is the dominant mode of transfer of energy in melts at 1600 C. Knowing and setting the Gaussian-type heat flux from the concentrated solar beam and the phase and temperature dependent thermal properties, the model predicts the dimensions and temperature profile of the melt. A validation of the model is presented in this paper through the experimental formation of a spherical cap melt realized by others. The Orbitec/PSI experimental setup uses an 3.6-cm diameter concentrated solar beam to create a hemispheric melt in a bed of lunar regolith simulant contained in a large pot. Upon cooling, the dimensions of the vitrified melt are measured to validate the thermal model. In a second phase, the model is augmented by multiphysics components to compute the passage of electrical currents between electrodes inserted in the molten regolith. The current through the melt generates Joule heating due to the high resistivity of the medium and this energy is transferred into the melt by conduction, convection and primarily by radiation. The model faces challenges in two major areas, the change of phase as temperature increases, and the dominance of radiative heat flux as heat transfer mechanism within the melt the change of phase concerns the regolith itself which is present in states ranging from a fine grain regolith with low thermal conductivity and low density to a vitrified melt with much higher thermal conductivity, and higher density. As the regolith is heated, it starts to soften around 1300 C the melt iS very viscous and evolving gas bubbles out in thick, lava-like fashion. By 1600 C the regolith is completely melted and the viscosity is low The second challenge resides in the proper modeling of the radiative heat flux requiring the addition of the computing-demanding radiative-heat-transfer function to the general heat transfer equation. The model Includes temperature-dependent properties (density, thermal conductivity, heat capacity, and viscosity, and absorption coefficients) and solves the radiative heat flux equation assuming gray (fine grains) and semi-transparent (melt) media and using an absorption coefficient spectral found in the literature for terrestrial minerals similar in composition to those of lunar regolith simulant

Dominguez, Jesus A.↗
NASA NTRS2024

Defining the Operational Envelope for Air Flows in the Miniature Arc-Jet Research Chamber (mARC II)

The second-generation 30 kW miniature Arc-jet Research Chamber (mARC II) at NASA Ames Research Center produces high enthalpy flows relevant for entry systems ground testing. The mARC II facility has recently undergone upgrades, including the installation of a new vacuum system to address the issues preventing it from maintaining underexpanded flow under test conditions. In this work, we present data obtained from an Integrated Systems Testing campaign and provide an initial assessment of arc-jet performance following the upgrades. Air is used as the working gas for the standard mARC II arc-heater configuration with two constrictor disks. Seven runs were investigated for five test conditions to assess the lowest achievable stagnation point heat fluxes for air flow rates of 0.15 or 0.25 g s −1 . The heat flux was measured using a water-cooled Gardon gauge (Ø4.76 mm, 3/16" hemispherical) at 70 mm from the nozzle exit plane. The new vacuum system produced test box pressures in the medium (fine) vacuum range (~0.03 torr, 4 Pa) prior to gas addition and successfully maintained underexpanded flow after gas addition. The upgrade yielded a ~4X reduction in heat flux relative to the previous system for the same set test conditions. We report the lowest heat fluxes measured in mARC II to date, ranging from 26 to 81 W cm −2 , for sonic flow enthalpies of 4–14 MJ kg −1 . Bulk enthalpies estimated using an energy balance method (EB 2 ) are reported for the first time using mARC II.Initial data suggests EB 2 generally estimates lower enthalpies than sonic flow methods for mARC II. Lastly, laminar axisymmetric Navier–Stokes simulations were performed using the NASA DPLR code. Numerical heat flux results show good agreement with experiments at low arc powers (3–15% difference at the minimum set arc current), but discrepancy increases with arc power (49% difference at the maximum set arc current).

arc-jet↗
DOE OSTI2024

Optimization of lithium vapor box divertor evaporator location on NSTX-U using SOLPS-ITER

Commercial fusion reactors will be faced with extremely high divertor target heat fluxes that will require mitigation. Simulations of detachment in an NSTX-U scenario projected to have 92 MW/m 2 unmitigated peak target heat flux are presented, which reaches sub-10 MW/m 2 target heat flux using a highly dissipating lithium vapor box divertor design. The lithium vapor box is a detached divertor design which employs lithium vapor evaporation and condensation to contain lithium below the X-point. Previous SOLPS modeling has indicated a lithium vapor box can reduce the heat flux down to 10 MW/m 2 via simultaneous evaporation from the Private Flux Region (PFR) and the Common Flux Region (CFR) sides of the vapor box. It is found here that PFR evaporation has improved access to the separatrix leading to significantly more efficient power dissipation than CFR evaporation. Simulations of target evaporation with an evaporation distribution that is self-consistent with the temperature of a Capillary Porous System with Fast flowing liquid lithium (CPSF) could reach n Li / n e ~ 0.025-0.030 at the Last Closed Flux Surface (LCFS) depending on the liquid metal flow speeds and lithium sputtering yield, while PFR-side evaporation can reach acceptable heat fluxes with n Li / n e ~ 0.038 at the LCFS. However, PFR evaporator performance can be improved if the target is allowed to be hot enough such that it reflects lithium, reaching n Li / n e ~ 0.028 and reducing required lithium evaporation. Ultimately PFR evaporation and target evaporation are found to have similar ability to produce acceptable heat flux solutions with minimal upstream concentration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Influence of Heat Flow from the Planetary Interior on Surface Habitability

an primordial and/or radiogenic heat flow from a rocky planet’s interior help maintain a surface environment amenable to life?On rocky Earth-like planets with shallow oceans (≤5 km depth) that are on the outer edge of the habitable zone, the additional energy provided might make the difference between a world supporting regional habitat space and a world uninhabitable by any known life-form. To test this hypothesis, we introduce heat flux from the planetary interior as a forcing to the NASA/GISSROCKE-3D exoplanet GCM, and explore whether it can enhance the habitability of marginal terrestrial rocky planets (e.g.,paleo Earth in a “snowball” climate state).GCM simulations of modern Earth have not included geothermal heat flux as a climate forcing, as it is more than an order of magnitude less than the 2.9 W/m heating caused by anthropogenic greenhouse gases [Flanner, GRL 2009]. Modern Earth measurements in the oceans [Downes et al., GRL 2019] and results from ocean GCMs [Barnes et al., Ocean Model. 2017]suggest that localized higher heat flux can have measurable thermodynamic and dynamic impacts, such as weakening deep ocean stratification and increased poleward ocean heat transports. A similar result with ROCKE-3D would be important, since our snowball Earth simulations without heat flux produce features like ocean stratification that should have hampered the survival of life. We have created a set of spatially variable geothermal heat fluxes at varying horizontal resolutions based on modern Earth heat flow patterns [Davies, Geochem. Geophys. Geosys. 2013] to test the sensitivity of the ROCKE-3D GCM initially in a modern Earth context. Preliminary results show that even modest heat fluxes have a measurable impact on deep ocean temperatures (upto +0.9°C locally) and that this heat can be transported at depth along ocean circulation paths. Estimates of heat flux impacts on vertical ocean mixing and surface ocean conditions require long simulations at the higher spatial resolutions needed to express heat flow at the scale of most ocean ridges. Ongoing simulations are focused on testing the impacts of horizontal resolution and flux scaling on surface environments, especially for “snowball” climate scenarios.

Heat flow↗
NASA NTRS1961

Base Flow Characteristics for Several Four-Clustered Rocket Configurations at Mach Numbers from 2.0 to 3.5

A generalized study of base flow phenomena has been conducted with four 500-pound-thrust JP-4 fuel-liquid-oxygen rocket motors installed in the base of a 12-inch-diameter cylindrical model. Data were obtained over a Mach number and nozzle pressure ratio range of 2.0 to 3.5 and 340 to 600, respectively. Base heat flux, gas temperature, and pressure were highest in the center of the cluster core and decreased in a radial direction. Although a maximum heat flux of 93 Btu per square foot per second was measured within the cluster core, peripheral heat fluxes were low, averaging about 5 Btu per square foot per second for all configurations. Generally base heat flux was found to be independent of Mach number over the range investigated. Base heat flux within the cluster core was decreased by increasing motor spacing, motor extension, a combination of increasing nozzle area ratio and decreasing exit angle and gimbaling the two side engines. Small amounts of nitrogen injected within the cluster core sharply reduced core heat flux.

Musial, Norman T.↗
NASA NTRS1994

Ocean response to surface heat anomalies

An ocean general circulation model (OGCM) is used to study the response of ocean heat and mass transport to positive and negative heat flux anomalies at the ocean surface. As expected, tropical and low-latitude mixed layers respond rapidly (e-folding time about 50-70 years) to external forcing, while the response of the high-latitude mixed layer, especially the Southern Ocean and northern North Atlantic, is very slow (e-folding time greater than 300 yr). The overall response is faster for negative than positive heat flux anomaly at the surface. The meridional heat transport changes by 15% in the first 50 yr in the southern high latitudes. Surprisingly, for the next 400-500 yr the change is very small. The analysis shows that the meridional mass transport intensifies in response to a negative surface heat flux anomaly but weakens in response to a positive heat flux anomaly. For example, at model year 100 the North Atlantic Deep Water (NADW) is reduced from about 18 Sv to about 10 Sv for the positive heat flux experiment but increased to about 26 Sv for the negative heat flux experiment.

Jiang, Xingjian↗
NASA NTRS2016

Flow Boiling and Condensation Experiment (FBCE) for the International Space Station

An effective means to reducing the size and weight of future space vehicles is to replace present mostly single-phase thermal management systems with two-phase counterparts. By capitalizing upon both latent and sensible heat of the coolant rather than sensible heat alone, two-phase thermal management systems can yield orders of magnitude enhancement in flow boiling and condensation heat transfer coefficients. Because the understanding of the influence of microgravity on two-phase flow and heat transfer is quite limited, there is an urgent need for a new experimental microgravity facility to enable investigators to perform long-duration flow boiling and condensation experiments in pursuit of reliable databases, correlations and models. This presentation will discuss recent progress in the development of the Flow Boiling and Condensation Experiment (FBCE) for the International Space Station (ISS) in collaboration between Purdue University and NASA Glenn Research Center. Emphasis will be placed on the design of the flow boiling module and on new flow boiling data that were measured in parabolic flight, along with extensive flow visualization of interfacial features at heat fluxes up to critical heat flux (CHF). Also discussed a theoretical model that will be shown to predict CHF with high accuracy.

flow boiling↗
NASA NTRS1997

Cabauw Experimental Results from the Project for Intercomparison of Land-Surface Parameterization Schemes

In the Project for Intercomparison of Land-Surface Parameterization Schemes phase 2a experiment, meteorological data for the year 1987 from Cabauw, the Netherlands, were used as inputs to 23 land-surface flux schemes designed for use in climate and weather models. Schemes were evaluated by comparing their outputs with long-term measurements of surface sensible heat fluxes into the atmosphere and the ground, and of upward longwave radiation and total net radiative fluxes, and also comparing them with latent heat fluxes derived from a surface energy balance. Tuning of schemes by use of the observed flux data was not permitted. On an annual basis, the predicted surface radiative temperature exhibits a range of 2 K across schemes, consistent with the range of about 10 W/m in predicted surface net radiation. Most modeled values of monthly net radiation differ from the observations by less than the estimated maximum monthly observational error (+/- 10 W/sq m). However, modeled radiative surface temperature appears to have a systematic positive bias in most schemes; this might be explained by an error in assumed emissivity and by models' neglect of canopy thermal heterogeneity. Annual means of sensible and latent heat fluxes, into which net radiation is partitioned, have ranges across schemes of 30 W/sq m and 25 W/sq m, respectively. Annual totals of evapotranspiration and runoff, into which the precipitation is partitioned, both have ranges of 315 mm. These ranges in annual heat and water fluxes were approximately halved upon exclusion of the three schemes that have no stomatal resistance under non-water-stressed conditions. Many schemes tend to underestimate latent heat flux and overestimate sensible heat flux in summer, with a reverse tendency in winter. For six schemes, root-mean-square deviations of predictions from monthly observations are less than the estimated upper bounds on observation errors (5 W/m for sensible beat flux and 10 W/m for latent heat flux). Actual runoff at the site is believed to be dominated by vertical drainage to ground water, but several schemes produced significant amounts of runoff as overland flow or interflow. There is a range across schemes of 184 mm (40% of total pore volume) in the simulated annual mean root-zone soil moisture. Unfortunately, no measurements of soil moisture were available for model evaluation. A theoretical analysis suggested that differences in boundary conditions used in various schemes are not sufficient to explain the large variance in soil moisture. However, many of the extreme values of soil moisture could be explained in terms of the particulars of experimental setup or excessive evapotranspiration.

Chen, Tian Hong↗
DOE OSTI2026

Simulation of Divertor Performance in ST40 Under Dynamic Double-Null Plasmas

A power fraction model was implemented for the simultaneous prediction of 3-D surface temperature evolution at all four divertor targets in near-double-null (DN) tokamak configurations, which is especially important for compact high-field devices that may not have the ability to dissipate large amounts of power on the high-field side. Evaluating the power-sharing between the four divertor strike points in a disconnected DN configuration is important for understanding the overall power balance, as well as for optimizing the power exhaust performance and prolonging the survivability of the plasma facing components (PFCs). This power-sharing is typically evaluated in terms of the separation between the primary and secondary separatrices at the outboard midplane, $\textit {dR}_{\text {sep}}$. The Heat flux Engineering Analysis Toolkit (HEAT) is coupled with Brunner’s power fraction model to simulate the deposited heat flux and resultant temperature change on 3-D divertor targets in a dynamic DN (DDN) pulse operation in ST40, a high-field spherical tokamak. The simulation results showed that with DDN operation, the operation time has significantly increased compared with single-null geometry configurations.

ST40↗
NASA NTRS1999

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.↗
DOE OSTI2025

How does ion temperature gradient turbulence depend on magnetic geometry? Insights from data and machine learning

Magnetic geometry has a significant effect on the level of turbulent transport in fusion plasmas. Here, we model and analyse this dependence using multiple machine learning methods and a dataset of >200 000 nonlinear gyrokinetic simulations of ion-temperature-gradient turbulence in diverse non-axisymmetric geometries. The dataset is generated using a large collection of both optimised and randomly generated stellarator equilibria. At fixed gradients and other input parameters, the turbulent heat flux varies between geometries by several orders of magnitude. Trends are apparent among the configurations with particularly high or particularly low heat flux. Regression and classification techniques from machine learning are then applied to extract patterns in the dataset. Due to a symmetry of the gyrokinetic equation, the heat flux and regressions thereof should be invariant to translations of the raw features in the parallel coordinate, similar to translation invariance in computer vision applications. Multiple regression models including convolutional neural networks (CNNs) and decision trees can achieve reasonable predictive power for the heat flux in held-out test configurations, with highest accuracy for the CNNs. Using Spearman correlation, sequential feature selection and Shapley values to measure feature importance, it is consistently found that the most important geometric lever on the heat flux is the flux surface compression in regions of bad curvature. The second most important geometric feature relates to the magnitude of geodesic curvature. These two features align remarkably with surrogates that have been proposed based on theory, while the methods here allow a natural extension to more features for increased accuracy. The dataset, released with this publication, may also be used to test other proposed surrogates, and we find that many previously published proxies do correlate well with both the heat flux and stability boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗