Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “unknowns”

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 1,063 records · Page 59

X-ray fluorescence and XANES spectroscopy revealed diverse potassium chemistries and colocalization with phosphorus in the ectomycorrhizal fungus Paxillus ammoniavirescens

Ectomycorrhizal (ECM) fungi play a major role in forest ecosystems and managed tree plantations. Particularly, they facilitate mineral weathering and nutrient transfer towards colonized roots. Among nutrients provided by these fungi, potassium (K) has been understudied compared to phosphorus (P) or nitrogen (N). The ECM fungus Paxillus ammoniavirescens is a generalist species that interacts with the root of many trees and can directly transfer K to them, including loblolly pine. However, the forms of K that ECM fungi can store is still unknown. Here, we used synchrotron potassium X-ray fluorescence (XRF) and K-edge X-ray Absorption Near Edge Structure (XANES) spectroscopy on P. ammoniavirescens growing in axenic conditions to investigate the K chemistries accumulating in the center and the edge of the mycelium. We observed that various K forms accumulated in different part of the mycelium, including K-nitrate (KNO 3 ), K-C-O compounds (such as K-tartrate K 2 (C 4 H 4 O 6 ) and K-oxalate (K 2 C 2 O 4 )), K-S and K-P compounds. Saprotrophic fungi have been shown to excrete carboxylic acids, which in turn play a role in soil mineral weathering. Our finding of several K counter-ions to carboxylic acids may suggest that, besides their direct transfer to colonized roots, K ions can also be involved in the production of compounds necessary for sourcing nutrients from their surrounding environment by ECM fungi. Additionally, this work reveals that XANES spectroscopy can be used to identify the various forms of K accumulating in biological systems.

Ectomycorrhizal symbiosis↗

Lignin-enhanced reduction of structural Fe(III) in nontronite: Dual roles of lignin as electron shuttle and donor

Lignin is a major component of plant-derived soil organic matter (SOM) in soils and sediments. Fe-bearing clay minerals are widely distributed in these environments and often co-exist with lignin. While previous studies have reported the electron shuttling and donating roles of certain redox-active SOM in the dissimilatory reduction of structural Fe(III) in Fe-bearing clay minerals, the role of lignin in this process remains unknown. Here we studied this role by incubating an Fe-rich smectite (nontronite NAu-2) with two types of lignin (soluble and insoluble) in the absence and presence of an Fe(III)-reducing bacterium Shewanella putrefaciens CN32 under anaerobic condition. Labile carbon source (i.e., lactate) was added in some experiments as an extra electron donor. The results demonstrated that both soluble and insoluble lignins abiotically reduced structural Fe(III) in NAu-2. The reduction extent was proportional to lignin concentration. After abiotic reaction, lignin served as either electron shuttle or electron donor in the presence of CN32: (1) When lactate was present, lignin served as an electron shuttle to enhance the rate of Fe(III) reduction; (2) When lactate was absent, lignin served as an electron donor for Fe(III) reduction. Although the biotic Fe(III) reduction extents were similar in the presence of either soluble or insoluble lignin, the reduction rates with soluble lignin were always higher than those with insoluble lignin, likely owing to their different electron transfer mechanisms. After interaction with NAu-2 and/or CN32, soluble lignin structure largely remained intact, but with some decreases of humic/fulvic acid-like and protein-like compounds, aromatic functional groups (e.g., C-H, C=O, COOH), and aliphatic/aromatic compounds. An increase of semiquinone-like organic radicals was observed after lignin interaction with NAu-2. These chemical changes of lignin were likely coupled with reduction of structural Fe(III) in nontronite. The nontronite structure remained largely unchanged without much dissolution and mineral transformation. The findings of this study provide insights into the role of lignin in promoting mineral-microbe interactions and have significant implications for coupled Fe and C biogeochemical cycles in soils and sediments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Co-precipitation induces changes to iron and carbon chemistry and spatial distribution at the nanometer scale

Association of organic matter (OM) with mineral phases via co-precipitation is expected to be a widespread process in environments with high OM input and frequent mineral dissolution and re-precipitation. In contrast to surface area-limited adsorption processes, co-precipitation may allow for greater carbon (C) accumulation. However, the potential sub-micrometer scale structural and compositional differences that affect the bioavailability of co-precipitated C are largely unknown. In this work, we used a combination of high-resolution analytical electron microscopy and bulk spectroscopy to probe interactions between a mineral phase (ferrihydrite, nominally Fe 2 O 3 •0.5H 2 O) and organic soil-derived water-extractable OM (WEOM). In co-precipitated WEOM-Fe, nanometer-scale scanning transmission electron microscopy with electron energy loss spectroscopy (STEM-EELS) revealed increased Fe(II) and less Fe aggregation relative to adsorbed WEOM-Fe. Spatially distinct lower- and higher-energy C regions were detected in both adsorbed and co-precipitated WEOM-Fe. In co-precipitates, lower-energy aromatic and/or substituted aromatic C was spatially associated with reduced Fe(II), but higher-energy oxidized C was enriched at the oxidized Fe(III) interface. Therefore, we show that co-precipitation does not constitute a non-specific physical encapsulation of C that only affects Fe chemistry and spatial distribution, but may cause a bi-directional set of reactions that lead to spatial separation and transformation of both Fe and C forms. In particular, we propose that abiotic redox reactions between Fe and C via substituted aromatic groups (e.g., hydroquinones) play a role in creating distinct co-precipitate composition, with potential implications for its mineralization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the dissolution of a borosilicate glass with the use of isotopic tracing – Insights into the mechanism for the long-term dissolution rate

An understanding of the mechanisms responsible for controlling the long-term corrosion rate of nuclear waste glass is paramount if reliable glass dissolution models are to be used to calculate the controlled release of radionuclides from nuclear waste glass under geologic disposal conditions. Additionally, understanding silicate glass dissolution is also important for natural glasses to understand the role from the dissolution of these glasses has on the composition of natural aquatic systems. Two general mechanisms appear to be responsible for the elemental release from glass to the surrounding biosphere - ion exchange (release principally of alkalis) and matrix dissolution (release of the structural components of the glass). Key unknowns related to these mechanisms are the impact of surface layers on the altered glass and how these layers participate in corrosion. To better understand the impact of these layers on glass corrosion, we made two glasses of the same elemental composition, but with selected elements enriched in specific isotopes, e.g. more 10 B than 11 B. We reacted these glasses for 1 y then the solutions from each were exchanged and the chemistry of the solutions and the solids were followed for another 2 y and in some cases 3 y. This allowed us to follow changes in the glasses with respect to these elements and how these elements were transported through the alteration layers with time and, for at least one element, into the glass. Results from these experiments demonstrate that the release behavior of different elements is strongly dependent on their structural role in the glass (e.g. network formers and modifiers), but, more importantly, the role of water transport and subsequent ion exchange has in the long-term dissolution of the glass studied here. In this article, we highlight the behavior of four elements: lithium, sodium, silicon, and boron. Lithium and sodium, network modifying elements, have similar chemistries. The behavior of Li could be tracked in more detail by following both 6 Li and 7 Li; Li penetrates through the gel layer in both directions without hinderance and into and released from the glass deeper than sodium. Silicon, a network forming element, reacts with the silica-rich alteration layer. Boron, a network forming element, does not accumulate in the gel or the pristine glass and has a very sharp elemental profile between the pristine glass and alteration layer. Boron, lithium, and sodium elemental profiles suggest that there is little if any transport control within the alteration layer. The correlation between the Li and Na profiles and that of a water species suggest that the limiting release is a steady state between the diffusion of a water species and matrix dissolution that results in a low release rate for alkali and this provides a steady driving force for the long-term dissolution of glass.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Reduction of selenite by bacterial exudates

Bacterial reduction of Se(IV) is a significant component of the global selenium cycle, and hence affects the fate and transport of selenium in both natural and contaminated environments. However, it is unknown whether bacterially-produced exudates are capable of reducing selenium independent of bacterial cells. In this study, we measured the rate and extent of Se(IV) reduction by exudates from three bacterial species, and we determined the importance of exudate sulfhydryl sites by conducting parallel experiments after treatment of the exudates with a sulfhydryl-specific blocking molecule. Here, we also conducted experiments with whole cell biomass for each of the three bacterial species to determine the importance of exudate-only reduction relative to whole biomass-promoted reduction. Under our experimental conditions, exudates from Bacillus subtilis and Pseudomonas putida remove Se(IV) from solution after an approximately 20-24 h lag period, but exudates from Shewanella oneidensis are ineffective at Se(IV) removal, except when the cells are grown in the presence of Se(IV). For both B. subtilis and P. putida , pretreatment blocking of the sulfhydryl sites on the exudate molecules dramatically decreases the rate and extent of Se(IV) removal, strongly suggesting that sulfhydryl groups on the exudate molecules play a key role in the Se(IV) reduction mechanism. The behavior of S. oneidensis exudates indicates that up-regulation of the Se(IV) reduction mechanism can occur in response to the Se content of the medium in which the cells grow. Our results demonstrate the capacity of some bacterial exudates to reduce Se(IV), and suggest that the activity of bacterial exudates should be accounted for when modeling selenium cycling in natural and engineered environments.

59 BASIC BIOLOGICAL SCIENCES↗

Comparison of the earliest NC and CC planetesimals: Evidence from ungrouped iron meteorites

Isotope anomalies in meteorites reveal a fundamental dichotomy between Non-Carbonaceous- (NC) and Carbonaceous-type (CC) planetary bodies. Until now, this dichotomy is established for the major meteorite groups, representing about 36 distinct parent bodies. Ungrouped meteorites represent an even larger number of additional parent bodies, but whether they conform to the overall NC-CC dichotomy is unknown. Here, the genetics and chronology of 26 ungrouped iron meteorites are considered through nucleosynthetic Mo and radiogenic W isotopic compositions. Secondary cosmic ray-induced modifications of these isotope compositions are corrected using Pt isotope measurements on the same samples. We find that all of the ungrouped irons have Mo isotope anomalies within the range of the major meteorite groups and confirm the NC-CC dichotomy for Mo, where NC and CC meteorites define two distinct, subparallel s-process mixing lines. All ungrouped NC irons fall on the NC-line, which is now precisely defined for 41 distinct parent bodies. The ungrouped CC irons show scatter around the CC-line indicative of small r-process Mo heterogeneities among these samples. These r-process Mo isotope variations correlate with O isotope anomalies, most likely reflecting mixing of CI chondrite-like matrix, chondrule precursors, and Ca-Al-rich inclusions. This implies that CC iron meteorite parent bodies accreted the same nebular components as the later-formed carbonaceous chondrites. The Hf-W model ages of core formation for the ungrouped irons overlap with those of the iron meteorite groups from each reservoir and reveal a narrow age peak at ∼3.3 Ma after Ca-Al-rich inclusions for the CC irons. By contrast, the NC irons display more variable ages, including younger ages indicative of impact-induced melting events, which seem absent among the CC irons. This is attributed to the more fragile and porous nature of the CC bodies, making impact-induced melting on their surfaces difficult. The chemical characteristics of all iron meteorites together reveal slightly more oxidizing conditions during core formation for CC compared to NC irons. More strikingly, strong depletions in moderately volatile elements, typical of many iron meteorite parent bodies, predominantly occur among CC irons, for reasons that remain unclear at present.

Chondrite components↗

The solubility of ErPO 4 and Er speciation in hydrothermal fluids at varying pH and salinity between 350 and 450 °C

The rare earth elements (REE) are important for the green-energy transition and can be incorporated into the REE phosphates, such as xenotime-(Y), which also hosts heavy REE (Tb– Lu). Xenotime-(Y) is a common accessory mineral in metamorphic rocks and a range of mineral deposits where it controls the mobility of heavy REE, however, the impact of high temperature aqueous fluids on the behavior of heavy REE is largely unknown. Thermodynamic modeling can be utilized as a tool to predict the mobility of REE in hydrothermal aqueous fluids, but must be supported by accurate experimental data. Here, we measured the solubility of endmember synthetic xenotime-structured ErPO4 in NaCl-HCl-NaOH-bearing aqueous solutions at 350 °C and water vapor saturation pressure, at 400 and at 450 °C and 500 bar using batch-type Inconel reactors. Erbium speciation was investigated as a function of pH from 2.8 to 8, where Er chloride species are predominant at acidic conditions (pH <3) and Er hydroxyl complexes are predominant at near- neutral to alkaline conditions (pH >3). At pH 7–9, the measured ErPO 4 solubility (-9.8 to -7.5 log m Er ) is up to 2.5 orders of magnitude lower than thermodynamic predictions (-9.4 to -6.7 log m Er ) using existing thermodynamic databases. At pH 2–3, the predicted ErPO 4 solubility is ~0.5 orders of magnitude higher at 350 °C and ~1 order of magnitude lower at 450 °C compared to experimentally measured Er concentrations. The thermodynamic properties of aqueous Er species were therefore revised in this study. The partial molal Gibbs energy of formation (Δ f G 0 T,P ) for aqueous Er hydroxyl and chloride species are optimized using GEMSFITS and the logarithmic formation constants (logβ n ( Cl,OH) ) were derived at each experimental temperature and pressure. The updated thermodynamic properties for Er hydroxyl species (Er(OH) +2 , Er(OH) 2 + , and Er(OH) 3 0 ) show that their stability shifts to more acidic conditions at and below 400 °C. The Er chloride species (ErCl +2 and ErCl 2 + ) show increased stability compared to Er hydroxyl species at temperatures of 450 °C and 0.01 mol/kg NaCl. The updated thermodynamic properties are implemented into the GEM-Selektor modeling package to investigate the mobility of Er in saline hydrothermal fluids in equilibrium with alkaline rocks. Importantly, the updated properties for Er hydroxyl species result in low Er solubility at rock equilibrated pH conditions due to an expanded hydroxyl predominance zone, but lower aqueous complex stability overall, whereas previous models suggest greater stability for aqueous Er species. Furthermore, ErPO 4 solubility increases with decreasing temperature due to the deprotonation of HCl, which increases the acidity of hydrothermal fluids and the availability of Cl - to complex with the REE. These simulations highlight how fluid-rock reaction and temperature affect the mobility of REE in hydrothermal ore-forming systems.

58 GEOSCIENCES↗

Multi-site evaluation of stratified and balanced sampling of soil organic carbon stocks in agricultural fields

Estimating soil organic carbon (SOC) stocks in agricultural fields is essential for environmental and agronomic research, management, and policy. Stratified sampling is a classic strategy for estimating mean soil properties, and has recently been codified in SOC monitoring protocols. However, for the specific task of estimating the SOC stock of an agricultural field, concrete guidance is needed for which covariates to stratify on and how much stratification can improve estimation efficiency. It is also unknown how stratified sampling of SOC stocks compares to modern alternatives, notably doubly balanced sampling. To address these gaps, we collected high-density (average of 7 samples ha -1 ) and deep (average of 75 cm) measurements of SOC stocks at eight commercial fields under maize-soybean production in two US Midwestern states. We combined these measurements with a Bayesian geostatistical model to evaluate stratified and balanced sampling strategies that use a set of readily-available geographic, topographic, spectroscopic, and soil survey data. We examined the number of samples needed to achieve a given level of SOC stock estimation accuracy. While stratified sampling using these variables enables an average sample size reduction of 17% (95% CI, 11% to 23%) compared to simple random sampling, doubly balanced sampling is consistently more efficient, reducing sample sizes by 32% (95% CI, 25% to 37%). The data most important to these efficiency gains are a remotely-sensed SOC index, SSURGO estimates of SOC stocks, and the topographic wetness index. We conclude that in order to meet the urgent challenge of climate change, SOC stocks in agricultural fields could be more efficiently estimated by taking advantage of this readily-available data, especially with doubly balanced sampling.

54 ENVIRONMENTAL SCIENCES↗

CO 2 storage site characterization using ensemble-based approaches with deep generative models

Estimating spatially distributed properties such as permeability from available sparse measurements is a great challenge in efficient subsurface CO 2 storage operations. In this paper, a deep generative model that can accurately capture complex subsurface structure is tested with an ensemble-based inversion method for accurate and accelerated characterization of CO 2 storage sites. We chose Wasserstein Generative Adversarial Network with Gradient Penalty (WGAN-GP) for its realistic reservoir property representation and Ensemble Smoother with Multiple Data Assimilation (ES-MDA) for its robust data fitting and uncertainty quantification capability. WGAN-GP are trained to generate high-dimensional permeability fields from a low-dimensional latent space and ES-MDA then updates the latent variables by assimilating available measurements. Several subsurface site characterization examples including Gaussian, channelized, and fractured reservoirs are used to evaluate the accuracy and computational efficiency of the proposed method and the main features of the unknown permeability fields are characterized accurately with reliable uncertainty quantification. Furthermore, the estimation performance is compared with a widely-used variational, i.e., optimization-based, inversion approach, and the proposed approach outperforms the variational inversion method in several benchmark cases. We explain such superior performance by visualizing the objective function in the latent space: because of nonlinear and aggressive dimension reduction via generative modeling, the objective function surface becomes extremely complex while the ensemble approximation can smooth out the multi-modal surface during the minimization. This suggests that the ensemble-based approach works well over the variational approach when combined with deep generative models at the cost of forward model runs unless convergence-ensuring modifications are implemented in the variational inversion.

42 ENGINEERING↗

Exploring density and strength variations in asteroid 16 Psyche’s composition with 3D hydrocode modeling of its deepest impact structure

Asteroid 16 Psyche is the largest metallic Main Belt Asteroid and is the subject of a forthcoming NASA mission. The composition of Psyche is still unknown and subject of recent debate. In particular, how much porosity is within Psyche, along with how much of Psyche consists of non-metallic versus metallic materials, are central questions to the issue of Psyche’s composition. If Psyche is indeed predominantly composed of metallic materials, it would need to have considerable porosity (~ 30%–50%) for a composition consistent with its expected bulk density (~ 3.7–4.1 g/cm). In this work, we vary the density and strength of Psyche by including uniform and layered fields of pseudo-microporosity, in addition to investigating the presence of macroscopic voids, i.e., spaces larger than the size of the simulation’s mesh cells, in rubble-pile configurations. Further, all configurations result in bulk densities within the uncertainties of measured values, however the strength of Psyche and the distribution of pseudo-pores are varied. Through 3D computational models of Psyche’s deepest impact structure, we show that Psyche’s composition is unlikely to contain only pseudo-microporosity. Rather, rubble pile structures, which include macroscopic voids, are shown to match the crater’s measured aspect ratio better than simulations of structures that included only pseudo-microporosity.

3D↗

Deep Learning Explicit Differentiable Predictive Control Laws for Buildings

We present a differentiable predictive control (DPC) methodology for learning constrained control laws for unknown nonlinear systems. DPC poses an approximate solution to multiparametric programming problems emerging from explicit nonlinear model predictive control (MPC). Contrary to approximate MPC, DPC does not require supervision by an expert controller. Instead, a system dynamics model is learned from a small dataset of recorded observations of the perturbed system's dynamics and the control law is optimized offline by interaction with the learned system model. The DPC method is based on two sequential steps, i) system identification using a constrained neural state-space model, and ii) optimization of an explicit control law parametrized by another neural network in closed-loop simulation with the identified neural state-space model. The combination of a differentiable closed-loop system and penalty methods for constraint handling of system outputs and inputs allows us to optimize the control law's parameters directly by backpropagating economic MPC loss through the learned system model. By incorporating domain knowledge and leveraging established techniques from optimal control, our method leverages deep neural networks as nonlinear function approximators for system identification and control while avoiding concomitant costs of intractably large datasets, and computationally expensive over-parametrized models. The scalability, data efficiency, and constrained optimal control capability of the proposed DPC method are demonstrated in simulation using a multi-zone building emulator.

Drgona, Jan↗

A review of active probing-based system identification techniques with applications in power systems

System identification is becoming a relevant research area for numerous applications in power grids due to the increasing complexity of the system. A paradigm shift in power system infrastructure driven by renewable energy resources, controllable loads, and new power electronics technologies have given rise to new challenges in power systems operation and control. This increased complexity of power systems and unavailability of physics-based models of most inverter-based resources requires that traditional modeling of power systems approaches be complemented by system identification-based black or grey box modeling techniques. In light of the importance of system identification and little attention paid to the applications of these techniques in transforming power systems, this paper provides a comprehensive review of active probing-based system identification methods in the context of power system applications. It reviews applications of both linear and nonlinear system identification with discussion on their potential and key takeaways. To motivate a further practical use in power systems, the paper provides an example of system identification, to develop a state space model of an unknown plant with step-by-step details. The paper highlights the advantages of using modern power electronics-based sources in the identification process and discusses the emerging research directions for future.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Physical solvents and techno-economic analysis for pre-combustion CO 2 capture: A review

In this manuscript we review the use, development, and economic performance of physical solvents for pre-combustion CO 2 capture from high pressure H 2 rich syngas streams. Commercially available physical solvents are presented, followed by an assessment of the ideal properties that are important for development of novel solvents for CO 2 capture from high-pressure syngas streams. To compare the technical and economic performance of traditional and novel physical solvents, a review of the methods, assumptions and models used in techno-economic analysis (TEA) studies was conducted. It was found that, although some novel solvents show promising technical performance in the laboratory (e.g., high CO 2 absorption capacity and low vapor pressure), other issues (e.g., solvent viscosity and cost) may limit their industrial applications. Process simulations were useful tools for modeling the technical performance of processes using traditional and novel solvents. However, model predictions are most reliable when the methods and correlations used to develop the process simulation are validated with representative experimental data, in particular highly accurate baseline models are required for fair comparison among physical solvents. The key inputs and assumptions in pre-combustion CO 2 capture TEAs have also been summarized. Some studies showed that the promising technical performance of novel physical solvents can be offset by the high and often unknown costs of these solvents. Future development of novel physical solvents for pre-combustion CO 2 capture will benefit from more studies that conduct in-depth techno-economic analysis, specifically with validated process simulations and transparent economic models.

42 ENGINEERING↗

Modeling the performance and faradaic efficiency of solid oxide electrolysis cells using doped barium zirconate perovskite electrolytes

Y-doped BaZrO 3 (BaZr 1 x Y x O 3 δ , or “BZY”), a proton-conducting ceramic featuring high bulk conductivity and good chemical stability, is a promising electrolyte material for solid oxide electrolysis cells. Further doping with Ce and/or Yb (creating materials “BCZY” and “BCZYYb”) can improve conductivity and sintering properties, but at significant penalty to cells’ faradaic efficiency (FE). Studies have proposed that reduction of lattice Ce can occur in the hydrogen electrode, which consumes some hydrogen produced by the hydrogen evolution reaction, leading to decreased FE. Despite studies suggesting this phenomenon, the mechanism is largely unknown. In this work, we developed a multiphysics model to study the transport of multiple defect species and the performance of BZY, BCZY, and BCZYYb, capturing the tradeoff between enhanced performance at the cost of FE for BCZY and BCZYYb electrolytes compared to BZY. We also found that increasing the water content of the anode gas supply lowers the current output of the cell but results in better FE. The model, which uses several parameters previously unavailable in the literature, was validated to experiments varying temperature, steam water content, and electrolyte material, as well as two performance metrics (performance curves and FE). Results verify and explain observed trends, informing future work on Ce-doped BZY electrolytes.

08 HYDROGEN↗

Plasma production and composition from hypervelocity impacts on solar cell cover glass

Satellites are subjected to detrimental effects from the space environment, including impacts by orbital debris and meteoroids that can cause mechanical or electrical damage. While mechanical damage is well studied, electrical effects resulting from hypervelocity impact (HVI) remain poorly understood. In HVIs, where the projectile speed exceeds the speed of sound in the target material, the impact energy ionizes the material near the surface, creating a plasma that can cause radio frequency emission. This emission is the main source of electrical anomalies in satellites, but the mechanism behind its creation is still unknown. Ground-based experiments have been used to characterize HVI plasma, including the empirical power law relations that describe the impact charge produced as a function of impactor mass and speed. Quantifying these power law relations gives us a better understanding of the species formed in the impact which can eventually lead to plasma evolution models and help us first, understand the mechanism behind radio frequency emission and second, ways to prevent its damage to the satellites. In this paper, we interpret the results of the HVI experiments on solar cell coverglass targets and identify the potential chemical or molecular reactions that form the impact ion species. This is done by looking simultaneously at multiple data sets for different species within the same impact configuration, and using a new technique for quantifying the power law relations for each emitted species. We use time of flight analysis (TOF) to detect the charged particle compositions and determine their initial speeds. The results suggest that the charge production power law is species specific, with different formation mechanisms involved in a single HVI event, some of which are associated with dissociation and some with ionization of a molecule. Furthermore, our results also show that electrons have different power laws that depend on the spacecraft surface charge. TOF measurements of ions from impact events on solar cell cover glass shows the presence of SiO + and SiO – 2 as the dominant species present in plasmas on positively and negatively charged targets, respectively, as well as the formation of FeO + .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the convexity of phase-field fracture formulations: Analytical study and comparison of various degradation functions

Efficient and accurate fracture modeling is of great importance in applications where catastrophic outcomes under extreme scenarios are possible. The phase-field (PF) approach to fracture received significant attention over the past decade, due to its capability to capture complicated fracture patterns (e.g., crack merging and branching). Specifically, crack initiation and propagation are modeled via minimization of the total energy functional, which is regularized with the aid of a phase field. Despite the promising results and modeling capabilities of the PF method in many applications, the solution of fracture problems remains computationally challenging mainly due to the non-convexity of the total energy functional with respect to the combined unknown (phase field and displacement) fields. Understanding the effects of their coupling on convexity is crucial in order to address frequently encountered hurdles in fracture modeling (e.g., inefficient solvers and non-physical crack nucleation). In this paper, we develop convexity criteria for a wide class of PF fracture formulations. For this class of formulations, the second variation of the total energy functional is expressed in terms of Hessian matrices (evaluated at individual material points). Depending on the choice of geometric crack functions and degradation functions, we classify the formulations into three categories and analytically study each one separately. To study the sign of the second variation, we derive inequalities which are satisfied at material points when the Hessian matrix is locally positive semi-definite. These inequalities provide objective criteria for comparing degradation functions. Finally, the applicability of the proposed convexity criteria is demonstrated in the context of a one-dimensional problem, solved using a conventional monolithic solver.

97 MATHEMATICS AND COMPUTING↗

Determining critical overlap concentration of polyethylene oxide to support excipient safety assessment of opioid products

Intravenous administration of abuse-deterrent opioid products poses high safety risks, in part due to the presence of high molecular weight polymeric excipients. Previous in vivo studies in animal models have shown that the higher molecular weight (Mw) polymeric excipients like polyethylene oxide (PEO) were directly linked to such adverse responses as intravenous hemolysis and kidney damage. PEO polymers have been widely used in abuse-deterrent formulations (ADF) of opioid products, adding to concerns over the general safety of the opioid category due to the unknown safety risk from abuse via unintended routes. Here, the current study focused on the determination of the critical overlap concentration (c*) at various PEO molecular weights to aid in explaining differences in observed adverse responses from previous animal studies on the intravenous administration of PEO solutions. Adverse in vivo responses may be related to the viscoelastic regime of the polymer solution, which depends not only on Mw but also on concentration. Having a localized polymer concentration in the blood above the c*, i.e., the transition from the dilute to semi-dilute entangled viscoelastic regime, may influence the flow behavior and interactions of cells in the blood. The relationship of c* to this combination of physical, chemical, and rheological effects is a possible driving force behind adverse in vivo responses.

60 APPLIED LIFE SCIENCES↗

Creep anisotropy modeling and uncertainty quantification of an additively manufactured Ni-based superalloy

The advantages offered by additive manufacturing over traditional processes has driven a great deal of industrial and academic interest in recent years. However, the process is relatively new and requires additional investigation to become sufficiently mature for wide scale industrial adoption. Electron beam melting powder bed fusion is one technology that has shown promise for fabricating high temperature resistant materials such as nickel based superalloys. The resulting microstructures typically exhibit a strong fiber texture in the build direction giving rise to anisotropic time-dependent deformation behavior. In order to accelerate the qualification of these materials for industrial adoption accurate numerical models are needed for simulating their behavior. In this work a crystal plasticity model including non-Schmid effects is presented for capturing creep anisotropy observed in additively manufactured IN738LC. The model is calibrated via a probabilistic framework where model parameters are treated as random variables. An iterative sequential design strategy is utilized to efficiently identify the probability density of the unknown model parameters. As a case study the model is utilized to investigate the behavior of randomly oriented equiaxed grain clusters sometimes observed embedded in the additively manufactured columnar structure. A synthetic realization is simulated and uncertainty is propagated through to the full-field response. Results indicate that these features are the source of significant creep relaxation and strain accumulation which partially explains observed grain boundary decohesion at these locations.

36 MATERIALS SCIENCE↗