Search NASA⌕ Search

SEARCH · Search NASA

Results for “Critical Velocity”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A critical look at the prediction of the temperature field around a laser-induced melt pool on metallic substrates

The study of microstructure evolution in additive manufacturing of metals would be aided by knowing the thermal history. Since temperature measurements beneath the surface are difficult, estimates are obtained from computational thermo-mechanical models calibrated against traces left in the sample revealed after etching, such as the trace of the melt pool boundary. Here we examine the question of how reliable thermal histories computed from a model that reproduces the melt pool trace are. To this end, we perform experiments in which one of two different laser beams moves with constant velocity and power over a substrate of 17-4PH SS or Ti-6Al-4V, with low enough power to avoid generating a keyhole. We find that thermal histories appear to be reliably computed provided that (a) the power density distribution of the laser beam over the substrate is well characterized, and (b) convective heat transport effects are accounted for. Poor control of the laser beam leads to potentially multiple three-dimensional melt pool shapes compatible with the melt pool trace, and therefore to multiple potential thermal histories. Ignoring convective effects leads to results that are inconsistent with experiments, even for the mild melt pools here.

42 ENGINEERING↗
DOE OSTI2025

Upper critical magnetic field and multiband superconductivity in artificial high-𝑇 𝑐 superlattices of nano quantum wells

Artificial high-T c superlattices (AHTS) composed of quantum building blocks with tunable superconducting critical temperature have been synthesized by engineering their nanoscale geometry using the Bianconi-Perali Valletta (BPV) two-gap superconductivity theory. These quantum heterostructures consist of quantum wells made of superconducting, modulation-doped Mott insulators (S), confined by a metallic (N) potential barrier. The lattice geometry has been carefully engineered to induce the predicted Fano-Feshbach shape resonance between the gaps, near a topological Lifshitz transition. Here, we validate the BPV theory by providing compelling experimental evidence that AHTS samples, at the peak of the superconducting dome, exhibit resonant two-band, two-gap superconductivity. This is demonstrated by measuring the temperature dependence of the upper critical magnetic field, μ 0 H c2 , in samples with superlattice periods 3.3 < d < 5.28 nm and L/d ratios close to the magic value 2/3 (where L is the thickness of the superconducting La 2 CuO 4 layer and d is the superlattice period). Here, the data reveal the predicted upward concavity in H c2 (T) and a characteristic kink in the coherence length as a function of temperature, confirming the predicted two-band superconductivity with Fermi velocity ratio ≈ 0.25 and significant pair-exchange term among the two condensates.

Heterostructures↗
DOE OSTI2021

Abrasion wear at high temperature in particle receiver type Concentrating Solar Power Systems

The use of solid particles as a heat transfer fluid in receivers for concentrating solar powerplants has garnered much attention in the past decade. Particles are not only stable at temperatures above 700 °C but are also able to absorb solar energy directly, which gives them a distinct advantage over conventional molten salts. The use of particles can, however, cause material degradation from the flow of hot or cold particles through discharge hoppers or along the inner receiver surfaces and other system components (e.g. tubes, valves etc.), depending on operating mode of the receiver. Stainless steel alloy 316, nickel-based alloy Inconel® 740H®, Haynes 230, and a 60% alumina, low-cement refractory castable from Allied Mineral Products were tested for wear from first-of-a-kind high-temperature low-velocity abrasion-erosion experiments. The tests are designed to simulate potential conditions in Generation 3 falling particle Concentrating Solar Power (CSP) systems. Simultaneous experiments were conducted at room temperature as well to draw comparison on the role of temperature on abrasion wear rates. It was noted that the abrasion wear rate in candidate materials was driven by the oxidation resistance of the material. The findings from these experiments will be used to determine life-time of critical components in particle receiver concentrated solar powerplants.

14 SOLAR ENERGY↗
DOE OSTI2022

Controlling pollutant emissions in a high-pressure combustor with fuel-diluent blending

This paper investigates the formation of nitric oxide (NO) and carbon monoxide (CO) emissions as a function of nitrogen (N 2 ) or carbon dioxide (CO 2 ) diluent content, added to a premixed reacting jet in crossflow. Here, reaction characteristics of a rich methane-air jet injected into a lean vitiated crossflow were analyzed at an elevated pressure of 5 atm. The jet was pre-heated and enriched with 0%, 15%, 30%, and 50% mass fraction to quantify the effect on pollutant emissions. Simulated results of the full chemistry Star-CCM+CFD model were verified with data taken in the experimental high-pressure combustion facility, which provides pressure, temperature and velocity profiles, as well as line-of-sight chemiluminescence and exit emission measurements. The analysis revealed the significant influence of the diluent to delay axial combustion; this mixing delay increased the axial flame lift-off and controlled any thermal ignition limitations. Increasing the diluent content increased the timescale for the flame to stabilize, which allowed for greater entrainment of crossflow oxidizer into the axial jet stream, and led to decreased pollutant emissions. Hence, crossflow entrainment is a critical driving force at high diluent content. Local nitric oxide (NO) emission formation in the axial stage was predicted numerically, showing the correlation between diluent addition, axial heat release, and the formation of nitric oxide pollution. The high diluent levels leaned out the jet mixture and delayed axial combustion, while reducing prompt and thermal NOx by mitigating flame hotspots and minimizing the timescale that the products remain at high temperatures.

03 NATURAL GAS↗
DOE OSTI2025

Numerical simulation of compressible fluid-dynamics in the chamber of inertial fusion energy systems

Here, this paper aims to establish new and innovative modeling capabilities for analyzing chambers in Inertial Fusion Energy (IFE) systems. IFE is emerging as a promising method to achieve fusion power production, but several challenges must be overcome to develop an IFE pilot plant or deploy commercial IFE systems. These challenges are both theoretical and technical, encompassing a deeper understanding of the underlying physical phenomena and the development of new technologies and materials. One of the needs is to develop mathematical models to describe IFE systems and numerical tools to simulate them. This paper contributes to this endeavor by presenting a new OpenFOAM solver for IFE systems, focusing on gas dynamics in their chambers. The analysis and development of chamber designs will play a significant role in the transition from single-shot experiments to high-repetition rates, as there is a need to protect the chamber walls from the intense radiation fields produced by fusion reactions. A promising design option, normally referred to as thick wall chamber design, consists in using lithium or molten salt jet arrays within the chamber. A critical phenomenon is the venting of high-pressure gases from the center to the external part of the chamber, passing through the blanket array. This process involves the propagation and attenuation of strong pressure waves, requiring suitable modeling approaches for compressible fluid-dynamics. The solver proposed in this work implements a multi-material hydrodynamics model tailored to accurately describe the non-linear propagation of pressure waves while avoiding numerical oscillation issues typical of high-velocity compressible simulation. This solver is verified against numerical test cases, validated against experimental data, and applied to the analysis of the High-Yield Lithium-Injection Fusion-Energy (HYLIFE-I) concept. The relevance of this paper is threefold. Firstly, it contributes to developing and testing modeling approaches for compressible fluid-dynamics phenomena, with specific focus on the new and unexplored topic of IFE thick-liquid-wall blanket modeling. Secondly, it marks one of the first applications of the OpenFOAM library in the research field of IFE systems. Finally, the investigated problem is of practical interest for IFE developers, as it provides useful indications about relevant phenomena in pressure wave propagation in the chamber of these systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Incoherent Thomson scattering system for PHAse space MApping (PHASMA) experiment

A new incoherent Thomson scattering system measures the evolution of electron velocity distribution functions perpendicular and parallel to the ambient magnetic field during kinking of a single flux rope and merging of two flux ropes through magnetic reconnection. The Thom-son scattering system provides sub-millimeter spatial resolution, sufficient to diagnose the several millimeters sized magnetic reconnection electron diffusion region in the PHAse Space MAppgin experiment. Due to the relatively modest plasma density ~ 10 19 m -3 and electron tem-perature ~1 eV, stray light suppression is critical for these measurements. Two volume Bragg gratings are used in series as a notch filter with a spectral bandwidth <0.1 nm in the collection branch. A CCD with a Gen III intensifier with peak quantum efficiency >47% is used as the detector in a 1.3 m spectrometer. Preliminary results of gun plasma electron temperature will be reported and compared with measurements obtained from a triple Langmuir probe.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Longitudinal coherence and short-term wind speed prediction based on a nacelle-mounted Doppler lidar

The spatial structure of turbulence in atmospheric boundary layer flows is highly relevant to wind energy. In particular, wind turbine control strategies based on inflow preview measurements require knowledge of the longitudinal evolution of turbulent flow as it approaches the rotor. These upstream measurements are usually obtained with nacelle-mounted wind lidars. In contrast to traditional in situ anemometry, lidars collect measurements within a probe volume which varies in size depending on the technology of the commercial system being used. Here, we address two issues related to the use of wind lidar to measure the incoming flow to a wind turbine: (i) whether existing longitudinal coherence models can be used to predict flow at the rotor, based on measurements performed at a distance away from the rotor; and (ii) what effect probe-volume averaging has on the inflow predictions. These two questions are critical to the design and implementation of robust wind turbine control strategies. To address these questions, we perform field measurements and large-eddy simulations to determine which incoming flow structures can be readily predicted with existing coherence models, and which require additional corrections to account for lidar volumetric averaging effects. Results reveal that the wind turbine induction zone has a negligible impact on the longitudinal coherence and first-order turbulence quantities, such as the standard deviation of velocity fluctuations. However, the phase of the signal, from which advection time periods of the turbulent structures are derived, is affected by the rotor blockage effect.

17 WIND ENERGY↗
DOE OSTI2020

Advanced Newton methods for geodynamical models of Stokes flow with viscoplastic rheologies

Strain localization and resulting plasticity and failure play an important role in the evolution of the lithosphere. These phenomena are commonly modeled by Stokes flows with viscoplastic rheologies. The nonlinearities of these rheologies make the numerical solution of the resulting systems challenging, and iterative methods often converge slowly or not at all. Yet accurate solutions are critical for representing the physics. Moreover, for some rheology laws, aspects of solvability are still unknown. We study a basic but representative viscoplastic rheology law. The law involves a yield stress that is independent of the dynamic pressure, referred to as von Mises yield criterion. Two commonly used variants, perfect/ideal and composite viscoplasticity, are compared. We derive both variants from energy minimization principles, and we use this perspective to argue when solutions are unique. We propose a new stress-velocity Newton solution algorithm that treats the stress as an independent variable during the Newton linearization but requires solution only of Stokes systems that are of the usual velocity-pressure form. To study different solution algorithms, we implement 2-D and 3-D finite element discretizations, and we generate Stokes problems with up to 7 orders of magnitude viscosity contrasts, in which compression or tension results in significant nonlinear localization effects. Comparing the performance of the proposed Newton method with the standard Newton method and the Picard fixed-point method, we observe a significant reduction in the number of iterations and improved stability with respect to problem nonlinearity, mesh refinement, and the polynomial order of the discretization.

Rudi, Johann↗
DOE OSTI2020

Dynamic Phase Alignment in Inertial Alfvén Turbulence

In weakly collisional plasma environments with sufficiently low electron beta, Alfvénic turbulence transforms into inertial Alfvénic turbulence at scales below the electron skin depth, k ⊥ d e ≳ 1. We argue that, in inertial Alfvénic turbulence, both energy and generalized kinetic helicity exhibit direct cascades. We demonstrate that the two cascades are compatible due to the existence of a strong scale dependence of the phase alignment angle between velocity and magnetic field fluctuations, with the phase alignment angle scaling as cos α k ∝ $k$$^{–1}_{⊥}$. The kinetic and magnetic energy spectra scale as ∝ $k$$^{–5/3}_{⊥}$ and ∝ $k$$^{–11/3}_{⊥}$, respectively. As a result of the dual direct cascade, the generalized helicity spectrum scales as ∝ $k$$^{–5/3}_{⊥}$, implying progressive balancing of the turbulence as the cascade proceeds to smaller scales in the k ⊥ d e >> 1 range. Turbulent eddies exhibit a phase-space anisotropy k ∥ ∝ $k$$^{5/3}_{⊥}$, consistent with critically balanced inertial Alfvén fluctuations. Furthermore, our results may be applicable to a variety of geophysical, space, and astrophysical environments, including the Earth’s magnetosheath and ionosphere, solar corona, and nonrelativistic pair plasmas, as well as to strongly rotating nonionized fluids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Optimization of a traveling wave superconducting radiofrequency cavity for upgrading the International Linear Collider

The standing wave TESLA niobium-based superconducting radio frequency structure is limited to an accelerating gradient of about 50 MV/m by the critical rf magnetic field. To break through this barrier, we explore the option of niobium-based traveling wave (TW) structures. Optimization of TW structures was done considering experimentally known limiting electric and magnetic fields. It is shown that a TW structure can have an accelerating gradient above 70 MeV/m that is about 1.5 times higher than contemporary standing wave structures with the same critical magnetic field. The other benefit of TW structures shown is R/Q about 2 times higher than the TESLA structure that reduces the dynamic heat load by a factor of 2. A method is proposed how to make TW structures multipactor-free. Some design proposals are offered to facilitate fabrication. Further increase of the real-estate gradient (equivalent to 80 MV/m active gradient) is also possible by increasing the length of the accelerating structure because of higher group velocity and cell-to-cell coupling. Realization of this work opens paths to International Linear Collider energy upgrades beyond 1 to 3 TeV in competition with CLIC. The paper will discuss corresponding opportunities and challenges.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Remarkably Weak Anisotropy in Thermal Conductivity of Two-Dimensional Hybrid Perovskite Butylammonium Lead Iodide Crystals

Two-dimensional (2D) hybrid organic-inorganic perovskites consisting of alternating organic and inorganic layers are a new class of layered structures. They have attracted increasing interest for photovoltaic, optoelectronic, and thermoelectric applications, where knowing their thermal transport properties is critical. In this study, we carry out both experimental and computational studies on thermal transport properties of 2 D butylammonium lead iodide crystals and find their thermal conductivity is ultralow (below 0.3 W m( -1 ) K -1 ) with very weak anisotropy (around 1.5) among layered crystals. Further analysis reveals that the unique structure with the preferential alignment of organic chains and complicated energy landscape leads to moderately smaller phonon lifetimes in the out-of-plane direction and comparable phonon group velocities in in-plane and out-of-plane directions. These new findings may guide the future design of novel hybrid materials with desired thermal conductivity for various applications.

2D hybrid perovskite↗
DOE OSTI2024

Phase-space entropy cascade and irreversibility of stochastic heating in nearly collisionless plasma turbulence

We consider a nearly collisionless plasma consisting of a species of “test particles” in one spatial and one velocity dimension, stirred by an externally imposed stochastic electric field—a kinetic analog of the Kraichnan model of passive advection. The mean effect on the particle distribution function is turbulent diffusion in velocity space—known as stochastic heating. Accompanying this heating is the generation of fine-scale structure in the distribution function, which we characterize with the collisionless (Casimir) invariant C 2 ∝ ∫ ∫ d x d v 〈 f 2 〉 —a quantity that here plays the role of (negative) entropy of the distribution function. We find that C 2 is transferred from large scales to small scales in both position and velocity space via a phase-space cascade enabled by both particle streaming and nonlinear interactions between particles and the stochastic electric field. We compute the steady-state fluxes and spectrum of C 2 in Fourier space, with k and s denoting spatial and velocity wave numbers, respectively. In our model, the nonlinearity in the evolution equation for the spectrum turns into a fractional Laplacian operator in k space, leading to anomalous diffusion. Whereas even the linear phase mixing alone would lead to a constant flux of C 2 to high s (towards the collisional dissipation range) at every k , the nonlinearity accelerates this cascade by intertwining velocity and position space so that the flux of C 2 is to both high k and high s simultaneously. Integrating over velocity (spatial) wave numbers, the k -space ( s -space) flux of C 2 is constant down to a dissipation length (velocity) scale that tends to zero as the collision frequency does, even though the rate of collisional dissipation remains finite. The resulting spectrum in the inertial range is a self-similar function in the ( k , s ) plane, with power-law asymptotics at large k and s . Our model is fully analytically solvable, but the asymptotic scalings of the spectrum can also be found via a simple phenomenological theory whose key assumption is that the cascade is governed by a “critical balance” in phase space between the linear and nonlinear timescales. We argue that stochastic heating is made irreversible by this entropy cascade and that, while collisional dissipation accessed via phase mixing occurs only at small spatial scales rather than at every scale as it would in a linear system, the cascade makes phase mixing even more effective overall in the nonlinear regime than in the linear one. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Space-Time Structured Plasma Waves

Electrostatic waves play a critical role in nearly every branch of plasma physics from fusion to advanced accelerators, to astro, solar, and ionospheric physics. The properties of planar electrostatic waves are fully determined by the plasma conditions, such as density, temperature, ionization state, or details of the distribution functions. Furthermore we demonstrate that electrostatic wave packets structured with space-time correlations can have properties that are independent of the plasma conditions. For instance, an appropriately structured electrostatic wave packet can travel at any group velocity, even backward with respect to its phase fronts, while maintaining a localized energy density. These linear, propagation-invariant wave packets can be constructed with or without orbital angular momentum by superposing natural modes of the plasma and can be ponderomotively excited by space-time structured laser pulses like the flying focus.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Electronic interactions in Dirac fluids visualized by nano-terahertz spacetime interference of electron-photon quasiparticles

Ultraclean graphene at charge neutrality hosts a quantum critical Dirac fluid of interacting electrons and holes. Interactions profoundly affect the charge dynamics of graphene, which is encoded in the properties of its electron-photon collective modes: surface plasmon polaritons (SPPs). Here, we show that polaritonic interference patterns are particularly well suited to unveil the interactions in Dirac fluids by tracking polaritonic interference in time at temporal scales commensurate with the electronic scattering. Spacetime SPP interference patterns recorded in terahertz (THz) frequency range provided unobstructed readouts of the group velocity and lifetime of polariton that can be directly mapped onto the electronic spectral weight and the relaxation rate. Our data uncovered prominent departures of the electron dynamics from the predictions of the conventional Fermi-liquid theory. The deviations are particularly strong when the densities of electrons and holes are approximately equal. The proposed spacetime imaging methodology can be broadly applied to probe the electrodynamics of quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Characterizing warm atmospheric boundary layer over land by combining Raman and Doppler lidar measurements

PBL plays a critical role in the atmosphere by transferring heat, moisture, and momentum. The warm PBL has a distinct diurnal cycle including daytime convective mixing layer (ML) and nighttime residual layer developments. Thus, for PBL characterization and process study, simultaneous determinations of PBL height (PBLH) and ML height (MLH) are necessary. Here, new approaches are developed to provide reliable PBLH and MLH to characterize warm PBL evolution. The approaches use Raman lidar (RL) water vapor mixing ratio (WVMR) and Doppler lidar (DL) vertical velocity measurements at the Southern Great Plains (SGP) atmospheric observatory, which was established by the Atmospheric Radiation Measurement (ARM) user facility. Compared with widely used lidar aerosol measurements for PBLH, WVMR is a better trace for PBL vertical mixing. For PBLH, the approach classifies PBL water vapor structures into a few general patterns, then uses a slope method and dynamic threshold method to determine PBLH. For MLH, wavelet analysis is used to re-construct 2-D variance from DL vertical wind velocity measurements according to the turbulence eddy size to minimize the impacts of gravity wave and eddy size on variance calculations; then, a dynamic threshold method is used to determine MLH. Remotely-sensed PBLHs and MLHs are compared with radiosonde measurements based on the Richardson number method. Good agreements between them confirm that the proposed new algorithms are reliable for PBLH and MLH characterization. The algorithms are applied to warm seasons’ RL and ML measurements at the SGP site for five years to study warm season PBL structure and processes. The weekly composited diurnal evolutions of PBLHs and MLHs in warm climate were provided to illustrate diurnal and seasonal PBL evolutions. This reliable PBLH and MLH dataset will be valuable for PBL process study, model evolution, and PBL parameterization improvement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Export of macroinvertebrate prey from tidal freshwater wetlands provides a significant energy subsidy for outmigrating juvenile salmon

Tidal freshwater wetlands linking terrestrial, riverine, and saline habitats are critical areas for material processing and exchange. Once historically widespread, herbaceous marsh and forested tidal freshwater wetlands especially are now highly degraded worldwide. Additionally, quantitative assessments of hydrology and material exchange from these systems are lacking compared to lotic and estuarine (saltmarsh) habitats. Here we investigate macroinvertebrate and energy export from tidal marsh and forested wetlands and consider potential benefits from this ecological process to endangered Pacific salmon in a large tidal freshwater system, the Columbia River (USA). Macroinvertebrate (salmon prey) concentration, water velocity, and discharge were measured at several wetland habitat types (forested swamp, emergent marsh, and restored marsh). We used these data to compute prey flux and transport metrics. Then, applying literature values to calculate prey energy equivalents and juvenile salmon metabolic requirements, we estimated the potential energy subsidy available to juvenile salmon. Numerically, larval stages of aquatic insects were the predominant type of prey exported from the wetlands, with Diptera chironomid fly abundance exceeding other groups. Energetically, however, non-chironomid dipterans and hemipteran prey comprised most of energy transport due to their higher energetic content (energy density × mean weight). We determined the prey energy transported from the sampled tidal channels was sufficient to meet energetic needs of tens to thousands of juvenile salmon per day, depending on prey production and hydrography. The prey taxonomic composition differed among organisms exiting forested swamp, emergent marsh, and restored marsh habitats with corresponding differences in energy transport, but all habitat types supported similar numbers of juvenile salmon. We conclude that macroinvertebrate prey exported from varied tidal freshwater wetlands likely provide significant benefits to juvenile salmon over a larger ecological footprint than the wetland area would suggest.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

ARM SGP PBLH and MLH datasets from Raman lidar and Doppler lidar

The planetary boundary layer (PBL) plays a critical role in the atmosphere by transferring heat, moisture, and momentum. The warm PBL has a distinct diurnal cycle including the daytime convective mixing layer (ML) and nighttime residual layer developments. Thus, simultaneous determinations of PBL height (PBLH) and ML height (MLH) are necessary for studying PBL characterization and processes. Here, new approaches are developed to provide reliable PBLH and MLH estimates to characterize warm PBL evolution. The approaches use Raman lidar (RL) water vapor mixing ratio (WVMR) and Doppler lidar (DL) vertical velocity measurements at the Southern Great Plains (SGP) atmospheric observatory, which was established by the Atmospheric Radiation Measurement (ARM) User Facility. Compared to widely used lidar aerosol measurements for PBLH, WVMR is a better tracer for PBL vertical mixing. For PBLH, the approach classifies PBL water vapor structures into a few general patterns, then uses a slope method and dynamic threshold method to determine PBLH. For MLH, wavelet analysis is used to reconstruct 2D variance from DL vertical wind velocity measurements according to the turbulence eddy size to minimize the impacts of gravity wave and eddy size on variance calculations; then, a dynamic threshold method is used to determine MLH. Remotely-sensed PBLHs and MLHs are compared with radiosonde measurements based on the Richardson number method. Good agreements between them confirm that the proposed new algorithms are reliable for PBLH and MLH characterization. The algorithms are applied to warm-season RL and ML measurements at the SGP site for five years to study warm-season PBL structure and processes. The weekly composited diurnal evolutions of PBLHs and MLHs in a warm climate were provided to illustrate diurnal and seasonal PBL evolutions. This reliable data set of PBLH and MLH values will be valuable for studying PBL processes, model evolution, and PBL parameterization improvements. The MLH dataset includes the MLH in values of km above ground level. The PBLH dataset includes the PBLH in values of km above ground level, along with a flag ("situation_PBLH") to determine the state of the PBL (1 = Cloudy Condition, 2 = Stable Layer, 3 = Multi-layer WVMR structure, 4 = Well-Mixed PBL, 5 = A de-coupled layer, 6 = Other).

mixing layer height↗
DOE OSTI2024

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

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

22 GENERAL STUDIES OF NUCLEAR REACTORS↗