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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Towards replacing physical testing of granular materials with a Topology-based Model

In the study of packed granular materials, the performance of a sample (e.g., the detonation of a high-energy explosive) often correlates to measurements of a fluid flowing through it. The “effective surface area,” the surface area accessible to the airflow, is typically measured using a permeametry apparatus that relates the flow conductance to the permeable surface area via the Carman-Kozeny equation. This equation allows calculating the flow rate of a fluid flowing through the granules packed in the sample for a given pressure drop. However, Carman-Kozeny makes inherent assumptions about tunnel shapes and flow paths that may not accurately hold in situations where the particles possess a wide distribution in shapes, sizes, and aspect ratios, as is true with many powdered systems of technological and commercial interest. To address this challenge, we replicate these measurements virtually on micro-CT images of the powdered material, introducing a new Pore Network Model based on the skeleton of the Morse-Smale complex. Pores are identified as basins of the complex, their incidence encodes adjacency, and the conductivity of the capillary between them is computed from the cross-section at their interface. We build and solve a resistive network to compute an approximate laminar fluid flow through the pore structure. Here, we provide two means of estimating flow-permeable surface area: (i) by direct computation of conductivity, and (ii) by identifying dead-ends in the flow coupled with isosurface extraction and the application of the Carman-Kozeny equation, with the aim of establishing consistency over a range of particle shapes, sizes, porosity levels, and void distribution patterns.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Seismic Characterization of the Blue Mountain Geothermal Field

Subsurface characterization is crucial for geothermal energy exploration and production. Yet hydrothermal reservoirs usually reside in highly fractured and faulted zones where accurate characterization is very challenging because of low signal-to-noise ratios of land seismic data and lack of coherent reflection signals. We perform an active-source seismic characterization for the Blue Mountain geothermal field in Nevada using active seismic data to reveal the elastic medium property complexity and fault distribution at this field. We first employ an unsupervised machine learning method to attenuate groundroll and near-surface guided-wave noise and enhance coherent reflection and scattering signals from noisy seismic data. We then build a smooth initial P-wave velocity model based on an existing magnetotellurics survey result, and use 3D first-arrival traveltime tomography to refine the initial velocity model. We then derive a set of elastic wave velocities and anisotropic parameters using elastic full-waveform inversion, and obtain PP and PS images using elastic reverse-time migration. We identify major faults by analyzing the variations of seismic velocities and anisotropy parameters, and reveal mid- to small-scale faults by applying a supervised machine learning method to the seismic migration images. Our characterization reveals complex velocity heterogeneities and anisotropies, as well as faults, with a high spatial resolution. These results can provide valuable information for optimal placement of future injection and production wells to increase geothermal energy production at the Blue Mountain geothermal power plant.

58 GEOSCIENCES↗
DOE OSTI2025

Dynamics and activation of membrane-bound B cell receptor assembly

B-cell receptor (BCR) complexes are expressed on the surface of a B-cell and are critical in antigen recognition and modulating the adaptive immune response. Even though the relevance of antibodies has been known for almost a hundred years, the antigen-dependent activation mechanism of B-cells has remained elusive. Several models have been proposed for BCR activation, including cross-linking, conformation-induced oligomerization, and dissociation activation models. Recently, the first cryo-EM structures of the human B-cell antigen receptor of the IgM and IgG isotypes have been published that validates the asymmetric organization of the BCR complex. Here, we carry out extensive molecular dynamics simulations to probe the conformational changes upon antigen binding and the influence of the membrane lipids. We identify two critical dynamical events that could be associated with antigen-dependent activation of BCR. First, antigen binding causes increased flexibility in regions distal to the antigen binding site. Second, antigen binding alters the rearrangement of IgM transmembrane helices, including the relative interaction of Igα/Igβ that mediates intracellular signaling. Furthermore, these transmembrane rearrangements lead to changes in localized lipid composition. Our work indirectly supports the conformational-change induced models of BCR activation and contributes to the understanding of the antigen-dependent activation mechanism of BCRs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Enhanced Menshutkin S N 2 Reactivity in Mesoporous Silica: The Influence of Surface Catalysis and Confinement

A significant enhancement in the Menshutkin S N 2 reaction between 1-methylimidazole (MeIm) and methyl thiocyanate (MeSCN) is observed when the reaction is confined in the nanoscale silica pores of MCM41 and SBA15. The experiments in the silica pores are conducted without surrounding bulk reaction mixture. The influences of temperature, pore radius, and surface chemistry on the kinetics of the confined reaction are analyzed with time-dependent infrared spectroscopy, molecular dynamics simulations, and ab initio calculations. The rate constant of the pseudo-first order reaction increases with decreasing pore size, and the activation energy is found to decrease by 5.6 kJ/mol in the smallest pore studied (2.8 nm) relative to the bulk reaction. The rate constant dependence on pore size is accurately described by a two-state model in which molecules within the 4.6 Å interfacial layer experience a 2.4-fold rate constant increase relative to those reacting at the bulk rate further away from the interface. The removal of polar silanol groups from the silica surface via passivation with trimethylsilyl chloride results in bulk-like kinetics despite a reduction in the pore diameter, demonstrating the role of silanols as catalytic sites. Here, electronic structure calculations of the energy profile on a model silica surface confirm that silanol groups, particularly those of the vicinal type, can reduce the activation energy and reaction endothermicity through the donation of hydrogen bonds to the reactant, transition state, and product complexes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

MCNP6.3 Unstructured Mesh Verification: GodivR and CANDU Models

A geometric cell of the Monte Carlo N-Particle (MCNP)1 transport code is traditionally created by using Boolean operators on defined surfaces. This constructive solid geometry (CSG) capability has been available in the MCNP code since its beginning. However, a CSG model approach is limited when it comes to constructing a representative geometry for a complex model in its ability to capture a correct model representation. Starting with the version 6.0, the MCNP code has the ability of embedding an unstructured mesh (UM) model into a CSG cell to create a hybrid geometry [1]. The MCNP UM feature provides the flexibility of defining very complex geometries because computer aided design (CAD) and mesh generation software packages can be utilized to construct UM models for MCNP simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

An Infrared Database of n/k Optical Constants for Calculating Aerosol Spectra for the PICARD Program

The objective of the PICARD Program is to develop fieldable sensing platforms for the rapid chemical identification of aerosol particles in plumes. Standoff detection involves interrogating the aerosol cloud from a distance using optical methods and probing the signal returned from direct backscattering from the aerosol particles or transmitted through the plume after reflection from a retroreflector or surface of opportunity (SOO). The identification of chemical species, however, in aerosols is complicated by their complex compositions and morphologies, chemical interferants, and non-uniform particle sizes. To advance standoff detection of aerosols, modelling the infrared transmittance, reflectance and scattering spectra of aerosolized liquid and solid chemical compounds is required, and then testing that model experimentally via laboratory and field experiments. To perform accurate modeling, the infrared optical constants (n/k), i.e., the complex refractive index, of the compounds of interest are required. Thus, PNNL was tasked to provide the optical constants for a set of analytes relevant to the PICARD program. This report describes the experimental techniques used to derive the optical constants of both liquid and solid compounds using established “gold-standard” protocols, how the experimental data are processed to produce the wavenumber dependent optical constant vectors, and how the data are used in the aerosol absorption spectra modeling.

complex refractive index↗
DOE OSTI2020

Experimental and theoretical study of energy transfer in a chromophore triad: What makes modeling dynamics successful?

Simulation of electronic dynamics in realistically large molecular systems is a demanding task that has not yet achieved the same level of quantitative prediction already realized for its static counterpart. This is particularly true for processes occurring beyond the Born–Oppenheimer regime. Non-adiabatic molecular dynamics (NAMD) simulations suffer from two convoluted sources of error: numerical algorithms for dynamics and electronic structure calculations. While the former has gained increasing attention, particularly addressing the validity of ad hoc methodologies, the effect of the latter remains relatively unexplored. Indeed, the required accuracy for electronic structure calculations to reach quantitative agreement with experiment in dynamics may be even more strict than that required for static simulations. In this study, we address this issue by modeling the electronic energy transfer in a donor–acceptor–donor (D–A–D) molecular light harvesting system using fewest switches surface hopping NAMD simulations. In the studied system, time-resolved experimental measurements deliver complete information on spectra and energy transfer rates. Subsequent modeling shows that the calculated electronic transition energies are “sufficiently good” to reproduce experimental spectra but produce over an order of magnitude error in simulated dynamical rates. We further perform simulations using artificially shifted energy gaps to investigate the complex relationship between transition energies and modeled dynamics to understand factors affecting non-radiative relaxation and energy transfer rates.

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

Groundwater characterization and monitoring at a complex industrial waste site using electrical resistivity imaging

A contaminated industrial waste site in Washington State (USA) containing buried, metallic-waste storage tanks, pipes, and wells, was evaluated to determine the feasibility of monitoring groundwater remediation activities associated with an underlying perched aquifer system using electrical resistivity tomography. The perched aquifer, located ~65 m below ground surface and ~10 m above the regional water table, contains high concentrations of nitrate, uranium, and other contaminants of concern from past tank leaks and intentional releases of wastes to surface disposal sites. The extent of the perched water aquifer is not well known, and the effectiveness of groundwater extraction for contaminant removal is uncertain, so supplemental characterization and monitoring technologies are being evaluated. Numerical simulations of subsurface flow and contaminant transport were performed with a highly resolved model of the hydrogeologic system and waste site infrastructure, and these simulations were used as the physical basis for electrical resistivity tomography modeling. The modeling explicitly accounted for metallic infrastructure at the site. The effectiveness of using surface electrodes versus surface and horizontal subsurface electrodes, for imaging groundwater extraction from the perched water aquifer, was investigated. Although directional drilling is a mature technology, its use for electrode emplacement in the deep subsurface under a complex industrial waste site via horizontal wells has not yet been demonstrated. Results from this study indicate that using horizontal subsurface electrode arrays could significantly improve the ability of electrical resistivity tomography to image deep subsurface features and monitor remediation activities under complex industrial waste sites.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Meshfree simulation and experimental validation of extreme thermomechanical conditions in friction stir extrusion

Abstract Friction stir extrusion (FSE) is a novel solid-phase processing technique that consolidates and extrudes metal powders, flakes, chips, or billets into high-performance parts by plastic deformation, which has the potential to save substantial processing time and energy. Currently, most studies on FSE are experimental and only a few numerical models have been developed to explain and predict the complex physics of the process. In this work, a meshfree simulation framework based on smoothed particle hydrodynamics (SPH) was developed for FSE. Unlike traditional grid-based methods, SPH is a Lagrangian particle-based method that can handle severe material deformations, capture moving interfaces and surfaces, and monitor the field variable histories explicitly without complicated tracking schemes. These aspects of SPH make it attractive for the FSE process, where in situ evolution of field variables is difficult to observe experimentally. To this end, a 3-D, fully thermomechanically coupled SPH model was developed to simulate the FSE of aluminum wires. The developed model was thoroughly validated by comparing the numerically predicted material flow, strain, temperature history, and extrusion force with experimental results for a certain set of process parameters. The validated SPH model can serve as an effective tool to predict and better understand the extreme thermomechanical conditions during the FSE process.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Model Package Report: Composite Analysis Solid Waste Release Model (CASWR Model)

This document describes the implementation of a solid waste form release model in GoldSim for the Hanford Site Composite Analysis (CA) Update. This Composite Analysis Solid Waste Release model (CASWR model) was designed to generate deterministic radionuclide release rates for Hanford’s Central Plateau solid waste disposal sites using single realizations of release model coefficients. Five generalized waste form types are used for the conceptual model of waste release: surplus reactor block, cement, soil-debris, grouted residual waste, and ancillary equipment. The surplus reactor block waste form consists of radionuclide leaching from surplus production graphite reactor core blocks via unspecified processes (White et al., 1984 as cited in PNNL-15965). The cement waste form represents solidified wastes whose permeability is much lower than that of the surrounding soil. The soil-debris waste form type is defined as unconsolidated waste mixed with soil material. Tanks and canyon complexes comprise the grouted residual waste form such that their solid waste will be grouted and capped with a surface barrier at the completion of their cleanup. Finally, the ancillary equipment waste form constitutes contaminant releases from ancillary and auxiliary waste form residues associated with tank farms at closure. Individual sub-models are implemented to numerically represent a respective waste form within the CASWR Model: Surplus Reactor Block Sub-model, Cement Sub-model, SoilDebris Sub-model, Grouted Residual Waste Sub-model, and Ancillary Equipment Submodel. Advection is assumed to be the primary transport process governing the release of radionuclides in the Ancillary Equipment and Soil-Debris Sub-models. Diffusion is assumed to be the primary release process in the Grouted Residual Waste and Cement Sub-models. An unspecified zero-order release process is considered in the Surplus Reactor Block Sub-model due to the lack of information regarding actual processes involved in irradiated graphite leaching. The Surplus Reactor Block, Cement, and SoilDebris Sub-models were compared against analytical solutions (PNNL-11800, Composite Analysis for Low-Level Waste Disposal in the 200 Area Plateau of the Hanford Site). The agreement between results of these analytical solutions and the corresponding waste form models verified their correct implementation in GoldSim. A 1-D modeling abstraction approach for the Grouted Residual Waste and Ancillary Equipment Sub-models was adopted from existing Performance Assessment (PA) models. Despite the simplifications made in these sub-models, they were found to be appropriate representations of the waste forms, similar to what was used in the Waste Management Area C PA model (RPP-ENV-58782, Performance Assessment of Waste Management Area C, Hanford Site, Washington, Rev. 0). A sensitivity analysis was conducted to identify the most influential parameters in each waste form sub-model. Suggestions for considering pH-dependent and redox-dependent release mechanisms are formulated through the development of a conditional constant approach in GoldSim.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

The Role of Evolving Interfacial Substrate Properties on Heterogeneous Cellulose Hydrolysis Kinetics

Interfacial enzyme reactions require formation of an enzyme–substrate complex at the surface of a heterogeneous substrate, but often multiple modes of enzyme binding and types of binding sites complicate analysis of their kinetics. Excess heterogeneous substrate is often used as a justification to model the substrate as unchanging, but using the study of the enzymatic hydrolysis of insoluble cellulose as an example, we argue that reaction rates are dependent on evolving substrate interfacial properties. We hypothesize that the relative abundance of binding sites on cellulose where hydrolysis can occur (productive binding sites) and binding sites where hydrolysis cannot be initiated or is inhibited (nonproductive binding sites) contribute to rate limitations. We show that the initial total number of productive binding sites (the productive binding capacity) determines the magnitude of the initial burst phase of cellulose hydrolysis, while productive binding site depletion explains overall hydrolysis kinetics. Furthermore, we show that irreversibly bound surface enzymes contribute to the depletion of productive binding sites. Our model shows that increasing the ratio of productive to nonproductive binding sites promotes hydrolysis, while maintaining an elevated productive binding capacity throughout conversion is key to preventing hydrolysis slowdown.

09 BIOMASS FUELS↗
DOE OSTI2022

A native IgE in complex with profilin provides insights into allergen recognition and cross-reactivity

Allergies have become a rising health problem, where plentiful substances can trigger IgE-mediated allergies in humans. While profilins are considered minor allergens, these ubiquitous proteins are primary molecules involved in cross-reactivity and pollen-food allergy syndrome. Here we report the first crystal structures of murine Fab/IgE, with its chains naturally paired, in complex with the allergen profilin from Hevea brasiliensis (Hev b 8). The crystallographic models revealed that the IgE’s six complementarity-determining regions (CDRs) interact with the allergen, comprising a rigid paratope-epitope surface of 926 Å 2 , which includes an extensive network of interactions. Interestingly, we also observed previously unreported flexibility at Fab/IgE’s elbow angle, which did not influence the shape of the paratope. The Fab/IgE exhibits a high affinity for Hev b 8, even when using 1 M NaCl in BLI experiments. Finally, based on the encouraging cross-reactivity assays using two mutants of the maize profilin (Zea m 12), this antibody could be a promising tool in IgE engineering for diagnosis and research applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Rapid simulations of hyperspectral near-field images of three-dimensional heterogeneous surfaces – part II

The modeling of the near-field interaction in the scattering-type scanning near-field optical microscope (s-SNOM) is rapidly advancing, although an accurate yet versatile modeling framework that can be easily adapted to various complex situations is still lacking. In this work, we propose a time-efficient numerical scheme in the quasi-electrostatic limit to capture the tip-sample interaction in the near field. This method considers an extended tip geometry, which is a significant advantage compared to the previously reported method based on the point-dipole approximation. Using this formalism, we investigate, among others, nontrivial questions such as uniaxial and biaxial anisotropy in the near-field interaction, the relationship between various experimental parameters (e.g. tip radius, tapping amplitude, etc.), and the tip-dependent spatial resolution. The demonstrated method further sheds light on the understanding of the contrast mechanism in s-SNOM imaging and spectroscopy, while also representing a valuable platform for future quantitative analysis of the experimental observations.

Chen, Xinzhong (ORCID:0000000261033848)↗
DOE OSTI2024

Surface Atmosphere Integrated Field Laboratory (SAIL) (Field Campaign Report)

Mountains are the natural water towers of the world, effectively turning water vapor into readily available fresh water through precipitation, snowpack, and runoff. They contribute disproportionately to precipitation over land, but are under-observed, leading to large gaps in the scientific understanding of convection, extreme precipitation and weather, and interactions between atmospheric circulation, radiation, and land-surface conditions. The mountain hydrometeorology community has repeatedly called for integrated atmospheric and land observations of water and energy budgets in complex terrain that span these scales to establish benchmarks against which scale-dependent models can be further developed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Probabilistic Evaluation of Geomechanical Risks in CO2 Storage: An Exploration of Caprock Integrity Metrics Using a Multilaminate Model

The probabilistic uncertainty assessment of geomechanical risk—specifically, caprock failure—attributable to CO2 injection, as presented in a simplified hypothetical geological model, was the focus of this study. Our approach amalgamates the implementation of a multilaminate model, the creation of a response surface model in conjunction with the Box–Behnken sampling design, the execution of associated numerical modeling experiments, and the utilization of Monte Carlo simulations. Probability distributions to encapsulate the inherent variability (elastic and mechanical properties of the caprock and reservoir) and uncertainty in prediction estimates (vertical displacement, total strain, and F value) were employed. Our findings reveal that the Young modulus of the caprock is a key factor controlling equivalent total strain but is insufficient as a stand-alone indicator of caprock integrity. It is confirmed that the caprock can accommodate significant deformation without failure, if it possesses a low Young’s modulus and high mechanical strength properties, such as the friction angle and uniaxial compressive strength. Similarly, vertical displacement was found to be an unreliable indicator for caprock integrity, as caprock failure can occur across a broad spectrum of vertical displacements, particularly when both the Young modulus and mechanical strength properties have wide ranges. This study introduces the F value as the most dependable indicator for caprock failure, although it is a theoretical attribute (the shortest distance between the Mohr circle and the nearest failure envelope used to measure the sensitivity to failure) and not physically measurable in the field. Deviatoric stress levels were found to vary based on stress regimes, with the maximum levels observed under extensive and compressive stress regimes. In conjunction with the use of the response surface method, this study demonstrates the efficacy of the multilaminate framework and the Mohr–Coulomb constitutive model in providing a simplified, yet effective, probabilistic model of the mechanical behavior of caprock failure, reducing mathematical and computational complexities.

Energy & Fuels↗
DOE OSTI2024

Hemispherically symmetric strategies for stratospheric aerosol injection

Stratospheric aerosol injection (SAI) comes with a wide range of possible design choices, such as the location and timing of the injection. Different stratospheric aerosol injection strategies can yield different climate responses; therefore, understanding the range of possible climate outcomes is crucial to making informed future decisions on SAI, along with the consideration of other factors. Yet, to date, there has been no systematic exploration of a broad range of SAI strategies. This limits the ability to determine which effects are robust across different strategies and which depend on specific injection choices. This study systematically explores how the choice of SAI strategy affects climate responses in one climate model. Here, we introduce four hemispherically symmetric injection strategies, all of which are designed to maintain the same global mean surface temperature: an annual injection at the Equator (EQ), an annual injection of equal amounts of SO 2 at 15° N and 15° S (15N+15S), an annual injection of equal amounts of SO 2 at 30° N and 30° S (30N+30S), and a polar injection strategy that injects equal amounts of SO 2 at 60° N and 60° S only during spring in each hemisphere (60N+60S). We compare these four hemispherically symmetric SAI strategies with a more complex injection strategy that injects different quantities of SO 2 at 30° N, 15° N, 15° S, and 30° S in order to maintain not only the global mean surface temperature but also its large-scale horizontal gradients. All five strategies are simulated using version 2 of the Community Earth System Model with the middle atmosphere version of the Whole Atmosphere Community Climate model, version 6, as the atmospheric component, CESM2(WACCM6-MA), with the global warming scenario, Shared Socioeconomic Pathway (SSP)2-4.5. We find that the choice of SAI strategy affects the spatial distribution of aerosol optical depths, injection efficiency, and various surface climate responses. In addition, injecting in the subtropics produces more global cooling per unit injection, with the EQ and the 60N+60S cases requiring, respectively, 59 % and 50 % more injection than the 30N+30S case to meet the same global mean temperature target. Injecting at higher latitudes results in larger Equator-to-pole temperature gradients. While all five strategies restore Arctic September sea ice, the high-latitude injection strategy is more effective due to the SAI-induced cooling occurring preferentially at higher latitudes. These results suggest trade-offs wherein different strategies appear better or worse, depending on which metrics are deemed important.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Computational Study on the Effect of Thermal Boundary Conditions and Axial Aspect Ratio on Catalytic Oxidative Coupling of Methane

Packed-bed reactor models with a large reactor-diameter-to-particle-diameter ratio (N > 10) are commonly used in catalysis research. Their use is underpinned by strong simplifying assumptions including unidirectional plug flow, negligible radial gradients, small axial diffusion, and isothermal operation. However, the complex exothermic–endothermic nature of reactions like oxidative coupling of methane (OCM) coupled with heat transfer from the catalyst surface and convective heat losses from the reactor wall provides an interesting case study to evaluate packed-bed assumptions and understand axial and radial temperature variations that may impact conversion and selectivity. This computational study investigates the effect of thermal boundary conditions on catalyst performance in laboratory OCM reactors (N = ~13) through computational fluid dynamics (CFD) simulation. Further, the 1-D model typically used for modeling plug flow reactors assumes negligible radial temperature and species concentration gradients. This assumption is evaluated for OCM reactors over a range of geometric axial aspect ratios, defined as the reactor length divided by the diameter, comparing the 1-D with 3-D CFD results. 1-D model predictions of catalyst performance start to deviate from the 3-D model predictions as the axial aspect ratio decreases, specifically as diameter increases at constant length. In the work, the 1-D model assumption is only found to be valid for axial aspect ratios higher than 1.6, a result considered to be relevant for all exothermic catalytic reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗