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 523 records · Page 29

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Strong Local Evaporative Cooling Over Land Due to Atmospheric Aerosols

Abstract Aerosols can enhance terrestrial productivity through increased absorption of solar radiation by the shaded portion of the plant canopy—the diffuse radiation fertilization effect. Although this process can, in principle, alter surface evaporation due to the coupling between plant water loss and carbon uptake, with the potential to change the surface temperature, aerosol‐climate interactions have been traditionally viewed in light of the radiative effects within the atmosphere. Here, we develop a modeling framework that combines global atmosphere and land model simulations with a conceptual diagnostic tool to investigate these interactions from a surface energy budget perspective. Aerosols increase the terrestrial evaporative fraction, or the portion of net incoming energy consumed by evaporation, by over 4% globally and as much as ∼40% regionally. The main mechanism for this is the increase in energy allocation from sensible to latent heat due to global dimming (reduction in global shortwave radiation) and slightly augmented by diffuse radiation fertilization. In regions with moderately dense vegetation (leaf area index >2), the local surface cooling response to aerosols is dominated by this evaporative pathway, not the reduction in incident radiation. Diffuse radiation fertilization alone has a stronger impact on gross primary productivity (+2.18 Pg C y −1 or +1.8%) than on land evaporation (+0.18 W m −2 or +0.48%) and surface temperature (−0.01 K). Our results suggest that it is important for land surface models to distinguish between quantity (change in total magnitude) and quality (change in diffuse fraction) of radiative forcing for properly simulating surface climate.

Chakraborty, TC↗

DOE OSTI2026

Xanthos-Lake Dataset

The Xanthos-Lake v1.0 dataset provides the input data, trained machine-learning models, and simulation outputs needed to characterize lake water balance, snow and ice conditions, and mixing-layer temperature within the Xanthos global hydrological modeling framework. The dataset supports lake representation across a wide range of lake sizes and hydroclimatic conditions by combining xLSIM, a basin-specific machine-learning emulator of lake snow, ice, ice-cover fraction, and mixing-layer temperature, with the Xanthos-Lake water-balance model. The archive contains NetCDF datasets used to train and evaluate xLSIM, trained model weights, processed meteorological and lake-property inputs, and basin- and lake-category-specific simulation outputs. These materials are organized into four primary data groups, described below. Snowice_model_inputs: Contains the NetCDF input data used to train xLSIM. The xLSIM machine-learning framework uses three lake-based datasets. The meteorological forcing dataset provides monthly relative humidity, specific humidity, surface wind speed, maximum and minimum air temperature, downward longwave and shortwave radiation, snowfall, surface air pressure, and total precipitation. Lake surface area is included as an additional static predictor. The target-state dataset provides lake ice thickness, snow depth, snow cover, and lake mixing-layer temperature, while a companion lake-surface dataset provides the lake ice-cover fraction. Before training, ice thickness and snow depth are converted from meters to centimeters, mixing-layer temperature is converted from kelvin to degrees Celsius and constrained to nonnegative values, and ice-cover fraction is converted from a fraction to a percentage. The predictor variables are normalized using statistics calculated across the selected lakes and time steps. Snowice_model_outputs: Contains the NetCDF outputs generated by xLSIM. For each basin, xLSIM produces a file containing observed and predicted lake-state variables for the training, validation, and testing periods. The modeled variables include lake ice thickness, snow depth, snow cover, mixing-layer temperature, and lake ice-cover fraction. For basins without a sufficiently persistent snow-and-ice signal, the emulator predicts only mixing-layer temperature. The outputs also include training and validation loss histories, the selected model configuration, identifiers of the lakes used in training, and SHAP-based feature-importance information at the global, lake, and seasonal-regime levels. The trained machine-learning model weights are provided separately within the dataset archive. Together, these files support model evaluation and subsequent coupling with the Xanthos-Lake water-balance framework. XanthosLAKES: Contains the NetCDF input data used by the Xanthos-Lake framework. Monthly meteorological inputs include relative and specific humidity, downward shortwave and longwave radiation, mean, maximum, and minimum air temperature, wind speed, precipitation, snowfall, and surface air pressure. Static lake-property datasets provide lake identifiers, geographic locations, surface area, volume, mean depth, elevation, drainage area, fetch, outlet-routing information, and associated Xanthos grid-cell attributes. Separate bathymetric datasets provide the coefficients of the area–depth and volume–depth relationships for each aggregated lake unit. GLEV-based records provide observed lake surface area and evaporation data used to initialize lake states, define reference conditions, and calibrate and evaluate the model. Xanthos-Lake Outputs: Contains the basin- and lake-category-specific NetCDF outputs generated by Xanthos-Lake. Monthly variables include lake surface area, storage volume, outlet discharge, evaporation rate, evaporation volume, lake–groundwater exchange, lake inflow, ice thickness, snow depth, snow-cover fraction, ice-cover fraction, and mixing-layer temperature. The files also contain lake-specific calibration and validation statistics, including normalized root-mean-square error, mean absolute error, Nash–Sutcliffe efficiency, Kling–Gupta efficiency, and percent bias. Stored calibrated and derived parameters include the weir discharge coefficient, fractional freeboard, groundwater exchange coefficient, reference water level, corresponding reference surface area and storage volume, weir-width adjustment factor, and the fraction of routed inflow entering the lake. Basin identifiers, lake category, simulation period, calibration and validation periods, and parameter-schema information are retained as NetCDF metadata.

Abeshu, Guta [Pacific Northwest National Laborator↗

DOE OSTI2017

Dynamics of canopy stomatal conductance, transpiration, and evaporation in a temperate deciduous forest, validated by carbonyl sulfide uptake

Abstract. Stomatal conductance influences both photosynthesis and transpiration, thereby coupling the carbon and water cycles and affecting surface–atmosphere energy exchange. The environmental response of stomatal conductance has been measured mainly on the leaf scale, and theoretical canopy models are relied on to upscale stomatal conductance for application in terrestrial ecosystem models and climate prediction. Here we estimate stomatal conductance and associated transpiration in a temperate deciduous forest directly on the canopy scale via two independent approaches: (i) from heat and water vapor exchange and (ii) from carbonyl sulfide (OCS) uptake. We use the eddy covariance method to measure the net ecosystem–atmosphere exchange of OCS, and we use a flux-gradient approach to separate canopy OCS uptake from soil OCS uptake. We find that the seasonal and diurnal patterns of canopy stomatal conductance obtained by the two approaches agree (to within ±6 % diurnally), validating both methods. Canopy stomatal conductance increases linearly with above-canopy light intensity (in contrast to the leaf scale, where stomatal conductance shows declining marginal increases) and otherwise depends only on the diffuse light fraction, the canopy-average leaf-to-air water vapor gradient, and the total leaf area. Based on stomatal conductance, we partition evapotranspiration (ET) and find that evaporation increases from 0 to 40 % of ET as the growing season progresses, driven primarily by rising soil temperature and secondarily by rainfall. Counterintuitively, evaporation peaks at the time of year when the soil is dry and the air is moist. Our method of ET partitioning avoids concerns about mismatched scales or measurement types because both ET and transpiration are derived from eddy covariance data. Neither of the two ecosystem models tested predicts the observed dynamics of evaporation or transpiration, indicating that ET partitioning such as that provided here is needed to further model development and improve our understanding of carbon and water cycling.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2018

On Which Microphysical Time Scales to Use in Studies of Entrainment-Mixing Mechanisms in Clouds

The commonly used time scales in entrainment-mixing studies are examined in this paper to seek the most appropriate one, based on aircraft observations of cumulus clouds from the RACORO campaign and numerical simulations with the Explicit Mixing Parcel Model. The time scales include: τ evap , the time for droplet complete evaporation; τ phase , the time for saturation ratio deficit (S) to reach 1/e of its initial value; τ satu , the time for S to reach -0.5%; τ react , the time for complete droplet evaporation or S to reach -0.5%. It is found that the proper time scale to use depends on the specific objectives of entrainment-mixing studies. First, if the focus is on the variations of liquid water content (LWC) and S, then τ react for saturation, τ satu and τ phase are almost equivalently appropriate, because they all represent the rate of dry air reaching saturation or of LWC decrease. Second, if one focuses on the variations of droplet size and number concentration, τ react for complete evaporation and τ evap are proper because they characterize how fast droplets evaporate and whether number concentration decreases. Moreover, τ react for complete evaporation and τ evap are always positively correlated with homogeneous mixing degree (ψ), thus the two time scales, especially τ evap , are recommended for developing parameterizations. However, ψ and the other time scales can be negatively, positively, or not correlated, depending on the dominant factors of the entrained air (i.e., relative humidity or aerosols). Third and finally, all time scales are proportional to each other under certain microphysical and thermodynamic conditions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Electrical and Spectroscopic Diagnostics as Real‐Time Metallization Indicators During Hydrogen Plasma Smelting Reduction

This study investigates the hydrogen (H 2 ) plasma reduction process of direct-reduced-iron-grade hematite ore at different arc currents (100–200 A) in an Ar–5% H 2 atmosphere at 0.9 bar. Iron ore samples (10 g) were exposed to a plasma arc, and the reduction/metallization kinetics were analyzed over fixed time intervals. Electrical diagnostics revealed that the arc voltage exhibited takeover-mode oscillations which were suppressed at higher currents due to stronger electromagnetic coupling. The voltage dropped significantly as metallization approached ∼95%, linked to increased electrical conductivity of the metallic iron (Fe) in the ore as well as Fe evaporation into the arc, lowering the arc resistance. A simplified Elenbaas–Heller model supported this explanation and confirmed that Fe vapor concentration enhances plasma conductivity. Optical emission spectroscopy focused on the plasma–metal interface revealed the plasma's optically thick nature, as the primary Fe I 526.95 nm line experienced self-absorption. However, weaker Fe I lines (404.58, 438.35 nm) normalized to Ar I 696.5 nm provided a reliable proxy for metallization. These diagnostics, electrical and spectroscopic, effectively track metallization in real-time during H 2 plasma smelting reduction.

08 HYDROGEN↗

DOE OSTI2025

Resistivity distribution and donor properties of antimony-doped n -type Czochralski silicon ingots

We investigate antimony (Sb)-doped Czochralski-grown silicon as an alternative n-type substrate for photovoltaic applications, and characterize their axial resistivity distribution, donor properties, and mechanical strength. We find that Sb-doped ingots can achieve a more uniform resistivity distribution along the axial direction compared to P-doped counterparts. Dopant concentration profiles in P-doped ingots can be accurately modelled using the standard Scheil's equation, accounting only for dopant segregation during solidification. In contrast, modelling Sb-doped ingots requires consideration of both dopant segregation and evaporation effects to fit the dopant distribution accurately. Using electron paramagnetic resonance spectroscopy at 9 K, we observe two hyperfine lines in P-doped samples, and six hyperfine lines for Sb 121 and eight for Sb 123 isotopes, with the number of hyperfine lines governed by the nuclear spins. We further identify two-atom Sb clustering in the Sb-doped wafers, confirmed through simulations of the additional weak electron paramagnetic resonance peaks. Finally, we find that 140 μm as-cut planar Sb-doped wafers exhibit slightly higher mechanical strength compared to P-doped wafers.

14 SOLAR ENERGY↗

DOE OSTI2017

Depth distribution of secondary phases in kesterite Cu 2 ZnSnS 4 by angle-resolved X-ray absorption spectroscopy

The depth distribution of secondary phases in the solar cell absorber material Cu 2 ZnSnS 4 (CZTS) is quantitatively investigated using X-ray Absorption Near Edge Structure (XANES) analysis at the K-edge of sulfur at varying incidence angles. Varying information depths from several nanometers up to the full thickness is achieved. A quantitative profile of the phase distribution is obtained by a self-consistent fit of a multilayer model to the XANES spectra for different angles. Single step co-evaporated CZTS thin-films are found to exhibit zinc and copper sulfide secondary phases preferentially at the front or back interfaces of the film.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Moving mirror model for quasithermal radiation fields

We analyze the flow of energy and entropy emitted by a class of moving mirror trajectories which provide models for important aspects of the radiation fields produced by black hole evaporation. The mirror radiation fields provide natural, concrete examples of processes that follow thermal distributions for long periods, accompanied by transients which are brief and carry little net energy, yet ultimately represent pure quantum states. In conclusion, a burst of negative energy flux is a generic feature of these fields, but it need not be prominent.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

DEVAP-EDDR-TES (Simulation framework for a desiccant assisted air conditioning system with heat pump regeneration and energy storage) [SWR-24-66]

This software is a simulation framework that models a load flexible air conditioner system. The system consists of an evaporatively cooled liquid desiccant air conditioner (eLD-AC) subsystem, an electrically driven desiccant regenerator (EDDR) subsystem, and a stratified liquid desiccant storage (SLDS) subsystem. The software can be used to 1) predict the steady-state performance of the system given user-specified convergence criteria; 2) predict the dynamic performance of the entire system over a typical drive cycle operation subjected to user-specified building thermal loads and desired electrical load profile; 3) evaluate the synergy of all three subsystems operating altogether and improve the energy storage control strategy.

Huang, Ransisi↗

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↗