Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Frit 625 Manufacturability Study

Based on the anticipated higher sodium concentration salt stream coming to the Defense Waste Processing Facility (DWPF) from the Salt Waste Processing Facility (SWPF), there may be a future need to reduce the alkali concentration in the frit depending on the composition of the sludge in Tank 40. SWPF is projected to start operating in late 2019. Thus, a lower alkali frit is beneficial for Sludge Batch (SB9) processing. Lower alkali frits may also be needed for Sludge Batch 10 (SB10) and later batches to ensure that the composition of the material in the Slurry Mix Evaporator (SME) is acceptable relative to the Product Composition Control System (PCCS) models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Preparation and Testing of Corrosion and Spallation-Resistant Coatings

This Energy & Environmental Research Center (EERC) project is designed to determine if plating APMT®, a specific highly oxidation-resistant oxide dispersion-strengthened FeCrAl alloy made by Kanthal, onto nickel-based superalloy turbine parts is a viable method for substantially improving the lifetimes and maximum use temperatures of the parts. The method for joining the APMT plate to the superalloys is called evaporative metal bonding and involves placing a thin foil of zinc between the plate and the superalloy, clamping them together, and heating in an atmosphere-controlled furnace. Upon heating, the zinc melts and dissolves the oxide skins of the alloys at the bond line, allowing the two alloys to diffuse into each other. The zinc then diffuses through the alloys and evaporates from their surfaces. During this annual reporting period, the finite element model was completed and used to design clamping jigs to hold the APMT plate to the larger blocks of superalloys during the bonding process. The clamping system was machined from titanium–zirconium–molybdenum and used to bond the APMT plate to the superalloy blocks. The bond between the APMT plate was weak for one of each of the superalloy blocks. We believe that this occurred because enough oxidation had occurred on the surface of the parts as a result of a 1-month time period between sandblasting to prepare the parts and the actual bonding process. The other blocks were, therefore, bonded within 1 day of preparing the parts for bonding, and their joints appear strong. Scanning electron microscopy analyses of representative joints showed that no zinc remained in the alloys after bonding. Also, phases rich in hafnium and tantalum had precipitated near the bond line in the APMT. Iron from the APMT had diffused into the superalloys during bonding, more extensively in the CM247LC than in the Rene 80. Nickel from the superalloys had diffused into the APMT, again more extensively in the joint with the CM247LC than with the Rene 80. One-inch-diameter buttons were machined from each of the bonded blocks and sent to Siemens for standard oxidation, spallation, and corrosion testing, which should be complete in the spring of 2016.

Hurley, John↗
DOE OSTI2023

Grounding our Understanding of the Impacts of Boreal Forest Expansion on Shallow Cumulus Clouds with a Simple Modeling Framework

Abstract The expansion of the boreal forest poleward is a potentially important driver of feedbacks between the land surface and Arctic climate. A growing body of work has highlighted the importance of differences in evaporative resistance between different possible future Arctic land covers, which in turn alters humidity and cloudiness in the boundary layer, for these feedbacks. While thus far this problem has been studied primarily with complex Earth system models, we turn to a locally focused, idealized model capable of diagnosing and testing the sensitivity of first-order processes connecting vegetation, the atmospheric boundary layer, and low clouds in this critical region. This allows us to benchmark the mechanisms and results at the center of predictions from larger-scale simulations. A surface dominated by broadleaf trees, characterized by higher albedo and lower surface evaporative resistance, drives cooling and moistening of the boundary layer relative to a surface of needleleaf trees, characterized by lower albedo and higher surface evaporative resistance. Differences in evaporative resistance between these hypothetical Arctic vegetation covers are of equal importance to changes in albedo for the initial response of the boundary layer to boreal expansion, even with our idealized approach. However, compensation between the elevation of the lifting condensation level (LCL) and more rapid growth of the mixed layer over higher evaporative resistance surfaces can minimize changes in the favorability of shallow clouds over different land cover types under some conditions. We then perform two tests on the sensitivity of this compensating effect, to changes in water availability, represented first by a reduction in boundary layer humidity and then by both a reduction in humidity and soil moisture available to our vegetation surface. Finally, given the importance of this potential LCL–mixed-layer height compensation in our idealized modeling results, we look to determine its relevance in observational data from a field campaign in boreal Finland. These observations do confirm that such a coupling plays an important role in cumulus-topped boundary layers over a needleleaf forest surface. While our results confirm some underlying mechanisms at the center of prior work with Earth system models, they also provide motivation for future work to constrain the impact of boreal forest expansion. This will include both large eddy simulations to examine the impact of processes and feedbacks not resolved by a mixed-layer model, as well as a more systematic evaluation and comparison of relevant observations at the site in Finland and sites from prior boreal field campaigns. Significance Statement Clouds and vegetation are both important components of the climate system that interact across a range of scales. These interactions are central to understanding how changes at the land surface feedback on climate. For example, if a forest expands or recedes, diagnosing how that will impact clouds will determine whether you predict warming or cooling temperatures from that shift in the forest area. These predictions are often made with complex Earth system models, but we look to a more idealized representation of the land–atmosphere system to diagnose how shallow clouds should respond to changes in surface properties with different scenarios of boreal forest expansion at a more foundational level. This both grounds our understanding of previous analysis and provides helpful direction for future studies of this relevant and impactful land cover change.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

Droplet Evaporation on Hot Micro-Structured Superhydrophobic Surfaces: Analysis of Evaporation from Droplet Cap and Base Surfaces

In this study, evaporation of sessile water droplets on hot micro-structured superhydrophobic surfaces is experimentally and theoretically investigated. Water droplets of 4 µL are placed on micro-pillared silicon substrates with the substrate temperature heated up to 120°C. A comprehensive thermal circuit model is developed to analyze the effects of substrate roughness and substrate temperature on the sessile droplet evaporation. For the first time, two components of heat and mass transfer, i.e., one from the droplet cap surface and the other from the droplet base surface, during droplet evaporation are distinguished and systematically studied. As such, the evaporation heat transfer rates from both the droplet cap surface and the interstitial liquid-vapor interface between micropillars at the droplet base are calculated in various conditions. For droplet evaporation on the heated substrates in the range of 40°C – 80°C, the predicted droplet cap temperature matches well with the experimental results. Furthermore, during the constant contact radius mode of droplet evaporation, the decrease of evaporation rate from the droplet base contributes most to the continuously decreasing overall evaporation heat transfer rate, whereas the decrease of evaporation rate from the droplet cap surface is dominant in the constant contact angle mode. The influence of internal fluid flow is considered for droplet evaporation on substrates heated above 100°C, and an effective thermal conductivity is adopted as a correction factor to account for the effect of convection heat transfer inside the droplet. Temperature differences between the droplet base and the substrate base are estimated to be about 2°C, 5°C, 8°C, 13°C and 18°C for droplet evaporation on substrates heated at 40°C, 60°C, 80°C, 100°C, and 120°C, respectively, elucidating the delayed or depressed boiling of water droplets on a heated rough surface due to evaporative cooling.

42 ENGINEERING↗
DOE OSTI2021

The Corrected Distortion model for Lagrangian spray simulation of transcritical fuel injection

In this work, we present a detailed implementation and validation of the droplet modeling framework proposed by Dahms and Oefelein (2016) into the engine commercial CFD software CONVERGE using the User Defined Function (UDF) interface. The model accounts for the nonlinear deformation and oscillation experienced by liquid spray droplet injected into high pressure and temperature. Lagrangian spray simulations of Engine Combustion Network (ECN) Spray A are performed. Model validation against standard experimental measurements of liquid velocity, vapor mixture fraction is conducted. To perform more rigorous model validation, new experimental measurements based on Diffused Back Illumination (DBI) are introduced. The new measurements are processed for Projected Liquid Volume (PLV), which offers as close as possible one-to-one model validation for liquid penetration while offering new insights into the spray physics. Comparison with a One-D model based on adiabatic mixing theory by Siebers (1999) and Desantes et al. (2007) are also conducted. Through these model validation exercises, it is shown that the new framework improves liquid-phase penetration predictions, following a tendency for enhanced evaporation, compared to the standard approach for both Reynolds Average Navier Stokes (RANS) and Large Eddy Simulation (LES). Additionally, at the liquid length, maximum mixture fraction values predicted by the new approach are in good agreement those of an adiabatic mixing model. Qualitative analysis of the spray behaviors during the early stage of the injection process reveals that the proposed framework predicts significant increase in droplet evaporation rate with lower drop drag compared to the current standard approach.

42 ENGINEERING↗
DOE OSTI2018

Fluid Dynamics Effects on Microstructure Prediction in Single Laser Tracks for Additive Manufacturing

The Laser Powder Bed Fusion Additive Manufacturing (LPBFAM) is one of the most important processes for the production of lightweight, cost-effective, complex, and high-performance enduse parts. At present, the cost and time associated with LPBFAM process development is very high due to a lack of fundamental process understanding. In this project, a multi-physics model was developed on a highly parallel open-source code, Truchas, with the ultimate goal of providing experimentally validated process maps for tailoring microstructure to achieve desired performance for LPBFAM. As a critical step towards fully LPBFAM modeling, modeling of single-track laser fusion (STLF) were conducted. Multi-physics simulations were conducted using Truchas for STLF by considering heat transfer, phase-change, fluid dynamics, surface tension phenomena, and evaporation. In order to assess the effect of the fluid dynamics on the solidification and ensuing microstructure, a heat transfer-and-solidification-only (HTS) model and a fully coupled heat-transfer-solidification and fluid-dynamics (HTS-FD) model were considered. In the HTS_FD model, the fluid flow was considered to be laminar while the molten alloy surface is assumed to be flat and non-deformable. This study is one the first attempt to understand the effect of the fluid flow on microstructure in STLF and LPBFAM. The results show that the fluid flow affects the solidification and ensuing microstructure in STLF. In order to validate Truchas for STLF and LPBFAM modeling, experimental data, which was obtained at GE Global Research (GEGR) for liquid pool shape and microstructure, were compared with those from numerical simulation results. For the keyhole regime case, numerical simulation results indicate that the melt-pool shape was typical to that of the conduction case with very large deviations from the measured melt-pool depths. Two sensitivity studies were conducted by varying the evaporation flux and the surface tension coefficient at fixed laser absorptivities of 0.5 and 0.3, respectively. The minimum value for the overall combined error, between the calculated values with the HTS_FD model and measured values for both the width and height of the melt-pool, was attained for a surface tension coefficient of -0.75e-4 N/m. An analytical solidification model for the columnar-to-equiaxed transition (CET) in rapid solidification was used to assess the microstructure variation within the melt-pool. From numerical simulation results, thermal gradient, G, and solidification velocity, V, were obtained in order to predict the microstructure type (e.g., dendritic, cellular). For the fluid-flow HTS_FD model, G was found to exhibit a maximum at the extremity of the solidified pool (i.e., at the free surface). By contrast, for HTS simulations, G was found to exhibit a maximum around the entire edge of the solidified pool. For the HTS_FD simulations, the minimum values of the cooling rate, GV, were found to be approximately half than their corresponding values for HTS simulations. By contrast with the min GV values, the maximum values for GV were found to occur for the HTS_FD simulations. Thus, HTS_FD simulations were found to exhibit a wider range of cooling rates than the HTS simulations. Concerning the solidification maps, the variation of the solidification velocity, V, as function of the thermal gradient, G, was obtained. It was found that the fluid flow model results (HTS_FD) exhibited an increased spread in the V(G) variation within the melt pool with respect to the HTS model results (without the fluid flow). A preliminary correlation model for the primary dendrite arm spacing (PDAS) based on a power law dependence on thermal gradient and solidification velocity was used to estimate PDAS from HTS and HTS_FD simulations. For low laser powers and low laser speeds, the PDAS obtained 6 with the fluid dynamics model (HTS_FD) was larger by more than 30% with respect to the PDAS calculated with the simple HTS model. In the second part of the project, fifty-seven simulations were conducted in order to obtain data on microstructure variables that can be used for process map development. In order to cover the entire processing space, thirty-three cases were selected STLF process simulations at six power levels and six laser scanning speeds. The 57 STLF process simulations were conducted as follows: 33 simulations with the heat-transfer-only (HTS) model and 24 simulations with the fluid flow model (HTS_FD). The solidification map data shows that for all simulations, a columnar dendritic microstructure would be expected. It was found that the minimum, average, and maximum thermal gradient exhibit exponential variations with respect to a process variable defined as the ratio between the power and square root of scan speed.

36 MATERIALS SCIENCE↗
DOE OSTI2018

The Role of Convective Gustiness in Reducing Seasonal Precipitation Biases in the Tropical West Pacific

Abstract Precipitation is an important climate quantity that is critically relevant to society. In spite of intense efforts, significant precipitation biases remain in most climate models. One pervasive and persistent bias found in many general circulation models occurs in the Tropical West Pacific where northern hemisphere summer‐time precipitation is often underestimated compared to observations. Using the DOE‐E3SM model, the inclusion of a missing process, convective gustiness, is shown to reduce those biases through a net increase in surface evaporation. Gustiness in surface wind fields is assumed to arise empirically in proportion to the intensity of convective precipitation. The increased evaporation can be treated as an increase in the moist static energy forcing into the atmosphere. A Normalized Gross Moist Stability (NGMS) framework (which characterizes the relationship between convective forcing and convective response) is used to explore the processes responsible for the precipitation bias, and the impact of the gustiness parameterization in reducing that bias. Because the NGMS of the Tropical West Pacific is less than unity in the E3SMv1 model, the increase in energy forcing amplifies the increase in precipitation to exceed that of the evaporative flux. Convective gustiness favors increased precipitation in regions where the resolved surface winds are weak and convection is present.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Dynamic domain kinematic modelling for predicting interflow over leaky impeding layers

Traditional Boussinesq or kinematic simulations of interflow (i.e., lateral subsurface flow) assume no leakage through the impeding layer and require a no-flow boundary condition at the ridge top. However, recent analyses of many interflow-producing landscapes indicate that leaky impeding layers are common, that most interflow percolates well before reaching the toe slope, and therefore, the downslope contributing length is shorter than the hillslope length. In watersheds characterised by perched interflow over a low conductivity layer through permeable topsoil, interflow with percolation may be modelled with a kinematic wave model using a mobile upslope boundary condition defining the hillslope portion contributing interflow to valleys. Here, we developed and applied a dynamic interflow model to simulate interflow using a downslope travel distance concept such that only the active contributing length is modelled at any time. The model defines a variable active area based on the depth of the perched layer, the topographic slope and the ratio of the hydraulic conductivity of topsoil to that of the impeding layer. It incorporates a two-layer soil moisture accounting water balance analysis, a pedo-transfer function, and percolation and evaporation routines to predict interflow rates in continuous and event-based scenarios. We tested the modelling concept on two sets of data (2-year dataset of rainfall observations for the continuous simulation and a multi-day irrigation experiment for the event simulation) from a 121-m-long open interflow collection trench on an experimental hillslope at the Savannah River Site, South Carolina. The continuous model simulation partially represented the observed interflow hydrograph and perched water depth in the experimental hillslope with correlation coefficients of 0.85 and 0.35, respectively. Model performance improved significantly at event-scale analysis. The modelling approach realistically represents interflow dynamics in hillslopes with leaky impeding layers and can be integrated into catchment-scale hydrology models for more detailed hillslope process modelling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

STAR CLUSTER FORMATION WITH STELLAR FEEDBACK AND LARGE-SCALE INFLOW

During star cluster formation, ongoing mass accretion is resisted by stellar feedback in the form of protostellar outflows from the low-mass stars and photo-ionization and radiation pressure feedback from the massive stars. We model the evolution of cluster-forming regions during a phase in which both accretion and feedback are present and use these models to investigate how star cluster formation might terminate. Protostellar outflows are the strongest form of feedback in low-mass regions, but these cannot stop cluster formation if matter continues to flow in. In more massive clusters, radiation pressure and photo-ionization rapidly clear the cluster-forming gas when its column density is too small. We assess the rates of dynamical mass ejection and of evaporation, while accounting for the important effect of dust opacity on photo-ionization. Our models are consistent with the census of protostellar outflows in NGC 1333 and Serpens South and with the dust temperatures observed in regions of massive star formation. Comparing observations of massive cluster-forming regions against our model parameter space, and against our expectations for accretion-driven evolution, we infer that massive-star feedback is a likely cause of gas disruption in regions with velocity dispersions less than a few kilometers per second, but that more massive and more turbulent regions are too strongly bound for stellar feedback to be disruptive.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Cost-constrained adaptive simulations of transient spray combustion in a gas turbine combustor

Predictive high-fidelity simulations of turbulent spray combustion must capture the combined effects of complex chemistry, multiphase evaporating flow and spray-flame interactions to achieve physical accuracy. Finite-rate chemistry (FRC) combined with a realistic chemical mechanism is a combustion model well-suited for this purpose, but has a high computational cost due to the large number and stiffness of transported chemical species. In contrast, flamelet-based models achieve lower cost by transporting a small number of quantities of reduced stiffness, but assumptions regarding local flame topology, boundary conditions and inter-phase coupling limit their physical accuracy. Recently, the Pareto-efficient combustion (PEC) framework was developed to dynamically assign combustion models based on local cost and accuracy metrics in gas-phase reacting flows. In this work, we extend this PEC framework to spray combustion through the rigorous analysis of the multiphase coupling terms in the governing equations. The derivation shows that spray evaporation causes errors in the prediction of species mass fractions for flamelet-based models due to the sensitivity of the local thermo-chemical state to changes in composition caused by fuel vaporization across combustion regimes present in practical spray combustion devices. Sub-model assignment is formulated as a multiple-choice knapsack problem, where computational cost is directly controlled through the fraction of the domain assigned to the FRC sub-model. The extended PEC formulation is applied to the simulation of a realistic rich-quench-lean gas turbine combustor at steady-state conditions, as well as transient operation resulting in lean blow-out (LBO). Analysis of transient simulations during LBO demonstrates the extended PEC formulation’s capacity to dynamically adapt to changing conditions within the combustor. Transient combustor dynamics are shown to approach convergence with limited increases in computational cost, while retaining substantial computational cost reduction compared to monolithic FRC simulations. Through PEC simulations with increasing fractions of the domain assigned to FRC, monolithic flamelet simulations are shown to over-predict flame stability during LBO. Finally, the extended PEC formulation is thus shown to overcome deficiencies of monolithic models by controlling modeling error for multiphase combustion modeling.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2011

Modelling and Simulation of the Advanced Plasma Source

Plasma ion assisted-deposition (PIAD) is a combination of conventional thermal evaporation deposition and plasma-beam surface modification; it serves as a well-established technology for the creation of high quality coatings on mirrors, lenses, and other optical devices. It is closely related to ion-assisted deposition to the extent that electrons preserve quasineutrality of the ion beam. This paper investigates the Advanced Plasma Source (APS), a plasma beam source employed for PIAD. A field enhanced glow discharge generates a radially expanding plasma flow with an ion energy of about 80-120 eV. Charge exchange collisions with the neutral background gas (pressure 0.1 Pa and below) produce a cold secondary plasma, which expands as well. A model is developed which describes the primary ions by a simplified Boltzmann equation, the secondary ions by the equations of continuity and momentum balance, and the electrons by the condition of Boltzmann equilibrium. Additionally, quasineutrality is assumed. The model can be reduced to a single nonlinear differential equation for the velocity of the secondary ions, which has several removable singularities and one essential singularity, identified as the Bohm singularity. Solving the model yields macroscopic plasma features, such as fluxes, densities, and the electrical field. An add-on Monte-Carlo simulation is employed to calculate the ion energy distribution function at the substrate. All results compare well to experiments conducted at a commercial APS system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Laboratory Evaporation Testing Of Hanford Waste Treatment Plant Low Activity Waste Off-Gas Condensate Simulant

The Hanford Waste Treatment and Immobilization Plant (WTP) Low Activity Waste (LAW) vitrification facility will generate an aqueous condensate recycle stream, LAW Off-Gas Condensate, from the off-gas system. The baseline plan for disposition of this stream is to send it to the WTP Pretreatment Facility, where it will be blended with LAW, concentrated by evaporation and recycled to the LAW vitrification facility again. Alternate disposition of this stream would eliminate recycling of problematic components, and would enable de-coupled operation of the LAW melter and the Pretreatment Facilities. Eliminating this stream from recycling within WTP would also decrease the LAW vitrification mission duration and quantity of canistered glass waste forms. This LAW Off-Gas Condensate stream contains components that are volatile at melter temperatures and are problematic for the glass waste form. Because this stream recycles within WTP, these components accumulate in the Condensate stream, exacerbating their impact on the number of LAW glass containers that must be produced. Approximately 32% of the sodium in Supplemental LAW comes from glass formers used to make the extra glass to dilute the halides to be within acceptable concentration ranges in the LAW glass. Diverting the stream reduces the halides in the recycled Condensate and is a key outcome of this work. Additionally, under possible scenarios where the LAW vitrification facility commences operation prior to the WTP Pretreatment facility, identifying a disposition path becomes vitally important. This task examines the impact of potential future disposition of this stream in the Hanford tank farms, and investigates auxiliary evaporation to enable another disposition path. Unless an auxiliary evaporator is used, returning the stream to the tank farms would require evaporation in the 242-A evaporator. This stream is expected to be unusual because it will be very high in corrosive species that are volatile in the melter (chloride, fluoride, sulfur), will have high ammonia, and will contain carryover particulates of glass-former chemicals. These species have potential to cause corrosion of tanks and equipment, precipitation of solids, release of ammonia gas vapors, and scale in the tank farm evaporator. Routing this stream to the tank farms does not permanently divert it from recycling into the WTP, only temporarily stores it prior to reprocessing. Testing is normally performed to demonstrate acceptable conditions and limits for these compounds in wastes sent to the tank farms. The primary parameter of this phase of the test program was measuring the formation of solids during evaporation in order to assess the compatibility of the stream with the evaporator and transfer and storage equipment. The origin of this LAW Off-Gas Condensate stream will be the liquids from the Submerged Bed Scrubber (SBS) and the Wet Electrostatic Precipitator (WESP) from the LAW facility melter offgas system. The stream is expected to be a dilute salt solution with near neutral pH, and will likely contain some insoluble solids from melter carryover. The soluble components are expected to be mostly sodium and ammonium salts of nitrate, chloride, and fluoride. This stream has not been generated yet, and, thus, the composition will not be available until the WTP begins operation, but a simulant has been produced based on models, calculations, and comparison with pilot-scale tests. This report discusses results of evaporation testing of the simulant. Two conditions were tested, one with the simulant at near neutral pH, and a second at alkaline pH. The neutral pH test is comparable to the conditions in the Hanford Effluent Treatment Facility (ETF) evaporator, although that evaporator operates at near atmospheric pressure and tests were done under vacuum. For the alkaline test, the target pH was based on the tank farm corrosion control program requirements, and the test protocol and equipment was comparable to that used for routine evaluation of feed compatibility studies for the 242-A evaporator. One of the radionuclides that is volatile in the melter and expected to be in high concentration in this LAW Off-Gas Condensate stream is Technetium-99 (99Tc). Technetium will not be removed from the aqueous waste in the Hanford WTP, and will primarily end up immobilized in the LAW glass by repeated recycle of the off-gas condensate into the LAW melter. Other radionuclides that are also expected to be in appreciable concentrations in the LAW Off-Gas Condensate are 129I, 90Sr, 137Cs, and 241Am. The concentrations of these radionuclides in this stream will be much lower than in the LAW, but they will still be higher than limits for some of the other disposition pathways currently available. At this time, these scoping tests did not evaluate the partitioning of the radionuclides to the evaporator condensate, since ample data are available separately from other experience in the DOE complex. Results from the evaporation testing show that the neutral SBS simulant first forms turbidity at ~7.5X concentration, while the alkaline-adjusted simulant became turbid at ~3X concentration. The major solid in both cases was Kogarkoite, Na3FSO4. Sodium and lithium fluorides were also detected. Minimal solids were formed in the evaporator bottoms until a substantial fraction of liquid was removed, indicating that evaporation could minimize storage volume issues. Achievable concentration factors without significant insoluble solids were 17X at alkaline pH, and 23X at neutral pH. In both runs, significant ammonia carried over and was captured in the condenser with the water condensate. Results also indicate that with low insoluble solids formation in the initial testing at neutral pH, the use of Reverse Osmosis is a potential alternate method for concentrating the solution, although an evaluation is needed to identify equipment that can tolerate insoluble solids. Most of the ammonia remains in the evaporator bottoms during the neutral pH evaporation, but partitions to the condensate during alkaline evaporation. Disposition of both streams needs to consider the management of ammonia vapor and its release. Since this is an initial phase of testing, additional tasks related to evaporation methods are expected to be identified for development. These tasks likely include evaluation and testing of composition variability testing and evaluations, corrosion and erosion testing, slurry storage and immobilization investigations, and evaporator condensate disposition.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

n-Propylbenzene Droplet Combustion: Experiments and Numerical Modeling

n-Propylbenzene (PB, C9H12, 432K boiling point) is important as a constituent in surrogates for some transportation fuels but little is understood about its droplet burning characteristics. We present an experimental and numerical study of a fuel droplet burning under conditions where the sole means of gas transport is fuel evaporation at the droplet surface thereby creating a one-dimensional gas transport that is amenable to detailed numerical modeling. The one-dimensional gas transport that results is amenable to detailed numerical modeling. Initial droplet diameters (Do) ranged from 0.3 mm to about 5 mm to better understand the mechanisms that control burning over such a wide range of Do. For Do < 1 mm droplets were burned while anchored to 14 m SiC fibers within a sealed container in free-fall during which time the droplets were ignited by spark discharge and the burning process recorded by video cameras. Larger and free-floating droplets were studied in the spaced-based platform of the orbiting International Space Station for Do > 1.0 mm which provided unlimited times for experimental observation. Measurements of droplet, flame and soot shell diameters were made and are discussed in the presentation. The experimental results showed formation of optically thick soot clouds developing during burning that required a robust image analysis algorithm to enable measurement of droplet diameter. The measured droplet diameters were well predicted by a numerical model that included soot formation, a PB kinetic mechanism comprising 576 species and 27369 reactions, unsteady gas and liquid transport, variable properties and radiative transport. Flame diameters showed more variability compared to the simulations, depending on how the flame was defined (peak gas temperature and OH concentration). Droplets burned slower as Do was increased and the PB sooting propensity increased with Do. Flame extinction was observed for Do > 3.3 mm indicating a strong effect of radiation as Do was increased compared to the smaller droplet diameters examined where extinction was not experimentally observed. For the largest droplets examined the possibility that the observed extinction is indicative of a transition to cool flame burning is still under investigation.

Guo, Songtao↗
DOE OSTI2024

Using an Isotope Enabled Mass Balance to Evaluate Existing Land Surface Models

Abstract Land surface models (LSMs) play a crucial role in elucidating water and carbon cycles by simulating processes such as plant transpiration and evaporation from bare soil, yet calibration often relies on comparing LSM outputs of landscape total evapotranspiration ( ET ) and discharge with measured bulk fluxes. Discrepancies in partitioning into component fluxes predicted by various LSMs have been noted, prompting the need for improved evaluation methods. Stable water isotopes serve as effective tracers of component hydrologic fluxes, but data and model integration challenges have hindered their widespread application. Leveraging National Ecological Observation Network measurements of water isotope ratios at 16 US sites over 3 years combined with LSM‐modeled fluxes, we employed an isotope‐enabled mass balance framework to simulate ET isotope values ( δET ) within three operational LSMs (Mosaic, Noah, and VIC) to evaluate their partitioning. Models simulating δET values consistent with observations were deemed more reflective of water cycling in these ecosystems. Mosaic exhibited the best overall performance (Kling‐Gupta Efficiency of 0.28). For both Mosaic and Noah there were robust correlations between bare soil evaporation fraction and error (negative) as well as transpiration fraction and error (positive). We found the point at which errors are smallest ( x ‐intercept of the multi‐site regression) is at a higher transpiration fraction than is currently specified in the models. Which means that transpiration fraction is underestimated on average. Stable isotope tracers offer an additional tool for model evaluation and identifying areas for improvement, potentially enhancing LSM simulations and our understanding of land‐surface hydrologic processes.

58 GEOSCIENCES↗
DOE OSTI2021

Towards a Techno-Economic Analysis of PCM Integrated Hybrid HVAC Systems

Thermal end uses dominate building energy consumption and are a major driver of peak demand. As heating is electrified, peak electrical power required will surge, prompting a need for innovative HVAC system designs and controls.These designs must incorporate novel technologies at the component level and new integration techniques at the system level. One such possibility involves the addition of thermal energy storage (TES) in heating and cooling equipment using a phase change material (PCM) heat exchanger. Here, thermal energy storage via phase change can be used to shift the HVAC system loads to times of lower electricity cost, reduced carbon intensity, and greater energy efficiency. Most of the current utilization of PCM in buildings involves passive components. By actively controlling when heat is stored and released from PCM, we can optimize the building HVAC system to cost-effectively meet consumer needs with the flexibility to draw on renewable energy resources when they are abundant and available. While this combination of technologies is promising in theory, simulation-based evaluation of a prototype can be difficult due to the modeling requirements at the component level and the large number of possible configurations and operating modes at the system level. To conduct this evaluation, we use the Modelica language for modeling and simulation because it enables users to represent the important physics of the problem, interchange and rearrange components in an efficient manner, and implement a range of control configurations. In this work, we considered three case studies: a portable building, a large commercial retail store, and a multifamily residential apartment unit. Each of these employs a different system design, ranging from a single package vertical unit incorporating PCM to a central plant with independent heat pump, evaporative cooling, and thermal energy storage components. This paper describes the technologies in question, presents modeling at the component and system levels, and demonstrates building energy and demand charge cost savings with local time-of-use tariffs in a hot-dry climate.

Helmns, Dre↗
DOE OSTI2024

Modeled estimates of permafrost hydrology and related fields for pan-Arctic region over the period 1980-2100

This dataset contains gridded estimates of total runoff (mm/day); subsurface runoff (mm/day); snow water equivalent (SWE, mm/day); vertical upward water flux of combined surface evaporation, evapotranspiration, and snow sublimation (mm/day); and soil temperature (C) across the pan-Arctic basin from simulations of the Permafrost Water Balance Model (PWBM v4). These data were generated to allow researchers to advance understanding of the changes unfolding in the Arctic's terrestrial water cycle under the influences of climate warming, hydrological cycle intensification, and permafrost thaw. The dataset will enable researchers to address important questions spanning multiple disciplines. For example, how will rising river flows exported to the coast impact critical zone ecosystems, near-shore ice dynamics, and ocean biogeochemical processes. Spatial resolution is 25x25 km with N=39926 grid cells spanning the domain. Three forcing datasets were used to drive simulations at an implicit daily time step over the period 1980-2100. Daily air temperature, precipitation, and wind speed data from the WATCH Forcing Data methodology applied to ERA5 reanalysis (W5E5) were used for meteorological forcings over the period 1980-2019. Simulations for the period 1980-2100 used meteorological data from two global climate models; IPSL-CM6A-LR and MPI-ESM1-2-HR. The PWBM hydrological outputs are organized by grid and month in gzipped ASCII text files within each archived tar file. Soil temperature is organized by grid cell and day of year from DOY 205 to DOY 260. No special software is required to work with these data. The package contains spatial domain and river basin information files. These data are most useful for analyses of the spatial and temporal variations in terrestrial hydrology and soil thermal dynamics. A manuscript describing the model simulations and associated analysis has been accepted for publication in The Cryosphere (Rawlins and Karmalkar, 2024).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

Utilization of municipal wastewater for cooling in thermoelectric power plants

A process simulation model has been developed using Aspen Plus® with the OLI (OLI System, Inc.) water chemistry model to predict water quality in the recirculating cooling loop utilizing secondary- and tertiary-treated municipal wastewater as the source of makeup water. Simulation results were compared with pilot-scale experimental data on makeup water alkalinity, loop pH, and ammonia evaporation. The effects of various parameters including makeup water quality, salt formation, NH 3 and CO 2 evaporation mass transfer coefficients, heat load, and operating temperatures were investigated. The results indicate that, although the simulation model can capture the general trends in the loop pH, experimental data on the rates of salt precipitation in the system are needed for more accurate prediction of the loop pH. It was also found that stripping of ammonia and carbon dioxide in the cooling tower can influence the cooling loop pH significantly. The effects of the NH 3 mass transfer coefficient on cooling loop pH appear to be more significant at lower values (e.g., k NH3 < 4×10 -3 m/s) when the makeup water alkalinity is low (e.g., <90 mg/L as CaCO 3 ). The effect of the CO2 mass transfer coefficient was found to be significant only at lower alkalinity values (e.g., k CO2 <4×10 -6 m/s).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

Analysis of 2H-Evaporator Acid Cleaning Samples

The 2H-Evaporator acid cleaning solution samples were analyzed by SRNL to determine a composition for the scale present in the evaporator before recent acid cleaning. Composite samples were formed from the solution samples from the two acid cleaning cycles. The solution composition was converted to a weight percent scale solids basis under an assumed chemical composition. The scale composition produced from the acid cleaning solution samples indicates a concentration of 6.85 wt% uranium. An upper bound, onesided 95% confidence interval on the weight percent uranium value may be given as 6.9 wt% + 1.645 × 0.596 wt% = 7.9 wt%. The comparison of the composition from the current acid cleaning solutions with the composition of recent scale samples along with the thermodynamic modeling results provides reasonable assurance that the sample results provide a good representation of the overall scale composition in the evaporator prior to acid cleaning. The small amount of scale solids dissolved in the 1.5 M nitric acid during the evaporator cleaning process likely produced only a small amount of precipitation based on modeling results and the visual appearance of the samples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗