Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “species dynamics”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Space life support engineering program

A comprehensive study to develop software to simulate the dynamic operation of water reclamation systems in long-term closed-loop life support systems is being carried out as part of an overall program for the design of systems for a moon station or a Mars voyage. This project is being done in parallel with a similar effort in the Department of Chemistry to develop durable accurate low-cost sensors for monitoring of trace chemical and biological species in recycled water supplies. Aspen-Plus software is being used on a group of high-performance work stations to develop the steady state descriptions for a number of existing technologies. Following completion, a dynamic simulation package will be developed for determining the response of such systems to changes in the metabolic needs of the crew and to upsets in system hardware performance.

Seagrave, Richard C.↗
DOE OSTI2023

Cost-constrained adaptive simulations of transient spray combustion in a gas turbine combustor

Predictive high-fidelity simulations of turbulent spray combustion must capture the combined effects of complex chemistry, multiphase evaporating flow and spray-flame interactions to achieve physical accuracy. Finite-rate chemistry (FRC) combined with a realistic chemical mechanism is a combustion model well-suited for this purpose, but has a high computational cost due to the large number and stiffness of transported chemical species. In contrast, flamelet-based models achieve lower cost by transporting a small number of quantities of reduced stiffness, but assumptions regarding local flame topology, boundary conditions and inter-phase coupling limit their physical accuracy. Recently, the Pareto-efficient combustion (PEC) framework was developed to dynamically assign combustion models based on local cost and accuracy metrics in gas-phase reacting flows. In this work, we extend this PEC framework to spray combustion through the rigorous analysis of the multiphase coupling terms in the governing equations. The derivation shows that spray evaporation causes errors in the prediction of species mass fractions for flamelet-based models due to the sensitivity of the local thermo-chemical state to changes in composition caused by fuel vaporization across combustion regimes present in practical spray combustion devices. Sub-model assignment is formulated as a multiple-choice knapsack problem, where computational cost is directly controlled through the fraction of the domain assigned to the FRC sub-model. The extended PEC formulation is applied to the simulation of a realistic rich-quench-lean gas turbine combustor at steady-state conditions, as well as transient operation resulting in lean blow-out (LBO). Analysis of transient simulations during LBO demonstrates the extended PEC formulation’s capacity to dynamically adapt to changing conditions within the combustor. Transient combustor dynamics are shown to approach convergence with limited increases in computational cost, while retaining substantial computational cost reduction compared to monolithic FRC simulations. Through PEC simulations with increasing fractions of the domain assigned to FRC, monolithic flamelet simulations are shown to over-predict flame stability during LBO. Finally, the extended PEC formulation is thus shown to overcome deficiencies of monolithic models by controlling modeling error for multiphase combustion modeling.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Computational fluid dynamics analysis of char conversion in Sandia’s pressurized entrained flow reactor

Design and analysis of practical reactors utilizing solid feedstocks rely on reaction rate parameters that are typically generated in lab-scale reactors. Evaluation of the reaction rate information often relies on assumptions of uniform temperature, velocity, and species distributions in the reactor, in lieu of detailed measurements that provide local information. This assumption might be a source of substantial error, since reactor designs can impose significant inhomogeneities, leading to data misinterpretation. Spatially resolved reactor simulations help understand the key processes within the reactor and support the identification of severe variations of temperature, velocity, and species distributions. In this work, Sandia’s pressurized entrained flow reactor is modeled to identify inhomogeneities in the reaction zone. Tracer particles are tracked through the reactor to estimate the residence times and burnout ratio of introduced coal char particles in gasifying environments. The results reveal a complex mixing environment for the cool gas and particles entering the reactor along the centerline and the main high-speed hot gas reactor flow. Furthermore, the computational fluid dynamics (CFD) results show that flow asymmetries are introduced through the use of a horizontal gas pre-heating section that connects to the vertical reactor tube. Computed particle temperatures and residence times in the reactor differ substantially from the idealized plug flow conditions typically evoked in interpreting experimental measurements. Furthermore, experimental measurements and CFD analysis of heat flow through porous refractory insulation suggest that for the investigated conditions (1350 °C, <20 atm), the thermal conductivity of the insulation does not increase substantially with increasing pressure.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Role of Surface Features on the Initial Dissolution of CH 3 NH 3 PbI 3 Perovskite in Liquid Water: An Ab Initio Molecular Dynamics Study

In this study, the degradation of CH 3 NH 3 PbI 3 (MAPbI 3 ) hybrid organic inorganic perovskite (HOIP) by water has been the major issue hampering its use in commercial perovskites solar cells (PSCs) as MAPbI 3 HOIP have been known to easily degrade in the presence of water molecules. However, even though there have been numerous studies investigating this phenomenon, there is still no consensus on the mechanisms of initial stages of dissolution. Here, we attempt to consolidate differing mechanistic interpretations previously reported in the literature through the use of the first-principles constrained ab-initio molecular dynamics (AIMD) to study the mechanisms, kinetics, and thermodynamics that accompany the degradation of MAPbI 3 HOIP in liquid water. We consider not only the dissolution of the species found on the pristine MAPbI 3 HOIP surfaces, but also the dissolution from defect sites to imitate the successive degradation steps, and propose a sequence of events in the initial phase of MAPbI 3 HOIP dissolution. By comparing the dissolution free energy barrier between surface species of different surficial types, we find that the dominant dissolution mechanisms of surface species varies widely based on the specific surface features adjacent to the dissolving ion. The high sensitivity of dissolution mechanism to surface features has contributed to the many dissolution mechanisms proposed in the literature. In contrast, the dissolution free energy barriers are mainly determined by the dissolving species rather than the type of surfaces, and the type of surfaces the ions are dissolving from are inconsequential toward the dissolution free energy barrier. However, the presence of surface defects such as vacancy sites are found to significantly lower the dissolution free energy barriers. Based on the estimated dissolution free energy barriers from different types of surfaces that we investigated in this study, we proposed that the dissolution of MAPbI 3 HOIP in liquid water originates from surface defect sites that propagate laterally along the surface layer of the MAPbI 3 HOIP crystal.

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

A comprehensive modeling method for proton exchange membrane electrolyzer development

Hydrogen attracts significant interests as an effective energy carrier that can be derived from renewable sources. Hydrogen production using a proton-exchange membrane (PEM) electrolyzer can efficiently convert renewable power via water splitting in wide scales—from large, centralized generation to on-site production. Mathematical models with multiple scales and fidelities facilitate the continuing improvements of PEM electrolyzer development to improve performance, cost, and reliability. The model scopes and methods are presented in this paper, which also introduces a comprehensive PEM electrolysis modeling tool based on computational fluid dynamics (CFD) software, ANSYS/Fluent. The modeling tool incorporates electrochemical model of a PEM electrolysis cell to simulate the performance of coupled thermal-fluid, species transport, and electrochemical processes in a product-scale cell or stack by leveraging the powerful meshing generation and CFD solver of ANSYS/Fluent. The thermal-fluid modeling includes liquid water/gas two-phase flow and simulates a PEM electrolysis cell by using Fluent user-defined functions as add-on modules accounting for PEM-specific species transport and electrochemical processes. The modeling outcomes expediate PEM electrolyzer scaling up from basic material development and laboratory testing.

08 HYDROGEN↗
DOE OSTI2025

Distinct Gas-Particle Partitioning and Viscosity Characteristics of Secondary Organic Aerosols Derived from α-Pinene versus Ocimene

Secondary organic aerosols (SOA) have complex, multicomponent composition that controls particle viscosity and gas-particle partitioning, key factors to their atmospheric evolution. This study investigates the chemical composition, volatility and viscosity of SOA formed by ozonolysis of cyclic α-pinene (PSOA) and acyclic ocimene (OSOA) monoterpenes. Using Temperature-Programmed Desorption combined with Direct Analysis in Real-Time ionization and High-Resolution Mass Spectrometry, we determined the molecular composition and saturation mass concentration of individual SOA constituents. These data enabled gas-particle partitioning and viscosity estimates under varied atmospheric conditions. PSOA, composed of higher molecular weight and less oxidized species, shows higher condensability and viscosity under high total organic mass (tOM) loadings. Here, in contrast, OSOA, consisting of more oxidized, lower molecular weight species, exhibits greater sensitivity to tOM, with viscosity increasing significantly upon dilution. Poke-flow experiments support this trend, indicating that OSOA undergoes more dynamic compositional and phase changes during atmospheric aging. These observations reveal distinct dynamic trends in the atmospheric transformations and reactivity of SOA from cyclic and acyclic monoterpenes, with the latter showing greater compositional changes during aging that alter viscosity and diffusion. This highlights the importance of incorporating such dynamic transformations into atmospheric models to improve predictions of SOA atmospheric loadings, lifetimes, and impacts.

cyclic and acyclic monoterpenes↗
DOE OSTI2025

Operando Unveiling of Hydrogen Spillover Mechanisms on Tungsten Oxide Surfaces

Hydrogen spillover is an important process in catalytic hydrogenation reactions, facilitating H 2 activation and modulating surface chemistry of reducible oxide catalysts. This study focuses on the operando unveiling of platinum-induced hydrogen spillover on monoclinic tungsten trioxide (γ-WO 3 ), employing ambient pressure X-ray photoelectron spectroscopy, density functional theory calculations and microkinetic modeling to investigate the dynamic evolution of surface states at varied temperatures. At room temperature, hydrogen spillover results in the formation of W 5+ and hydrogen intermediates (hydroxyl species and adsorbed water), facilitated by Pt metal clusters. With increasing temperature, water desorption, reverse hydrogen spillover and surface-to-bulk diffusion of hydrogen atoms compete with each other, leading initially to reoxidation and then further reduction of W atoms in the near-surface. The combined experimental results and simulations provide a comprehensive understanding of the mechanisms underlying hydrogen interaction with reducible metal oxides, lending insights of relevance to the design of enhanced hydrogenation catalysts.

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

Unexpected electrochemical behavior of an anolyte redoxmer in flow battery electrolytes: solvating cations help to fight against the thermodynamic–kinetic dilemma

Redoxmers are redox-active molecules that can store energy in electrolytes for redox flow batteries (RFBs), and their electrochemical properties are significantly affected by the choice of supporting electrolytes. Herein, we use 2,1,3-benzothiadiazole (BzNSN) as a model system to scrutinize the supporting electrolyte impact. By systemically varying the components of supporting salts, BzNSN not only shows substantial redox potential shifts but also exhibits varying electrochemical stabilities. Specifically, changing the size of cations can effectively alter the coordination between the supporting salt and BzNSN species. From Li + , Na + , K + , to NEt 4 + , the redox potential of BzNSN shifts negatively, from -1.63 V to -1.82 V vs. Ag/Ag + . Molecular dynamics and density functional theory simulations revealed that smaller cations, like Li + , are closer to the charged BzNSN when coordinated, implying stronger coordination, while larger cations, like K + and NEt 4 + , are farther away. Interestingly, the large cation electrolytes also lead to much improved electrochemical stability, evidenced by the extraordinarily enhanced kinetic lifetime from electron paramagnetic resonance measurement. This study demonstrates the first example of tuning an anolyte redoxmer toward a concurrent improvement of lowered redox potentials AND enhanced calendar lives via solvation means, which is usually constrained by the thermodynamic–kinetic relation.

2,1,3-benzothiadiazole↗
DOE OSTI2025

Genetic control of morphological transitions in a coacervating protein template

Nature routinely exploits liquid–liquid phase separation (LLPS) of proteins to control the assembly and mineralization of hybrid materials. Here, we show that fusion of the Car9 silica-binding peptide to an elastin-like polypeptide (ELP) yields temperature- and sequence-programmable soft matter templates for the synthesis of silicified architectures ranging in size from nanometers to micrometers. Specifically, we demonstrate unprecedented control over the diameter of silica nanoparticles (SiNP) in the 30–60 nm range with 4 nm precision, show that a single arginine residue (R4) in the Car9 sequence underpins the transition from micelles to proteinosomes, and find that substitutions in other basic residues modulate electrostatic repulsion and solvation to enable access to kinetically trapped species. These structures, which include interconnected micelles, small (∼200 nm) and large (>5 µm) vesicles, are readily visualized by SEM imaging following silicification. Molecular dynamics (MD) simulations and AlphaFold predictions reveal that mutations in positively charged residues alter interfacial packing, hydration, and conformational freedom of the silica-binding segments. Overall, our results establish sequence and thermal energy as synergistic levers for morphological control across length scales using solid-binding ELPs and establish mineralization as a powerful tool to visualize the structure of dynamic soft matter assemblies.

hierarchy↗
DOE OSTI2022

Mitigation of Molten Salt Corrosion

The desire to increase the peak operational temperature of Concentrated Solar Power (CSP) plants to lower the levelized cost of CSP-derived electricity, has led to consideration of an earth-abundant, low-cost, thermally-stable, molten MgCl 2 -KCl-NaCl (40-40-20 mole%) salt for high-temperature thermal energy storage and heat transfer. The reliable, long-term containment of such a molten salt is challenging, however, due to the corrosion of cost-effective, high-temperature structural metal alloys in this salt (particularly if this salt is contaminated with oxygen-bearing species). This project has been focused on: i) development of equipment for evaluating corrosion kinetics in molten chlorides under well-controlled thermal, fluid dynamic, and atmospheric conditions at CSP-relevant temperatures, ii) use of such equipment to evaluate the influence of forced convection on corrosion kinetics (to determine the role of liquid phase diffusion on such corrosion), and iii) identifying a corrosion-mitigation strategy that can be industrially applied to achieve a projected annual corrosion ≤30 μm at CSP-relevant temperatures. The overall goals of this project have been: i) to understand the rate-limiting step(s) and associated kinetic laws governing corrosion in molten chlorides at 750°C, ii) to determine key operational parameters affecting the corrosion kinetics in molten salts, and iii) to demonstrate the tailoring of such parameters to appreciably reduce such corrosion.

14 SOLAR ENERGY↗
DOE OSTI2025

Advancing Multiscale Simulation of Plasma-Surface Interfaces

We report the development of an atomistic-informed, surface-state-dependent predictive model for particle exchange in a carbon-tungsten plasma-surface interface. The predictive model uses machine learning (ML) techniques to learn the energy and angular distributions for particle exchange and rate functions for surface state evolution from molecular dynamics simulations of cumulative bombardment of tungsten by energetic carbon ions. Each predictive component is sensitive to the energy and trajectory of incident plasma species and the surface state. The surface state is represented by a set of surface state descriptors, which were derived from the atomistic surface state for each independent carbon bombardment event. These descriptors are representative of the composition and degree of amorphization of the outermost angstrom of surface material and were chosen to optimize predictive performance for particle exchange at the interface. The distributions for particle exchange (reflection/sputtering) are demonstrated to vary with each surface state descriptor, motivating the development of surface-state-dependent particle exchange models for plasma simulations. The performance of various ML methods was compared, including polynomial quantile regression, artificial neural networks, k-nearest neighbors, and random forest algorithms, with polynomial regression performing the best for interpolation and extrapolation of learned relationships. In addition to the particle exchange model, a neutral network was developed and used to identify data sufficiency throughout surface descriptor space, which will enable real-time feedback during future data production to ensure data is produced where it is most needed, and we provide commentary on improvements to the data production workflow for future endeavors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Library Screening, In Vivo Confirmation, and Structural and Bioinformatic Analysis of Pentapeptide Sequences as Substrates for Protein Farnesyltransferase

Protein farnesylation is a post-translational modification where a 15-carbon farnesyl isoprenoid is appended to the C-terminal end of a protein by farnesyltransferase (FTase). This process often causes proteins to associate with the membrane and participate in signal transduction pathways. The most common substrates of FTase are proteins that have C-terminal tetrapeptide CaaX box sequences where the cysteine is the site of modification. However, recent work has shown that five amino acid sequences can also be recognized, including the pentapeptides CMIIM and CSLMQ. In this work, peptide libraries were initially used to systematically vary the residues in those two parental sequences using an assay based on Matrix Assisted Laser Desorption Ionization–Mass Spectrometry (MALDI-MS). In addition, 192 pentapeptide sequences from the human proteome were screened using that assay to discover additional extended CaaaX-box motifs. Selected hits from that screening effort were rescreened using an in vivo yeast reporter protein assay. The X-ray crystal structure of CMIIM bound to FTase was also solved, showing that the C-terminal tripeptide of that sequence interacted with the enzyme in a similar manner as the C-terminal tripeptide of CVVM, suggesting that the tripeptide comprises a common structural element for substrate recognition in both tetrapeptide and pentapeptide sequences. Molecular dynamics simulation of CMIIM bound to FTase further shed light on the molecular interactions involved, showing that a putative catalytically competent Zn(II)-thiolate species was able to form. Bioinformatic predictions of tetrapeptide (CaaX-box) reactivity correlated well with the reactivity of pentapeptides obtained from in vivo analysis, reinforcing the importance of the C-terminal tripeptide motif. This analysis provides a structural framework for understanding the reactivity of extended CaaaX-box motifs and a method that may be useful for predicting the reactivity of additional FTase substrates bearing CaaaX-box sequences.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1993

Axisymmetric computational fluid dynamics analysis of a film/dump-cooled rocket nozzle plume

Prediction of convective base heating rates for a new launch vehicle presents significant challenges to analysts concerned with base environments. The present effort seeks to augment classical base heating scaling techniques via a detailed investigation of the exhaust plume shear layer of a single H2/O2 Space Transportation Main Engine (STME). Use of fuel-rich turbine exhaust to cool the STME nozzle presented concerns regarding potential recirculation of these gases to the base region with attendant increase in the base heating rate. A pressure-based full Navier-Stokes computational fluid dynamics (CFD) code with finite rate chemistry is used to predict plumes for vehicle altitudes of 10 kft and 50 kft. Levels of combustible species within the plume shear layers are calculated in order to assess assumptions made in the base heating analysis.

Tucker, P. K.↗
NASA NTRS2002

Aeroheating Environments for a Mars Smart Lander

A proposed Mars Smart Lander is designed to reach the surface via lifting-body atmospheric entry (alpha = 16 deg) to within 10 km of the target site. CFD (computational fluid dynamics) predictions of the forebody aeroheating environments are given for a direct entry from a 2005 launch. The solutions were obtained using an 8-species gas in thermal and chemical nonequilibrium with a radiative-equilibrium wall temperature boundary condition. Select wind tunnel data are presented from tests at NASA Langley Research Center. Turbulence effects are included to account for both smooth body transition and turbulence due to heatshield penetrations. Natural transition is based on a momentum-thickness Reynolds number value of 200. The effects of heatshield penetrations on turbulence are estimated from wind tunnel tests of various cavity sizes and locations. Both natural transition and heatshield penetrations are predicted to cause turbulence prior to the nominal trajectory peak heating time. Laminar and turbulent CFD predictions along the trajectory are used to estimate heat rates and loads. The predicted peak turbulent heat rate of 63 W/sq cm on the heatshield leeward flank is 70% higher than the laminar peak. The maximum integrated heat load for a fully turbulent heat pulse is 38% higher than the laminar load on the heatshield nose. The predicted aeroheating environments with uncertainty factors will be used to design a thermal protection system.

Edquist, Karl T.↗
DOE OSTI2022

Utilizing CellBox to Describe Varied Dynamical Systems

CellBox is a hybrid modeling framework that integrates machine-learning approaches with explicit mathematical modeling to determine the behavior of dynamical systems, intended for determining the complex interactions of molecule and phenotype within a cellular environment using perturbation/response-based data. The original work utilized analysis of a melanoma cell line, SK-Mel-133, which we sought to expand to and evaluate effectiveness in describing other nonlinear systems both within and outside of biology. We evaluated with two models. One model is based upon the free-fall of an object under a velocity-dependent drag force; this produced inherent perturbation/response data. The other model is based upon a form of mass-action molecular kinetics; for this, we ran the simulation to steady state in each species for each perturbation and select the steady state as input for CellBox. We determined that CellBox, outside of its original use case, can still strongly predict even non-biological dynamical systems in the small-perturbation regime, with r = 0.99977 to 0.85123, p < 0.05 for freefall. However, CellBox cannot accurately predict behavior within large-perturbation regimes, p >> 0.05 for freefall. In addition, we performed a naïve comparison of CellBox network structure to known network structure in the mass-action model to assess accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Long-term, large-scale experiment reveals the effects of seed limitation, climate, and anthropogenic disturbance on restoration of plant communities in a biodiversity hotspot

Ecological restoration is essential for maintaining biodiversity in the face of dynamic, global changes in climate, human land use, and disturbance regimes. Effective restoration requires understanding bottlenecks in plant community recovery that exist today, while recognizing that these bottlenecks may relate to complex histories of environmental change. Such understanding has been a challenge because few long-term, well-replicated experiments exist to decipher the demographic processes influencing recovery for numerous species against the backdrop of multiyear variation in climate and management. Here we address this challenge through a long-term and geographically expansive experiment in longleaf pine savannas, an imperiled ecosystem and biodiversity hotspot in the southeastern United States. Using 48 sites at three locations spanning 480 km, the 8-y experiment manipulated initial seed arrival for 24 herbaceous plant species and presence of competitors to evaluate the impacts of climate variability and management actions (e.g., prescribed burning) on plant establishment and persistence. Adding seeds increased plant establishment of many species. Cool and wet climatic conditions, low tree density, and reduced litter depth also promoted establishment. Once established, most species persisted for the duration of the 8-y experiment. Plant traits were most predictive when tightly coupled to the process of establishment. Our results illustrate how seed additions can restore plant diversity and how interannual climatic variation affects the dynamics of plant communities across a large region. The significant effects of temperature and precipitation inform how future climate may affect restoration and conservation via large-scale changes in the fundamental processes of establishment and persistence.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Simulation of the dynamics of gas mixtures during plasma processing in the C75 Cavity

Plasma processing using a mixture of noble gas and oxygen is a technique that is currently being used to reduce field emission and multipacting in accelerating cavities. Plasma is created inside the cavity when the gas mixture is exposed to an electromagnetic field that is generated by applying RF power through the fundamental power or higher-order mode couplers. Oxygen ions and atomic oxygen are created in the plasma which breaks down the hydrocarbons on the surface of the cavity and the residuals from this process are removed as part of the process gas flow. Removal of hydrocarbons from the surface increases the work function and reduces the secondary emission coefficient. This work describes the initial results of plasma simulation, which provides insight into the ignition process, distribution of different species, and interactions of free oxygen and oxygen ions with the cavity surfaces. The simulations have been done with an Ar/?2 plasma using COMSOL® multiphysics. These simulations help in understanding the dynamics and control of plasma inside the cavity and the exploration of different gas mixtures.

Dhakal, Pashupati↗