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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Studies of Inviscid Flux Schemes for Acoustics and Turbulence Problems

The last two decades have witnessed tremendous growth in computational power, the development of computational fluid dynamics (CFD) codes which scale well over thousands of processors, and the refinement of unstructured grid-generation tools which facilitate rapid surface and volume gridding of complex geometries. Thus, engineering calculations of 10(exp 7) - 10(exp 8) finite-volume cells have become routine for some types of problems. Although the Reynolds Averaged Navier Stokes (RANS) approach to modeling turbulence is still in extensive and wide use, increasingly large-eddy simulation (LES) and hybrid RANS-LES approaches are being applied to resolve the largest scales of turbulence in many engineering problems. However, it has also become evident that LES places different requirements on the numerical approaches for both the spatial and temporal discretization of the Navier Stokes equations than does RANS. In particular, LES requires high time accuracy and minimal intrinsic numerical dispersion and dissipation over a wide spectral range. In this paper, the performance of both central-difference and upwind-biased spatial discretizations is examined for a one-dimensional acoustic standing wave problem, the Taylor-Green vortex problem, and the turbulent channel ow problem.

Morris, Christopher I.↗
NASA NTRSDate not supplied

Degradation of Carbon Fiber Microstructures Due to Oxidative Etch Pitting

In order to understand the influence of etch pits on the material properties of carbon fiber microstructures, we use the Porous Microstructure Analysis (PuMA) software to create realistic pitted structures from virgin fiber structures and then analyze their degradation as the degree of pitting increases. For that purpose, we developed a pitting module within PuMA that creates pits with user-defined geometry (shape, size) and distribution (defect density) on the surface of arbitrary microstructures. In this work, we consider FiberForm microstructures, which are commonly used as the base material for ablative thermal protection systems (TPS). Using the new pitting module, we mimic the evolution of etch pits on FiberForm as it occurs from oxidation in ablative heating environments, by first importing the virgin FiberForm structure from a micrograph scan and then imposing different pitting conditions by varying pit radii and defect densities. We then leverage PuMA’s capabilities to calculate material properties for each structure. The results clearly show the degradation of thermal conductivity and structural integrity as etch pits grow, and furthermore highlight the complex evolution of the surface topology, which results in significant changes of the flowgeometry and chemical gas-surface interactions. Furthermore, a comparison of these changes with the case of shrinking fibers, which has heretofore been the standard for modelling ablation of carbon structures, yields significant differences. For example, we observe the formation of small chunks of fibers as pits grow in size, which could be released into the surrounding flow environment as spallation particles. This motivates us to continue our efforts of modelling pitting as an integral part of the structural degradation of carbon fiber structures in oxidation environments, which will ultimately improve our understanding of failure mechanisms in those materials.

PuMA↗
NASA NTRSDate not supplied

Degradation of Carbon Fiber Microstructures Due to Oxidative Etch Pitting

In order to understand the influence of etch pits on the material properties of carbon fiber microstructures, we use the Porous Microstructure Analysis (PuMA) software to create realistic pitted structures from virgin fiber structures and then analyze their degradation as the degree of pitting increases. For that purpose, we developed a pitting module within PuMA that creates pits with user-defined geometry (shape, size) and distribution (defect density) on the surface of arbitrary microstructures. In this work, we consider FiberForm microstructures, which are commonly used as the base material for ablative thermal protection systems (TPS). Using the new pitting module, we mimic the evolution of etch pits on FiberForm as it occurs from oxidation in ablative heating environments, by first importing the virgin FiberForm structure from a micrograph scan and then imposing different pitting conditions by varying pit radii and defect densities. We then leverage PuMA’s capabilities to calculate material properties for each structure. The results clearly show the degradation of thermal conductivity and structural integrity as etch pits grow, and furthermore highlight the complex evolution of the surface topology, which results in significant changes of the flowgeometry and chemical gas-surface interactions. Furthermore, a comparison of these changes with the case of shrinking fibers, which has heretofore been the standard for modelling ablation of carbon structures, yields significant differences. For example, we observe the formation of small chunks of fibers as pits grow in size, which could be released into the surrounding flow environment as spallation particles. This motivates us to continue our efforts of modelling pitting as an integral part of the structural degradation of carbon fiber structures in oxidation environments, which will ultimately improve our understanding of failure mechanisms in those materials.

PuMA↗
NASA NTRS1971

Hypersonic shock wave interaction and impingement

An experimental investigation was conducted on space shuttle type, body-wing configurations. The purpose of the investigation was to determine the effects of body and wing geometry on the hypersonic shock structure about these vehicles and on the resulting surface impingement of interior flow field shock and expansion waves. Schlieren photographs and thermographic phosphor paint data were obtained on three body cross sections with three wing planforms at 40, 50 and 60 degree angles of attack. Specific configuration data were obtained at 0 and 30 degree angles of attack to develop trends. These data were obtained at a nominal Mach number of 13.5 and a freestream unit Reynolds number of 0.7 million per foot. For comparison with these straight wing configurations, data were also obtained on a model of a point design, high cross-range, delta wing orbiter at 40, 50 and 60 degree angles of attack. As expected, the data on this delta wing orbiter indicated that the shock intersection/impingement phenomena associated with straight wing vehicles are considerably more complex than, and result in both windward and leeward surface heating regions not present on, the delta configuration.

Kessler, W. C.↗
DOE OSTI2026

Scalable Hot-Water-Repellent Superhydrophobicity via Thermal Insulation

Superhydrophobic surfaces, which rely on a combination of surface texture and chemistry, often lose their repellent behavior when contacted by hot water (≳40 °C) because the impinging hot water replaces the requisite air layer within the surface texture via evaporation and recondensation. In contrast to previous approaches targeting this condensation-induced failure mode that rely on intricately tailored surface structures or complex chemical treatments, we present a scalable approach based on thermal design: the multilayered insulated superhydrophobic (MISH) coating mitigates condensation-induced failure by preventing heat transfer. Superhydrophobicity is retained at impinging water temperatures up to 90 °C, with durability demonstrated via long-term (>1 million impacts) droplet impingement experiments. We explain the mechanism for this approach with a detailed thermal model; the model reveals that the underlying physical behavior is self-similar across coating parameters and impinging fluid temperatures. Finally, the MISH coating accommodates curved geometries and large surfaces, and it is over 4 orders of magnitude less expensive than cleanroomnanofabricated alternatives, indicating promise for practical use in the energy sector, chemical processing, and the food and medical industries.

coating materials↗
DOE OSTI2021

Interface creation on a mixed-terminated perovskite surface

In the field of complex oxide heterostructures, understanding of the initial substrate surface can be critical to fundamental studies regarding the development of emergent properties at the film–substrate interface. For this reason, a considerable amount of effort has gone into the development of techniques to achieve surfaces with single termination for a variety of perovskite single crystals. However, a decisive understanding of how an interface is created when epitaxial growth occurs on a mixed terminated surface remains lacking. Employing in situ synchrotron X-ray scattering during thin film growth by molecular beam epitaxy, we investigate the initial stages of growth on (LaAlO 3 )(Sr 2 AlTaO 6 ) (001) substrates with mixed termination. Using LaNiO 3 as a model system, we find that the surface layers of the substrate restructure during deposition such that while a NiO 2 monolayer is weakly bound, a LaO monolayer bonds strongly, effectively incorporating with (Al, Ta)O 2 from the surface and forming a La(Al, Ta)O 3 ultrathin film. <!--/abstract content--><!--fulltext content--> With regard to the synthesis of high-quality complex oxide heterostructures, it is important that the substrate meet an array of criteria related to the crystal structure, lattice parameter, thermodynamic properties, electrical and optical properties, and so on. In particular, the surface of the substrate should be atomically smooth and exhibit the step-terrace structure ideal for epitaxial thin film growth. Perovskite oxides with the formula ABO 3 have two possible (001) surfaces: AO (A-site termination) and BO 2 (B-site termination). The ability to consistently achieve one of these surface terminations on SrTiO 3 , i.e., TiO 2 , has helped SrTiO 3 become the predominant oxide in the field.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Enabling Interoperability in Earth System Digital Twins (ESDT): Integrating Observations, Models, and AI for Actionable Insights Through NASA'S Intelligent Systems Technology Program

NASA’s Intelligent Systems Technology Program (IST) is driving a paradigm shift in Earth science through the development of Earth System Digital Twins (ESDT). These integrated information systems create a dynamic "digital replica" of the Earth by harmonizing continuous, multi-source observations with high-fidelity models and state-of-the-art artificial intelligence (AI) that enable “What now?”, “What next?”, and “What if?” scenario building. These scenarios are reflected in NASA IST’s series of ESDTs, from the Coastal Zone Digital Twin that integrates complex data on the current state of the Chesapeake Bay to the Terrestrial Environmental Rapid-Replication and Assimilation Hydrometeorological (TerraHydro) AI-based ESDT that forecasts water movement across Earth’s surface, to the Agriculture Land Information System (AgLIS) which can be used to assess optimal planting dates and crop yield estimates. By bridging the gap between vast data archives and actionable insights, these projects enable a system-of-systems approach to understanding complex, interacting Earth processes. This poster will highlight recent innovations and future directions from NASA’s ESDT initiatives: Continuous Data Assimilation & Multi-Source Fusion. A core requirement of the ESDT work is the transition from static models to dynamic "living" replicas. This involves creating frameworks for the continual assimilation of near-real-time data from uncoordinated, heterogeneous sources, including satellite observations and airborne assets, and ground-based Internet of Things (IoT) sensors. These systems link design, operational status, and environmental data, ensuring the digital twin accurately reflects the current state of the physical Earth system. High-Fidelity Hybrid Modeling & Computational Acceleration to enable interactive "what-if" explorations, programs are moving beyond traditional, slow physical solvers by developing fast surrogate machine learning models and Deep Generative Models (DGMs). These hybrid approaches use neural networks to emulate complex physics, such as cloud feedback or ocean dynamics, at a fraction of the original computing cost, often leveraging advanced hardware like Graphics Processing Units (GPUs) to achieve the necessary scale. Federated Ecosystems & Interoperable Frameworks rather than building isolated tools, NASA IST is moving toward federated ESDTs and reusable analytic collaborative frameworks. This theme focuses on interoperability standards and common ontologies that allow specialized digital twins to interact and share data. This system-of-systems architecture supports multi-discipline investigations, such as analyzing how upstream watershed changes impact downstream urban flooding or how wildfire emissions affect regional air quality. By leveraging these advancements, ESDTs empower researchers and decision-makers to conduct real-time analysis and run complex hypothetical scenarios, ultimately improving our understanding of Earth’s evolving systems and informing critical real-world applications.

Earth System↗
NASA NTRS1991

Design and test of three active flutter suppression controllers

Three flutter suppression control law design techniques are presented. Each uses multiple control surfaces and/or sensors. The first uses linear combinations of several accelerometer signals together with dynamic compensation to synthesize the modal rate of the critical mode for feedback to distributed control surfaces. The second uses traditional tools (pole/zero loci and Nyquist diagrams) to develop a good understanding of the flutter mechanism and produce a controller with minimal complexity and good robustness to plant uncertainty. The third starts with a minimum energy Linear Quadratic Gaussian controller, applies controller order reduction, and then modifies weight and noise covariance matrices to improve multi-variable robustness. The resulting designs were implemented digitally and tested subsonically on the Active Flexible Wing (AFW) wind tunnel model. Test results presented here include plant characteristics, maximum attained closed-loop dynamic pressure, and Root Mean Square control surface activity. A key result is that simultaneous symmetric and antisymmetric flutter suppression was achieved by the second control law, with a 24 percent increase in attainable dynamic pressure.

Christhilf, David M.↗
DOE OSTI2020

Plasmon coupling in topological insulator multilayers

Topological insulators (TIs) house two-dimensional Dirac plasmons on their surfaces. In TI thin films, plasmons on the top and bottom surfaces couple electrostatically, giving rise to an acoustic and an optical plasmon mode. By extension, a superlattice comprising TI layers and trivially-insulating layers could house multiple complex plasmon modes which could be used to create a new type of Dirac metamaterial. In this paper, we synthesize TI superlattices, fabricate them into stripe arrays, characterize their optical and plasmonic properties, and model the samples using transfer matrices. Furthermore, we excite plasmon modes that couple to the phonons in the superlattice, resulting in hybrid plasmon phonon polaritons. These modes are modeled using an analytical Fano resonance model and the extracted resonant positions are reproduced by transfer matrix modeling. We also excite an epsilon-near-zero mode in the top dielectric material (Bi 0.5 In 0.5 ) 2 Se 3 . Understanding the behavior of the polariton modes in this complex system will lend insight into the many-body interaction in low-dimensional systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Calculation of the Heat Transfer Coefficient in the Outer Body for a Rotational Detonation

Unsteady heat transfer characterization on the combustion surfaces of Rotational Detonation Engines (RDE) is not well understood. It is generally thought that the complex nature of the unsteady, reacting, compressible fluid flow inside the combustion anulus of the RDE causes the convective heat transfer coefficient to be significantly higher than it is in other applications. Empirical models that have been used to analyze this strictly apply to steady flows where dimensionless groups can be employed. The reacting flow fields that are characteristic of RDEs are inherently compressible, three dimensional, unsteady, and turbulent, having properties that change by orders of magnitude throughout the flow field. They will therefore contain multiple length scales and time scales operating everywhere in the flow during all times. Because of this it is not likely that the RDE flow fields will lend themselves to explanation using simple dimensionless parameters. The dimensionless groups have meaning only in situations where length scales and time scales are singular and well defined. In spite of this it may be possible to get a relatively good idea of what the convection heat transfer coefficient is. In this work a numerical study is performed where the inside wall surface temperature distribution in the RDE outer body is systematically changed over a given range that would be characteristic of the start-up flows inside an RDE. For each case, temperature distributions inside the outer containment wall of the RDE was calculated and compared with experimental data. The closet match can then be used to directly calculate the convection heat transfer coefficient on the inside surface of the RDE.

VanOsdol, John↗
NASA NTRS2005

Photo-realistic Terrain Modeling and Visualization for Mars Exploration Rover Science Operations

Modern NASA planetary exploration missions employ complex systems of hardware and software managed by large teams of. engineers and scientists in order to study remote environments. The most complex and successful of these recent projects is the Mars Exploration Rover mission. The Computational Sciences Division at NASA Ames Research Center delivered a 30 visualization program, Viz, to the MER mission that provides an immersive, interactive environment for science analysis of the remote planetary surface. In addition, Ames provided the Athena Science Team with high-quality terrain reconstructions generated with the Ames Stereo-pipeline. The on-site support team for these software systems responded to unanticipated opportunities to generate 30 terrain models during the primary MER mission. This paper describes Viz, the Stereo-pipeline, and the experiences of the on-site team supporting the scientists at JPL during the primary MER mission.

Edwards, Laurence↗
DOE OSTI2025

Closing the Loop between In Situ Stress Complexity and EGS Fracture Complexity

We present an agent-guided approach to CAD geometry decomposition that automates hex/hybrid meshing with graph neural networks (GNNs) to accelerate next-generation ModSim workflows. Our end-to-end pipeline (i) reduces 3D boundary-representation (B-Rep) models to a 2D chordal axis skeleton (CAT) and then to a 1D bipartite graph of surface and curve nodes, (ii) assigns per node labels as Cubit® WebCut actions, (iii) trains a multi-action GNN under supervised learning, and (iv) predicts five surface-node and three curve-node actions on out-of-distribution test geometries. Each graph node carries geometric, topological, and meshing attributes drawn from the B-Rep “skin” and CAT “skeleton,” with two-way mappings across 3D↔2D↔1D representations to maintain traceability back to 3D CAD. The supervised learning model exhibits stable convergence of the binary cross-entropy loss and achieves 98.7% accuracy on unseen lattice models. To operationalize decision-making, we rank predicted commands by geometric significance and prototyped the agent-guided workflow through the Cubit® Meshing PowerTool GUI. As a stretch goal, we explore reinforcement learning (RL) to reduce or remove label requirements and to learn policies for action sequences that maximize total reward (e.g., size of hex-meshable regions and resulting hex mesh quality). When all-hex meshing is not feasible, the agent assists in producing hybrid meshes—prioritizing hex in critical regions and transitioning to tetrahedral elements (tets) elsewhere—maintaining fidelity while ensuring robustness. The overarching objective is to replace manual, heuristics-based decomposition with data-driven, reproducible automation, cutting meshing turnaround time by orders of magnitude. We anticipate direct impact on simulation workflows through intelligent, scalable decomposition of complex CAD models into hex-meshable subdomains.

42 ENGINEERING↗
DOE OSTI2023

Mechanisms of Ethylene Epoxidation over Silver from Machine Learning-Accelerated First-Principles Modeling and Microkinetic Simulations

In this work we employed machine learning-augmented density functional theory (DFT) thermodynamic calculations to assess the stability of different AgO x structures under catalytic ethylene epoxidation reaction conditions. We found that there are multiple AgO x surface motifs that could co-exist under the relevant conditions. These included Ag surface oxides (e.g., AgO_p(4 × 4) and Ag 1.83 O) and atomic oxygen-covered Ag(111) surfaces. Furthermore, we employed DFT calculations to evaluate the energetics of different reaction mechanisms by which ethylene and oxygen can react on these surfaces. These studies revealed several energetically viable reaction pathways for ethylene epoxidation. Microkinetic modeling analyses, based on the DFT-calculated reaction pathways, showed that ethylene epoxidation can proceed on all surfaces and that multiple pathways, including those involving Langmuir–Hinshelwood and Eley–Rideal mechanisms, could be involved in selective and unselective reactions. The diversity of mechanisms that we discovered in the context of the relatively simple ethylene epoxidation reaction on Ag suggests that the richness and complexity of surface chemistry are most likely a rule rather than an exception in heterogeneous catalytic chemical transformations on metal surfaces and that the concept of a single or even a dominant mechanism and reaction intermediates might need to be revisited for many reactions.

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

Predicting Pulsed-Laser Deposition SrTiO 3 Homoepitaxy Growth Dynamics Using High-Speed Reflection High-Energy Electron Diffraction

Pulsed-laser deposition (PLD) is a powerful technique for growing complex oxides with controlled stoichiometry. To understand growth dynamics therein, it is common to leverage in situ spectroscopies, such as reflection high-energy electron diffraction (RHEED), to monitor surface crystallinity. Most commercial systems rely on video-rate cameras operating at 60-120 Hz that lack sufficient temporal resolution to capture growth dynamics at practical deposition frequencies. Here, a high-speed platform to record in situ dynamics via RHEED at >500 Hz is implemented. An open-source analysis package is designed to fit diffraction spots to 2D Gaussians, allowing single-pulse surface reconstruction kinetics extraction. Using homoepitaxially deposited (001)-oriented SrTiO 3 as a model system, we demonstrate how high-speed RHEED can provide real-time insight into growth processes obscured by slower acquisition systems. By fitting the single-pulse intensity to a set of exponential functions, we observe changes in the characteristic decay time and mechanism correlated to the substrate step width and surface termination. We observe distinct surface effects, with diffraction intensity decaying on lower-energy TiO 2 -terminated surfaces and stabilizing on SrO- or mixed-terminated surfaces. Similarly, using an exponential model, the extracted characteristic time of adatom deposition decreases with increased density of bonding sites associated with mixed termination and narrower step widths. Ultimately, this work shows how increasing RHEED temporal resolution can uncover new insights into growth processes, with practical implications for the design and control of PLD processes. This experimental platform provides new capabilities to enable data-driven machine learning analysis and autonomous control systems to enhance the complexity and fecundity of PLD.

(SrO)↗
DOE OSTI2020

Soil Respiration Response to Rainfall Modulated by Plant Phenology in a Montane Meadow, East River, Colorado, USA

Soil respiration is a primary component of the terrestrial carbon cycle. However, predicting the response of soil respiration to climate change remains a challenge due to the complex interactions between environmental drivers, especially plant phenology, temperature, and soil moisture. In this study, we use a 1-D diffusion-reaction model to calculate depth-resolved CO 2 production rates from soil CO 2 concentrations and surface efflux observations in a subalpine meadow in the East River watershed, CO. Modeled rates are compared to in situ soil temperature and moisture conditions and MODIS satellite enhanced vegetation index (EVI) representing plant phenology across three hydrologically distinct growing seasons from 2016-2018. While soil respiration correlated with temperature on diel timescales (p<0.05), seasonal variability was dominated by soil moisture and plant phenology (p<0.05). We observed significant respiration increases in response to precipitation events; however, magnitude and duration were significantly higher in 2017 than 2016 despite similar wetting characteristics. Based on MODIS EVI, we suggest the respiration response to rainfall is controlled by plant phenology, which in turn reflects the capacity of plants to respond to precipitation via increased photosynthesis and autotrophic respiration, behavior that is not captured in typical soil respiration pulse models. Projected changes in montane climate such as earlier snowmelt and prolonged fore-summer drought may decrease soil respiration fluxes by decreasing the overlap between peak productivity and the summer monsoon. Lastly, we observed significant late-season CO 2 fluxes from the deep sub-soil (>165 cm) that support growing evidence for the importance of sub-soil processes in driving integrated respiration fluxes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Anolyte Enhances Catalyst Utilization and Ion Transport Inside a CO 2 Electrolyzer Cathode

Electrochemical CO 2 reduction is a promising technology to capture and convert CO 2 to valuable chemicals. High Faradaic efficiencies of CO 2 reduction products are achieved with zero-gap alkaline CO 2 electrolyzers with a supporting electrolyte at the anode (anolyte). Herein, we investigate the effect of anolyte on the electrode properties such as catalyst utilization, ionic accessibility etc. of a CO 2 reduction cathode using electrochemical techniques and cell configurations that avoid the complexities related to co-electrolysis. Using 1M KOH as the anolyte and a Cu gas-diffusion-electrode with low Nafion content as the model CO 2 reduction electrode, we find that electrode capacitance (proxy for electrochemically active surface area) and ionic conductivity inside the cathode increase approximately 4 and 447 times, respectively, in presence of KOH. Liquid anolyte wets the electrode’s pore structure more efficiently than capillary condensation of feed water vapor. The ionomer coverage is very low, and its distribution inside the electrode is highly fragmented. Surface ion conduction mechanisms inside the electrode are orders of magnitude lower than the bulk ion conduction in presence of anolyte. This study shows that when an anolyte (e.g., KOH) is used, catalyst utilization and ionic accessibility inside the electrode increase significantly.

30 DIRECT ENERGY CONVERSION↗
NASA NTRS1974

Six-color photometry of Iapetus, Titan, Rhea, Dione and Tethys

The extensive photometric observations of Titan, Iapetus, Rhea, Dione and Tethys have made it possible to separate the solar phase and orbital phase contributions to the observed light variations of these satellites. For Titan, the wavelength dependence of its solar phase coefficient has been obtained. This dependence should prove useful in constructing future model atmospheres. The other four satellites show a surprising array of different photometric behaviors. Despite these differences, all four satellites have similar spectral reflectivities. Clearly Iapetus, Rhea, Dione and Tethys are complex objects, varying substantially from one another in surface structure and/or composition.

Noland, M.↗