SEARCH · Engineering Papers

Results for “interfacial engineering”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Interfacial thermal conductance between multi-layer graphene sheets and solid/liquid octadecane: A molecular dynamics study

Mixtures of paraffin and carbon nanofillers have promising potential for thermal storage, as paraffin (the matrix) possesses high latent heat and the nanofillers compensate for the low thermal conductivity (TC) of paraffin. Understanding thermal transport in these materials is essential for practical applications, as weak thermal transport hinders fast charge/discharge of thermal energy. Here, we use non-equilibrium molecular dynamics (NEMD) simulations to study the interfacial thermal conductance (ITC) between graphene sheets and octadecane (C18H38) matrix under the limiting conditions of the sheets being parallel or perpendicular to the direction of the imposed heat flux. The findings show that the systems containing thin graphene layers exhibit higher values of ITC. This study captures the asymptotic saturation of thermal conductance for the liquid phase of the perpendicular structure. Besides, given the greater number of structured layers of paraffin upon phase change, the ITC for the solid paraffin-graphene system is higher than the conductance of the liquid paraffin-graphene interface. We use the Pearson correlation coefficients of the vibrational power spectrum (VPS) of interfacing materials to explain the orders of magnitude variations of the observed ITC.

25 ENERGY STORAGE↗

DOE OSTI2017

A highly scalable spray coating technique for electrode infiltration: Barium carbonate infiltrated La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3-δ perovskite structured electrocatalyst with demonstrated long term durability

Barium carbonate (BaCO 3 ) nanoparticles have demonstrated excellent catalytic activity for the oxygen reduction reaction by spray coating onto La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3-δ (LSCF) and LSCF-SDC (Sm 0.2 Ce 0.8 O 2-δ ) cathode materials. In this work, barium acetate solutions were modified by a surfactant to lower the surface tension and decrease the contact angle on LSCF which is a benefit for the infiltration process. In the LSCF electrode, BaCO 3 nano-particles exhibited significant interfacial contact with LSCF particles by the spray coating technique. As a result, the polarization resistance of BaCO 3 infiltrated LSCF was reduced from 2.5 to 1.2 Ωcm 2 at 700 °C. In addition, commercial full cell SOFCs with BaCO 3 infiltrated LSCF-SDC cathodes also demonstrated higher performance due to the reduced cathode resistance. At 750 °C, the electrode overpotential of the BaCO3 infiltrated cell was much lower than that of baseline cell during long term testing (500 hours). The polarization resistance of the BaCO 3 infiltrated LSCF-SDC electrode only increased by 1.6% after 500 hours. This work demonstrated a highly scalable spray coating process of a promising electrocatalyst (BaCO 3 ) with excellent long-term stability for the oxygen reduction reaction.

42 ENGINEERING↗

DOE OSTI2026

Deciphering the Dynamic Balance Between Solvation Strength and Polysulfides Reaction Heterogeneity in Practical Lithium‐Sulfur Batteries

Achieving stable interfacial chemistry in lithium–sulfur batteries under practical conditions remains a key barrier to commercialization. Here, we demonstrate that interfacial dynamics can be effectively regulated by coupling solvation-power control with intrinsic heterogeneity of sulfur redox chemistry through the introduction of a weakly solvating fluorinated cosolvent, LIB 1200ET (1200ET). Compared with conventional fluorinated ethers, 1200ET efficiently shifts Li + solvation environment toward a more non-coordinated configuration at low volume fractions, enabling substantial solvation modulation without significantly impairing sulfur redox kinetics. This solvation transition weakens Li + –solvent interactions while strengthening Li + –anion and Li + –lithium polysulfide (LPS) coordination, suppressing LPS solubility and promoting reconstruction of solid–electrolyte interphase (SEI). Regulated LPS chemistry, together with 1200ET, leads to formation of a S 4+ -rich, LiF-reinforced SEI with enhanced ionic conductivity and mechanical robustness. Spatially resolved sulfur K-edge X-ray absorption spectroscopy on pouch cells reveals pronounced current-density-dependent chemical heterogeneity, distinguishing kinetically dominated and solvation-controlled regions. Under practical conditions (3.7 mg cm −2 sulfur loading, E/S = 6 µL mg −1 ), a single-layer pouch cell delivers 527 mAh g −1 over 200 cycles at C/3, while an Ah-level multilayer pouch cell achieves an energy density of 358 Wh kg −1 . These results establish non-coordinating cosolvent-driven solvation engineering as a scalable strategy for practical Li–S batteries.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Enhanced Electrocatalytic and Cathode‐Electrolyte Interfacial Properties With a Pr‐Based Simple Perovskite/Ruddlesden‐Popper Nanocomposite Cathode in Protonic Ceramic Fuel Cells

The sluggish kinetics and poor stability of the oxygen reduction reaction (ORR) remain the primary bottleneck for achieving high performance in protonic ceramic fuel cells (PCFCs) at intermediate temperatures (400–650°C). In this work, a Pr-based nanocomposite cathode comprised of simple perovskite phase (PrNi 0.7 Co 0.3 O 3-δ ) and Ruddlesden-Popper phase (Co-doped Pr 4 Ni 3 O 10+δ ) is developed. Although PrNi 0.7 Co 0.3 O 3-δ solely stands as a good cathode with facile proton transfer, combining the superior catalytic activity against oxygen on the Ruddlesden-Popper phase boosts the ORR performance further. The designed nanocomposite cathode outperforms the simple perovskite cathode, attributed to enhanced oxygen absorption and surface diffusion with the Ruddlesden-Popper phase. A precursor-based cathode deposition technique is also developed to achieve cathode grain sizes of ∼100 nm. A single cell with the nanocomposite cathode delivers a peak power density of 1.38 W cm −2 at 650°C, among the highest in reported PCFCs with Pr-based cathodes, with a small degradation rate of 0.145 mV h −1 during 250 h stability test. Further investigation of cathode-electrolyte interface revealed interfacial PrO 2 phase formation, promoted by abundant Pr 6 O 11 in the nanocomposite precursor powder, thereby improving both ohmic resistance and stability. These findings highlight the effectiveness of the nanocomposite cathode and underscore its advantages on interfacial properties.

08 - HYDROGEN↗

DOE OSTI2024

Magnet-in-ferroelectric crystals exhibiting photomultiferroicity

Growing crystallographically incommensurate and dissimilar organic materials is fundamentally intriguing but challenging for the prominent cross-correlation phenomenon enabling unique magnetic, electronic, and optical functionalities. Here, we report the growth of molecular layered magnet-in-ferroelectric crystals, demonstrating photomanipulation of interfacial ferroic coupling. The heterocrystals exhibit striking photomagnetization and magnetoelectricity, resulting in photomultiferroic coupling and complete change of their color while inheriting ferroelectricity and magnetism from the parent phases. Under a light illumination, ferromagnetic resonance shifts of 910 Oe are observed in heterocrystals while showing a magnetization change of 0.015 emu/g. In addition, a noticeable magnetization change (8% of magnetization at a 1,000 Oe external field) in the vicinity of ferro-to-paraelectric transition is observed. The mechanistic electric-field-dependent studies suggest the photoinduced ferroelectric field effect responsible for the tailoring of photo-piezo-magnetism. The crystallographic analyses further evidence the lattice coupling of a magnet-in-ferroelectric heterocrystal system.

42 ENGINEERING↗

DOE OSTI2018

Molecular Dynamics-based Simulations of Bulk/Interfacial Structures and Diffusion Behaviors in Nuclear Waste Glasses

This NEUP Project aimed to generate accurate atomic structural models of nuclear waste glasses by using large-scale molecular dynamics-based computer simulations and to use these models to investigate self-diffusion behaviors, interfacial structures, and hydrated gel structures formed during dissolution of these glasses. The goal was to obtain realistic and accurate short and medium range structures of these complex oxide glasses, to provide a mechanistic understanding of the dissolution behaviors, and to generate reliable information with predictive power in designing nuclear waste glasses for long-term geological storage. Looking back of the research accomplishments of this project, most of the scientific goals initially proposed have been achieved through intensive research in the three and a half year period of the project. This project has also generated a wealth of scientific data and vibrant discussions with various groups through collaborations within and outside of this project. Throughout the project one book chapter and 14 peer reviewed journal publications have been generated (including one under review) and 16 presentations (including 8 invited talks) have been made to disseminate the results of this project in national and international conference. Furthermore, this project has trained several outstanding graduate students and young researchers for future workforce in nuclear related field, especially on nuclear waste immobilization. One postdoc and four PhD students have been fully or partially supported through the project with intensive training in the field material science and engineering with expertise on glass science and nuclear waste disposal

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS2004

Electron-Spin Filters Based on the Rashba Effect

Semiconductor electron-spin filters of a proposed type would be based on the Rashba effect, which is described briefly below. Electron-spin filters more precisely, sources of spin-polarized electron currents have been sought for research on, and development of, the emerging technological discipline of spintronics (spin-based electronics). There have been a number of successful demonstrations of injection of spin-polarized electrons from diluted magnetic semiconductors and from ferromagnetic metals into nonmagnetic semiconductors. In contrast, a device according to the proposal would be made from nonmagnetic semiconductor materials and would function without an applied magnetic field. The Rashba effect, named after one of its discoverers, is an energy splitting, of what would otherwise be degenerate quantum states, caused by a spin-orbit interaction in conjunction with a structural-inversion asymmetry in the presence of interfacial electric fields in a semiconductor heterostructure. The magnitude of the energy split is proportional to the electron wave number. The present proposal evolved from recent theoretical studies that suggested the possibility of devices in which electron energy states would be split by the Rashba effect and spin-polarized currents would be extracted by resonant quantum-mechanical tunneling. Accordingly, a device according to the proposal would be denoted an asymmetric resonant interband tunneling diode [a-RITD]. An a-RITD could be implemented in a variety of forms, the form favored in the proposal being a double-barrier heterostructure containing an asymmetric quantum well. It is envisioned that a-RITDs would be designed and fabricated in the InAs/GaSb/AlSb material system for several reasons: Heterostructures in this material system are strong candidates for pronounced Rashba spin splitting because InAs and GaSb exhibit large spin-orbit interactions and because both InAs and GaSb would be available for the construction of highly asymmetric quantum wells. This mate-rial system affords a variety of energy-band alignments that can be exploited to obtain resonant tunneling and other desired effects. The no-common-atom InAs/GaSb and InAs/AlSb interfaces would present opportunities for engineering interface potentials for optimizing Rashba spin splitting.

Ting, David Z.-Y.↗

NASA NTRS2015

Effects of Atomic-Scale Structure on the Fracture Properties of Amorphous Carbon - Carbon Nanotube Composites

The fracture of carbon materials is a complex process, the understanding of which is critical to the development of next generation high performance materials. While quantum mechanical (QM) calculations are the most accurate way to model fracture, the fracture behavior of many carbon-based composite engineering materials, such as carbon nanotube (CNT) composites, is a multi-scale process that occurs on time and length scales beyond the practical limitations of QM methods. The Reax Force Field (ReaxFF) is capable of predicting mechanical properties involving strong deformation, bond breaking and bond formation in the classical molecular dynamics framework. This has been achieved by adding to the potential energy function a bond-order term that varies continuously with distance. The use of an empirical bond order potential, such as ReaxFF, enables the simulation of failure in molecular systems that are several orders of magnitude larger than would be possible in QM techniques. In this work, the fracture behavior of an amorphous carbon (AC) matrix reinforced with CNTs was modeled using molecular dynamics with the ReaxFF reactive forcefield. Care was taken to select the appropriate simulation parameters, which can be different from those required when using traditional fixed-bond force fields. The effect of CNT arrangement was investigated with three systems: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. For each arrangement, covalent bonds are added between the CNTs and AC, with crosslink fractions ranging from 0-25% of the interfacial CNT atoms. The SWNT and MWNT array systems represent ideal cases with evenly spaced CNTs; the SWNT bundle system represents a more realistic case because, in practice, van der Waals interactions lead to the agglomeration of CNTs into bundles. The simulation results will serve as guidance in setting experimental processing conditions to optimize the mechanical properties of CNT composites.

Jensen, Benjamin D.↗

DOE OSTI2021

Design Principles for Efficient and Stable Water Splitting Photoelectrocatalysts

Photoelectrochemical water splitting is a promising avenue for sustainable production of hydrogen used in the chemical industry and hydrogen fuel cells. The basic components of most photoelectrochemical water splitting systems are semiconductor light absorbers coupled to electrocatalysts, which perform the desired chemical reactions. A critical challenge for the design of these systems is the lack of stability for the majority of desired semiconductors under operating water splitting conditions. One strategy to address this issue is to protect the semiconductor by covering it with a stabilizing insulator layer, creating a metal–insulator–semiconductor (MIS) architecture, which has demonstrated improved stability. In addition to enhanced stability, the insulator layer may significantly affect the electron and hole transfer, which governs the recombination rates. Furthermore, the insertion of an insulator layer leads to the introduction of additional insulator/electrocatalyst and insulator/semiconductor interfaces. These interfaces can impact the system’s performance significantly, and they need to be carefully engineered to optimize the efficiencies of MIS systems. In this Account, we describe our recent progress in shedding light on the critical role of the insulator and the interfaces on the performance of MIS systems. We discuss our findings by focusing on the concrete example of planar n-type Si protected by a HfO 2 insulator layer and coupled to a Ni or Ir electrocatalyst that performs the oxygen evolution reaction, one of the water splitting half-reactions. To improve our fundamental understanding of the insulator layer, we precisely control the HfO 2 insulator thickness using atomic layer deposition (ALD), and we perform a series of rigorous electrochemical experiments coupled with theory and modeling. We demonstrate that by tuning the insulator thickness, we can control the flux and recombination of photogenerated electrons and holes to optimize the generated photovoltage. Despite optimizing the thickness, we find that the maximum generated photovoltage in MIS systems is often significantly lower than the upper performance limit, i.e., there are additional losses in the system that could not be addressed by optimizing the insulator thickness. We identify the sources of these losses and describe strategies to minimize them by a combination of improving the semiconductor light absorption, removing nonidealities associated with interfacial defects, and finding alternative insulators with improved charge carrier selectivity. Finally, we quantify the improvements that can be obtained by implementing these specific strategies. Our collective work outlines strategies to analyze MIS systems, identify the sources of efficiency losses, and optimize the design to approach the fundamental performance limits. Furthermore, these general approaches are broadly applicable to photoelectrochemical materials that utilize sunlight to produce value-added chemicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Dynamic structure of active sites in ceria-supported Pt catalysts for the water gas shift reaction

Oxide-supported noble metal catalysts have been extensively studied for decades for the water gas shift (WGS) reaction, a catalytic transformation central to a host of large volume processes that variously utilize or produce hydrogen. There remains considerable uncertainty as to how the specific features of the active metal-support interfacial bonding—perhaps most importantly the temporal dynamic changes occurring therein—serve to enable high activity and selectivity. Here we report the dynamic characteristics of a Pt/CeO 2 system at the atomic level for the WGS reaction and specifically reveal the synergistic effects of metal-support bonding at the perimeter region. We find that the perimeter Pt 0 -O vacancy -Ce 3+ sites are formed in the active structure, transformed at working temperatures and their appearance regulates the adsorbate behaviors. We find that the dynamic nature of this site is a key mechanistic step for the WGS reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Mathematical Modelling of Gas–Liquid, Two-Phase Flows in a Ladle Shroud

Differential and macroscopic models of argon–steel flows in ladle shroud have been developed. In this, argon–steel, two-phase flow phenomena have been formulated via a transient, three dimensional, turbulent flow model, based on the volume of fluid (VOF) calculation procedure. While realizable k–ε turbulence model has been applied to map turbulence, commercial, CFD software ANSYS-Fluent™ (Version 18), has been applied to carry out numerical calculations. Predictions from the model have been directly assessed against experimental measurements across the range of shroud dimensions and volumetric flow rates typically practiced in the industry. It is demonstrated that the two-phase turbulent flow model captures the general features of gas–liquid flows in ladle shroud providing estimates of free jet length and threshold gas flow rates (required to halt air ingression) which are in agreement with corresponding experimental measurements. In the absence of differential solutions, a macroscopic model has been worked out through dimensional analysis embodying multiple non-linear regression. It is shown that dimensionless free jet length in bloom and slab casting shrouds can be estimated reasonably accurately from the following correlation (in SI unit), viz., ((L{sub jet})/(D{sub CN}))=2.8×10{sup −2}(((Q{sub G})/(Q{sub L}))){sup 1.14}(((gD{sub CN}{sup 5})/(Q{sub L}{sup 2}))){sup 0.8}(((σD{sub CN}{sup 3})/(ρ{sub L}Q{sub L}{sup 2}))){sup −0.9}(((D{sub sh})/(D{sub CN}))){sup 2.0}(((ρ{sub G})/(ρ{sub L}))){sup −0.30}in which, L{sub jet} is the free liquid jet length (m), Q{sub G} is the gas flow rate (m{sup 3}/s), Q{sub L} is the liquid flow rate (m{sup 3}/s), D{sub sh} is shroud diameter (m), D{sub CN} is the collector nozzle diameter (m), σ is the interfacial tension (N/m), and ρ{sub G} as well as ρ{sub L} are respectively density of gas and liquid (kg/m{sup 3}). It is demonstrated that the proposed correlation is consistent with the laws of physical modeling and leads to estimates that are in good agreement with predictions from the differential models, for both air-water as well as argon–steel systems. Numerical simulations as well as macroscopic modeling have indicated that thermo-physical properties of the gas–liquid system are important and exert some influences on the gas–liquid, two-phase, flow in ladle shrouds, albeit not to a large extent. Despite dissimilar thermo-physical properties, full scale water modeling appears to be sufficiently predictive and provides reasonable macroscopic descriptions of the two-phase flow phenomena in industrial ladle shroud systems.

42 ENGINEERING↗

DOE OSTI2025

High-Nickel Cathodes with Mechanical and Interfacial Robustness via Tailored Concentration Gradients for Stable Li-Ion Batteries

Here, we have developed a versatile mathematical framework integrated with an automated reactor system to design and reify highly customizable full concentration gradient (FCG) in high-nickel cathodes for advanced Li-ion batteries. This method provides precise and independent control of the average composition, slope, and curvature of FCGs, enabling the optimization of structural and mechanical properties of the cathode materials. We have showcased this method with Ni 0.8 Co 0.1 Mn 0.1 (OH) 2 precursors of controlled FCGs, which unlocked an optimized cathode with excellent cycling stability without crack formation after repeated cycles. This work opens up new possibilities for the design and manufacturing of advanced cathode materials, enabling safer, high-performance batteries.

25 ENERGY STORAGE↗

DOE OSTI2017

Numerical Modeling Describing the Effects of Heterogeneous Distributions of Asperities on the Quasi-static Evolution of Frictional Slip

A better understanding of how slip accumulates along faults and its relation to the breakdown of shear stress is beneficial to many engineering disciplines, such as, hydraulic fracture and understanding induced seismicity (among others). Asperities forming along a preexisting fault resist the relative motion of the two sides of the interface and occur due to the interaction of the surface topographies. We employ a finite element model to simulate circular partial slip asperities along a nominally flat frictional interface. Shear behavior of our partial slip asperity model closely matched the theory described by Cattaneo. The asperity model was employed to simulate a small section of an experimental fault formed between two bodies of polymethyl methacrylate, which consisted of multiple asperities whose location and sizes were directly measured using a pressure sensitive film. The quasi-static shear behavior of the interface was modeled for cyclical loading conditions, and the frictional dissipation (hysteresis) was normal stress dependent. We further our understanding by synthetically modeling lognormal size distributions of asperities that were randomly distributed in space. Synthetic distributions conserved the real contact area and aspects of the size distributions from the experimental case, allowing us to compare the constitutive behaviors based solely on spacing effects. Traction-slip behavior of the experimental interface appears to be considerably affected by spatial clustering of asperities that was not present in the randomly spaced, synthetic asperity distributions. Finally, estimates of bulk interfacial shear stiffness were determined from the constitutive traction-slip behavior and were comparable to the theoretical estimates of multi-contact interfaces with non-interacting asperities.

58 GEOSCIENCES↗

DOE OSTI2022

Waste to worth: A high-temperature water-gas shift magnetite catalyst with encapsulated core-shell structure from coal fly ash

In this paper, we demonstrate the synthesis of magnetite-based, encapsulated catalyst with a core-shell structure from coal fly ash, an industrial solid waste, and its application to industrially important high temperature water-gas shift reaction (HT-WGSR). The catalyst precursor was first synthesized through a leaching-precipitation-hydrothermal process, leading to the formation of a nanosheet-like hematite precursor consisting of Fe 3+ as primary element and minor cations including Al 3+ , Mg 2+ and other cations (e.g. Ca 2+ and Ti 4+ ). During the WGSR, the nanosheets broke into smaller particles with a core of magnetite with a diameter of ~20 nm, whereas the inherent impurity elements were activated and migrated to the outer surface forming a shell with a thickness of ~4 nm. Due to the protection of shell, the as-synthesized catalyst exhibited improved activity, thermal stability and long-term durability for HT-WGSR at 400 °C and 450 °C. Compared to the pure hematite, this catalyst exhibited much higher activity and stability, i.e., increased CO conversion and decreased deactivation rate constant, by a factor of four and three, respectively. Additionally, through a bottom-up reverse synthesis, we have confirmed that Al 3+ functioned as a textual promoter enhancing the interfacial area of the catalyst. On the other hand, the presence of minor impurity elements, in particular 2.58 wt% Mg 2+ and 0.18 wt% Ti 4+ , were found to increase the basicity, and promote the adsorption rate of CO, the redox of catalysts, and the alteration of the reaction mechanism from regenerative to associative with the formation of a strong intermediate formate. Finally, this study not only opens a new direction on the valorization of solid waste into high-value catalysts, promoting the sustainability of solid waste management, but also provides an alternative, cost-effective and scalable approach for the design and fine-tuning of encapsulated catalysts.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2020

The Role of Ligand–Ligand Interactions in Multimodal Ligand Conformational Equilibria and Surface Pattern Formation

Multimodal chromatography uses multiple modes of interaction such as charge, hydrophobic, or hydrogen bonding to separate proteins. Recently, we used molecular dynamics (MD) simulations to show that ligands immobilized on surfaces can interact and associate with neighboring ligands to form hydrophobic and charge patches, which may have important implications for the nature of protein–surface interactions. Here, we study interfacial systems of increasing complexity—from a single immobilized multimodal ligand to high density surfaces—to better understand how ligand behavior is affected by the presence of a surface and the presence of other ligands in the vicinity, and how this behavior scales to larger systems. Furthermore, we find that tethering a ligand to a surface restricts its conformations to a subset of those observed in free solution, yet the ligand maintains flexibility in the plane of the surface and can form contacts with neighboring ligands. We find that although the formation of a contact between two neighboring ligands is slightly unfavorable, three neighboring ligands exhibit a preference for the formation of a fully connected cluster. To explore how these trends in ligand association extend to a larger surface with high density of ligands, we performed coarse-grained Monte Carlo (MC) simulations of a 132-ligand surface using ligand interactions parametrized based on free energies obtained from the three-ligand MD simulations. Despite their simplicity, the coarse-grained simulations qualitatively capture the cluster size distribution of ligands observed in detailed MD simulations. Quantitative differences between the two suggest opportunities for improvements in the coarse-grained energy function for efficient predictions of cluster and pattern formations. Our approach presents a promising route to the engineering of multimodal patterns for future chromatographic resin design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

An Eulerian multimaterial framework for simulating high-explosive aquarium tests

Aquarium tests of cylindrical high-explosive charges provide optical data of the detonation front velocity and shape, propagation of the shock wave in the surrounding water, and expansion rates of the detonation products behind the front. Data from aquarium experiments is often used for calibration of reactive burn models based on phenomenological equations of state (EOS) and reaction rate laws. This paper presents a multimaterial numerical modeling framework to solve the 2D axisymmetric reactive Euler equations for high-explosive aquarium tests, in particular for ammonium nitrate - fuel oil (ANFO) explosives. An extension of the Ghost Fluid Method (GFM) is used to handle the dynamic material interfaces for the ANFO explosion products, the charge-confining material (polymethyl methacrylate PMMA), and the surrounding water. This study analyzes the sensitivity of calculations (both computational efficiency and numerical accuracy) to different algorithms for the material interface models including the original GFM versus Riemann solver-based strategies. A novel method for defining the left and right states in the interfacial Riemann problem eliminates the need for sorting or nodal interpolation during the projection along the material interface. Numerical tests indicate that populating the interface node values using the Riemann solution mitigate the overheating error observed in steady-state calculations. Solution convergence and computational efficiency are explored as a function of the spatial and temporal order of the schemes. Results from the computational model with analytical equations of state and fitted reaction rate parameters show very good quantitative agreement with experimentally observed detonation front velocity, reaction products expansion, and shock wave propagation in the surrounding water for a cylindrical ANFO charge encased in PMMA. Finally, the proposed modeling framework, in conjunction with experimental tests, provides a reliable tool to assess equations of state and reaction rate expressions for reactive burn models of confined high explosives.

42 ENGINEERING↗

DOE OSTI2017

Non-Darcy interfacial dynamics of air-water two-phase flow in rough fractures under drainage conditions

Two-phase flow interfacial dynamics in rough fractures is fundamental to understanding fluid transport in fractured media. The Haines jump of non-Darcy flow in porous media has been investigated at pore scales, but its fundamental processes in rough fractures remain unclear. In this study, the micron-scale Haines jump of the air-water interface in rough fractures was investigated under drainage conditions, with the air-water interface tracked using dyed water and an imaging system. The results indicate that the interfacial velocities represent significant Haines jumps when the meniscus passes from a narrow "throat" to a wide "body", with jump velocities as high as five times the bulk drainage velocity. Locally, each velocity jump corresponds to a fracture aperture variation; statistically, the velocity variations follow an exponential function of the aperture variations at a length scale of ~100 µm to ~100 mm. This spatial-scale-invariant correlation may indicate that the high-speed local velocities during the Haines jump would not average out spatially for a bulk system. The results may help in understanding the origin of interface instabilities and the resulting non-uniform phase distribution, as well as the micron-scale essence of the spatial and temporal instability of two-phase flow in fractured media at the macroscopic scale.

42 ENGINEERING↗

NASA NTRS1993

Engine materials characterization and damage monitoring by using x ray technologies

X ray attenuation measurement systems that are capable of characterizing density variations in monolithic ceramics and damage due to processing and/or mechanical testing in ceramic and intermetallic matrix composites are developed and applied. Noninvasive monitoring of damage accumulation and failure sequences in ceramic matrix composites is used during room-temperature tensile testing. This work resulted in the development of a point-scan digital radiography system and an in situ x ray material testing system. The former is used to characterize silicon carbide and silicon nitride specimens, and the latter is used to image the failure behavior of silicon-carbide-fiber-reinforced, reaction-bonded silicon nitride matrix composites. State-of-the-art x ray computed tomography is investigated to determine its capabilities and limitations in characterizing density variations of subscale engine components (e.g., a silicon carbide rotor, a silicon nitride blade, and a silicon-carbide-fiber-reinforced beta titanium matrix rod, rotor, and ring). Microfocus radiography, conventional radiography, scanning acoustic microscopy, and metallography are used to substantiate the x ray computed tomography findings. Point-scan digital radiography is a viable technique for characterizing density variations in monolithic ceramic specimens. But it is very limited and time consuming in characterizing ceramic matrix composites. Precise x ray attenuation measurements, reflecting minute density variations, are achieved by photon counting and by using microcollimators at the source and the detector. X ray computed tomography is found to be a unique x ray attenuation measurement technique capable of providing cross-sectional spatial density information in monolithic ceramics and metal matrix composites. X ray computed tomography is proven to accelerate generic composite component development. Radiographic evaluation before, during, and after loading shows the effect of preexisting volume flaws on the fracture behavior of composites. Results from one-, three-, five-, and eight-ply ceramic composite specimens show that x ray film radiography can monitor damage accumulation during tensile loading. Matrix cracking, fiber-matrix debonding, fiber bridging, and fiber pullout are imaged throughout the tensile loading of the specimens. In situ film radiography is found to be a practical technique for estimating interfacial shear strength between the silicon carbide fibers and the reaction-bonded silicon nitride matrix. It is concluded that pretest, in situ, and post-test x ray imaging can provide greater understanding of ceramic matrix composite mechanical behavior.

Baaklini, George Y.↗