SEARCH · Engineering Papers
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.
Illustration of hot hole transfer through interfacial S–Pb interactions, and a peak efficiency of 27.30% was achieved under 5.9 suns via ultrafast hot hole extraction.
First-principles derived computational fluid dynamics (CFD) simulations have been proposed as a fundamental tool for investigating solvent-based CO 2 absorption in packed columns due to their ability to accurately represent the underlying nonlinear, multiscale dynamics. Numerous studies have previously utilized such CFD simulations to investigate hydrodynamics of columns with structured and random packings by assessing the key hydrodynamic metrics such as the interfacial and wetted areas. While mapping such metrics for different conditions is essential to the optimization of absorption columns, it is not sufficient, as the CO 2 capture rate depends also on the coupled, nonlinear dynamics from the underlying chemical reaction kinetics, thermodynamics, and heat-transfer rates. In this work, we present detailed CFD simulation results augmented by incorporating the effects of interfacial physical mass transfer of CO 2 , heat release from chemical reaction kinetics, and thermophysical property variations from resulting temperature gradients. We demonstrate the applicability of the proposed approach in numerically assessing the performance of packed columns by evaluating key hydrodynamic quantities, CO 2 absorption rates, and temperature rise in a reference column with packings that are structurally similar to the Sulzer Mellapak™ 250.Y packing, for different solvent inflow velocities and temperatures. Predictions from simulation results are found to be consistent with the trends in experimental observations from the literature, suggesting that the predictive capabilities of the simulation framework can be leveraged to guide the future development of absorber-column designs and optimized process flowsheets.
An experimental investigation which seeks the fundamental relationship between the interfacial condensation rate and the parameters which control it when the liquid side is turbulent is discussed. The scaling laws for free-surface condensation are discussed for this case. It is argued that the condensation of cryogenic liquids can, in principle, be simulated in experiments using steam and water. Data are presented for the condensation rate in terms of the dimensionless scaling parameters which involve the fluid properties and the liquid-side turbulence velocity and length scales.
Plasma-based water treatment (PWT) is a promising technology that can degrade various emerging contaminants. However, PWT application on an industrially viable scale is hindered by the lack of an efficient reactor design that combines enhanced plasma-liquid contact with high liquid throughput. This work investigates the applicability of a bubble column gas-liquid contactor to PWT. A pulsed plasma bubble column reactor (PBCR) with a concentric rod-cylinder electrode configuration was used to correlate contaminant removal performance with the gas-liquid contact parameters of the bubble column. A surfactant, rhodamine B dye, and a nonsurfactant, caffeine, were used as model contaminants at µM concentration levels. The bubble column characteristics, i.e., gas holdup, bubble size distribution, and gas-liquid area, were measured as a function of superficial gas velocity using image-based methods. Degradation rates of both contaminants increased with gas flowrate. For caffeine, the increase was attributed to intensified bulk liquid mixing, while dye degradation increased due to the increased gas-liquid area. Ultimately, we show that bubble column contactors significantly improve the utilization of plasma-generated reactive species toward contaminant degradation by distributing them over a large contact area. As a result, a better match between the plasma species interfacial flux and the interfacial contaminant concentration leads to improved treatment energy efficiency. Furthermore, typical degradation energy efficiencies were ~10 g/kWh for caffeine and ~60 g/kWh for rhodamine B.
Neutron reflectometry is a technique for measuring structure near planar interfaces that has been previously used to non-destructively characterize the polymer density of hydrated, dilute, and soft materials. Previous investigations have conducted neutron reflectometry measurements of liquids, gels, emulsion, and polymer solutions at rest, in compression, and subject to shear stress. However, correlating structure with tribological properties of soft materials presents significant experimental challenges for prior instruments due to wall slip, sample thickness, and structural heterogeneity (e.g., depth-wise gradients). A linear reciprocating tribometer offers several advantages for in situ neutron reflectometry studies, including uniform velocity profiles, constant shear stress over large regions of interest, and independent control of normal force and sliding velocity during measurements. This work outlines basic considerations for the design of a custom linear reciprocating tribometer that operates in a neutron beamline and includes commissioning measurements. The tribometer is designed to compress soft and hydrated materials against linearly reciprocating silicon disks. The three key design considerations for this tribometer are (1) safety, (2) neutron transmission, and (3) sample positioning. This instrument design will enable in situ studies of soft matter and illuminate the role of interfacial structure on tribological phenomena.
The major focus of this investigation is to study the fundamentals of layer spreading mechanisms during growth of doped Ge (a facet forming material), and to determine the conditions for morphological instability of vicinal solid-liquid interfaces. The investigation will also lead to the determination of the effect of dopants on the layer growth kinetics, step free energy, and dopant capture by the advancing ledges. The theoretical treatment of growth of faceted interfaces indicates that the kinetics of a step on a growing vicinal interface considerably depends on its angle of inclination, the melt concentration, and characteristics of flow currents in the melt. The morphological stability of the interface also depends on these parameters, as well as on the density and spreading velocity of the steps. However, the treatment of the instability of the interface by the layer growth mechanism is rather difficult because it requires exact knowledge of the thermal and solutal fields, hydrodynamics of the melt, and supercooling at the interface. The results of recent space experiments of the principal investigator involving directional solidification of faceted Bi-Sn alloys have shown that the morphological stability of various crystallographic orientations is significantly affected by the anistropy in interfacial properties of the faceted alloy in general, and the interface kinetics in particular. These findings have also raised many important and fundamental questions, particularly with respect to the behavior of interfacial steps, which need to be addressed via additional groundbased and microgravity experiments. For the present investigation we will use a novel crystal growth technique which provides axial heat flux close to the solid-liquid boundary. The Axial Heat Processing (AHP) technique allows for precise control and determination of the heat and mass transfer close to the crystallization front, and the establishment of a planar interface over the entire cross-section of the growing crystal.
Here, in this work, the design and execution of an experiment with the goal of demonstrating control over the evolution of a copper jet is described. Simulations show that when using simple multi-material buffers placed between a copper target with a conical cavity and a cylinder of high-explosive, a variety of jetting behaviors occur based on material placement, including both jet velocity augmentation and mitigation. A parameter sweep was performed to determine optimal buffer designs in two configurations. Experiments using the optimal buffer designs verified the effectiveness of the buffers at altering jet velocities. Similar trends were shown between the experimental results and the modeling.
The signal feedback process at the edge of the dendrite tip is investigated, and the global instability mechanism of the system is derived. A mechanism is developed to describe a discrete set of unstable global modes for the system. Called WEASR, the method considers the wave emission at the turning point and signal reflections between the turning point and the front edge of the tip. It is shown that the method can obtain the asymptotic solutions for the unstable global modes and the quantum condition for the corresponding eigenvalues. A turning point called the pattern formation condition is shown to be crucial in the formation of dendritic structure and the choice of the final tip velocity. The wave emission process is outlined, and the importance of a signal feedback process at the edge of the dendrite tip is demonstrated. Parameters such as stability and surface tension can be analyzed in terms of their effects on WEASR modes.
Numerous crystal- and microstructural-level mechanisms are at play in the formation of hotspots, which are known to govern high explosives initiation behavior. Most of these mechanisms, including pore collapse, interfacial friction, and shear banding, involve both compressive and shear work done within the material and have thus far remained difficult to separate. Here, we assess hotspots formed at shocked crystal–crystal interfaces using quasi-1D molecular dynamics simulations that isolate effects due to compression and shear. Two high explosive materials are considered (TATB and PETN) that exhibit distinctly different levels of molecular conformational flexibility and crystal packing anisotropy. Temperature and intramolecular strain energy localization in the hotspot are assessed through parametric variation of the crystal orientation and two velocity components that respectively modulate compression and shear work. The resulting hotspots are found to be highly localized to a region within 5–20 nm of the crystal–crystal interface. Compressive work plays a considerably larger role in localizing temperature and intramolecular strain energy for both materials and all crystal orientations considered. Shear induces a moderate increase in energy localization relative to unsheared cases only for relatively weak compressive shock pressures of approximately 10 GPa. These results help isolate and rank the relative importance of hotspot generation mechanisms and are anticipated to guide the treatment of crystal–crystal interfaces in coarse-grained models of polycrystalline high-explosive materials.
Liquid-phase turbulence plays a vital role in determining various gas-liquid two-phase flow parameters, such as void fraction distribution, bubble morphology, bubble-bubble interactions, and interfacial area concentration. In this study, a two-phase flow database including both the gas- and liquid-phases measurements was developed focusing on three bubbly flow conditions in an air-water two-phase flow loop with a vertical one-inch diameter circular pipe test section. A particle image velocimetry (PIV) system integrating an optical phase separation method, i.e., the planar laser-induced fluorescence (PLIF) technique using fluorescent particles and optical filtration, was applied to measure the liquid-phase turbulence information, including the time-averaged velocity, Reynolds stress, and turbulent kinetic energy for the liquid phase. The PIV measurements were taken at three ports along the test section at 14.5, 51.5, and 88.5 pipe inner diameters downstream of a bubble injector. In addition, a double-sensor conductivity probe was used to measure radial distributions of the local time-averaged void fraction and gas velocity. The measured liquid-phase turbulence was used to benchmark Sato's turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. Lastly, in the two bubbly flows tested that have low void fractions being less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction.
Even though slow convective flow is difficult to detect in solidifying metals, it can readily be observed in transparent materials by observing the motion of small neutrally buoyant particles. Succinonitrile, which solidifies with an unfaceted solid/liquid interface and has well characterized physical properties, is considered an excellent material for such studies. For studies of solute-induced convection, ethanol is a useful addition to succinonitrile since it has a lower density and a somewhat similar molecular structure. Samples of high purity and ethanol-doped succinonitrile are unidirectionally solidified in a vertical temperature gradient. Latex mimcrospheres 2 microns in diameter are suspended in the liquid to reveal the convective flow. Convective and morphological stability is observed as a function of solute concentration and growth velocity. These measurements are compared with theoretical calculations that predict the transition from stability to instability as a function of solidification conditions. The predicted transitions occur at low concentrations and solidification velocities; for this reason, extreme care must be taken in order to eliminate the effects of impurities or thermally induced convection.
Wave Energy Converters (WECs) rely on hydraulic Power Take-Off (PTO) systems to transform the kinetic energy of ocean waves into electricity, where hydraulic seals play a critical role in preventing fluid leakage and ensuring watertight operation under harsh marine conditions. However, seal degradation due to wear, fatigue, and corrosion remains a key challenge limiting system reliability and increasing maintenance costs. In this study, the tribological performance of elastomeric materials for hydraulic sealing applications that are based on commercial polyurethane and custom-compounded Ethylene Propylene Diene Monomer (EPDM) is investigated. Using a CSM tribometer, reciprocating wear tests were conducted in accordance with ASTM G133 under dry sliding conditions, with normal loads ranging from 3 to 10 N, and sliding speeds of 10 to 30 mm/s. Coefficient of friction (COF), specific wear rate, and worn surface morphology are evaluated to understand how load, velocity, and material composition affect seal durability. Improved wear stability and high Shore A hardness (93–94) were demonstrated by the commercial thermoplastic polyurethane (Duralast 4203 and 4758). EPDM-based compounds showed different filler-dependent behaviors, formulations filled with carbon black (E3) had the lowest wear rate (1.45 × 10?4 mm³ N?¹ m?¹) and increased load-bearing capacity, while samples modified with silica and plasticizer (E4, E5) had better conformability and damping but less hardness. Compounded EPDMs showed wider deformation zones because of their lower modulus and filler-matrix interfacial fluctuations, while commercial grades exhibited shallower wear scars (= 380 µm), according to optical profilometry. The Archard model-based lifetime estimation showed that, under 10 N and 30 mm/s, operational endurance ranged from 3.1 to 6.2 years. Because of its low wear rate and steady interfacial behavior, E3 had the longest anticipated life. For medium-load hydraulic sealing applications, these results demonstrate the possibility of customized EPDM formulations as affordable substitutes for commercial polyurethanes. They also offer a quantitative foundation for material selection and lifespan prediction.
Bubble column bioreactors are being actively considered for gas fermentation applications, specifically for CO2 utilization, and sugars to fuels conversion. Their main advantages include good mass transfer without any moving parts and low-cost of operation and maintenance. However, the design and scale-up of such reactors is challenging specifically for carbon capture applications where a mixture of gases (e.g. CO2/CO/H2) with variable solubilities is used. The overall performance of scaled-up bioreactors (e.g., mass transfer rate) is largely affected by gas holdup, bubble size distribution (BSD), and multiphase hydrodynamics. We investigate the effect of gas sparger designs on the performance of these large-scale bioreactors using computational fluid dynamics simulations in this work, so as to improve CO2 conversion at scale. The gas distribution systems in bubble column reactors not only determines operational regime, but also affects the evolution of the BSD, which in turn influences interfacial mass transfer and ultimately the efficiency of the gas-liquid exchange process. In addition to the BSD, uniformity in gas sparging affects gas holdup and bubble residence time which constitute important metrics of performance in gas-liquid systems. In this work, we use computational models to simulate high fidelity representations of different sparger designs and their effect on the operation of a bubble column reactor. Four different types of spargers have been selected for the computational study (Fig. 1): ladder, multi-ring, single-ring and toroidal. Their effect on superficial velocity, gas holdup mixing efficiency, and BSD will be evaluated in this work. The model uses a multiphase Eulerian framework similar to [1] and include a composition of mixtures of H2/CO/CO2 gases, common in fermentation applications.
Lowest symmetric and lowest antisymmetric plate wave modes were excited and identified in broad-band acousto-ultrasonic (AU) signals collected from various high temperature composite materials. Group velocities have been determined for these nearly nondispersive modes. An algorithm has been developed and applied to determine phase velocities and hence dispersion curves for the frequency ranges of the broad-band pulses. It is demonstrated that these data are sensitive to changes in the various stiffness moduli of the materials, in agreement by analogy, with the theoretical and experimental results of Tang and Henneke on fiber reinforced polymers. Diffuse field decay rates have been determined in the same specimen geometries and AU configuration as for the plate wave measurements. These decay rates are of value in assessing degradation such as matrix cracking in ceramic matrix composites. In addition, we verify that diffuse field decay rates respond to fiber/matrix interfacial shear strength and density in ceramic matrix composites. This work shows that velocity/stiffness and decay rate measurements can be obtained in the same set of AU experiments for characterizing materials and in specimens with geometries useful for mechanical measurements.
The molecular distribution at the liquid-vapor interface and evolution of the hydrogen bond interactions in mixtures of glycerol and choline chloride are investigated using X-ray photoelectron spectroscopy (XPS) and Raman spectroscopy. Nanoscale depth profiles of supersaturated deep eutectic solvent (DES) mixtures up to ~2 nm measured by ambient-pressure XPS show the enhancement of choline cation (Ch + ) concentration by a factor of 2 at the liquid-vapor interface compared to the bulk. In addition, Raman spectral analysis of a wide range of DES mixtures reveals the conversion of gauche-conformer Ch + into the anti-conformer in relatively lower ChCl concentrations. Finally, the depletion of Ch + from the interface (probing depth = 0.4 nm) is demonstrated by aerosol-based velocity map imaging XPS measurements of glyceline and water mixtures. The nanostructure of liquid-vapor interfaces and structural rearrangement by hydration can provide critical insight into the molecular origin of the deep eutectic behavior and gas-capturing application of DESs.
Back-contacting p-type CdTe has been identified as one of the major areas of loss in CdTe photovoltaic (PV) power conversion efficiency (PCE). In research settings, Au is a common contact material due to its ease of use and decent performance. This work provides a detailed investigation into using gold for back contacting As-doped, CdCl 2 -treated, polycrystalline CdTe that has been exposed to air after absorber processing, another routine practice. First, X-ray photoemission spectroscopy (XPS) is used to determine the native oxide to be 1.6 nm of CdTeO 3 by using a combination of angle-resolved XPS and the cadmium-modified Auger parameter. During gold metallization of CdTe, oxygen and oxidized tellurium are eliminated from the thin CdTeO 3 native oxide. The fates of the released oxygen and possibly cadmium and tellurium are not known, but these reaction byproducts can enter the absorber bulk or grain boundaries, stay at the interface, or dissolve in the Au. Interfacial hole barriers between CdTe and Au are measured for samples with and without the native oxide present prior to metallization. Results show that the thin CdTeO 3 alleviates the downward band bending by 40 meV from 470 to 430 meV even though it is consumed during interface formation. The implications of these chemical reactions on the device are assessed through photoluminescence (PL) spectroscopy, which shows losses in internal open circuit voltage (iVoc) from 820 to 795 meV, carrier lifetime from 123 to 45 ns, and PL quantum yield from 2.9 × 10 −5 to 1.2 × 10 −5 . Modeling timeresolved PL lifetimes demonstrates the back surface recombination velocity due to metallization reduces minority carrier lifetimes. These results identify the native oxide and show that it plays an important role in mediating downward band bending along with how the back interface reaction can negatively impact device-scale parameters and reduce PV PCE.
Thermal insulation materials with ultra-low effective thermal conductivity are crucial for a multitude of applications. Over the years, a significant amount of experimental work has been dedicated to creating new insulation materials with low thermal conductivity. Similarly, substantial efforts have been made to enhance the theoretical understanding of thermal transport mechanisms in thermal insulation materials and to push the boundaries of lower thermal conductivity. However, ample room remains for enhancing the thermal resistivity of closed-cell foam insulations. To aid in the development of ultra-low effective thermal conductivity foam insulations, Oak Ridge National Laboratory has introduced a unique online tool, ThermoPI. This tool calculates not only the effective thermal conductivity but also the thermal conductivity components for porous materials, including gases, solids, and radiative, based on the materials’ structural information (e.g., porosity, pore size, gas species, pressure, solid species, pore geometry, and temperature).This tool collected and improved the existing theoretical models for thermal insulation materials’ gas, solid, and radiative thermal conductivity. These improved models have been validated with experimental data. Effective medium theory and sound velocity softening effects are considered for solid thermal conductivity. The built-in solid materials include polystyrene, polyurethane, polyethylene, and silica. Users are allowed to input new solid materials that are not built inside the tool. For gas thermal conductivity, the subcontinuum Knudsen effect is considered. In addition to several built-in gases, users can input new gases. For solid and radiative thermal conductivity calculations, several models are available to select, and the tool can determine the best model to choose based on the materials that users input. Users are also allowed to change selections manually. In addition, the tool can also calculate the effect of interfacial resistance on the overall thermal conductivity of layered materials. This work will elaborate on the tool and discuss how it can guide the development of new insulation products.
We study the coarsening of the interfacial pattern separating two coexisting phases after the pattern becomes morphologically unstable. Shortly after the instability, the scale of the structure is given by the most unstable wave number in the linear regime. At later times, nonlinear interactions cause the structure to coarsen. Coarsening is studied by monitoring the time dependence of the two-dimensional power spectrum of the pattern, especially for wave vectors transverse to the direction of growth. Characteristic length scales of the pattern obtained from moments of the power spectrum are asymptotically linear in time. Furthermore, the power spectrum is seen to satisfy a scaling relation, in agreement with previous studies. A normal-velocity autocorrelation function is calculated and found to decay substantially over length scales that are of the order of the scale of the pattern. The issue of spatial anisotropy in the correlation functions is also discussed.