Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 523 records · Page 29

Lunar rock surfaces as detectors of solar processes

Lunar rock surfaces exposed at or just below the lunar surface are considered as detectors of the solar wind, solar flares and solar-derived magnetic fields through their interactions with galactic cosmic rays. The degradation of the solar detector capabilities of lunar surface rocks by meteoroid impact erosion, accreta deposition, loose dust, and sputtering, amorphous layer formation and accelerated diffusion due to solar particles and illumination is discussed, and it is noted that the complex interactions of factors affecting the outer micron of exposed surface material has so far prevented the development of a satisfactory model for a particle detector on the submicron scale. Methods for the determination of surface exposure ages based on the accumulation of light solar wind noble gases, Fe and Mg, impact craters, solar flare tracks, and cosmogenic Kr isotopes are examined, and the systematic variations in the ages determined by the various clocks are discussed. It is concluded that a means of obtaining satisfactory quantitative rate or flux data has not yet been established.

Hartung, J. B.↗

Effects of vegetation on passive microwave estimates of soil moisture

Data collected over plots using a microwave radiometer with a wavelength of 21 cm were used to evaluate a theoretical model of the effect of vegetative cover on the emission of soils. The theoretical model was used as a basis of a simplified model which predicts the emissivity of a soil-vegetation complex with a full canopy cover. This model, when transformed to predict soil moisture, requires measurements of the microwave brightness temperature, surface temperature and vegetation biomass. Since surface temperature can be determined using thermal infrared measurements and vegetation biomass could be estimated with visible and near infrared measurements, the simplified model could be applied using only remotely sensed data.

Jackson, T. J.↗

Rolling DICE to Advance Knowledge of Land-Atmosphere Interactions

The Diurnal Land-Atmosphere Coupling Experiment (DICE) aims to explore the complex interactions between the land surface and atmospheric boundary layer, which are generally not well understood and difficult to isolate in models. The project involves over 10 different models, combining expertise from both land surface and atmospheric boundary-layer modelling groups. A simple three-stage methodology is designed to assess land-atmosphere feedbacks. Stage 1: the individual components are assessed in isolation, driven and evaluated against observational data; stage 2: the impact of coupling is investigated; stage 3: the sensitivity of the stand-alone models to variations in driving data is explored. For this initial study, a three-day clear-sky period in the mid-west USA over, an assumed simple, predominantly grass surface was simulated using data from the CASES-99 field campaign. Key conclusions from the study include: (1) the memory of vegetation state within Land Surface Models (LSMs) needs attention; (2) the height of atmospheric forcing for LSMs is important, particularly for the nocturnal boundary layer, and this has implications for both observations and vertical resolution for atmospheric models; (3) land-atmosphere feedbacks reduce errors in simulated surface fluxes at the expense of the accuracy of the variables that the models are designed to simulate (e.g., temperature, humidity and wind speed); (4) problems remain in representing the stable boundary layer in atmospheric models; (5) the mixing of temperature and humidity within the boundary layer may need to be represented separately; (6) differences in daytime profiles of heat, moisture and momentum between models are mainly due to the way the models erode the inversion at the top of the boundary layer, rather than differences in the surface fluxes. Resultant variations in modelled boundary layer heights have a substantial impact on relative humidity and could partially explain variations in coupling strength between models in the GLACE experiment.

Martin Best↗

Mesoflow: An Open-Source Reacting Flow Solver for Catalysis at Mesoscale

We present the capabilities and software performance metrics of our open-source continuum solver for catalysis, Mesoflow, developed specifically for modeling transport and chemistry at the mesoscale. Our solver utilizes Cartesian block-structured adaptive mesh refinement to resolve complex catalyst surface morphologies directly obtained from X-ray tomography data. An immersed boundary based formulation enables rapid representation of complex geometries prevalent in most mesoporous catalyst interfaces. The solver is developed on top of open-source performance portable library, AMReX, providing parallel execution capabilities on current and upcoming high-performance-computing (HPC) architectures. Our flexible software framework enables integration of complex chemical mechanisms at heterogenous interfaces and time-split algorithms for circumventing highly disparate reaction and flow time-scales. Our current studies indicate a ten-fold performance gain by using graphics-processing-units (GPUs) compared to a single processor for representative problem sizes (2 million cell mesh). We will also present a brief introduction on how to build and use this software for application problems pertaining to catalytic upgrading and gas transport within porous catalyst particles.

adaptive meshing↗

DIPS-Plus: The enhanced database of interacting protein structures for interface prediction

Abstract In this work, we expand on a dataset recently introduced for protein interface prediction (PIP), the Database of Interacting Protein Structures (DIPS), to present DIPS-Plus, an enhanced, feature-rich dataset of 42,112 complexes for machine learning of protein interfaces. While the original DIPS dataset contains only the Cartesian coordinates for atoms contained in the protein complex along with their types, DIPS-Plus contains multiple residue-level features including surface proximities, half-sphere amino acid compositions, and new profile hidden Markov model (HMM)-based sequence features for each amino acid, providing researchers a curated feature bank for training protein interface prediction methods. We demonstrate through rigorous benchmarks that training an existing state-of-the-art (SOTA) model for PIP on DIPS-Plus yields new SOTA results, surpassing the performance of some of the latest models trained on residue-level and atom-level encodings of protein complexes to date.

59 BASIC BIOLOGICAL SCIENCES↗

Automated calculation and convergence of defect transport tensors

Defect diffusion is a key process in materials science and catalysis, but as migration mechanisms are often too complex to enumerate a priori, calculation of transport tensors typically have no measure of convergence and require significant end-user intervention. These two bottlenecks prevent high-throughput implementations essential to propagate model-form uncertainty from interatomic interactions to predictive simulations. In order to address these issues, we extend a massively parallel accelerated sampling scheme, autonomously controlled by Bayesian estimators of statewide sampling completeness, to build atomistic kinetic Monte Carlo models on a state-space irreducible under exchange and space group symmetries. Focusing on isolated defects, we derive analytic expressions for drift and diffusion coefficients, providing a convergence metric by calculating the Kullback–Leibler divergence across the ensemble of diffusion processes consistent with the sampling uncertainty. The autonomy and efficacy of the method is demonstrated on surface trimers in tungsten and Hexa-interstitials in magnesium oxide, both of which exhibit complex, correlated migration mechanisms.

36 MATERIALS SCIENCE↗

Advances in Land Data Assimilation Systems

Assimilation of remotely sensed land surface observations into regional to global scale numerical models have the potential to significantly advance our ability, to assess, understand, and predict surface water, energy, and carbon cycles. This session seeks to assess the state-of-the-art in data assimilation methods for integrating land surface remote sensing and modeling, with a focus on practical applications and techniques. Assimilated land surface variables of interest include (but are not limited to, soil moisture, surface temperature, snowpack, streamflow, vegetation dynamics, and carbon storage. Contributions describing the development of practical land surface data assimilation methods, multivariate land surface data assimilation strategies, evaluation of the required accuracy and resolution of remote sensing observations, the effects of scale, process complexity, and uncertainty on data assimilation, and the optimal treatment of model and observation errors are encouraged.

Houser, Paul R.↗

A tentative unified sea model for scattering and emission

This work is an attempt to relate active and passive observations of the sea surface to each other, considering that more information can be obtained from a set of instruments when they are not just considered as a juxtaposition of separate and independent boxes, but as a single complex instrument. A careful review of several approaches to the solution of the scattering problem for a random rough surface has been performed. Inherent simplifications, and inconsistencies, are pointed out. To reduce the computation time to reasonable values, an approximate model is selected and some numerical results will be presented in order to illustrate the benefit that can be gained by a unified approach.

Guissard, A.↗

Theoretical evaluation of a V/STOL fighter model utilizing the PAN AIR code

The PAN AIR computer code was investigated as a tool for predicting closely coupled aerodynamic and propulsive flowfields of arbitrary configurations. The NASA/Ames V/STOL fighter model, a configuration of complex geometry, was analyzed with the PAN AIR code. A successful solution for this configuration was obtained when the nozzle exit was treated as an impermeable surface and no wakes were included around the nozzle exit. When separated flow was simulated from the end of the nacelle, requiring the use of wake networks emanating from the nozzle exit, a number of problems were encountered. A circular body nacelle model was used to investigate various techniques for simulating the exhaust plume in PAN AIR. Several approaches were tested and eliminated because they could not correctly simulate the interference effects. Only one plume modeling technique gave good results. A PAN AIR computation that used a plume shape and inflow velocities obtained from the Navier-Stokes solution for the plume produced results for the effects of power that compared well with experimental data.

Howell, G. A.↗

Beyond Idealized Models of Nanoscale Metal Hydrides for Hydrogen Storage

Metal hydrides are attractive for compact, low-pressure hydrogen storage, yet a foundational understanding of factors governing their thermodynamics and kinetics is still lacking. Predictive modeling from the atomic to the microstructural scale plays a critical role in addressing these gaps, particularly for nanoscale materials, which promise improved performance but are difficult to probe. In this paper, we summarize strategies being developed within the Hydrogen Materials—Advanced Research Consortium (HyMARC) for going beyond conventional models to incorporate more complex physics, more realistic structures, and better approximation of operation conditions in simulations of nanoscale metal hydrides. We highlight four beyond-ideal factors that influence predicted performance: (1) surface anharmonic dynamics, (2) interface and surface energy penalties, (3) mechanical stress under confinement, and (4) the presence of native surface oxide. Approaches for addressing these factors are demonstrated on model materials representative of high-capacity hydrogen storage systems, and implications for understanding performance under operating conditions are discussed.

08 HYDROGEN↗

Radiation exchange in large space structures and frames

The analysis of radiation heat transfer is determined by the need to determine the geometrical quantity, F(ij), which is commonly called the view factor, and which represents the fraction of energy leaving surface i which strikes surface j. The present investigation shows that the use of a rotating plane to model a one-dimensional radiation element is both computationally efficient and accurate. Because of the reduced computational cost of this element, radiation view factors for complex structures can be calculated directly. Under some conditions, the characteristics of the one-dimensional element may lead to a loss in accuracy when it is used to model obstructing cylinders. In cases in which the element is too close to the radiating surfaces, curved surfaces must be represented by a collection of flat planes. However, it appears that as long as the one-dimensional element is more than 5 diameters from the other surfaces, its use is valid.

Emery, A. F.↗

Revisiting the Anisotropic Complex Refractive Indices of Sodium Nitrate for Interpretation of the Reflectance Spectrum of Pressed Pellets

Reflectance spectroscopy is notoriously confounding in that the measured spectral response is highly dependent upon the morphology of the sample. Fortunately, all such perturbations are neatly encoded by the complex refractive index of the sample. Herein, we seek to quantitatively model the measured infrared reflectance spectrum of a specularly flat pressed pellet sample of the birefringent compound, sodium nitrate. Single crystals of sodium nitrate were synthesized via a slow evaporation process and spectroscopically analyzed using polarization-dependent infrared single-angle reflectance spectroscopy. The anisotropic complex refractive index was measured from 7500 to 300 cm-1 (1.33 to 33.33 µm). The deduced anisotropic optical constants were found to be consistent with those previously reported. Once measured and validated, the optical constants were applied to model the pressed pellet reflectance spectrum. It was evident that an average of the anisotropic refractive indices was insufficient to account for the measured pellet reflectance. In order to account for contributions of all possible microcrystalline orientations within the pellet, the Python package PYELLI was used to calculate a basis set of orientation-dependent reflection spectra spanning the distinct ? and ? Euler rotations of the uniaxial crystal. When the population of orientations was allowed to vary freely in a spectral fit analysis, the fit-deduced orientations were tightly clustered along f = 45º, hinting at residual anisotropy in the pressed pellet sample. Conversely, an equally valid spectral fit (with marginally worse fit metric) was obtained when the population was constrained to an isotropic distribution of orientations. Subsequent non-zero cross-polarization reflectance measurements likewise suggested anisotropy in the pellet. However, both grazing incidence wide-angle x-ray scattering and scanning electron microscopy measurements revealed that the microcrystal orientations at the surface of the pressed pellet sample were isotropically distributed (and that the average crystallite size was larger than ?/10). Application of the measured complex refractive indices for modeling the reflectance spectrum of the pressed pellet, and rectification of these seemingly contradictory observations will be discussed.

Wilhelm, Michael J.↗

Retrieval of Surface Snow Grain Size and Melt Water from AVIRIS Spectra

The Earth's energy balance and hydrology are affected by the distribution and characteristics of snow cover on the surface. Snow grain size and snow melt influence surface albedo and hydrology. A model of snow reflectance that depends on both grain size and surface melt water was developed to derive these parameters from remote spectral measurements. This reflectance model is based on a discrete ordinate radiative transfer approach that uses Mie calculations of snow optical properties, which are based on the complex refractive index of ice and water. This snow model was linked to an atmospheric radiative transfer code and a nonlinear least squares fitting algorithm. The resulting combined algorithm was applied to an AVIRIS snow data set acquired over Mammoth Mountain, California. Maps of grain size and surface snow melt were generated that are consistent with the expected ranges and distributions for conditions at the site.

Green Robert O.↗

Estimating Watershed Subsurface Permeability From Stream Discharge Data Using Deep Neural Networks

Subsurface permeability is a key parameter in watershed models that controls the contribution from the subsurface flow to stream flows. Since the permeability is difficult and expensive to measure directly at the spatial extent and resolution required by fully distributed watershed models, estimation through inverse modeling has had a long history in subsurface hydrology. The wide availability of stream surface flow data, compared to groundwater monitoring data, provides a new data source to infer soil and geologic properties using integrated surface and subsurface hydrologic models. As most of the existing methods have shown difficulty in dealing with highly nonlinear inverse problems, we explore the use of deep neural networks for inversion owing to their successes in mapping complex, highly nonlinear relationships. We train various deep neural network (DNN) models with different architectures to predict subsurface permeability from stream discharge hydrograph at the watershed outlet. The training data are obtained from ensemble simulations of hydrographs corresponding to an permeability ensemble using a fully-distributed, integrated surface-subsurface hydrologic model. The trained model is then applied to estimate the permeability of the real watershed using its observed hydrograph at the outlet. Our study demonstrates that the permeabilities of the soil and geologic facies that make significant contributions to the outlet discharge can be more accurately estimated from the discharge data. Their estimations are also more robust with observation errors. Compared to the traditional ensemble smoother method, DNNs show stronger performance in capturing the nonlinear relationship between permeability and stream hydrograph to accurately estimate permeability. Our study sheds new light on the value of the emerging deep learning methods in assisting integrated watershed modeling by improving parameter estimation, which will eventually reduce the uncertainty in predictive watershed models.

54 ENVIRONMENTAL SCIENCES↗

Mechanistic and Kinetic Analysis of Complete Methane Oxidation on a Practical PtPd/Al 2 O 3 Catalyst

A PtPd/Al 2 O 3 catalyst developed for the complete oxidation of methane from the ventilation air of underground coal mines is compared against a model PdO/Al 2 O 3 catalyst. Although the PtPd/Al 2 O 3 catalyst is substantially more active and stable than the model catalyst, the nature of active sites between the two catalysts is deemed to be fundamentally the same based on their response to different feed gas compositions and the evolution of surface CO adsorption complexes during time-resolved CO adsorption DRIFTS experiment. For both catalysts, coordinatively unsaturated Pd sites are considered the active centers for methane activation and the subsequent oxidation reaction. H 2 O competes with CH 4 for the same active sites, resulting in severe inhibition. Additionally, the CH 4 oxidation reaction also causes self-inhibition. Taking both inhibition effects into consideration, a relatively simple kinetic model is developed. The model provides a good fit of the 72 sets of kinetic data collected on the PtPd/Al 2 O 3 catalyst under practically relevant reaction conditions with CH 4 concentration in the range of 0.05–0.4%, H 2 O concentration of 1.0–5.0%, and reaction temperatures of 450–700 °C. Kinetic parameters based on the model suggest that the CH 4 activation energy on the PtPd/Al 2 O 3 catalyst is 96.7 kJ/mol, and the H 2 O adsorption energy is –31.0 kJ/mol. Both values are consistent with the parameters reported in the literature. The model can be used to develop catalyst sizing guidelines and be incorporated into the control algorithm of the catalytic system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulation studies for surfaces and materials strength

Computer simulation studies were conducted to investigate materials properties. All the calculations were carried out using atomic level simulation techniques which are based on semiempirical or model functions. During this project, functions with varying degrees of complexity were derived and employed in simulations. Investigations include covalently bonded materials as well as systems involving metal atoms. For small clusters calculated results provided information on various energy- and structure-related properties along with vibrational characteristics. Also, energy barriers for configurational transitions were calculated for selected cases. Simulation calculations for surfaces produced new results in areas related to surface energies, reconstructions and relaxations, surface defects, surface stresses as well as adsorption and nucleation processes. Simulations were also performed on bulk materials. Calculated results provided an atomic level understanding on energetics and structures of point defects, crystal stability, elastic properties, and materials strength for various systems. Calculations involving polymeric materials include studies of polar polymer melts and polymer/solid interfaces. Simulations employing specially developed codes provided significant information about energetics and conformational characteristics of different polymeric chains.

Halicioglu, Timur↗

Prediction of hemiwicking dynamics in micropillar arrays

Dynamic hemiwicking behavior is observable in both nature and a wide range of industrial applications ranging from biomedical devices to thermal management. We present a semi-analytical modeling framework (without empirical fitting coefficients) to predict transient capillary-driven hemiwicking behavior of a liquid through a nano/microstructured surface, specifically a micropillar array. In our model framework, the liquid domain is discretized into micropillar unit cells to enable the time marching of the hemiwicking front. A simplified linear pressure drop is assumed along the hemiwicking length such that the local meniscus curvature, contact angle, and effective liquid height are determined at each time step in our transient model. This semi-analytical model is validated with experimental data from our own experiments and from published literature for different fluids. Our model predicts hemiwicking dynamics with <20% error over a broad range of micropillar geometries with height-to-pitch ratio ranging between ≈0.34 and 6.7 and diameter-to-pitch ratio in the range of ≈0.25–0.7 and without any fitting parameters. For lower diameter-to-pitch ratio data points related to sparse micropillar array arrangements, we suggest modifications to the semi-analytical model. This work sheds light on complex and dynamic solid–liquid–vapor interfacial interactions which could serve as a guide for the design of textured surfaces for wicking enhancement in multi-phase thermal and mass transport technologies and applications.

Mechanics↗

Theoretical insights into the surface physics and chemistry of redox-active oxides

Redox active oxides are ubiquitous in materials applications, including catalysis, photovoltaics, self-cleaning glasses, chemical sensors, and electronic components. Their utility derives from their unique ability to access multiple metal charge states within a finite energy window. However, this property also confounds our ability to study reducible oxides, because it leads to structural, compositional, and electronic complexities that elude simplistic models of materials structure and function. Oxygen vacancies play a critical role in shaping the functional properties of such oxides; most notably, they lead to mobile charge imbalances that impact surface processes at significant distances from the originating defect. Atomistic simulations are inherently equipped to illuminate these phenomena at a fundamental level; however, reducible oxides pose great challenges due to the high level of electron correlation needed to correctly describe them. Knowing how defects form, couple, propagate, agglomerate, or repel each other and influence the surface properties of reducible oxides is only now coming into the grasp of modern theory and simulation capabilities. This knowledge is also key to discovering and controlling emergent materials processes at nanometer scale and beyond. The manuscript was written with support of all three authors from U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. Work by RR and VG was performed at Pacific Northwest National Laboratory (PNNL), which is a multi-program laboratory operated by Battelle for DOE under Contract DE AC05 76RL01830. AS was supported under Award DE-SC0007347.

Rousseau, Roger J.↗