Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Wetness fraction”

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 37 records · Page 2

Predicting the Mineral Composition of Dust Aerosols: Model Evaluation and Identification of Key Processes with Observations - Part 2

A global compilation of nearly sixty measurement studies is used to evaluate two methods of simulating the mineral composition of dust aerosols in an Earth system model. Both methods are based upon a Mean Mineralogical Table (MMT) that relates the soil mineral fractions to a global atlas of arid soil type. The Soil Mineral Fraction (SMF) method assumes that the aerosol mineral fractions match the fractions of the soil. The MMT is based upon soil measurements after wet sieving, a process that destroys aggregates of soil particles that would have been emitted from the original, undisturbed soil. The second method approximately reconstructs the emitted aggregates. This model is referred to as the Aerosol Mineral Fraction (AMF) method because the mineral fractions of the aerosols differ from those of the wet-sieved parent soil, partly due to reaggregation. The AMF method remedies some of the deficiencies of the SMF method in comparison to observations. Only the AMF method exhibits phyllosilicate mass at silt sizes, where they are abundant according to observations. In addition, the AMF quartz fraction of silt particles is in better agreement with measured values, in contrast to the overestimated SMF fraction. Measurements at distinct clay and silt particle sizes are shown to be more useful for evaluation of the models, in contrast to the sum over all particles sizes that is susceptible to compensating errors, as illustrated by the SMF experiment. Model errors suggest that allocation of the emitted silt fraction of each mineral into the corresponding transported size categories is an important remaining source of uncertainty. Evaluation of both models and the MMT is hindered by the limited number of size-resolved measurements of mineral content that sparsely sample aerosols from the major dust sources. The importance of climate processes dependent upon aerosol mineral composition shows the need for global and routine mineral measurements.

aerosols

Investigation of Critical Heat Flux in Reduced Gravity Using Photomicrographic Techniques

Experiments were performed to examine the effects of body force on flow boiling critical heat flux (CHF). FC-72 was boiled along one wall of a transparent rectangular flow channel that permitted photographic study of the vapor-liquid interface just prior to CHF. High-speed video imaging techniques were used to identify dominant CHF mechanisms corresponding to different flow orientations and liquid velocities. Six different CHF regimes were identified: Wavy Vapor Layer, Pool Boiling, Stratification, Vapor Counterflow, Vapor Stagnation, and Separated Concurrent Vapor Flow. CHF showed significant sensitivity to orientation for flow velocities below 0.2 m/s, where extremely low CHF values where measured, especially with downward-facing heated wall and downflow orientations. High flow velocities dampened the effects of orientation considerably. The CHF data were used to assess the suitability of previous CHF models and correlations. It is shown the Interfacial Lift-off Model is very effective at predicting CHF for high velocities at all orientations. The flooding limit, on the other hand, is useful at estimating CHF at low velocities and for downflow orientations. A new method consisting of three dimensionless criteria is developed for determining the minimum flow velocity required to overcome body force effects on near-saturated flow boiling CHF. Vertical upflow boiling experiments were performed in pursuit of identifying the trigger mechanism for subcooled flow boiling CHF. While virtually all prior studies on flow boiling CHF concern the prediction or measurement of conditions that lead to CHF, this study was focused on events that take place during the CHF transient. High-speed video imaging and photomicrographic techniques were used to record the transient behavior of interfacial features from the last steady-state power level before CHF until the moment of power cut-off following CHF. The video records show the development of a wavy vapor layer which propagates along the heated wall, permitting cooling prior to CHF only in wetting fronts corresponding to the wave troughs. Image analysis software was developed to estimate void fraction from the individual video images. The void fraction records for subcooled flow boiling show the CHF transient is accompanied by gradual lift-off of wetting fronts culminating in some maximum vapor layer mean thickness, following which the vapor layer begins to thin down as the transition to film boiling ensues. This study proves the Interfacial Lift-off Model, which has been validated for near-saturated flow boiling CHF, is equally valid for subcooled conditions.

Mudawar, Issam

Modeling the effect of water on mantle rheology

To study the thermal history of the Earth we use a parameterized model of mantle convection. This model includes a mathematical description of de- and regassing processes of water from the Earth's mantle. The rates of this processes are considered to be directly proportional to the seafloor spreading rate. The kinematic viscosity of the mantle depends on the temperature/pressure as well as on the volatile content. Dissolved volatiles such as water weaken the minerals by reducing their activation energy for solid state creep. Karato and Toriumi showed a power law dependence between creep rate and water fugacity derived from experimental results. Therefore, we use such flow parameters of diffusion creep in olivine under wet and dry conditions to calculate the mantle viscosity as a function of the water content. Because the creep rate is proportional to the concentration of water-related point deflects we assume that the water fugacity is proportional to the water weight fraction. An equation for the steady-state strain rate under wet conditions is established. To assess the unknown constant K in this equation, we use flow law parameters given by Karato and Wu as well as the results of McGovern and Schubert.

Bounama, CH.

Predicting the Mineral Composition of Dust Aerosols: Representing Key Processes - Part 1

Soil dust aerosols created by wind erosion are typically assigned globally uniform physical and chemical properties within Earth system models, despite known regional variations in the mineral content of the parent soil. Mineral composition of the aerosol particles is important to their interaction with climate, including shortwave absorption and radiative forcing, nucleation of cloud droplets and ice crystals, heterogeneous formation of sulfates and nitrates, and atmospheric processing of iron into bioavailable forms that increase the productivity of marine phytoplankton. Here, aerosol mineral composition is derived by extending a method that provides the composition of a wet-sieved soil. The extension accounts for measurements showing significant differences between the mineral fractions of the wetsieved soil and the emitted aerosol concentration. For example, some phyllosilicate aerosols are more prevalent at silt sizes, even though they are nearly absent at these diameters in a soil whose aggregates are dispersed by wet sieving. We calculate the emitted mass of each mineral with respect to size by accounting for the disintegration of soil aggregates during wet sieving. These aggregates are emitted during mobilization and fragmentation of the original undispersed soil that is subject to wind erosion. The emitted aggregates are carried far downwind from their parent soil. The soil mineral fractions used to calculate the aggregates also include larger particles that are suspended only in the vicinity of the source. We calculate the emitted size distribution of these particles using a normalized distribution derived from aerosol measurements. In addition, a method is proposed for mixing minerals with small impurities composed of iron oxides. These mixtures are important for transporting iron far from the dust source, because pure iron oxides are more dense and vulnerable to gravitational removal than most minerals comprising dust aerosols. A limited comparison to measurements from North Africa shows that the model extensions result in better agreement, consistent with a more extensive comparison to global observations as well as measurements of elemental composition downwind of the Sahara, as described in companion articles.

aerosols

Climatic imprint on interfacially-controlled platinum-palladium resources

Abstract Iron oxide-rich laterites, soils, and regolith formed from the weathering of ultramafic rocks represent untapped unconventional resources for the critical minerals platinum and palladium, but the fundamental surficial geochemistry of these elements remains poorly understood. Depletion of Pd relative to Pt occurs in some weathering zones in semi-arid climates. The accepted model attributes this platinum-palladium chemical fractionation to preferential complexation of Pd by dissolved chloride. However, similar fractionation is not observed in laterites of humid equatorial regions despite substantial wet deposition of chloride. The established mechanistic model for Pt and Pd behavior during weathering thus inaccurately predicts the distribution of these critical minerals in many settings, hindering global resource assessment. We show through mineral-fluid partitioning experiments coupled to element-specific spectroscopy that this canonical explanation for platinum-palladium fractionation is invalid: chloride complexation does not differentially mobilize Pd versus Pt. Instead, mineral-specific interfacial reactions control Pd and Pt accumulation. Modeling of platinum-palladium fractionation in representative weathering zone profiles demonstrates sub-equal retention in goethite-rich settings and Pd depletion in hematite-rich zones, accurately predicting trends observed in soils and laterites. Iron oxide mineralogy, reflecting modern and past regional climate conditions, is likely the primary determinant of Pt and Pd endowment in weathering zone resources. This new model for Pt and Pd mobilization and accumulation behavior provides a mechanistic foundation for exploration and recovery of platinum group elements from novel ultramafic regolith deposits.

58 GEOSCIENCES

Importance of Tree- and Species-Level Interactions With Wildfire, Climate, and Soils in Interior Alaska: Implications for Forest Change Under A Warming Climate

The boreal zone of Alaska is dominated by interactions between disturbances, vegetation, and soils. These interactions are likely to change in the future through increasing permafrost thaw, more frequent and intense wildfires, and vegetation change from drought and competition. We utilize an individual tree-based vegetation model, the University of Virginia Forest Model Enhanced (UVAFME), to estimate current and future forest conditions across sites within interior Alaska. We updated UVAFME for application within interior Alaska, including improved simulation of permafrost dynamics, litter decay, nutrient dynamics, fire mortality, and post-fire regrowth. Following these updates, UVAFME output on species-specific biomass and stem density was comparable to inventory measurements at various forest types within interior Alaska. We then simulated forest response to climate change at specific inventory locations and across the Tanana Valley River Basin on a 2 × 2 km2 grid. We derived projected temperature and precipitation from a five-model average taken from the CMIP5 archive under the RCP 4.5 and 8.5 scenarios. Results suggest that climate change and the concomitant impacts on wildfire and permafrost dynamics will result in overall decreases in biomass (particularly for spruce (Picea spp.)) within the interior Tanana Valley, despite increases in quaking aspen (Populus tremuloides) biomass, and a resulting shift towards higher deciduous fraction. Simulation results also predict increases in biomass at cold, wet locations and at high elevations, and decreases in biomass in dry locations, under both moderate (RCP 4.5) and extreme (RCP 8.5) climate change scenarios. These simulations demonstrate that a highly detailed, species interactive model can be used across a large region within Alaska to investigate interactions between vegetation, climate, wildfire, and permafrost. The vegetation changes predicted here have the capacity to feed back to broader scale climate-forest interactions in the North American boreal forest, a region which contributes significantly to the global carbon and energy budgets.

forest modeling

Modeling Dust Mineralogical Composition: Sensitivity to Soil Mineralogy Atlases and Their Expected Climate Impacts

Soil dust aerosols are a key component of the climate system, as they interact with short- and long-wave radiation, alter cloud formation processes, affect atmospheric chemistry and play a role in biogeochemical cycles by providing nutrient inputs such as iron and phosphorus. The influence of dust on these processes depends on its physicochemical properties, which, far from being homogeneous, are shaped by its regionally varying mineral composition. The relative amount of minerals in dust depends on the source region and shows a large geographical variability. However, many state-of-the-art Earth system models (ESMs), upon which climate analyses and projections rely, still consider dust mineralogy to be invariant. The explicit representation of minerals in ESMs is more hindered by our limited knowledge of the global soil composition along with the resulting size-resolved airborne mineralogy than by computational constraints. In this work we introduce an explicit mineralogy representation within the state-of-the-art Multiscale Online Nonhydrostatic AtmospheRe CHemistry (MONARCH) model. We review and compare two existing soil mineralogy datasets, which remain a source of uncertainty for dust mineralogy modeling and provide an evaluation of multiannual simulations against available mineralogy observations. Soil mineralogy datasets are based on measurements performed after wet sieving, which breaks the aggregates found in the parent soil. Our model predicts the emitted particle size distribution (PSD) in terms of its constituent minerals based on brittle fragmentation theory (BFT), which reconstructs the emitted mineral aggregates destroyed by wet sieving. Our simulations broadly reproduce the most abundant mineral fractions independently of the soil composition data used. Feldspars and calcite are highly sensitive to the soil mineralogy map, mainly due to the different assumptions made in each soil dataset to extrapolate a handful of soil measurements to arid and semi-arid regions worldwide. For the least abundant or more difficult-to-determine minerals, such as iron oxides, uncertainties in soil mineralogy yield differences in annual mean aerosol mass fractions of up to ∼ 100 %. Although BFT restores coarse aggregates including phyllosilicates that usually break during soil analysis, we still identify an overestimation of coarse quartz mass fractions (above 2 µm in diameter). In a dedicated experiment, we estimate the fraction of dust with undetermined composition as given by a soil map, which makes up ∼ 10 % of the emitted dust mass at the global scale and can be regionally larger. Changes in the underlying soil mineralogy impact our estimates of climate-relevant variables, particularly affecting the regional variability of the single-scattering albedo at solar wavelengths or the total iron deposited over oceans. All in all, this assessment represents a baseline for future model experiments including new mineralogical maps constrained by high-quality spaceborne hyperspectral measurements, such as those arising from the NASA Earth Surface Mineral Dust Source Investigation (EMIT) mission.

Soil dust

A cloudiness transition in a marine boundary layer

In situ aircraft data and lidar data are used to analyze a transition in the boundary layer thermodynamic structure from a clear boundary layer through small cumulus and broken stratocumulus to a deck of solid stratocumulus. The data was collected in conjunction with a Landsat overpass on July 7, 1987 off the coast of southern California. A steady progression in mixing line stability is seen associated with the change in cloudiness. The (empirically based) stability threshold for the breakup of this stratocumulus is that the slope of the mixing line is 0.66 + or - 0.04 of the slope of the wet virtual adiabat. A simple linear parameterization is proposed for cloud fraction in terms of mixing line stability. Surface flux measurements are consistent with bulk aerodynamic estimates.

Boers, Reinout

Characteristics of Pool Boiling on Graphite-Copper Composite Surfaces

Nucleate pool boiling performance of different liquids on graphite-copper composite (Gr-Cu) surfaces has been experimentally studied and modeled. Both highly wetting fluids, such as freon-113 and pentane, and a moderately wetting fluid (water) were tested on the Gr-Cu surfaces with different graphite-fiber volume fractions to reveal the enhancement effects of the composite surfaces on the nucleate pool boiling. Results of the experiments show that the graphite-fiber volume fraction has an optimum value. The Gr-Cu composite surface with 25 percent graphite-fiber volume (f=0.25) has a maximum enhancement effect on the nucleate boiling heat transfer comparing to the pure copper surface. For the highly wetting fluid, the nucleate boiling heat transfer is generally enhanced on the Gr- Cu composite surfaces by 3 to 6 times shown. In the low heat flux region, the enhancement is over 6 times, but in the high heat flux region, the enhancement is reduced to about 40%. For the moderately wetting fluid (water), stronger enhancement of nucleate boiling heat transfer is achieved on the composite surface. It shown the experimental results in which one observes the nucleate boiling heat transfer enhancement of 5 to 10 times in the low heat flux region and an enhancement of 3 to 5 times in the high heat flux region. Photographs of bubble departure during the initial stage of nucleate boiling indicate that the bubbles detached from the composite surface are much smaller in diameter than those detached from the pure copper surface. Typical photographs are presented.It shows that the bubbles departed from the composite surface have diameters of only O(0.1) mm, while those departed from the pure copper surface have diameters of O(1) mm. It is also found that the bubbles depart from the composite surface at a much higher frequency, thus forming vapor columns. These two phenomena combined with high thermal conductivity of the graphite fiber are considered the mechanisms for such a significant augmentation in nucleate boiling heat transfer on the composite surfaces. A physical model is developed to describe the phenomenon of bubble departure from the composite surface: The preferred site of bubble nucleation is the fiber tip because of higher tip temperature than the surrounding copper base and poor wettability of the graphite tip compared with that of the base material (copper). The high evaporation rate near the contact line produces the vapor cutback due to the vapor recoil pushing the three-phase line outwards from the fiber tip, and so a neck of the bubble is formed near the bubble bottom. Evaporation and surface tension accelerate the necking process and finally result in the bubble departure while a new small bubble is formed at the tip when the surface tension pushes the three-phase line back to the tip. The process is schematically shown. The proposed model is based on and confirmed by experimental results.

Zhang, Nengli

Reconciling LCROSS and Orbital Neutron Water Abundance Estimates in Cabeus Crater

The Lunar Prospector Neutron Spectrometer (LPNS) first revealed Cabeus crater (84.9 deg S, 35.5degW) as having the highest inferred hydrogen on the Moon. Because of the broad LPNS footprint (approximately 40 km FWHM), the apparent peak water-equivalent hydrogen (WEH) concentration is only approximately 0.25 wt%, but could be much higher in smaller areas than the spectrometer footprint. Earlier image reconstruction work suggested that areas within permanent shadow have abundances approximately 1 wt% WEH. However, the LCROSS impact yielded total water estimates, ice plus vapor, of between 3 and 10 wt%. The large disagreement between LCROSS and apparent orbital values imply that either the ice is buried, by perhaps as much as 50 to 100 cm; or the ice distribution within Cabeus is spatially inhomogeneous, or both. Modeling reveals that the areal extent of a "shallow permafrost zone" is far greater than the area of permanent shadow. Ice can be virtually stable for billions of years within a few tens of centimeters of the surface in these areas. However, the LCROSS impact took place in an area of permanent shadow. If stably-trapped volatiles can be found in locales that receive occasional, oblique sunlight, landed missions may target these sites and eventual resource exploitation may be done more easily. Are orbital neutron data consistent with areally-extensive, volatile-rich cold traps? Orbital epithermal neutron data over the northern half of Cabeus (near the LCROSS impact site) are consistent with 0.2 wt% WEH or less in the "permafrost zone" near the crater. On the other hand, pixon reconstructions that confine the hydrogen enhancements to permanent shadow result in higher abundance estimates -- around 1 wt% if homogeneously mixed. But if the PSR abundance is increased to 10 wt%, consistent with the sum of all H-bearing compounds seen by LCROSS, a much larger-than-observed reduction in neutron count rate would be seen from orbit. It is likely that volatiles are inhomogeneously distributed, due to both impact processes and emplacement history. Two possibilities may bring consistency to the orbital and LCROSS measurements. Inhomogeneous lateral distribution: Consider the extreme case of a bimodal distribution within the crater -- dry and wet. In this case the epithermal leakage flux seen from orbit is a mixture of two different values, weighted according to fractional areas. Two possible outcomes, depending on whether the inferred leakage flux for the PSR or "permafrost" areas are considered. In the first case, approximately 40% of the PSR may be "wet", the remainder dry (and LCROSS was slightly lucky). However, if the whole area of permafrost is considered, then as little as 20% of the area will be as "wet" as the LCROSS results (and LCROSS was quite lucky). Inhomogeneous depth distribution: The leakage flux of thermal and epithermal neutrons depends on depth of burial of an icy layer beneath dry ferroan anorthosite soil (FAn). For the Cabeus PSR, the pixon reconstruction values for the epithermal flux allows a range of abundance and burial depth, while that of the thermal+epi detector constrains this range. (Uncertainties in iron abundance in the FAn can have significant impact on thermal neutron leakage flux estimates.) Between 20% and 40% of the Cabeus floor may be "wet", or alternatively a 5-10 wt% "wet" layer exists between 50 and 100 cm beneath a layer of dry regolith within the PSR. But volatile abundances of 5 wt% or more, distributed uniformly and homogeneously throughout the Cabeus PSR do not agree with orbital measurements

Elphic, Richard

Free-spray characteristics and spray-wall interactions of methanol on a gasoline direct injector under flash-boiling and non-flash-boiling conditions

Methanol is considered a promising alternative fuel for internal combustion engines (ICEs) due to its high-octane number, fast laminar flame speed, and elevated latent heat of vaporization, all of which support higher compression ratios and improved thermal efficiency. However, its substantial latent heat of vaporization also poses cold-start challenges, such as misfire and fuel film deposition. This study aims to investigate methanol spray morphology and spray-wall interaction using the Spray M injector from the Engine Combustion Network within a constant-pressure flow vessel. A recently developed unified numerical framework capable of modeling both flash and non-flash boiling sprays is validated against experimental liquid volume fraction data acquired via 3-D computed tomography. Here, the results reveal that flash boiling significantly alters the spray morphology, leading to smaller droplets and spray collapse due to enhanced air-entrainment-induced turbulence. Quantitative agreement between experiments and simulations confirms this behavior. Coupled 0-D equilibrium and 3-D computational fluid dynamics analyses show that flash boiling accelerates evaporation and reduces fuel residence time, while non-flash conditions maintain a persistent liquid core more susceptible to wall wetting. Wall temperature diagnostics reveal that spray collapse alters heat transfer patterns by shifting cooling effects. Mixture fraction analysis indicates that evaporation is primarily governed by shear-layer turbulence, though deviations from adiabatic equilibrium mixing emerge under low-turbulence conditions. Finally, increasing fuel, ambient, and wall temperatures reduces wall wetting and film thickness, mitigating cold-start risks. These findings enhance the understanding of methanol sprays’ behavior and support its adoption as a viable, alternative fuel for ICEs.

Engine Combustion Network

IM3 Projected US Data Center Locations

IM3 Projected US Data Center Locations This dataset contains model projections of new data center facilities in the contiguous United States (CONUS) through 2035 using the CERF – Data Centers model. Data center locations are modeled across four data center electricity demand growth scenarios (low, moderate, high, higher) and five market gravity scenarios (0%, 25%, 50%, 75%, 100%). Projected locations are intended to be regional representations of feasible siting locations in the future to assess potential grid and water stress impacts. The data center load growth scenarios correspond with the rates outlined in EPRI (2024) and include 3.71%, 5%, 10%, and 15% annual growth of electricity demand for data centers from 2023 values in 37 states across the CONUS. Market gravity scenarios correspond to the relative importance of proximity to data center markets or high population areas compared to locational cost in the siting algorithm. 0% market gravity means that siting decisions were entirely determined by the locational cost in each feasible location. 100% market gravity means that only market proximity was considered when siting. Other scenarios have weight placed on both components where total weight always equals 100%. Locational cost is dependent on facility cooling type and corresponding electricity cost, taxes, and other factors. Facility cooling type is spatially determined where high water stress and/or areas with high summer wet bulb temperatures are assumed to operate with mechanical cooling for a higher fraction of the year rather than evaporative cooling. Feasible data center siting areas are based on geospatial suitability raster data developed with open-source information. The following areas are excluded from siting: Areas within 300 m of a federal airport runway Waterbodies Areas with slope >16% Areas susceptible to sinkholes High coastal or inland flood risk areas Local, state, and federal parks, leisure areas, and cemeteries Areas >2 km away from electric substations Areas >5 km away from a municipal water supplier service area Areas >2 km away from high-speed fiber provider service territory Protected Areas Database of the United States (PAD-US) areas Railroads, major roadways, and minor roadways Military areas and training grounds NLCD developed lands Areas >0.8 km (0.5 miles) from NLCD developed lands Because we use open-source information, proprietary information that can influence siting decisions such as individual tax agreements with cities, detailed fiber line connectivity, electric grid power capacity agreements, and others, are not currently accounted for in the modeling process. Using specific building locations and footprints in the dataset for local planning purposes is not advised. Technical Information Geospatial data is provided in geojson format using the Albers Equal Area Conic (ESRI:102003) coordinate reference system. The datasets contain the following parameters: id - unique identification number within given scenario file growth_scenario – data center demand growth scenario market_gravity_weight – market gravity weight scenario (%) region – name of region (i.e., US State) total_cost_million_usd – locational siting cost ($million) campus_size_square_ft – total land acquired for data center facility (square ft) data_center_it_power_mw – IT power of data center facility (MW) mechanical_cooling_frac – fraction of year when data center uses mechanical cooling system water_cooling_frac– fraction of year when data center uses evaporative cooling system cooling_energy_demand_mwh – total annual facility energy demand for cooling (MWh) cooling_water_demand_mgy – total annual facility water demand for cooling (MG) cooling_water_consumption_mgy – total annual facility water consumed (MG) normalized_locational_cost – normalized total locational cost score for location normalized_gravity_score – normalized market gravity score for location weighted_siting_score – total weighted siting score of locational cost and gravity score geometry – polygon geometry of facility Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program. License This data is made available under a CCBY4.0 License Disclaimer This material was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor the United States Department of Energy, nor the Contractor, nor any or their employees, nor any jurisdiction or organization that has cooperated in the development of these materials, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness or any information, apparatus, product, software, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or Battelle Memorial Institute. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. PACIFIC NORTHWEST NATIONAL LABORATORYoperated byBATTELLEfor theUNITED STATES DEPARTMENT OF ENERGYunder Contract DE-AC05-76RL01830

Mongird, Kendall (ORCID:0000000328077088)

IM3 Projected US Data Center Locations

IM3 Projected US Data Center Locations This dataset contains model projections of new data center facilities in the contiguous United States (CONUS) through 2035 using the CERF – Data Centers model. Data center locations are modeled across four data center electricity demand growth scenarios (low, moderate, high, higher) and five market gravity scenarios (0%, 25%, 50%, 75%, 100%). Projected locations are intended to be regional representations of feasible siting locations in the future to assess potential grid and water stress impacts. The data center load growth scenarios correspond with the rates outlined in EPRI (2024) and include 3.71%, 5%, 10%, and 15% annual growth of electricity demand for data centers from 2023 values in 37 states across the CONUS. Market gravity scenarios correspond to the relative importance of proximity to data center markets or high population areas compared to locational cost in the siting algorithm. 0% market gravity means that siting decisions were entirely determined by the locational cost in each feasible location. 100% market gravity means that only market proximity was considered when siting. Other scenarios have weight placed on both components where total weight always equals 100%. Locational cost is dependent on facility cooling type and corresponding electricity cost, taxes, and other factors. Facility cooling type is spatially determined where high water stress and/or areas with high summer wet bulb temperatures are assumed to operate with mechanical cooling for a higher fraction of the year rather than evaporative cooling. Feasible data center siting areas are based on geospatial suitability raster data developed with open-source information. The following areas are excluded from siting: Areas within 300 m of a federal airport runway or within an airport area boundary Waterbodies Areas with slope >16% Areas susceptible to sinkholes High coastal or inland flood risk areas Local, state, and federal parks, leisure areas, and cemeteries Areas >2 km away from electric substations Areas >5 km away from a municipal water supplier service area Areas >2 km away from high-speed fiber provider service territory USGS Protected Areas Database of the United States (PAD-US) GAP status 1, 2, or 3 areas US National Parks Wetlands USFWS critical habitats BIA land areas Railroads, major roadways, and minor roadways Military areas and training grounds NLCD developed lands Areas >0.8 km (0.5 miles) from NLCD developed lands Because we use open-source information, proprietary information that can influence siting decisions such as individual tax agreements with cities, detailed fiber line connectivity, electric grid power capacity agreements, and others, are not currently accounted for in the modeling process. Using specific building locations and footprints in the dataset for local planning purposes is not advised. Technical Information Geospatial data is provided in geojson format using the Albers Equal Area Conic (ESRI:102003) coordinate reference system. The datasets contain the following parameters: id - unique identification number within given scenario file growth_scenario – data center demand growth scenario market_gravity_weight – market gravity weight scenario (%) region – name of region (i.e., US State) total_cost_million_usd – locational siting cost ($million) campus_size_square_ft – total land acquired for data center facility (square ft) data_center_it_power_mw – IT power of data center facility (MW) mechanical_cooling_frac – fraction of year when data center uses mechanical cooling system water_cooling_frac– fraction of year when data center uses evaporative cooling system cooling_energy_demand_mwh – total annual facility energy demand for cooling (MWh) cooling_water_demand_mgy – total annual facility water demand for cooling (MG) cooling_water_consumption_mgy – total annual facility water consumed (MG) normalized_locational_cost – normalized total locational cost score for location normalized_gravity_score – normalized market gravity score for location weighted_siting_score – total weighted siting score of locational cost and gravity score geometry – polygon geometry of facility Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program. License This data is made available under a CCBY4.0 License Disclaimer This material was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor the United States Department of Energy, nor the Contractor, nor any or their employees, nor any jurisdiction or organization that has cooperated in the development of these materials, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness or any information, apparatus, product, software, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or Battelle Memorial Institute. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. PACIFIC NORTHWEST NATIONAL LABORATORYoperated byBATTELLEfor theUNITED STATES DEPARTMENT OF ENERGYunder Contract DE-AC05-76RL01830

Mongird, Kendall (ORCID:0000000328077088)

Close-up analysis of aircraft ice accretion

Various types of ice formation have been studied by analysis of high magnification video observations. All testing was conducted in the NASA Lewis Icing Research Tunnel (IRT). A faired 8.9 cm (3.5 in.) diameter metal-clad cylinder and a 5.1 (2 in.) aluminum cylinder were observed by close-up and overview video cameras for several wind tunnel conditions. These included close-up grazing angle, close-up side view, as well as overhead and side overview cameras. Still photographs were taken at the end of each spray along with tracings of the subsequent ice shape. While in earlier tests only the stagnation region was observed, the entire area from the stagnation line to the horn region of glaze ice shapes was observed in this test. The modes or horn formation have been identified within the range of conditions observed. In the horn region, Horn Type A ice is formed by 'dry' feather growth into the flow direction and Horn Type B is formed by a 'wet' growth normal to the surface. The feather growth occurs when the freezing fraction is near unity and roughness elements exist to provide an initial growth site.

Hansman, R. John

Close-up analysis of aircraft ice accretion

Various types of ice formation have been studied by analysis of high magnification video observations. All testing was conducted in the NASA Lewis Icing Research Tunnel (IRT). A faired 8.9 cm (3.5 in.) diameter metal-clad cylinder and a 5.1 (2 in.) aluminum cylinder were observed by close-up and overview video cameras for several wind tunnel conditions. These included close-up grazing angle, close-up side view, as well as overhead and side overview cameras. Still photographs were taken at the end of each spray along with tracings of the subsequent ice shape. While in earlier tests only the stagnation region was observed, the entire area from the stagnation line to the horn region of glaze ice shapes was observed in this test. The modes or horn formation have been identified within the range of conditions observed. In the horn region, Horn Type A ice is formed by 'dry' feather growth into the flow direction and Horn Type B is formed by a 'wet' growth normal to the surface. The feather growth occurs when the freezing fraction is near unity and roughness elements exists to provide an initial growth site.

Hansman, R. J., Jr.

VARTM Model Development and Verification

In this investigation, a comprehensive Vacuum Assisted Resin Transfer Molding (VARTM) process simulation model was developed and verified. The model incorporates resin flow through the preform, compaction and relaxation of the preform, and viscosity and cure kinetics of the resin. The computer model can be used to analyze the resin flow details, track the thickness change of the preform, predict the total infiltration time and final fiber volume fraction of the parts, and determine whether the resin could completely infiltrate and uniformly wet out the preform.

Cano, Roberto J.

El Nino-Southern Oscillation Correlated Aerosol Angstrom Exponent Anomaly Over the Tropical Pacific Discovered in Satellite Measurements

El Nino.Southern Oscillation (ENSO) is the dominant mode of interannual variability in the tropical atmosphere. ENSO could potentially impact local and global aerosol properties through atmospheric circulation anomalies and teleconnections. By analyzing aerosol properties, including aerosol optical depth (AOD) and Angstrom exponent (AE; often used as a qualitative indicator of aerosol particle size) from the Moderate Resolution Imaging Spectrometer, the Multiangle Imaging Spectroradiometer and the Sea ]viewing Wide Field ]of ]view Sensor for the period 2000.2011, we find a strong correlation between the AE data and the multivariate ENSO index (MEI) over the tropical Pacific. Over the western tropical Pacific (WTP), AE increases during El Nino events and decreases during La Nina events, while the opposite is true over the eastern tropical Pacific (ETP). The difference between AE anomalies in the WTP and ETP has a higher correlation coefficient (>0.7) with the MEI than the individual time series and could be considered another type of ENSO index. As no significant ENSO correlation is found in AOD over the same region, the change in AE (and hence aerosol size) is likely to be associated with aerosol composition changes due to anomalous meteorological conditions induced by the ENSO. Several physical parameters or mechanisms that might be responsible for the correlation are discussed. Preliminary analysis indicates surface wind anomaly might be the major contributor, as it reduces sea ]salt production and aerosol transport during El Nino events. Precipitation and cloud fraction are also found to be correlated with tropical Pacific AE. Possible mechanisms, including wet removal and cloud shielding effects, are considered. Variations in relative humidity, tropospheric ozone concentration, and ocean color during El Nino have been ruled out. Further investigation is needed to fully understand this AE ]ENSO covariability and the underlying physical processes responsible for it.

Li, Jing

A Robotic High-Throughput Grid-Search Platform for Mapping Phase Behavior in Triblock Copolymer–Homopolymer Blends

We present a high-throughput experimental investigation of the phase behavior in triblock copolymers (PS-b-PB-b-PS and PS-b-PI-b-PS) and polystyrene (PS) homopolymer blends as a function of homopolymer molecular weight (MW) and blend ratio. Using a robotic thin-film processing platform (NOVA) integrated with Grazing Incidence Small-Angle X-ray Scattering (GISAXS) and Atomic Force Microscopy (AFM), we systematically mapped the order–disorder transition (ODT) boundaries and domain spacing evolution across a broad MW range (4.0–101.3 kDa) with varying homopolymer loadings (10% to 90%). The results reveal three distinct regimes: low-MW homopolymers, corresponding to the wet-brush regime produced only gradual domain swelling before disordering at high blend ratios (weight fraction); medium-MW homopolymers, corresponding to thedry-brush regime induced significant domain spacing increase up to 80% followed by earlier disordering, while high-MW homopolymers led to macrophase separation with minimal changes in domain spacing. Additionally, coarse-grained molecular dynamics simulations confirmed our experimental finding that in the low-MW region, the PS homopolymer uniformly distributed in the PS domain. These findings demonstrate that homopolymer molecular weight critically governs both the extent of domain swelling and the onset of disorder in triblock copolymer systems. This high-throughput platform enables the rapid mapping of composition–morphology relationships and can be integrated with AI/ML tools for designing next-generation nanostructured polymers.

36 MATERIALS SCIENCE