Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation 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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Subsetted model output for surface-related variables from a LASSO-CACTI WRF simulation

Surface-related model output from a Weather Research and Forecasting (WRF) simulation from the Large-Eddy Simulation (LES) Atmospheric Radiation Measurement (ARM) Symbiotic Simulation and Observation (LASSO) deep-convection scenario for the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign. The LASSO-CACTI simulations span grid spacings from 7.5 km to 100 m through the use of four nested domains, labeled d1 through d4. The simulations are of convection near the Sierras de Córdoba mountain range, roughly centered on the ARM Mobile Facility. More information can be found at https://www.arm.gov/capabilities/modeling/lasso. This version of the output is a collection of variables subsetted from the raw WRF output to reduce the file overhead for users not needing the full raw dataset. The subset files are in netCDF format with the variables in these files at the surface and sub-surface levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Three-Dimensional Surface Downwelling Longwave Radiation Clear-Sky Effects in the Upper Colorado River Basin

In complex terrain, non-parallel surfaces receive emitted radiation from adjacent surfaces. Qualitatively, where surface skin temperatures and lower tropospheric temperature and humidity are not uniform, the downwelling longwave radiation (DLR) will be determined not just by radiation from the atmosphere above a given location, but also by adjacent surface temperatures. We quantify this three-dimensional longwave radiative effect over the Upper Colorado River Basin in clear-sky conditions by calculating surface DLR with observed land-surface temperatures from ECOSTRESS. We find that this effect is due to terrain-subtended sky-view and represents ~22% of the surface longwave flux, rising to ~28% and ~24% in the East and Southeast of the Basin, respectively, and can be >50% in extreme cases. The common omission of this effect in atmospheric radiation models leads to an underestimation of DLR in complex terrain, especially at higher elevations, which has significant implications for mountainous ecohydrology simulations.

3-dimensional effects↗
DOE OSTI2022

Tritium Transport Modeling at the Pore Scale in Ceramic Breeder Materials Using TMAP8

Fusion reactors depend on the blanket material to breed and release tritium at the same rate or faster than it is consumed by the fusion reaction. Cellular ceramic breeders (CCBs) are dense materials that can maintain a high tritium breeding ratio while promoting tritium release because of highly connected pores. Assessing the tritium breeding capabilities of these materials requires a combination of extensive experimental and modeling efforts. In this work, we develop and calibrate a multiphysics model of tritium transport. This novel model accounts for ceramic and pore diffusion, trapping and detrapping, and several surface reactions at the pore surface. We perform a sensitivity analysis and calibrate the model by comparing its predictions against experimental measurements of deuterium absorption. The calibrated model is then used to model tritium absorption in samples with different pore microstructures to investigate the effect of pore interconnectivity on tritium absorption. The model is part of the development of the multiscale, multiphysics framework for tritium transport [i.e., the Tritium Migration Analysis Program (TMAP8)], which is itself built on top of the finite-element multiphysics framework multiphysics object-oriented simulation environment (MOOSE). This study demonstrates some of TMAP8’s capabilities and is the first step toward assessing the tritium breeding capabilities of ceramic breeder material designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A coarse-grained model of clay colloidal aggregation and consolidation with explicit representation of the electrical double layer

The aggregation of clay minerals in liquid water exemplifies colloidal self-assembly in nature. These negatively charged aluminosilicate platelets interact through multiple mechanisms with different sensitivities to particle shape, surface charge, aqueous chemistry, and interparticle distance and exhibit complex aggregation structures. Experiments have difficulty resolving the associated colloidal assemblages at the scale of individual particles. Conversely, all-atom molecular dynamics (MD) simulations provide detailed insight on clay colloidal interaction mechanisms, but they are limited to systems containing a few particles. We develop a new coarse-grained (CG) model capable of representing assemblages of hundreds of clay particles with accuracy approaching that of MD simulations, at a fraction of the computational cost. Our CG model is parameterized based on MD simulations of a pair of smectite clay particles in liquid water. A distinctive feature of our model is that it explicitly represents the electrical double layer (EDL), i.e., the cloud of charge-compensating cations that surrounds the clay particles. Our model captures the simultaneous importance of long-range colloidal interactions (i.e., interactions consistent with simplified analytical models, already included in extant clay CG models) and short-range interactions such as ion correlation and surface and ion hydration effects. The resulting simulations correctly predict, at low solid-water ratios, the existence of ordered arrangements of parallel particles separated by water films with a thickness up to ~10 nm and, at high solid-water ratios, the coexistence of crystalline and osmotic swelling states, in agreement with experimental observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Surface Quantitative Precipitation Estimates (SQUIRE) of Snow Water Equivalent from the Surface Atmospheric Integrated Field Laboratory

The upper Colorado River basin is the primary source of water for 40 million people. With declining snowpack in the basin, forecasting hydrological budgets in the Southwest United States is more important than ever. However, due in part, to a lack of reliable observations of precipitation in complex terrain, hydrological models struggle to assess and forecast snowpack snow water equivalent (SWE) in the upper Colorado River basin (UCRB). Therefore, the need for more reliable SWE forecasts in the UCRB motivated the U.S. Department of Energy Atmospheric Radiation Measurement Facility’s Surface Atmospheric Integrated Field Laboratory (SAIL) that occurred from June 2021 to June 2023. During SAIL, the X-band precipitation radar from Colorado State University conducted volume scans sampling the precipitation properties over the UCRB. The ARM facility developed a gridded Surface Quantitative Precipitation Estimates (SQUIRE) product from the radar observations. To do this, various daily SWE estimates from radar using the radar reflectivity factor Z e and specific differential phase K dp were compared against ground-based precipitation gauges. SWE in precipitation calculated from Wolfe and Snider’s S–Z e estimator was in best agreement with the rain gauges for the days when SWE < 12 mm. For days with SWE > 12 mm, the WSR-88D Intermountain West relationship had the best agreement with the precipitation gauges. Airborne snow depth observations show that SQUIRE captures regions of orographic enhancement in the mountains to the west and northwest of the SAIL study area, indicating that the scientific community should focus on understanding and ultimately simulating orographic atmospheric precipitation processes to improve UCRB snowpack SWE assessment and forecasting.

Hydrology↗
DOE OSTI2025

Complex surface topographic changes from explosive experiments at the Dry Alluvium Geology site, southern Nevada, United States

Understanding the surface topographic change that results from underground explosions is important for global security. Current techniques to relate the surface change to underground explosion characteristics usually involve assuming the earth has homogenous properties, leading to highly variable interpretations. Here we use an unoccupied aerial platform and a digital single lens reflex camera along with 200+ ground control points surveyed with a real-time kinematic global navigation satellite system to measure the surface topographic change resulting from two underground explosions at the Dry Alluvium Geology site in Yucca Flat, Nevada National Security Site, southern Nevada, United States. We find areas of 5–7 cm of subsidence that are not directly above the explosion source but rather 200–300 m away. For experiment DAG2, this zone is located south and west of the explosion, while for DAG4, there is a zone of subsidence located northeast of the explosion. In addition, late-time measurements show as much as 5 cm of horizontal change without measurable associated vertical change in the weeks following DAG4 but not DAG2. These indicate that the deformation resulting from underground chemical explosions can be very complex and bear little to no resemblance to predictions using half-space models. It is likely the tectonic environment plays a significant role in controlling the surface change, but the details are not fully understood.

58 GEOSCIENCES↗
DOE OSTI2025

JAX-CanVeg: A Differentiable Land Surface Model

Land surface models consider the exchange of water, energy, and carbon along the soil-canopy-atmosphere continuum, which is challenging to model due to their complex interdependency and associated challenges in representing and parameterizing them. Differentiable modeling provides a new opportunity to capture these complex interactions by seamlessly hybridizing process-based models with deep neural networks (DNNs), benefiting both worlds, that is, the physical interpretation of process-based models and the learning power of DNNs. Here, we developed a differentiable land model, JAX-CanVeg. The new model builds on the legacy CanVeg by incorporating advanced functionalities through JAX in the graphic processing unit support, automatic differentiation, and integration with DNNs. We demonstrated JAX-CanVeg's hybrid modeling capability by applying the model at four flux tower sites with varying aridity. To this end, we developed a hybrid version of the Ball-Berry equation that emulates the water stress impact on stomatal closure to explore the capability of the hybrid model in (a) improving the simulations of latent heat fluxes (LE) and net ecosystem exchange (NEE), (b) improving the optimization trade-off when learning observations of both LE and NEE, and (c) benefiting a multi-layer canopy model setup. Our results show that the proposed hybrid model improved the simulations of LE and NEE at all sites, with an improved optimization trade-off over the process-based model. Additionally, the multi-layer canopy set benefited hybrid modeling at some sites. Anchored in differentiable modeling, our study provides a new avenue for modeling land-atmosphere interactions by leveraging the benefits of both data-driven learning and process-based modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Dynamic clay microstructures emerge via ion complexation waves

Clays control carbon, water and nutrient transport in the lithosphere, promote cloud formation5 and lubricate fault slip through interactions among hydrated mineral interfaces. Clay mineral properties are difficult to model because their structures are disordered, curved and dynamic. Consequently, interactions at the clay mineral-aqueous interface have been approximated using electric double layer models based on single crystals of mica and atomistic simulations. We discover that waves of complexation dipoles at dynamically curving interfaces create an emergent long-range force that drives exfoliation and restacking over time- and length-scales that are not captured in existing models. Curvature delocalizes electrostatic interactions in ways that fundamentally differ from planar surfaces, altering the ratio of ions bound to the convex and concave sides of a layer. Multiple-scattering reconstruction of low-dose energy-filtered cryo electron tomography enabled direct imaging of ion complexes and electrolyte distributions at hydrated and curved mineral interfaces with {\aa}ngstrom resolution over micron length scales. Layers exfoliate and restack abruptly and repeatedly over timescales that depend strongly on the counterion identity, demonstrating that the strong coupling between elastic, electrostatic and hydration forces in clays promote collective reorganization previously thought to be a feature only of active matter.

Whittaker, Michael L↗
DOE OSTI2024

Assessing the Atmospheric Response to Subgrid Surface Heterogeneity in the Single-Column Community Earth System Model, Version 2 (CESM2)

Land-atmosphere interactions are central to the evolution of the atmospheric boundary layer and the subsequent formation of clouds and precipitation. Existing global climate models represent these connections with bulk approximations on coarse spatial scales, but observations suggest that small-scale variations in surface characteristics and co-located turbulent and momentum fluxes can significantly impact the atmosphere. Recent model development efforts have attempted to capture this phenomenon by coupling existing representations of subgrid-scale (SGS) heterogeneity between land and atmosphere models. Such approaches are in their infancy and it is not yet clear if they can produce a realistic atmospheric response to surface heterogeneity. Here, we implement a parameterization to capture the effects of SGS heterogeneity in the Community Earth System Model (CESM2), and compare single-column simulations against high-resolution Weather Research and Forecasting (WRF) large-eddy simulations (LESs), which we use as a proxy for observations. The CESM2 experiments increase the temperature and humidity variances in the lowest atmospheric levels, but the response is weaker than in WRF-LES. In part, this is attributed to an underestimate of surface heterogeneity in the land model due to a lack of SGS meteorology, a separation between deep and shallow convection schemes in the atmosphere, and a lack of explicitly represented mesoscale secondary circulations. These results highlight the complex processes involved in capturing the effects of SGS heterogeneity and suggest the need for parameterizations that communicate their influence not only at the surface but also vertically.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Measurement of atomic oxygen densities using TALIF on a dielectric barrier discharge: insights into the volume above a micro cavity plasma array

Dielectric barrier discharges, particularly micro cavity plasma arrays, offer significant potential for plasma-catalytic research due to their ability to ignite plasma in direct contact with a catalytic surface, enabling the observation of plasma-surface interactions. A key factor in their application is the generation of reactive species, such as atomic oxygen, within the cavities. These species can interact with both the surface (e.g. for activation or cleaning) and the gas being treated (e.g. for oxidation). Given the central role of oxygen atoms in plasma catalysis and their use as a model for more complex species, this work investigates the transport of these atoms out of the cavities. Two-photon absorption laser-induced fluorescence spectroscopy with picosecond laser excitation is performed in the volume above the cavities. The results are compared with a basic diffusion model. The reactor operates with a He/O 2 mixture at a flow rate of 1 slm and atmospheric pressure. Densities of up to 10 16 cm -3 are measured near the surface. Time-dependent measurements show that, at a distance of 350 µm from the surface, a density equilibrium is reached within less than 3 ms of reactor operation. Decay times due to ozone formation after the reactor is turned off are on a similar scale. Spatially resolved measurements show that the oxygen density decreases exponentially from the surface but remains detectable up to approximately 1 mm above the surface, indicating significant application potential. Variations in the O 2 admixture show a density maximum at 0.4%, confirming previous helium state enhanced actinometry measurements within the cavities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

From material properties to multiscale modeling to improve lithium-ion energy storage safety

Energy storage using lithium-ion cells dominates consumer electronics and is rapidly becoming predominant in electric vehicles and grid-scale energy storage, but the high energy densities attained lead to the potential for release of this stored chemical energy. This article introduces some of the paths by which this energy might be unintentionally released, relating cell material properties to the physical processes associated with this potential release. The selected paths focus on the anode–electrolyte and cathode–electrolyte interactions that are of typical concern for current and near-future systems. Relevant material processes include bulk phase transformations, bulk diffusion, surface reactions, transport limitations across insulating passivation layers, and the potential for more complex material structures to enhance safety. We also discuss the development, parameterization, and application of predictive models for this energy release and give examples of the application of these models to gain further insight into the development of safer energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2022

Learning Optimal Aerodynamic Designs

This project created a framework for efficient, accurate, and scalable deep neural network representations of design optimization problem solutions. The inputs to these DNN representations are the vector of design requirement parameters, the outputs are the optimal design variables, and the goal is to learn the map from inputs to outputs (i.e., inverse design). The team addressed the problem of the optimal shape design of aerodynamic lifting surfaces—in particular aircraft wings—using a Reynolds-Average Navier Stokes model to govern the CFD-based aerodynamic shape optimization. The inverse design map for such problems is very complex and high-dimensional, involving inputs and outputs on the order of 1000s. To approximate this inverse design map, the team developed algorithms to construct parsimonious DNN architectures, which automatically identify low-dimensional manifolds in which design requirements affect optimal shape parameters, and trained these architectures with multifidelity optimization methods. The resulting methodology accurately and automatically designs optimal aerodynamic lifting surfaces with very high accuracy (99%) at interactive speeds, of the order of milliseconds, resulting in factors of one million or more speedup relative to CFD-based design optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Plasma Surface Interactions: Predicting the Performance and Impact of Dynamic PFC Surfaces

The objective of this project is to develop, and integrate, high-performance simulation tools capable of predicting plasma-facing component (PFC) operating lifetime and the impact of the evolving surface morphology of tungsten-based PFCs on plasma contamination, including the dynamic recycling of fuel species and tritium retention, in future magnetic fusion devices. Establishing a fundamental physical understanding and developing predictive capabilities of plasma-surface interactions (PSI) requires simultaneously addressing complex and diverse physics occurring over a wide range of length (Angstroms to meters) and time (femtoseconds to years) scales, as well as integrating extensive physical processes across the plasma–surface interface. This requires development of not only detailed physics models and computational strategies at each scale, but also algorithms and methods to couple them effectively in a way that can be robustly validated. Deploying these tools requires the continued development and coupling of leadership-scale computational codes to describe the boundary plasma and the evolving PFC surface, as well as a host of simulations that bridge disparate scales to address complex physical and computational issues at the plasma–surface interface in multi-component materials systems for magnetic fusion energy development beyond ITER.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Competitive formation of intercalated versus supported metal nanoclusters during deposition on layered materials with surface point defects

Intercalated metal nanoclusters (NCs) can be formed under the surface of graphite after sputtering to generate surface “portal” defects that allow deposited atoms to reach the subsurface gallery. Although, there is a competition between formation of supported NCs on top of the surface and intercalated NCs under the surface, the latter only dominating at sufficiently high temperature. A stochastic model incorporating appropriate system thermodynamics and kinetics is developed to capture this complex and competitive nucleation and growth process. Kinetic Monte Carlo simulation shows that the model captures experimental trends observed for Cu and other metals and reveals that higher temperatures are needed to facilitate detachment of atoms from supported NCs enabling them to reach the gallery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Characterizing Spatial and Temporal Variability of California’s Coastal Upwelling Using a Satellite-Derived CUTI Index

Ocean upwelling off California is a wind-driven coastal process with complex interactions between spatial wind patterns and sea surface temperatures. Large-scale spatial and temporal upwelling characterizations are traditionally based on numerical models, while quantifying upwelling with direct observations is not straightforward. To complement model-based upwelling characterizations, here we use satellite estimates of winds, sea surface temperatures, and ocean currents along the California coast to calculate the coastal upwelling transport index (CUTI), originally calculated from a numerical ocean model. The CUTI framework accounts for both Ekman and geostrophic offshore-directed ocean transport, the latter of which is often neglected in satellite-based upwelling assessments. The approach yields the CUTI over a 30-yr period at a daily temporal and 0.25° spatial resolution to analyze upwelling patterns, with focus on potential Morro Bay and Humboldt wind energy areas. Morro Bay (CUTI median of 0.63 m 2 s −1 ) exhibits stronger, more consistent upwelling throughout the year, whereas Humboldt (CUTI median of 0.45 m 2 s −1 ) exhibits a pronounced seasonal cycle with strongest upwelling in the summer and occasional downwelling in winter. The cluster analysis of 1993–2022 satellite maps of CUTI identifies four characteristic spatial upwelling patterns along the coast, corresponding to distinct spatial distributions of sea surface temperature anomalies, Ekman transport, and geostrophic ocean water transport. This dataset provides an observation-based foundation for characterizing natural coastal upwelling off the California coast, showing good agreement with the original model CUTI while being based on satellite data and offering increased spatial resolution.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Temperature, Humidity, and Time-Lapse Video Data from Yosemite National Park, Water Year 2024

This dataset contains time-lapse imagery and distributed measurements of air temperature, relative humidity, dew point, and soil temperature across Yosemite National Park from October 2023 to September 2024. Instruments were deployed at 10 sites in two cross-valley transects as part of the DOE Grant: Seasonal Cycles Unravel Mysteries of Missing Mountain Water organized by Jessica Lundquist (University of Washington), Rosemary Carroll (Desert Research Institute), and Ethan Gutmann (National Center for Atmospheric Research). The data are intended to support hydrologic modeling efforts to better resolve the fate of mountain water, and are published to support studies of surface climate or hydrologic processes in complex terrain. Measurements were collected with low-cost data loggers installed 2 m high on evergreen trees or buried just below the soil surface. A time-lapse camera at one site captures valley-scale seasonal snow cover variability.Dataset files are organized by site and variable (air measurements, ground measurements, or time-lapse video). Air and ground measurements are packaged in LoggerData.zip, and time-lapse imagery is compiled into a short video stored in TimelapseVideos.zip. File-level metadata contains details for each file included in the dataset. A data dictionary provides units and descriptions for column or row names in all files. The locations metadata file describes site characteristics, locations, and associated GPS methods.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Seismic Observations of the OSIRIS-REx Sample Return Capsule Reentry: Deployment, Signal Characteristics, and Wavefield Phenomenology

Controlled spacecraft reentries from interplanetary trajectories provide rare, well-characterized hypersonic sources for advancing seismoacoustic observation techniques. Here we present seismic observations of the OSIRIS-REx sample return capsule (SRC) reentry on 24 September 2023, recorded by 16 three-component nodal seismometers deployed near Eureka, Nevada, at ground distances of 7–20 km from the capsule trajectory. Air-to-ground coupled signals are detected at all stations, exhibiting impulsive onsets consistent with ballistic shock arrivals from the descending Mach cone. We characterize the seismic wavefield through signal amplitude, period, waveform cross-correlation, and array processing. Signal periods decrease systematically with increasing distance from the trajectory within the airport array, indicating that higher-frequency content becomes more prominent at greater offsets, opposite to expectations from geometric spreading and atmospheric absorption. Seismic array processing identifies frequency-dependent back-azimuth variations whose origin remains unresolved; possible contributing factors include source geometry, scattering by fine-scale layered structure in the stratosphere, and near-surface effects. These observations document a spatially complex seismic wavefield from a well-characterized hypersonic line source and provide constraints for future modeling of atmospheric propagation and air-to-ground coupling.

OSIRIS-REx↗
DOE OSTI2020

Electrode Mesoscale as a Collection of Particles: Coupled Electrochemical and Mechanical Analysis of NMC Cathodes

Battery electrodes are composed of polydisperse particles and a porous, composite binder domain. These materials are arranged into a complex mesostructure whose morphology impacts both electrochemical performance and mechanical response. We present image-based, particle-resolved, mesoscale finite element model simulations of coupled electrochemical-mechanical performance on a representative NMC electrode domain. Beyond predicting macroscale quantities such as half-cell voltage and evolving electrical conductivity, studying behaviors on a per-particle and per-surface basis enables performance and material design insights previously unachievable. Voltage losses are primarily attributable to a complex interplay between interfacial charge transfer kinetics, lithium diffusion, and, locally, electrical conductivity. Mesoscale heterogeneities arise from particle polydispersity and lead to material underutilization at high current densities. Particle-particle contacts, however, reduce heterogeneities by enabling lithium diffusion between connected particle groups. While the porous composite binder domain (CBD) may have slower ionic transport and less available area for electrochemical reactions, its high electrical conductivity makes it the preferred reaction site late in electrode discharge. Mesoscale results are favorably compared to both experimental data and macrohomogeneous models. This work enables improvements in materials design by providing a tool for optimization of particle sizes, CBD morphology, and manufacturing conditions.

42 ENGINEERING↗