SEARCH · Engineering Papers

Results for “EVAPORATION MODEL”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Computational Study of Droplet Dynamics in Dispersed Flow Film Boiling Regime in PWR Geometries

Dispersed flow film boiling (DFFB) regime is exhibited during low-flooding rate conditions following the postulated Loss-of-Coolant Accident (LOCA) in nuclear reactor flow channels. LOCA is a type of design-basis accident scenario, which is crucial for reactor design and licensing. Accurate models to predict the heat transfer in DFFB regime is critical for the reactor to maintain safety limits. The nominal operating linear heat rate is also set based on the maximum linear heat rate determined for a LOCA transient. It is understood that the behavior of entrained droplets in DFFB regime greatly affects the heat transfer by evaporating and desuperheating the surrounding steam. Furthermore, interaction with flow structures also aids in increasing the evaporation rate of entrained droplets. In order to quantify this enhancement in heat transfer, the size of droplets downstream of flow structures (Spacer Grids and Mixing Vanes) is of crucial importance. Thermal hydraulic analysis codes like COBRA-TF, CATHARE-3 and adopt a three-field equation model that allows them to model the droplet field separately. However, to ensure the accuracy of these field equations, various correlations and physical models are required. This paper proposes a method for characterizing the mean Sauter Diameter of droplets across spacer-grids and mixing vanes in PWR sub-channels using PHASTA - a high fidelity flow solver with multiphase capabilities using the level-set method. The PHASTA code is a finite element based computational fluid dynamics code designed for solving compressible and incompressible turbulent flows. It has been demonstrated to be massively scalable (up to 786,432 cores), which is crucial for simulating large realistic domains. The code incorporates the level-set method for interface tracking of the gas-liquid interface.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2021

Boron nitride nanotube precursor formation during high-temperature synthesis: kinetic and thermodynamic modelling

In this work, we performed integrated modelling of the chemical pathways of formation for boron nitride nanotube (BNNT) precursors during high-temperature synthesis in a B/N 2 mixture. Integrated modelling includes quantum chemistry, Quantum–classical molecular dynamics, thermodynamic modelling, and kinetic approaches. We demonstrate that BN compounds are formed via the interaction of molecular nitrogen with small boron clusters, rather than through interactions with less reactive liquid boron. (This process can also be described as N 2 molecule fixation.) Liquid boron evaporates to produce these boron clusters (B m with m ≤ 5), which are subsequently converted into B m N n chains. The production of such chains is crucial to the growth of BNNTs because these chains form the building blocks of bigger and longer BN chains and rings, which are in turn the building blocks of fullborenes and BNNTs. Additionally, kinetic modelling revealed that B 4 N 4 and B 5 N 4 species in particular play a major role in the N 2 molecule fixation process. Furthermore, the formation of these species via reactions with B 4 and B 5 clusters is not adequately described under the assumption of thermodynamic equilibrium, as is demonstrated in our kinetic modelling. Thus, the accumulation of both B 4 N 4 and B 5 N 4 depends on the background gas pressure and the gas cooling rate. Long BN chains and rings, which are precursors of the fullborene and BNNT growth, form via self-assembly of components B 4 N 4 and B 5 N 4 . Our modelling results—particularly the increased densities of B 4 N 4 and B 5 N 4 species at higher gas pressures—explain the experimentally observed effect of gas pressure on the yield of high-quality BNNTs. The catalytic role of hydrogen was also studied; it is shown that HBNH molecules can be the main precursor of BNNT synthesis in the presence of hydrogen.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2018

Impacts of Broad-Scale Surface Freshening of the Southern Ocean in a Coupled Climate Model

The Southern Ocean surface has freshened in recent decades, increasing water column stability and reducing upwelling of warmer subsurface waters. The majority of CMIP5 models underestimate or fail to capture this historical surface freshening, yet little is known about the impact of this model bias on regional ocean circulation and hydrography. Here experiments are performed using a global coupled climate model with additional freshwater applied to the Southern Ocean to assess the influence of recent surface freshening. The simulations explore the impact of persistent and long-term broad-scale freshening as a result of processes including precipitation minus evaporation changes. Thus, unlike previous studies, the freshening is applied as far north as 55°S, beyond the Antarctic ice margin. It is found that imposing a large-scale surface freshening causes a surface cooling and sea ice increase under preindustrial conditions, because of a reduction in ocean convection and weakened entrainment of warm subsurface waters into the surface ocean. This is consistent with intermodel relationships between CMIP5 models and the simulations, suggesting that models with larger surface freshening also exhibit stronger surface cooling and increased sea ice. Additional experiments are conducted with surface salinity restoration applied to capture observed regional salinity trends. Remarkably, without any mechanical wind trend forcing, these simulations accurately represent the spatial pattern of observed surface temperature and sea ice trends around Antarctica. This study highlights the importance of accurately simulating changes in Southern Ocean salinity to capture changes in ocean circulation, sea surface temperature, and sea ice.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Surface heat flux and its association with the MJO in the tropical western Pacific using ARM observations

The Madden-Julian oscillation (MJO) is a major atmospheric phenomenon in the tropics that moves eastward every 20 to 100 days. It brings heavy rain and strong winds, influencing extreme weather events far beyond the tropics – including flooding, hurricanes, tornadoes and heavy snow in the United States. Therefore, an improved understanding of the MJO is critical for enhancing weather forecasts and supporting better decision-making for communities, emergency managers, and the private sector. Past studies, based on short-term observations or long-term model data, have emphasized the dominant role of atmospheric humidity in driving the MJO. Our research, using long-term observations (2000-2014) from three U.S. Department of Energy’s Atmospheric Radiation Measurement (ARM) sites in the tropical western Pacific (Manus, Nauru, and Darwin) confirms those past findings. However, we also show that atmospheric temperature, especially in the mid and upper troposphere, also plays a key role in the MJO’s evolution from its quiet (suppressed) to convective (active) phase. These findings provide valuable insights into how the MJO evolves and may help evaluate and improve weather and climate models. Our study also examined how the exchange of heat between the Earth's surface and the atmosphere, called surface heat flux, interacts with the MJO. This flux includes components such as latent heat (related to evaporation), sensible heat (related to temperature difference between surface and atmosphere), and radiation. We found that modulation in MJO convection is well connected to surface heat flux over the tropical western Pacific. One surface heat flux component often overlooked is the sensible heat flux due to precipitation (Q P ). Since falling raindrops are typically cooler than the land or ocean surface, they can cool the surface and affect local weather. To test the impact of Q P on convection during MJO, we incorporated it into a weather model, and ran simulations for two periods: April 2009 (when the MJO was active) and June 2006 (when it was not). Including Q P reduced errors in simulating the daily peak of rainfall—by 83% during the strong MJO and 23% during the inactive phase. It also improved the timing of peak daily rainfall and reduced the overall precipitation error by about 10%. These improvements were especially noticeable during periods of heavy rainfall. Our results suggest that even small heat exchanges from rainfall can play a significant role in shaping local weather. Accounting for these effects, particularly in regions like the tropical islands, can lead to more accurate simulations of the MJO-associated precipitation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2018

CAUSES: Diagnosis of the Summertime Warm Bias in CMIP5 Climate Models at the ARM Southern Great Plains Site

Abstract All the weather and climate models participating in the Clouds Above the United States and Errors at the Surface project show a summertime surface air temperature (T2 m) warm bias in the region of the central United States. To understand the warm bias in long‐term climate simulations, we assess the Atmospheric Model Intercomparison Project simulations from the Coupled Model Intercomparison Project Phase 5, with long‐term observations mainly from the Atmospheric Radiation Measurement program Southern Great Plains site. Quantities related to the surface energy and water budget, and large‐scale circulation are analyzed to identify possible factors and plausible links involved in the warm bias. The systematic warm season bias is characterized by an overestimation of T2 m and underestimation of surface humidity, precipitation, and precipitable water. Accompanying the warm bias is an overestimation of absorbed solar radiation at the surface, which is due to a combination of insufficient cloud reflection and clear‐sky shortwave absorption by water vapor and an underestimation in surface albedo. The bias in cloud is shown to contribute most to the radiation bias. The surface layer soil moisture impacts T2 m through its control on evaporative fraction. The error in evaporative fraction is another important contributor to T2 m. Similar sources of error are found in hindcast from other Clouds Above the United States and Errors at the Surface studies. In Atmospheric Model Intercomparison Project simulations, biases in meridional wind velocity associated with the low‐level jet and the 500 hPa vertical velocity may also relate to T2 m bias through their control on the surface energy and water budget.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2018

Does capillary evaporation limit the accessibility of nonaqueous electrolytes to the ultrasmall pores of carbon electrodes?

Porous carbons have been widely utilized as electrode materials for capacitive energy storage. Whereas the significance of pore size and geometry on the device performance has been well recognized, little guidance is available for identification of carbon materials with ideal porous structures. Here, we study the phase behavior of ionic fluids in slit pores using the classical density functional theory. Within the framework of the restricted primitive model for nonaqueous electrolytes, we demonstrate that the accessibility of micropores depends not only on the ionic diameters (or desolvation) but also on their wetting behavior intrinsically related to the vapor-liquid or liquid-liquid phase separation of the bulk ionic systems. Narrowing the pore size from several tens of nanometers to subnanometers may lead to a drastic reduction in the capacitance due to capillary evaporation. The wettability of micropores deteriorates as the pore size is reduced but can be noticeably improved by raising the surface electrical potential. The theoretical findings provide fresh insights into the properties of confined ionic systems beyond electric double layer models commonly employed for rational design/selection of electrolytes and electrode materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Development and validation of a model for predicting direct heat transfer from the fuel to droplets in the post dryout regime

Following dryout after a large loss of coolant accident, but before eventual final quenching of the fuel, the overheated fuel is cooled by a flow of superheated vapor with entrained saturated liquid droplets. The main mechanism for cooling the cladding is convective heat transfer to the superheated vapor, with this vapor in turn cooled by evaporation into it of the saturated droplets. An additional mechanism that is not normally considered explicitly is the cooling of the metal by direct droplet impingement of non-wetting droplets. It is difficult to measure and to analyze, and it is believed probably to contribute relatively little to the cooling, but it could be helpful additional cooling mechanism. In this paper is published validation of a simple analytical model of the heat removed from the wall during the few milliseconds that a non-wetting droplet is in 'contact' with the wall. It is found that the model's estimates of the heat removal are consistent with such limited experimental evidence as is available, paving the way for incorporation of the model into computations of post dryout heat transfer. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2014

Eddy Fluxes and Sensitivity of the Water Cycle to Spatial Resolution in Idealized Regional Aquaplanet Model Simulations

A multi-scale moisture budget analysis is used to identify the mechanisms responsible for the sensitivity of the water cycle to spatial resolution using idealized regional aquaplanet simulations. In the higher resolution simulations, moisture transport by eddies fluxes dry the boundary layer enhancing evaporation and precipitation. This effect of eddies, which is underestimated by the physics parameterizations in the low-resolution simulations, is found to be responsible for the sensitivity of the water cycle both directly, and through its upscale effect, on the mean circulation. Correlations among moisture transport by eddies at adjacent ranges of scales provides the potential for reducing this sensitivity by representing the unresolved eddies by their marginally resolved counterparts.

cycle↗

DOE OSTI2021

Evaluation of hybrid evaporative-vapor compression air conditioners for different global climates

Currently, a majority of global residential air-conditioning requirements are met using cooling systems based on the vapor compression cycle (VCC). To meet the future demand for space cooling and to reduce the corresponding environmental impact, there is a need for alternative cooling systems which require less resources and exhibit a higher coefficient of performance (COP). One of proposed alternatives is called the hybrid evaporative-vapor compression (HEVC) cycle, which utilizes adiabatic latent cooling in combination with the VCC. In this paper, we conducted a feasibility analysis of a HEVC system by comparing its performance to that of VCC for various global climates. To accomplish this, we developed and validated models to simulate the performance of residential VCC and HEVC systems. ASHRAE weather data and design conditions were then used to compare the performance of the two systems. In general, HEVC systems are best suited for hot arid climates with energy reduction greater than 20%. Conversely, humid climates are not suitable for HEVC adoption due to degraded performance of evaporative coolers in these climates. In addition to comparing the energy savings from HEVC systems, it is also critical to analyze their water consumption. Generally, climates that benefit the most from HEVC technology tend to also experience water scarcity issues. Thus, both the energy savings and water consumption of HEVC systems must be analyzed to guide the discussion on their adoption.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2019

Simulating a Mesoscale Convective System Using WRF with a New Spectral Bin Microphysics: Part 1: Hail vs Graupel

A modified Fast Spectral Bin Microphysics scheme (FSBM-2) embedded into the Weather Research and Forecasting Model (WRF) is used to simulate a mesoscale deep convective system observed during the Midlatitude Continental Convective Clouds Experiment (MC3E). FSBM-2 uses modified source codes as compared to the current FSBM (FSBM-1). In contrast to FSBM-1, FSBM-2 can simulate hail of several cm in diameter and includes additional processes such as spontaneous breakup of raindrops and aerosol regeneration by drop evaporation. It is shown that allowing large hail particles of diameters exceeding about 1 cm substantially increases the agreement between the simulated and observed squall-line structures both in the convective and stratiform regions. In contrast, if graupel particles are used to represent high-density hydrometeors in convective areas, the ratio of convective-to-stratiform areas diverges from the ratio seen in observations and maximum radar reflectivities are substantially underestimated. Analysis of snow size distributions in the stratiform area shows an important link between the ice crystals formed by homogeneous freezing in the convective area to ice particle number concentration in the stratiform region. Simulated raindrop size distributions in the stratiform area below the melting level from FSBM-2 show a good agreement with observations. The regeneration of CCN by droplet evaporation and detrainment of these CCN to the stratiform region increases the concentration of CCN there up to 8-9 km altitude; the additional CCN penetrate clouds and produce new droplets, leading to some convection intensification. Key words: Mesoscale convective systems (MCS), squall lines, spectral bin microphysics, large hail, in-cloud nucleation, size distributions of hydrometeors in MCS.

Shpund, Kobby↗

DOE OSTI2011

Role of the coupling between neck and radial degrees of freedom in evolution from dinucleus to mononucleus

The coupled Langevin equations in two-dimensional collective space are used to study the dynamics of nuclear neck growth. Special attention is paid to the effects of coupling between neck and radial degrees of freedom on the evolution from dinucleus to mononuleus. The dynamic model is applied for the study of neck evolution of the mass asymmetric system {sup 50}Ti + {sup 249}Cf. In order to estimate the effects of the coupling, we use the equations in the coupled and uncoupled cases. Our results show that the coupling between neck and radial motions reduces the neck growth velocity and delays the transition from dinuclear to mononuclear regimes. In addition, by solving these dynamic equations we get the probability distributions of radial degree of freedom at the injection point s{sub in} in the asymmetric fission valley. In this way, s{sub in} is no longer an adjustable parameter in the fusion-by-diffusion model. The distributions obtained are located at positions very close to s=0. Moreover, the coupling significantly reduces the fluctuation in the s{sub in} space. Our results show that the quasifissionlike events are unlikely to take place during the transition period for the mass asymmetric systems. Based on the s{sub in} distribution obtained, the evaporation residue (ER) cross sections for 3n and 4n evaporation channels in the {sup 50}Ti + {sup 249}Cf reaction leading to formation of {sup 296}120 and {sup 295}120 isotopes are calculated. The maximum ER cross sections in 3n and 4n channels with the s{sub in} distributions are equal to 0.1 and 0.065 pb, respectively, which are more than two times larger than those obtained by the fusion-by-diffusion model with the s{sub in}=0 assumption.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Dual-season thermal energy storage-integrated heat pump system for load shifting: experiments, modeling and nationwide analysis

Here, this paper proposes and validates a dual-season thermal energy storage-integrated heat pump (TES-HP) system that shifts building electric load in both cooling and heating modes using a single phase-change TES. The design uses the TES as a heat sink for the condenser in summer and a heat source for the evaporator in winter, thereby reducing on-peak electric demand while ensuring compatible with existing air-distribution systems. A 14-kW TES-HP prototype with six operating modes and novel refrigerant charge management strategies was experimentally tested under varied conditions. Subsequently, data-driven polynomial performance curves were developed and validated against experiments, then coupled with a phase change material (PCM) model and a calibrated building model within a rule-based controller. Single-day case studies demonstrated effective on-peak demand reduction while maintaining thermal comfort. During cooling, hourly power was lowered by up to 1.5 kW; in heating, average hourly power decreased from over 3.5 kW to 1.8 kW, compared to an HP-only system, achieving 64.5 % load shifting, and the need for resistance heating was eliminated. Seasonal simulations showed typical on-peak electricity savings of 25–35 % in cooling and 40–65 % in heating, with the largest benefits on the hottest and coldest days. Extended response-surface analysis and nationwide mapping quantified load shifting as functions of ambient and TES temperatures, revealing a seasonal trade-off in TES phase-change temperature. These results demonstrate the TES-HP system as a practical and scalable solution for grid-interactive buildings that reduces on-peak demand and electricity use while maintaining thermal comfort and enhancing grid flexibility and reliability.

Grid-interactive buildings↗

DOE OSTI2022

The Dependence of Mean Climate State on Shortwave Absorption by Water Vapor

Abstract State-of-the-art climate models exhibit significant spread in the climatological value of atmospheric shortwave absorption (SWA). This study investigates both the possible causes and climatic impacts of this SWA intermodel spread. The intermodel spread of global-mean SWA largely originates from the intermodel difference in water vapor shortwave absorptivity. Hence, we alter the water vapor shortwave absorptivity in the Community Earth System Model, version 1, with the Community Atmosphere Model, version 4 (CESM1-CAM4). Increasing the water vapor shortwave absorptivity leads to a reduction in global-mean precipitation and a La Niña–like cooling over the tropical Pacific. The global-mean atmospheric energy budget suggests that the precipitation is suppressed as a way to compensate for the increased SWA. The precipitation reduction is driven by the weakened surface winds, stabilized planetary boundary layer, and surface cooling. The La Niña–like cooling over the tropical Pacific is attributed to the zonal asymmetry of climatological evaporative damping efficiency and the low cloud enhancement over the eastern basin. Complementary fixed SSTs simulations suggest that the latter is more fundamental and that it primarily arises from atmospheric processes. Consistent with our experiments, the CMIP5/6 models with a higher global-mean SWA tend to produce tropical Pacific toward a more La Niña–like mean state, highlighting the possible role of water vapor shortwave absorptivity for shaping the mean-state climate patterns.

Kim, Hanjun↗

DOE OSTI2025

Sensitivities of Large Eddy Simulations of Aerosol Plume Transport and Cloud Response

Abstract Cloud responses to surface‐based sources of aerosol perturbation partially depend on how turbulent transport of the aerosol to cloud base affects the spatial and temporal distribution of aerosol. Here, scenarios of plume injection below a marine stratocumulus cloud are modeled using large eddy simulations coupled to a prognostic bulk aerosol and cloud microphysics scheme. Both passive plumes, consisting of an inert tracer, and active plumes are investigated, where the latter are representative of saltwater droplet plumes such as have been proposed for marine cloud brightening. Passive plume scenarios show higher in‐plume cloud brightness (relative to out‐of‐plume) due to the predominant transport of the passive plume tracer from the near‐surface to the cloud layer within updrafts. These updrafts rise into brighter areas within the cloud deck, even in the absence of an aerosol perturbation associated with an active plume. Comparing albedo at in‐plume to out‐of‐plume locations associates the inert plume with the brightest cloud locations, without any causal effect of the plume on the cloud. Numerical sensitivities are first assessed to establish a suitable model configuration. Then sensitivity to particle injection rate is investigated. Trade‐offs are identified between the number of injected particles and the suppressive effect of droplet evaporation on plume loft and spread. Furthermore, as the near‐field in‐plume brightening effect does not depend significantly on injection rate given a suitable definition of perturbed versus unperturbed regions of the flow, plume area is a key controlling factor on the overall cloud brightening effect of an aerosol perturbation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2017

Optical analysis of nanoparticle packing after drying in microdroplets

The optical spectra of nanostructured patterns with different degrees of order, obtained by dehydration self-assembly of nanoparticles in solution microdroplets placed on a planar substrate, have been analyzed. The reflectance spectrum of the best-ordered pattern contains a characteristic peak, the position of which is determined by the average lattice period in the layer. It is shown that the height of this peak depends on the degree of order of the pattern structure. Patterns with different degrees of order have been simulated using a software for modeling the self-assembly process, which takes into account the properties of particles, substrate, and solvent, as well as the contact line dynamics during solvent evaporation. A peak occurs in the reflectance spectrum due to the resonance Rayleigh scattering in the ordered layer. The described technique for optical diagnostics can be used to analyze the degree of particle order in a pattern.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Volatility Basis Set Distributions and Viscosity of Organic Aerosol Mixtures: Insights from Chemical Characterization Using Temperature-Programmed Desorption–Direct Analysis in Real-Time High-Resolution Mass Spectrometry

Quantitative assessment of gas-particle partitioning of individual components within complex atmospheric organic aerosols (OA) mixtures is critical for predicting and comprehending the formation and evolution of OA particles in the atmosphere. This investigation leverages previously documented data obtained through a temperature programmed desorption - direct analysis in real time – high resolution mass spectrometry (TPD-DART-HRMS) platform. This methodology facilitates the bottom-up construction of volatility basis set (VBS) distributions for constituents found in three biogenic secondary organic aerosol (SOA) mixtures produced through the ozonolysis of a-pinene, limonene, and ocimene. The apparent enthalpies (ΔH*, kJ mol -1 ) and saturated vapor mass concentrations (C T *, µg∙m -3 ) of individual SOA components, determined as a function of temperature (T, K), facilitated an assessment of changes in VBS distributions and gas-particle partitioning with respect to T and atmospheric total organic mass loadings (tOM, µg∙m -3 ). Further, the VBS distributions reveal distinct differences in volatilities among monomers, dimers, and trimers, enabling their categorization into separate volatility bins. At the ambient temperature of T = 298 K, only monomers efficiently partition between gas and particle phases across a broad range of atmospherically relevant total organic mass loadings (tOM) values of 1–100 µg∙m -3 . Partitioning of dimers and trimers becomes notable only at T > 360 K and T > 420 K, respectively. The viscosity of SOA mixtures is assessed using a bottom-up calculation approach, incorporating the input of elemental formulas, ΔH*, C T *, and particle-phase mass fractions of the SOA components. Through this approach, we are able to accurately estimate the variations in SOA viscosity that result from the evaporation of its components. These variations are, in turn, influenced by atmospherically relevant changes in tOM and T. Comparison of the calculated SOA viscosity and diffusivity values with literature reported experimental results shows close agreement, thereby validating the employed calculation approach. These findings underscore the significant potential for TPD-DART-HRMS measurements in enabling the untargeted analysis of organic molecules within OA mixtures. This approach facilitates quantitative assessment of their gas-particle partitioning and allows for the estimation of their viscosity and condensed-phase diffusion, thereby contributing valuable insights to atmospheric models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Complementarity and the unitarity of the black hole S -matrix

Recently, Akers et al. proposed a non-isometric holographic map from the interior of a black hole to its exterior. Within this model, we study properties of the black hole S-matrix, which are in principle accessible to observers who stay outside the black hole. Specifically, we investigate a scenario in which an infalling agent interacts with radiation both outside and inside the black hole. Because the holographic map involves postselection, the unitarity of the S-matrix is not guaranteed in this scenario, but we find that unitarity is satisfied to very high precision if suitable conditions are met. If the internal black hole dynamics is described by a pseudorandom unitary transformation, and if the operations performed by the infaller have computational complexity scaling polynomially with the black hole entropy, then the S-matrix is unitary up to corrections that are superpolynomially small in the black hole entropy. Furthermore, while in principle quantum computation assisted by postselection can be very powerful, we find under similar assumptions that the S-matrix of an evaporating black hole has polynomial computational complexity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2019

A simplified lithium vapor box divertor

A detached divertor plasma is predicted to be necessary to mitigate the heat fluxes and ion energy, and so sputtering, at the divertor plate of a fusion power plant. It has proven difficult in current experiments, however, to prevent the detachment front, once formed, from running up to the main plasma resulting in deterioration of pedestal performance. The lithium vapor box divertor localizes a dense cloud of lithium vapor away from the main plasma, in order to induce detachment at a stable, distant location. The vapor localization is created by local evaporation and nearby condensation, a configuration inaccessible for non-condensing seed impurities. This paper provides simulations of lithium vapor flow using the SPARTA direct simulation Monte Carlo (DSMC) code, in which the lithium vapor dynamics are governed by Li–Li collisions that should dominate in regions with little plasma. Inertia is included in this code, permitting the formation of shocks. The plasma model implemented in SPARTA is based on UEDGE results for the location of ionization and recombination. We find that a simplified lithium vapor box configuration, without baffles, provides robust stabilization of the detachment front and acceptable lithium vapor flow to the main chamber. Lithium is evaporated from a capillary porous surface on the private-flux side of the divertor, distant from the main plasma, and is condensed on the other side walls of the divertor chamber. The condensed lithium flows back to the evaporator through 2 cm internal diameter pipes. The capillary pressure differential at the surface of the evaporator is capable of overcoming MHD back-force in the piping, with a great deal of margin if a sandwich flow channel insert is implemented. Very little lithium should be accumulated on the first wall of the main chamber if it is maintained at 600 °C, as expected for a power-producing fusion system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗