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.
Enhanced TcO 4 − reduction by metallic Fe 0 in the presence of particulate and structural Si. Rhythmical precipitation of dissolved iron leads to formation of layered structures related to geological phenomena such as orbicular rocks and Liesegang rings.
Diethylenetriaminepentaacetic acid (DTPA) is a frequently used chelator in the nuclear and medical industries, especially for the complexation of trivalent actinides. However, structural data on these complexes in the solid-state have long remained elusive. Herein, a detailed structural analysis of the presented crystal structures of [C(NH 2 ) 3 ] 4 [Nd(DTPA)] 2 · n H 2 O and [C(NH 2 ) 3 ] 4 [Am(DTPA)] 2 · n H 2 O, where [C(NH 2 ) 3 ] + is guanidinium, details the subtle differences in the Lewis acidity between a lanthanide/actinide pair of similar ionic sizes. Contractions in nitrogen–metal bond lengths between neodymium(III) and americium(III) were observed, while the metal–oxygen bonds remained relatively consistent, highlighting the marginal favorability for actinide complexation over the lanthanides with moderately soft N-donors. Spectroscopic analysis shows significant splitting of many transitions and relatively strong electronic interactions with traditionally low-intensity transitions in the americium complex, as is demonstrated in the 7 F 0 → 7 F 5 transitions. Pressure-induced spectroscopic analysis showed surprisingly little effect on the americium complex, with 5 f →5 f transitions either not shifting or marginally shifting from 2 to 3 nm at 11.93 ± 0.06 GPa─atypical of a soft, N-donor americium complex under pressure. Finally, large voids occupied by water molecules in between the complexes within the crystal structure may be responsible for the lack of pressure response in the 5 f →5 f transitions.
Observations collected during cold-season precipitation periods at Utquagvik, Alaska and at the multidisciplinary drifting observatory for the study of Arctic climate (MOSAiC) are used to statistically analyze the relations among the atmospheric water cycle parameters including the columnar supercooled liquid and ice amounts (expressed as liquid-water and ice-water paths, i.e., LWP and IWP), the integrated water vapor (IWV) and the near-surface snowfall rate. Data come from radar and radiometer-based retrievals and from optical precipitation sensors. While the correlation between snowfall rate and LWP is rather weak, correlation coefficients between radar-derived snowfall rate and IWP are high (~0.8), which is explained, in part, by the generally low LWP/IWP ratios during significant precipitation. Correlation coefficients between snowfall rate and IWV are moderate (~0.45). Correlations are generally weaker if snowfall is estimated by optical sensors, which is, in part, due to blowing snow. Correlation coefficients between near-surface temperature and snowfall rates are low (r < 0.3). The results from the Alaska and MOSAiC sites are generally similar. These results are not very sensitive to the amount of time averaging (e.g., 15 min averaging versus daily averages). Observationally based relations among the water cycle parameters are informative about atmospheric moisture conversion processes and can be used for model evaluations.
We study the initial stages of homogeneous melting of a hexagonal ice crystal at coexistence and at moderate superheating. Our trajectory-based computer simulation approach provides a comprehensive picture of the events that lead to melting, from the initial accumulation of 5+7 defects, via the formation of L-D and interstitial-vacancy pairs, to the formation of a liquid nucleus. Of the different types of defects that we observe to be involved in melting, a particular kind of 5+7 type defect (type 5) plays a prominent role as it often forms prior to the formation of the initial liquid nucleus and close to the site where the nucleus forms. Hence, like other solids, ice homogeneously melts via the prior accumulation of defects.
Hyperstoichiometric UO 2.15 was characterized by neutron total scattering at high temperature in the single-phase UO 2+x region of the U/O phase diagram. The diffraction data confirmed a single-phase fluorite structure at high temperature. Analysis of the short-range data showed that the same structural model does not fit the pair distribution functions well. Instead, structural models containing specific configurations of oxygen defect clusters best represent the local atomic arrangement. In conclusion, prevalent defect clusters previously proposed were fit to the experimental data, and moderately distorted oxygen cuboctahedra hypothesized by recent molecular dynamics simulations fit the data most accurately.
Clathrate hydrates are crystalline solids characterized by their ability to accommodate large quantities of guest molecules. Although CH 4 and CO 2 are the traditional guests found in natural systems, incorporating smaller molecules (e.g., H 2 ) is challenging due to the need to apply higher pressures to stabilize the hydrogen-bonded network. Another critical limitation of hydrates is the slow nucleation and growth kinetics. Here, we show that specially designed activated carbon materials can surpass these obstacles by acting as nanoreactors promoting the nucleation and growth of H 2 hydrates. The confinement effects in the inner cavities promote the massive growth of hydrogen hydrates at moderate temperatures, using pure water, with extremely fast kinetics and much lower pressures than the bulk system.
Hybrid density functional theory (DFT) remains intractable for large periodic systems due to the demanding computational cost of exact exchange. We apply the tensor hypercontraction (THC) (or interpolative separable density fitting) approximation to periodic hybrid DFT calculations with Gaussian-type orbitals using the Gaussian plane wave approach. This is done to lower the computational scaling with respect to the number of basis functions (N) and k-points (N k ) at a fixed system size. Additionally, we propose an algorithm to fit only occupied orbital products via THC (i.e., a set of points, N ISDF ) to further reduce computation time and memory usage. This algorithm has linear scaling cost with k-points, no explicit dependence of N ISDF on basis set size, and overall cubic scaling with unit cell size. Significant speedups and reduced memory usage may be obtained for moderately sized k-point meshes, with additional gains for large k-point meshes. Adequate accuracy can be obtained using THC-oo-K for self-consistent calculations. In conclusion, we perform illustrative hybrid density function theory calculations on the benzene crystal in the basis set and thermodynamic limits to highlight the utility of this algorithm.
Detailed ground-based observations of snow are scarce in remote regions, such as the Arctic. Here, Multi-Angle Snowflake Camera measurements of over 55,000 solid hydrometeors—obtained during a two-year period from August 2016 to August 2018 at Oliktok Point, Alaska—are analyzed and compared to similar measurements from an earlier experiment at Alta, Utah. In general, distributions of hydrometeor fall speed, fall orientation, aspect ratio, flatness, and complexity (i.e., riming degree) were observed to be very similar between the two locations, except that Arctic hydrometeors tended to be smaller. In total, the slope parameter defining a negative exponential of the size distribution was approximately 50% steeper in the Arctic as at Alta. Sixty-six percent of particles were observed to be rimed or moderately rimed with some suggestion that riming is favored by weak boundary layer stability. On average, the fall speed of rimed particles was not notably different from aggregates. However, graupel density and fall speed increase as cloud temperatures approach the melting point.
Ion-containing block copolymers (BCPs) are promising solid-state electrolytes as they can combine ion conduction and mechanical robustness via multiple microphases. However, salt-doped BCPs can have low ion conduction and transference number. One strategy to increase transference number is to tether anions to the polymer backbone (creating a single-ion conductor), though this may impact cation mobility. Leveraging coarse-grained molecular dynamics simulations, we compare the effects of ion concentration and polymer dielectric constant on BCPs with and without bonding anions to the chain (without changing any other chemical properties). We find that single-ion BCPs have higher cation conductivity than analogous salt-doped BCPs, especially at moderate ion concentrations. Furthermore, we also compare several structural features of these systems; notably, tethering anions causes the ions to spread more evenly throughout the conducting phase, which is related to improved cation conduction.
Fin-and-tube heat exchangers (FTHEs) are widely used for high-temperature flue-gas heat recovery, but their performance is often limited by wake regions and non-uniform fin-surface temperatures. This study proposes and numerically evaluates four bio-inspired longitudinal vortex generator (VG) configurations in a high-temperature FTHE with flue-gas inlet temperature ∼1230 K: double-delta, curved double-delta, alula, and a new curved-alula geometry. The reference fin is not hydraulically plain; it already incorporates leading-edge separation columns and convex protrusions, so the alula-type winglets are assessed as downstream add-ons acting on a strongly disturbed flow. In a second step, perforations (one, two and three circular holes) are introduced into the curved-alula VGs to further tailor the flow field. Three-dimensional simulations with the Shear Stress Transpor (SST) $k - ω$ model, temperature-dependent flue-gas properties and conjugate conduction are carried out for gas-side Reynolds numbers $Re_g ≈ 8.0$ x $10^2 - 3.6$ x $10^3$ (mass flow rates 0.5 – 2.5 g/s), and the designs are compared in terms of surface heat flux, Nusselt number, friction factor and hydrothermal performance factor (HTPF). For this already-promoted fin, the additional downstream winglets provide moderate, incremental hydrothermal gains. At the highest Reynolds number, the best non-perforated design (curved-alula) increases surface heat flux from 1630.9 to 1794.7 kW/m² (∼ 10 % gain) and the Nusselt number from 227.6 to 242.6 (∼ 7 % gain), while the friction factor rises from 0.26 to about 0.30, yielding HTPF values close to unity (∼ 0.9 – 1.0). Introducing circular perforations into the curved-alula winglets acts mainly as a wake-bleeding refinement: the three-hole configuration provides a heat flux of 1824.7 kW/m² and a pressure drop of 127.9 Pa, with HTPF in the range ∼ 1.03 – 1.14 and a small (∼ 1 – 3 %) improvement over the solid curved-alula design. Flow-field analysis shows that the perforated curved-alula VGs shrink tube-wake regions, thin the thermal boundary layer and homogenize the fin-surface temperature (outlet-gas temperature ∼ 510 – 520 K and fin-surface temperature ∼ 420 – 421 K for the three-hole case). An optimal flue-gas mass flow rate of ∼ 1 g/s ($Re_g ≈ 1.5$ x $10^3$) is identified, beyond which additional heat-transfer gains are offset by rapidly increasing pressure losses. Overall, the results highlight that initial fin geometry and VG placement are as important as VG shape: alula-based winglets are expected to yield larger relative gains on simpler flat-fin layouts or when positioned closer to the fin leading edge and tube
Abstract Due to its moderate strength (≈700 MPa) and impressive extensibility before breaking (≈60–80%), orb‐weaving spider aciniform (AC) prey‐wrapping silks are actually the toughest of the spider silks but are remarkably understudied. The previous results indicate that native AC silk fibers are an α‐helix rich coiled‐coil/β‐sheet hybrid nanofiber, and that conversion of disordered or helical domains to β‐sheet aggregates is surprisingly minimal and overall β‐sheet content is low (≈15%). In this work, it is demonstrated through scanning electron microscopy that native AC silk fibers undergo matted cross‐linking upon exposure to moisture that increases silk stiffness. The unique molecular mechanism of water‐induced cross‐linking is revealed with solid‐state NMR (SSNMR) methods; water‐induced morphological changes are correlated with an increase in AC silk protein β‐sheet content, and additionally a minor unfolding of coiled‐coil regions is observed. Continued and increased β‐sheet cross‐linking is observed upon application of mechanical shear. The size of these β‐sheet domains to be 4–6 nm using Wide‐Line Separation SSNMR is determined. The observation that merely water treatment can be used to convert a protein‐based material from a flexible/extensible α‐helix‐rich fiber to a rigid crossed‐linked β‐sheet mat is a novel observation that should provide new avenues in bioinspired materials design.
Since solids are only sometimes seen en masse in a pure bulk form, and for liquids other than water almost never, a capability to model reflectance spectra from analytes deposited on various substrates would be highly advantageous. If available, the real, n ( ν ) , and imaginary, k ( ν ) , components of the complex refractive index, n n + i k , can be used to simulate infrared spectra, accounting for reflection, refraction, and absorption phenomena as a function of wavelength. We focus on using the Pacific Northwest National Laboratory (PNNL) derived n / k vectors for solid and liquid analytes deposited as thin layers on different types of substrates including conductors, such as aluminum, and inorganic dielectrics, such as glass. The model is an adaptation of the Monte Carlo ray trace modeling program, TracePro, extended through the use of its macrolanguage. The model is tested using thin films of organic liquids including silicone oil and no. 2 diesel fuel, as well as organic solids such as caffeine and acetaminophen on aluminum and glass. The predicted spectra for the solid films were compared to experimental hemispherical reflectance data measured using a Fourier transform spectrometer with an integrating sphere. The thickness of the calculated layer is a parameter for predicting the (transflectance) spectra and is obtained using the areal density measured from gravimetric methods to generate the thin-layer samples. Comparison of the calculated spectra with experimental hemispherical reflectance data shows excellent agreement, indicating promise for the use of measured n / k data to synthesize reference spectral data. In particular, accounting for the inhomogeneity of the deposits greatly improved the match with experimental data. Finally, the theoretical modeling shows that for thicker layers (ca. 20 to 100 μm) of typical organics possessing moderately strong k values, the longwave infrared features are often saturated and better spectral contrast is obtained from the overtone/combination bands in the shortwave infrared.
The formation of ice, which plays an important role in diverse contexts ranging from cryopreservation to atmospheric science, is often mediated by solid surfaces. Although surfaces that interact favorably with ice (relative to liquid water) can facilitate ice formation by lowering nucleation barriers, the molecular characteristics that confer ice-philicity to a surface are complex and incompletely understood. To address this challenge, here we introduce a robust and computationally efficient method for characterizing surface ice-philicity that combines molecular simulations and enhanced sampling techniques to quantify the free energetic cost of increasing surface–ice contact at the expense of surface–water contact. Using this method to characterize the ice-philicity of a family of model surfaces that are lattice matched with ice but vary in their polarity, we find that the nonpolar surfaces are moderately ice-phobic, whereas the polar surfaces are highly ice-philic. In contrast, for surfaces that display no complementarity to the ice lattice, we find that ice-philicity is independent of surface polarity and that both nonpolar and polar surfaces are moderately ice-phobic. Furthermore, our work thus provides a prescription for quantitatively characterizing surface ice-philicity and sheds light on how ice-philicity is influenced by lattice matching and polarity.
The chalcogenide perovskite family has been steadily gaining increasing attention from the research community due to its optoelectronic properties and potential for diverse applications. While BaZrS 3 and BaTiS 3 have been the most extensively studied, other promising compounds in this family, such as Sr x TiS 3 (1.05 < x < 1.22), are now being explored for various optical, optoelectronic, and energy storage applications. However, challenges remain in achieving the low-temperature synthesis of Sr x TiS 3 . In this study, we report, for the first time, the synthesis of Sr x TiS 3 nanocrystals at temperatures below 400 °C. The synthesized nanocrystals exhibit a rod-like morphology. Additionally, we have developed solution-processing routes to synthesize phase-pure Sr x TiS 3 thin films, marking the first reported instance of such films, at temperatures below 600 °C. We also demonstrate the solid-state synthesis of Sr x TiS 3 powder below 600 °C. Our work paves the way for new and exciting application avenues for Sr x TiS 3 .
In amorphous solids, soft vibrational modes derived from normal mode analysis can be utilized to identify the soft particles that are prone to irreversible rearrangements. However, the normal mode analysis approach cannot explain why the spatial distributions of clustered soft particles do not change over time. We define a softness parameter based on the vibrational density of states calculated directly from molecular dynamics simulations with both the harmonic vibrations and anharmonic relaxations being properly captured at finite temperatures. This parameter spontaneously correlates with the real space atomic motions and the dynamics heterogeneity. Using the softness parameter, we show that the softest particles are confined within rigid cages. These particles keep rearranging reversibly without long-range diffusion. The moderately soft particles rearrange irreversibly, and the hard particles mainly participate in vibrations without rearrangement. Here, we also show that the soft particles form locally disordered structures, while the hard particles present strong ordering. These findings confirm the defective nature of soft particles, and provide insights on the nature of softness as the ability to rearrange, but not necessarily irreversibly.
The necessity to reduce global energy consumption calls for innovative strategies in building thermal management. Passive thermal regulation, particularly through bio-inspired designs, offers a promising avenue by mimicking nature's efficient control of optical properties. This research introduces a novel, climate-responsive coating that integrates optimized bio-inspired surface features with a solid-solid phase change material (SS-PCM) to dynamically manage solar absorptivity without adding additional thickness, enabling both heating and cooling as needed. Drawing on the photonic architectures of the Saharan silver ant and Morpho Didius butterfly, we employed a modeling and multi-objective optimization framework to tailor these surface features. Simulations reveal that surface texture, rather than the intrinsic phase transition of the SS PCM, dominates optical control. Relative to a flat SS PCM coating, optimized isotropic random roughness and broader range features yielded the highest passive heating power increase of about 144 % and 319 % respectively suitable for cold climates. Saharan ant-inspired features enhanced passive cooling for hot climates, achieving a 21.8 % improvement. For moderate climates, Butterfly-wing-inspired surface features provided a balanced enhancement of 19 % for heating and 7 % for cooling. Across all cases, the optimized surface features reduced combined heating and cooling energy demand more effectively than the baseline coating, while preserving material thickness. These findings demonstrate that climate-adaptive, optimized bio-inspired surface features can unlock the full potential of SS PCM coatings, providing a versatile pathway to significant energy savings in buildings and other applications. The methodology establishes a framework for designing next-generation adaptive envelopes that leverage natural photonic principles for high-impact, low-cost thermal regulation.
The use of oxide cathodes in Mg batteries would unlock a potential energy storage system that delivers high-energy density. However, poor kinetics of Mg diffusion in known solid oxide lattices strongly limits reversible intercalation, which motivates the sustained exploration of new candidates. Herein, nanocrystals of a few-layer phyllomanganate, reminiscent of the mineral vernadite, were shown to have considerable electrochemical activity toward Mg intercalation at room temperature, where it delivered similar to 190 mAh g- 1 at similar to 1.9 V (vs Mg/Mg 2+ ) in batteries paired with a Mg metal anode. Multimodal characterization confirmed the notable degree of reversible intercalation by probing the structural, compositional, and redox changes undertaken by the oxide. Distinct levels of Mg activity were also observed while varying the content of small amounts of lattice water and the temperature of the reaction. The results reaffirm the prospects for operational Mg batteries using oxide cathode in moderate conditions, overcoming current limits of performance of this prospective technology.
The structures adopted by solids under pressure are often assumed to reflect thermodynamic equilibrium, yet in many materials phase selection is strongly influenced by kinetic pathways and microstructural inheritance. Elemental cerium (Ce) exemplifies this challenge, with decades of conflicting reports describing different high-pressure crystal structures emerging under nominally identical conditions. Here we use neutron diffraction from large ( ~ 60 mm 3 ) sample volumes to follow the structural evolution in ultra-high-purity Ce during controlled pressure-temperature cycling between 85 and 295 K and up to 8 GPa. We find that the crystal structure formed at high pressure depends on the compression pathway: slow compression ( ~ 0.25 GPa hr −1 ) at room temperature favors an orthorhombic phase (α'), whereas slow ( ~ 0.25 GPa hr −1 ) and also moderately faster ( ~ 0.5 GPa hr −1 ) compression at low temperature stabilizes a pure monoclinic phase (α"). The low-temperature phase persists metastably over a wide temperature range but transforms irreversibly upon heating above ~ 280 K, or modest pressure cycling. Remarkably, the lower-pressure γ phase remains trapped far beyond the expected stability field, persisting to the highest pressures studied. These observations show that phase selection in Ce is governed by kinetics and microstructural memory rather than equilibrium thermodynamics or sample morphology alone, establishing path dependence as a defining feature of its high-pressure behavior.