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 451 records · Page 25

Evaluating Probability of Containment Effectiveness at a GCS Sites using integrated assessment modeling approach with Bayesian decision Networks

Improved scientific and engineering understanding of the behavior of geologic CO2 storage together with established regulatory framework and incentive structures raise the prospects for accelerated, large-scale deployment of this greenhouse gas emissions reduction approach. Incentive structures call for the establishment of appropriate verification and accounting approaches to support claims of the integrity of a geologic storage complex and to justify taking credit for long-term storage. In this study, we present a framework for assessing the probability of containment effectiveness over the lifetime of a geologic carbon storage site (e.g., after 70 years of injection and post-injection site performance) using forward stochastic model realizations based on site characterization data and using a monitoring-informed Bayesian network based on hypothetical detectability from surface seismic surveys over the site injection and post-injection phases. The National Risk Assessment Partnership’s open-source Integrated Assessment Model (NRAP-Open-IAM) was utilized to develop an ensemble of 10,000 a priori stochastic forecasts of CO2 containment. Those simulations were used to train the Bayesian network model to estimate the prior probabilities of the CO2 leakage mass into overlying, monitorable aquifers considering the uncertainties in the reservoir properties, permeability of potentially leaky wells and the overlying aquifers. The conditional probabilities in the Bayesian network were either learned from the NRAP-Open-IAM simulations or derived from the predefined detection thresholds for the monitoring method. Observations obtained from monitoring, over time during the site operation phases were then used to generate updated posterior probabilities of containment (and any loss from containment) in the Bayesian network by propagating the prior probabilities through the conditional probabilities. We demonstrate how to construct and use the Bayesian network for verifying the long-term storage complex effectiveness informed by monitoring based on the NRAP-Open-IAM simulations previously developed for the FutureGen 2.0 site. This approach may have relevance for stake holders to demonstrate secure geologic storage, provide a defensible, probabilistic approach to claim credit for geologic storage, and to estimate the likelihood that any fraction of the claimed credit may need to be refunded to the creditor based on available monitoring information.

Bayesian network, Risk assessment, Monitoring, car↗

Molecular-level understanding of metal ion retention in clay-rich materials

Clay minerals retain or adsorb metal ions in the Earth’s critical zone. Rocks, sediments and soils rich in clay minerals can concentrate rare earth elements (REEs) in ion adsorption-type deposits (IADs) and are similarly effective at metallic contaminant remediation. However, the molecular-scale chemical and physical mechanisms of metal ion retention remain only partly understood. Here, we describe the nature, location and energy requirements of metal retention at clay mineral surfaces. Retention originates mainly from electrostatic interactions during cation exchange at low pH and chemical bonding in surface complexation and precipitation at neutral and high pH. Surface complexation can induce surface redox reactions and precipitation mechanisms including neoformation of clay mineral layered structures. In IADs, outer-sphere adsorption is the major retention mechanism of REE ions. By contrast, the use of clay minerals in pollution control relies on various mechanisms that can coexist, including cation exchange, surface complexation and nucleation growth. To more effectively leverage clay mineral–metal interactions in resource recovery and contaminant remediation, complex mechanisms such as surface precipitation and redox reactions must be better understood; for instance, by utilizing advances in quantum mechanical calculations, close combination between synchrotron and simulation techniques, and upscaling of molecular-level information in macroscopic thermokinetic predictive models.

54 ENVIRONMENTAL SCIENCES↗

Performance of Closed Loop Venturi Aspirated Aeration System: Experimental Study and Numerical Analysis with Discrete Bubble Model

In wastewater treatment plants, aeration plays a significant role as it encourages aerobic respiration of microbes, which are necessary to break down carbonaceous matter in the waste stream. This process can account for the majority of energy use in wastewater treatment plants. The aeration process is also necessary in odor control in lagoons and in the aquaculture industry. Generally, the aeration process is accomplished with compressors or blowers which can be of low efficiency due to ideal gas laws. This study introduces the idea of increasing aeration efficiency by looping water from a reservoir through a piping network which includes a venturi aspirator at its inlet. For this purpose, an experimental study has been conducted in a laboratory setup with a pump which pulls water from a tank, passes it through a Venturi aspirator and a helical piping network intended to increase bubble residence time, before depositing it back into the bulk fluid tank. This same process is modeled computationally, using a discrete bubble method (DBM), with good agreement with experiments. The overall purpose here is to determine the optimal configurations for standard aeration efficiency (SAE) and the standard oxygen transfer rate (SOTR). A parametric study has been implemented using the DBM based on different hydraulic and flow parameters. The model is also used to predict the SAE of a hypothetical aeration system. Results indicate that it is possible to achieve SAE values in the range of surface aerators or submersed jet aerators using the proposed aeration system with less complex components, thereby decreasing overall costs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Computational insights into hydrogen adsorption energies on medium-entropy oxides

High entropy oxides (HEOs) have emerged as promising catalysts for several important chemical transformations including alkane activation. Hydrogen adsorption energy (HAE) has been used as a key descriptor for many reactions including methane C–H activation and hydrogen evolution reactions. Hence, understanding the relationship between HAEs and the surface chemistry of HEO surfaces could lay the foundation for meaningful correlations among methane C–H activation, HAE, and the complex, local environment of HEO surfaces. Here, we used a medium-entropy oxide as a prototypical system – Mg 0.25 Ni 0.25 Cu 0.25 Zn 0.25 O with a rock-salt structure – to interrogate these relationships. We sampled 2000 different surfaces of its (100) plane and calculated the HAEs at randomly chosen surface O sites using density functional theory (DFT). Our analysis of the 2000 data points reveals that the HAEs at the surface O sites are significantly influenced by the local environment around the adsorption sites, particularly the nature of the metal atom directly below the surface O site where H adsorbs. After comparing several popular graph-neural-network-based machine learning models, we found that the DimeNet++ model performed best achieving satisfactory accuracy in predicting HAEs for both Mg 0.25 Ni 0.25 Cu 0.25 Zn 0.25 O and slightly varied compositions. Our work underscores the promise of such models and the need for further refinement to address the complexity of HEOs.

Song, Haohong [Vanderbilt Univ., Nashville, TN (Un↗

Revolutionizing Materials Design: The Intersection of Quantum Mechanics and Data Modeling

The field of materials design is currently experiencing a notable evolution, driven by the convergence of sophisticated computational methodologies based on first principles and data-driven modeling approaches. I will review our recent endeavors employing AI/ML to expedite first-principles simulations and mitigate traditional methods' temporal and spatial limitations. Central to our efforts is developing and utilizing ML interatomic potentials (MLPs) across a diverse spectrum of materials. We show that MLPs serve as invaluable tools for navigating the complexities of the simulations, such as understanding the behavior of MgO at extreme environments of ~1 terapascal and temperatures >10,000 Kelvin. Moreover, we show that MLPs can provide precise details of the intricate dynamics governing the oxidation processes of binary alloy systems due to the competition between surface segregation and reconstruction tendencies. In summation, advancements in MLPs open the door to fresh possibilities in material modeling and, ultimately, discovery.

Saidi, Wissam↗

Atomic-Scale View of Redox Induced Changes for Monolayer MoO x on α-TiO 2 (110) with Chemical-State Sensitivity

Supported molybdenum oxide (MoO x ) plays an important role in catalytic transformations from alcohol dehydrogenation to transesterification. During these reactions, Mo and oxygen surface species undergo structural and chemical changes. A detailed, chemical-state specific, atomic-scale structural analysis of the catalyst under redox conditions is important for improving catalytic properties. In this study, monolayer of MoO x grown on α-TiO 2 (110) by atomic-layer deposition is analyzed by X-ray standing wave (XSW) excited X-ray photoelectron spectroscopy (XPS). The chemical shifts for Mo 2p 3/2 and O 1s peaks are used to distinguish Mo 6+ from Mo 4+ and surface O from bulk O. Excitation of XPS by XSW allows pin-pointing the location of these surface species relative to the underlying substrate lattice. Finally, measured 3D composite atomic density maps for the oxidized and reduced interfaces compare well with our density functional theory models, and collectively create a unique view of the redoxdriven dynamics for this complex catalytic structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydrobiogeochemical Controls on the Delivery of Dissolved Organic Matter to Boreal Headwater Streams

Understanding controls on the delivery of dissolved organic matter (DOM) from terrestrial to aquatic systems is key for constraining carbon balances and fluxes in boreal headwater catchments. The largest export of DOM generally occurs during intense periods of water delivery to the landscape (e.g., rain events); however, the timing and magnitude of DOM is complicated by geomorphological, hydrometeorological, and biogeochemical variability. To better assess mechanisms controlling the delivery of DOM to boreal streams, DOM was investigated in two morphologically distinct catchments of an experimental forest watershed in western Newfoundland, one dominated by low relief wetlands and the other by steep hillslopes. The DOM was compared during two fall storms of similar magnitude but contrasting moisture conditions during the autumn transition. Variable concentrations and optical character highlighted antecedent conditions are important for the delivery of DOM to boreal streams. Concentration-discharge relationships with stream hydrograph separation analysis suggest preferential flowpaths through shallow mineral horizons as a key pathway for delivery of DOM in the hillslope dominated catchment compared to rapid input from near-stream wetlands in the low relief catchment. Loss of DOM via preferential pathways suggests an additional mechanism for surface soil carbon loss not directly informed by lateral flow and soil hydrology, which has important implications for our conceptualization and modeling of carbon fluxes in boreal systems. These results suggest watershed scale responses of DOM in boreal headwater catchments are complex, and better linkages between catchment scale hydrology and landscape hydrobiogeochemistry are required to constrain terrestrial to aquatic carbon fluxes in boreal landscapes.

54 ENVIRONMENTAL SCIENCES↗

Emergent actin flows explain distinct modes of gliding motility

During host infection, Toxoplasma gondii and related unicellular parasites move using gliding, which differs fundamentally from other known mechanisms of eukaryotic cell motility. Gliding is thought to be powered by a thin layer of flowing filamentous (F)-actin sandwiched between the plasma membrane and a myosin-covered inner membrane complex. How this surface actin layer drives the various gliding modes observed in experiments—helical, circular, twirling and patch, pendulum or rolling—is unclear. Here we suggest that F-actin flows arise through self-organization and develop a continuum model of emergent F-actin flow within the confines provided by Toxoplasma geometry. In the presence of F-actin turnover, our model predicts the emergence of a steady-state mode in which actin transport is largely directed rearward. Removing F-actin turnover leads to actin patches that recirculate up and down the cell, which we observe experimentally for drug-stabilized actin bundles in live Toxoplasma gondii parasites. These distinct self-organized actin states can account for observed gliding modes, illustrating how different forms of gliding motility can emerge as an intrinsic consequence of the self-organizing properties of F-actin flow in a confined geometry.

59 BASIC BIOLOGICAL SCIENCES↗

Machine Learning Inference of Random Medium Properties

Earth materials are heterogeneous across a range of spatial scales, but the resolvability of small structures is limited by sparse data coverage, noise, bandlimitedness, and other difficulties. In practice, heterogeneities below a certain size cannot be recovered from seismic data except through statistical medium descriptions, which even then can be difficult to uniquely determine. To improve the characterization of such heterogeneities, we develop a novel supervised machine learning (ML) model that provides insight about the recoverability of statistical medium properties from elastic waveform data and succeeds despite cycle-skipping and other challenges well known from elastic waveform inversion. We demonstrate the approach using random media generated by superimposing self-affine random variations on homogeneous and layered background structures. After training on sparsely-recorded, high-frequency waveforms from hundreds of different random medium realizations, we show the ability of our ML model to recover correlation lengths and other statistical properties of interest to near-surface and crustal seismology, among other fields. For frequency passbands and spatial offsets encountered in seismology, Gaussian correlation lengths and the amplitude of the random variations relative to the background model are recovered even in challenging scenarios involving unknown medium parameters, complex crustal structures, and low signal-to-noise ratio. In comparison, von Kármán correlation lengths, which are related to larger-wavelength variations of the medium than Gaussian correlation lengths, are not as well recovered. These results provide one of the first and most systematic investigations of the recoverability of statistical properties of heterogeneities below the resolution limit of deterministic seismic tomography, and suggest practical ML strategies for high-frequency waveform seismology.

58 GEOSCIENCES↗

Atmospheric pressure plasmas interacting with wet and dry microchannels: reverse surface ionization waves

Atmospheric pressure plasma jets (APPJs) are increasingly being used to functionalize polymers and dielectric materials for biomedical and biotechnology applications. Once such application is microfluidic labs-on-a-chip consisting of dielectric slabs with microchannel grooves hundreds of microns in width and depth. The periodic channels, an example of a complex surface, present challenges in terms of directly and uniformly exposing the surface to the plasma. In this paper, we discuss results from computational and experimental investigations of negative APPJs sustained in Ar/N 2 mixtures flowing into ambient air and incident onto a series of microchannels. Results from two-dimensional plasma hydrodynamics modeling are compared to experimental measurements of electric field and fast-camera imaging. The propagation of the plasma across dry microchannels largely consists of a sequence of surface ionization waves (SIWs) on the top ridges of the channels and bulk ionization waves (IWs) crossing over the channels. The IWs are directed into electric field enhanced vertices of the next ridge. The charging of these ridges produce reverse IWs responsible for the majority of the ionization. The propagation of the plasma across water filled microchannels evolve into hopping SIWs between the leading edges of the water channels, regions of electric enhancement due to polarization of the water. Positive, reverse IWs follow the pre-ionized path of the initial negative waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A level-set immersed boundary method for reactive transport in complex topologies with moving interfaces

A simulation framework based on the level-set and the immersed boundary methods (LS-IBM) has been developed for reactive transport problems in porous media involving a moving solid-fluid interface. The interface movement due to surface reactions is tracked by the level-set method, while the immersed boundary method captures the momentum and mass transport at the interface. The proposed method is capable of accurately modeling transport near evolving boundaries in Cartesian grids. The framework formulation guarantees second order accuracy in space. Since the interface velocity is only defined at the moving boundary, an interface velocity propagation method is also proposed. The method can be applied to other moving interface problems of the “Stefan” type. Here, we validate the proposed LS-IBM both for flow and transport close to an immersed object with reactive boundaries as well as for crystal growth. Lastly, the proposed method provides a powerful tool to model more realistic problems involving moving reactive interfaces in complex domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Microstructure Generation via Generative Adversarial Network for Heterogeneous, Topologically Complex 3D Materials

Using a large-scale, experimentally captured 3D microstructure dataset, we implement the generative adversarial network (GAN) framework to learn and generate 3D microstructures of solid oxide fuel cell electrodes. The generated microstructures are visually, statistically, and topologically realistic, with distributions of microstructural parameters, including volume fraction, particle size, surface area, tortuosity, and triple phase boundary density, being highly similar to those of the original microstructure. These results are compared and contrasted with those from an established, grain-based generation algorithm (DREAM.3D). Importantly, simulations of electrochemical performance, using a locally resolved finite element model, demonstrate that the GAN generated microstructures closely match the performance distribution of the original, while DREAM.3D leads to significant differences. Finally, the ability of the generative machine learning model to recreate microstructures with high fidelity suggests that the essence of complex microstructures may be captured and represented in a compact and manipulatable form.

36 MATERIALS SCIENCE↗

Projecting Antarctica's contribution to future sea level rise from basal ice shelf melt using linear response functions of 16 ice sheet models (LARMIP-2)

Abstract. The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty of Antarctica's future contribution to global sea level rise that arises from large uncertainty in the oceanic forcing and the associated ice shelf melting. Ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean. However, by computing only the sea level contribution in response to ice shelf melting, our study is neglecting a number of processes such as surface-mass-balance-related contributions. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any self-dampening or self-amplifying processes. This is particularly relevant in situations in which an instability is dominating the ice loss. The results obtained here are thus relevant, in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong ocean warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014) but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP5 models in relation to the global mean surface warming), and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 of the ice loss due to basal ice shelf melting is 10.2 mm, with a likely range between 5.2 and 21.3 mm. For the same period the Antarctic ice sheet lost mass equivalent to 7.4 mm of global sea level rise, with a standard deviation of 3.7 mm (Shepherd et al., 2018) including all processes, especially surface-mass-balance changes. For the unabated warming path, Representative Concentration Pathway 8.5 (RCP8.5), we obtain a median contribution of the Antarctic ice sheet to global mean sea level rise from basal ice shelf melting within the 21st century of 17 cm, with a likely range (66th percentile around the mean) between 9 and 36 cm and a very likely range (90th percentile around the mean) between 6 and 58 cm. For the RCP2.6 warming path, which will keep the global mean temperature below 2 ∘C of global warming and is thus consistent with the Paris Climate Agreement, the procedure yields a median of 13 cm of global mean sea level contribution. The likely range for the RCP2.6 scenario is between 7 and 24 cm, and the very likely range is between 4 and 37 cm. The structural uncertainties in the method do not allow for an interpretation of any higher uncertainty percentiles. We provide projections for the five Antarctic regions and for each model and each scenario separately. The rate of sea level contribution is highest under the RCP8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade, with a likely range between 2 and 9 cm per decade and a very likely range between 1 and 14 cm per decade.

54 ENVIRONMENTAL SCIENCES↗

Goethite and Hematite Nucleation and Growth from Ferrihydrite: Effects of Oxyanion Surface Complexes

The presence of oxyanions, such as nitrate (NO 3 – ) and phosphate (PO 4 3– ), regulates the nucleation and growth of goethite (Gt) and hematite (Hm) during the transformation of ferrihydrite (Fh). Our previous studies showed that oxyanion surface complexes control the rate and pathway of Fh transformation to Gt and Hm. However, how oxyanion surface complexes control the mechanism of Gt and Hm nucleation and growth during the Fh transformation is still unclear. Here, we used synchrotron scattering methods and cryogenic transmission electron microscopy to investigate the effects of NO 3 – outer-sphere complexes and PO 4 3– inner-sphere complexes on the mechanism of Gt and Hm formation from Fh. Our TEM results indicated that Gt particles form through a two-step model in which Fh particles first transform to Gt nanoparticles and then crystallographically align and grow to larger particles by oriented attachment (OA). In contrast, for the formation of Hm, imaging shows that Fh particles first aggregate and then transform to Hm through interface nucleation. This is consistent with our X-ray scattering results, which demonstrate that NO 3 – outer-sphere and PO 4 3– inner-sphere complexes promote the formation of Gt and Hm, respectively. These results have implications for understanding the coupled interactions of oxyanions and iron oxy-hydroxides in Earth-surface environments.

54 ENVIRONMENTAL SCIENCES↗

Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry

A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.

Mechanics↗

Ca-dimers, solvent layering, and dominant electrochemically active species in Ca(BH4)2 in THF

Abstract Divalent ions (Mg, Ca, and Zn) are being considered as competitive, safe, and earth-abundant alternatives to Li-ion electrochemistry, but present challenges for stable cycling due to undesirable interfacial phenomena. We explore the formation of electroactive species in the electrolyte Ca(BH 4 ) 2 ∣THF using molecular dynamics coupled with a continuum model of bulk and interfacial speciation. Free-energy analysis and unsupervised learning indicate a majority population of neutral Ca dimers and monomers with diverse molecular conformations and an order of magnitude lower concentration of the primary electroactive charged species – the monocation, $${\rm{CaBH}}_{4}^{+}$$ CaBH 4 + – produced via disproportionation of neutral complexes. Dense layering of THF molecules within ~1 nm of the electrode surface strongly modulates local electrolyte species populations. A dramatic increase in monocation population in this interfacial zone is induced at negative bias. We see no evidence for electrochemical activity of fully-solvated Ca 2+ . The consequences for performance are discussed in light of this molecular-scale insight.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transport and Assembly of Magnetic Surface Rotors**

This minireview focuses on recent advances with surface magnetic rotors, namely field-responsive spherical or anisotropic microparticles that translate close to, or are embedded in a confining surface. The application of external magnetic modulations allows these microscopic wheels to be remotely spun and steered while also tuning their interactions and inducing assembly from a collection of disordered, moving units. With optical microscopy one can observe and characterize the complex collective phenomena that emerge in dissipative colloidal systems driven far from equilibrium by external fields. From a technological point of view, magnetic surface rotors envisage implementation into microfluidic devices, and can be used for drug delivery, mixing as well as a model system for biological active matter.

36 MATERIALS SCIENCE↗

The black hole interior from non-isometric codes and complexity

Quantum error correction has given us a natural language for the emergence of spacetime, but the black hole interior poses a challenge for this framework: at late times the apparent number of interior degrees of freedom in effective field theory can vastly exceed the true number of fundamental degrees of freedom, so there can be no isometric (i.e. inner-product preserving) encoding of the former into the latter. In this paper we explain how quantum error correction nonetheless can be used to explain the emergence of the black hole interior, via the idea of “non-isometric codes protected by computational complexity”. We show that many previous ideas, such as the existence of a large number of “null states”, a breakdown of effective field theory for operations of exponential complexity, the quantum extremal surface calculation of the Page curve, post-selection, “state-dependent/state-specific” operator reconstruction, and the “simple entropy” approach to complexity coarse-graining, all fit naturally into this framework, and we illustrate all of these phenomena simultaneously in a soluble model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗