Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “layered solid element”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The role of alkali metal exchange in zeolite-templated carbon synthesis

Despite great efforts to achieve ideal atomistic packing of carbon in the pore networks of even the largest pore zeolites, templating inaccuracies are ubiquitous in all presently synthesized zeolite-templated carbon (ZTC) materials. Idealized models suggest that the long-sought schwarzite-like structures could in principle be obtained by zeolite templating if the appropriate zeolite template were chosen (e.g., faujasite for D surface schwarzites) and if perfect template fidelity (insertion of a pristine layer of pure carbon directly on the surface of the zeolite) could be achieved. A requirement to achieve such structures is increased carbon density within the zeolite. We report the investigation of a series of alkali metal cation-exchanged zeolites to determine how the periodic trends in the group 1 elements influence zeolite templating, with a specific focus on the metric of structural packing density (SPD) as resolved by ex situ thermogravimetry. In a survey based on controlled synthesis temperature, time, and flow conditions, an increasing SPD was observed with decreasing cation size, an effect that is consistent with the increasing strength of cation-π interactions. Furthermore, this effect could be promising for future work to increase the SPD of ZTCs for the synthesis of closed-tube, schwarzite-like carbonaceous solids.

03 NATURAL GAS↗
NASA NTRS2018

New Formation Models for the Kepler-36 System

Formation of the planets in the Kepler-36 system is modeled by detailed numerical simulations according to the core-nucleated accretion scenario. The standard model is updated to include the dissolution of accreting rocky planetesimals in the gaseous envelope of the planet, leading to substantial enrichment of the envelope mass in heavy elements and a non-uniform composition with depth. For Kepler-36 c, models involving in situ formation and models involving orbital migration are considered. The results are compared with standard formation models. The calculations include the formation (accretion) phase, as well as the subsequent cooling phase, up to the age of Kepler-36 (7 Gyr). During the latter phase, mass loss induced by stellar XUV radiation is included. In all cases, the results fit the measured mass, 7.84 M, and radius, 3.68 R, of Kepler-36 c. Two parameters are varied to obtain these fits: the disk solid surface density at the formation location, and the "efficiency" factor in the XUV mass loss rate. The updated models are hotter and therefore less dense in the silicate portion of the planet and in the overlying layers of H/He, as compared with standard models. The lower densities mean that only about half as much H/He is needed to be accreted to fit the present-day mass and radius constraints. For Kepler-36 b, an updated in situ calculation shows that the entire H/He envelope is lost, early in the cooling phase, in agreement with observation.

planets and satellites:formation↗
DOE OSTI2025

The Origin of Improved Performance in Boron‐Alloyed Silicon Nanoparticle‐Based Anodes for Lithium‐Ion Batteries

Stabilizing the solid electrolyte interphase (SEI) remains a key challenge for silicon‐based lithium‐ion battery anodes. Alloying silicon with secondary elements like boron has emerged as a promising strategy to improve the cycle life of silicon anodes, yet the underlying mechanism remains unclear. To address this knowledge gap, how boron concentration influences battery performance is systematically investigated. These results show a near‐monotonic increase in cycle lifetime with higher boron content, with boron‐rich electrodes significantly outperforming pure silicon. Additionally, silicon‐boron alloy anodes exhibit nearly three times longer calendar life than pure silicon. Through detailed mechanistic analysis, alternative contributing factors are systematically ruled out, and it is proposed that improved passivation arises from a strong permanent dipole at the nanoparticle surface. This dipole, formed by undercoordinated and highly Lewis acidic boron, creates a static, ion‐dense layer that stabilizes the electrochemical interface, reducing parasitic electrolyte decomposition and enhancing long‐term stability. These findings suggest that, within the SEI framework, the electric double layer is an important consideration in surface passivation. This insight provides an underexplored parameter space for optimizing silicon anodes in next‐generation lithium‐ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2025

Anodized Aluminum Oxide Membrane Ionic Memristors

Memory effect in ion transport (IT) at the solid–solution interface is uniquely attractive in that the conductance depends on or “memorizes” the previous states. Hysteretic and rectified transport properties offer exciting potential to developing advanced iontronics and neuromorphic functions, improving the efficiency of energy conversion and electrochemical processes, and overcoming the selectivity-throughput bottleneck in the enrichment of low abundant species for environment- and energy-friendly separations, among others. Herein, memory effects are discovered in the rectified electrokinetic IT through anodized aluminum oxide (AAO) membranes containing densely packed highly ordered nanochannels (10 10 per cm 2 ). Characteristic memristor responses of pinched current–potential loops are resolved in voltammetric experiments and successfully reproduced through finite element simulation. Excitatory and inhibitory conductance states are shown to arise from the enrichment and depletion of mobile charge carriers. Structurewise, the transport symmetry is broken by the barrier oxide layer (BOL) on the one end of the cylindrical nanochannels across the AAO membranes. Charge selectivity is attributed to the gradient(s) of the space charge density across the BOL characterized by depth profiling via X-ray photoelectron spectroscopy analysis. The space charge gradient(s) overcomes the fundamental limitation of widely exploited surface charge effects to enable intense rectification and hysteresis prevailing at very high ionic concentrations up to 1–2 M. A new strategy is developed for controlling the preferential IT direction and selectivity via counterion intercalation and extraction/exchange. Mechanistic understanding is further confirmed through parameter variations such as potential scan rate and ionic strength, which also demonstrates convenient controls of the related functions.

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

An finite element analysis surrogate model with boundary oriented graph embedding approach for rapid design

Abstract In this work, we present a boundary oriented graph embedding (BOGE) approach for the graph neural network to assist in rapid design and digital prototyping. The cantilever beam problem has been solved as an example to validate its potential of providing physical field results and optimized designs using only 10 ms. Providing shortcuts for both boundary elements and local neighbor elements, the BOGE approach can embed unstructured mesh elements into the graph and performs an efficient regression on large-scale triangular-mesh-based finite element analysis (FEA) results, which cannot be realized by other machine-learning-based surrogate methods. It has the potential to serve as a surrogate model for other boundary value problems. Focusing on the cantilever beam problem, the BOGE approach with 3-layer DeepGCN model achieves the regression with mean square error (MSE) of 0.011 706 (2.41% mean absolute percentage error) for stress field prediction and 0.002 735 MSE (with 1.58% elements having error larger than 0.01) for topological optimization. The overall concept of the BOGE approach paves the way for a general and efficient deep-learning-based FEA simulator that will benefit both industry and Computer Aided Design (CAD) design-related areas.

42 ENGINEERING↗
DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗
NASA NTRS1990

Vertical bloch line memory

Vertical Bloch Line (VBL) memory is a recently conceived, integrated, solid-state, block-access, VLSI memory which offers the potential of 1Gbit/sq cm real storage density, gigabit per second data rates, and sub-millisecond average access times simultaneously at relatively low mass, volume, and power values when compared to alternative technologies. VBL's are micromagnetic structures within magnetic domain walls which can be manipulated using magnetic fields from integrated conductors. The presence or absence of VBL pairs are used to store binary information. At present, efforts are being directed at developing a single-chip memory using 25Mbit/sq cm technology in magnetic garnet material which integrates, at a single operating point, the writing, storage, reading, and amplification functions needed in a memory. This paper describes the current design architecture, functional elements, and supercomputer simulation results which are used to assist the design process. The current design architecture uses three metal layers, two ion implantation steps for modulating the thickness of the magnetic layer, one ion implantation step for assisting propagation in the major line track, one NiFe soft magnetic layer, one CoPt hard magnetic layer, and one reflective Cr layer for facilitating magneto-optic observation of magnetic structure. Data are stored in a series of elongated magnetic domains, called stripes, which serve as storage sites for arrays of VBL pairs. The ends of these stripes are placed near conductors which serve as VBL read/write gates. A major line track is present to provide a source and propagation path for magnetic bubbles. Writing and reading, respectively, are achieved by converting magnetic bubbles to VBL's and vice versa. The output function is effected by stretching a magnetic bubble and detecting it magnetoresistively. Experimental results from the past design cycle created four design goals for the current design cycle. First, the bias field ranges for the stripes and the major line needed to be matched. Second, the magnetic field barrier between the stripe and the read/write gates needed to be reduced. Third, current conductor routing needed to be improved to reduce occurrences of open-circuiting, short-circuiting, and eddy-current shielding. Fourth, a modified Co-alloy was needed with an increased coercivity and controlled magnetization to allow VBL stabilization to occur without affecting stripe stability.

Katti, R.↗
NASA NTRS1984

Surface sensitive photoemission EXAFS: A new structural tool

A new technique for surface EXAFS has been developed. In this technique the absorption fine structure is observed by monitoring the photoelectrons emitted from core levels as a function of photon energy. The surface sensitivity is achieved by detecting elastic photoelectrons that escape to the vacuum from a few atomic layers near the surface. EXAFS from Mn 3p and F2s outer core levels in MnF2 were measured, the results were compared with transmission EXAFS from Mn ls core level and good agreement was found. This development is significant because: (1) almost all elements have suitable core levels of low binding energy; (2) since photopeaks from adjacent core levels do not cross as the photon energy is varied, the limitations imposed on other surface EXAFS techniques by this interference is removed. Auger electrons, which have fixed kinetic energies and so can cross photopeaks, will generally have low energies if they originate from outer core levels and be unlikely to appear in the EXAFS region of kinetic energies, i.e., in excess of about 50 eV; (3) VUV light is used, which greatly extends the usable range of photon energies; (4) more than one element in a solid may be studied with the same monochromator. The experiments were carried out on the VUV beam line U14A at NSLS.

Choudhary, K. M.↗
NASA NTRS2004

A Critical Review of the Transport and Decay of Wake Vortices in Ground Effect

This slide presentation reviews the transport and decay of wake vortices in ground effect and cites a need for a physics-based parametric model. The encounter of a vortex with a solid body is always a complex event involving turbulence enhancement, unsteadiness, and very large gradients of velocity and pressure. Wake counter in ground effect is the most dangerous of them all. The interaction of diverging, area-varying, and decaying aircraft wake vortices with the ground is very complex because both the vortices and the flow field generated by them are altered to accommodate the presence of the ground (where there is very little room to maneuver) and the background turbulent flow. Previous research regarding vortex models, wake vortex decay mechanisms, time evolution within in ground effect of a wake vortex pair, laminar flow in ground effect, and the interaction of the existing boundary layer with a convected vortex are reviewed. Additionally, numerical simulations, 3-dimensional large-eddy simulations, a probabilistic 2-phase wake vortex decay and transport model and a vortex element method are discussed. The devising of physics-based, parametric models for the prediction of (operational) real-time response, mindful of the highly three-dimensional and unsteady structure of vortices, boundary layers, atmospheric thermodynamics, and weather convective phenomena is required. In creating a model, LES and field data will be the most powerful tools.

Sarpkaya, T.↗
DOE OSTI2023

Asymmetric Lithium Extraction and Insertion in High Voltage Spinel at Fast Rate

Spinel-structured ordered-LiNi 0.5 Mn 1.5 O 4 (o-LNMO) has experienced a resurgence of interest in the context of reducing scarce elements such as cobalt from lithium-ion battery. O-LNMO undergoes two two-phase reaction at slow rates. However, it is not known if such phenomenon also applies at fast rates. Herein, we investigate the rate-dependent phase transition behavior of o-LNMO through in operando time-resolved X-ray diffraction. The results indicate that a narrow region of solid solution reaction exists for charge and discharge at both slow and fast rate. The overall phase transition is highly asymmetric at fast rates. During fast charge, it is a particle-by-particle mechanism resulting from asynchronized reaction among the particles. Finally, during fast discharge, it is likely a core-shell mechanism involving transition from Li 0+x Ni 0.5 Mn 1.5 O 4 to Li 1+x Ni 0.5 Mn 1.5 O 4 in the outer layer of particles. The Li 0.5 Ni 0.5 Mn 1.5 O 4 phase is suppressed during fast discharge and only appears through Li redistribution upon relaxation.

25 ENERGY STORAGE↗
NASA NTRS2001

The CE/SE Method: a CFD Framework for the Challenges of the New Millennium

The space-time conservation element and solution element (CE/SE) method, which was originated and is continuously being developed at NASA Glenn Research Center, is a high-resolution, genuinely multidimensional and unstructured-mesh compatible numerical method for solving conservation laws. Since its inception in 1991, the CE/SE method has been used to obtain highly accurate numerical solutions for 1D, 2D and 3D flow problems involving shocks, contact discontinuities, acoustic waves, vortices, shock/acoustic waves/vortices interactions, shock/boundary layers interactions and chemical reactions. Without the aid of preconditioning or other special techniques, it has been applied to both steady and unsteady flows with speeds ranging from Mach number = 0.00288 to 10. In addition, the method has unique features that allow for (i) the use of very simple non-reflecting boundary conditions, and (ii) a unified wall boundary treatment for viscous and inviscid flows. The CE/SE method was developed with the conviction that, with a solid foundation in physics, a robust, coherent and accurate numerical framework can be built without involving overly complex mathematics. As a result, the method was constructed using a set of design principles that facilitate simplicity, robustness and accuracy. The most important among them are: (i) enforcing both local and global flux conservation in space and time, with flux evaluation at an interface being an integral part of the solution procedure and requiring no interpolation or extrapolation; (ii) unifying space and time and treating them as a single entity; and (iii) requiring that a numerical scheme be built from a nondissipative core scheme such that the numerical dissipation can be effectively controlled and, as a result, will not overwhelm the physical dissipation. Part I of the workshop will be devoted to a discussion of these principles along with a description of how the ID, 2D and 3D CE/SE schemes are constructed. In Part II, various applications of the CE/SE method, particularly those involving chemical reactions and acoustics, will be presented. The workshop will be concluded with a sketch of the future research directions.

Chang, Sin-Chung↗
DOE OSTI2025

Crystal structure of Cu2Zn(GexSi1−x)Se4 solid solution: the kesterite to wurtz–kesterite structural phase transition

Developing low-cost, sustainable, and environmentally friendly top absorber layers for tandem solar cells is essential to advancing photovoltaic technologies and accelerating the transition to renewable energy. In this work, we explore the potential of tetravalent (Cu2Zn(GexSi1−x)Se4) cation mutations in chalcogenide compound semiconductors with the aim of finding a material with increased band gap and reduced structural disorder. A combination of high-resolution synchrotron powder diffraction and neutron powder diffraction was used to determine the atomic positions and monoclinic angles in monoclinic wurtz–kesterite type Cu2Zn(GexSi1−x)Se4 mixed crystals as well as to determine the cation distribution in the crystal structure of Ge-rich kesterite-type and Si-rich wurtz–kesterite type mixed crystals. These investigations enabled us to deduce the structural transition scenario within the Cu2Zn(GexSi1−x)Se4 series. The transition occurs via a region where two phases with different crystal structures, tetragonal and monoclinic and thus a different distortion of the coordination tetrahedra, but the same cation distribution within the element specific cation sites co-exist. Thus, the structural transition between the kesterite and the wurtz–kesterite structure within the Cu2Zn(GexSi1−x)Se4 series is a distortion driven transition. The study identifies cation mutation in quaternary chalcogenides as a promising strategy beyond chalcopyrites and kesterites for low cost and environmentally friendly top absorbers in tandem solar cells.

Gurieva, Galina [Helmholtz Center Berlin for Mater↗
DOE OSTI2024

Reactive flash sintering and characterization of bulk high entropy nitrides

Over the past decade, numerous high-entropy ceramics have been synthesized, often displaying attractive properties. However, the study on facile preparation of bulk high entropy nitrides (HEN) are limited, despite its broad potential applications. This research demonstrates for the first time rapid fabrication (within ∼6 min) of bulk high-entropy nitrides, especially (Al 0.17 Nb 0.17 Ta 0.17 Ti 0.32 Zr 0.17 )N, from binary nitride powder mixtures using a highly efficient reactive flash sintering (RFS) technique. X-ray diffraction (XRD) shows the HENs from RFS are near single-phase solid solutions with a rock salt crystal structure, while in situ synchrotron study carried out during RFS captured in real time the formation of HEN, which was preserved upon cooling, suggesting thermodynamic stability of the HEN phase, even up to extreme pressure (∼35.6 GPa). Microscopic analyses using SEM, STEM, and EDS reveal decent uniformity for HEN with no obvious segregation of elements, even to submicron scale. Some properties of the obtained bulk HENs are consistent with expectations. For example, their hardness and bulk modulus are close to estimates based on rule-of-mixture (ROM) values from the constituent binary nitrides. Meanwhile, some other measured properties seem to show surprises. For example, the fracture toughness for the HENs (e.g., 7.81 ± 1.40 MPa•m 1/2 or higher) turns out to be more than double of the expected ROM estimates. The significantly improved fracture toughness is attributed to the observed nano-layered structure of the HENs, despite the HEN’s cubic crystal structure and high hardness. In addition, the oxidation resistance shows improvement up till ∼800°C, possibly due to Ta doping that suppress oxygen vacancy formation in the oxide shell, while the 5-metal HEN of (Al 0.17 Nb 0.17 Ta 0.17 Ti 0.32 Zr 0.17 )N displays superconductivity (T c of ∼5–7 K from magnetism and resistivity measurements, slightly lower than ROM estimate), despite insulating property of starting AlN. Furthermore, future study combining experimental investigation using larger samples to confirm the observed increase in fracture toughness and oxidation resistance, theoretical modeling at different length scale, and more detailed structural/chemical characterization, especially at the atomic scale, are all needed to fully understand the inter-relationships between composition, processing, structure, and novel properties for these HENs and the development of related new materials for different applications.

Flash sintering↗
DOE OSTI2022

Metal fluxes at the sediment–water interface in a free water surface constructed wetland

Metal biogeochemistry in the sediment, water, and the sediment–water interface (SWI) was studied in a free water surface constructed wetland. Concentrations of labile copper (Cu), zinc (Zn), sulfate, chloride, and dissolved organic carbon (DOC) were measured with the diffusive gradients in thin film (DGT) and peeper. A good agreement between peeper- and DGT-measured metals was observed for Cu (regression r 2 = 0.3, 95% CI of the slopes > 0), but not for Zn (95% CI of the slopes overlapped with 0), which was attributed to the different complexed compounds between Cu and Zn in porewater. The depth profile of labile metals in sediment porewater varied with time and was consistent with the solid-phase metal deposition, showing higher concentrations in the surface layer (3 to –3 cm) than in the bottom layer (–4 to –13 cm). The depth-averaged labile Cu and Zn concentrations measured by DGT were 1.0 and 3.1 µg/L, and labile sulfate, chloride, and DOC concentrations measured by peeper were 1.8, 3.6, and 2.1 mg/L, respectively. A sharp decrease in sulfate occurred in September when sulfate concentrations became the lowest among sampling months. This was caused by the seasonal sulfur cycles in the wetland, where the dominant sulfur reaction is sulfate reduction in warm seasons and sulfide oxidation in cold seasons. Different metal-removal mechanisms were observed in the two wetland cells; sulfur dynamics controlled the removal processes in the cell without frequent disturbance but failed to influence metal removal in the cell with frequent disturbance due to the interruption of anoxic layers. The flux ratios that compare labile element concentrations between the water column and the SWI (R-Cu, R-Zn, R-DOC, R-sulfate, and R-chloride) were generated to determine metal diffusive fluxes at the interface. Labile Zn was mostly transported from the water to the SWI during all seasons (R-Zn < 1 for all months except January). Labile Cu moved from the SWI to the water during the warm months (R-Cu < 1), which was explained by the bioturbation-induced transport of organic matter based on the positive correlations between R-Cu and R-DOC. In general, sediment can serve either as a sink or a source depending on the environmental conditions, metal speciation, and presence of living organisms. Here, metal flux at the SWI is a key component in the biogeochemical cycling of a constructed wetland.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Elucidating the phase transformations and grain growth behavior of O3-type sodium-ion layered oxide cathode materials during high temperature synthesis

Understanding the formation mechanism of layered oxide cathodes via solid-state synthesis is imperative to achieving controllability over their materials properties and electrochemical behaviors. In this work, we investigate the phase and microstructure evolution during the synthesis of NaNi 1/3 Fe 1/3 Mn 1/3 O 2 , a model sodium-ion layered oxide cathode, using a combination of imaging, diffraction, and spectroscopic techniques. We unravel the synthetic mechanistic pathways involved in the high-temperature calcination reaction, as well as elaborate the synthesis-microstructure-performance relationship of this material. The formation of the final layered oxide phase involves a gradual transformation through a sodiated oxyhydroxide intermediate. During the reaction, the precursor dehydration reaction dominates at 250–550 °C, while the major sodiation reaction occurs at 550–850 °C. Alongside multiple stages of phase transformations, the final grain structure formation occurs through the continuous growth of the (003) and (104) facets. During the reaction, Mn acts as the charge-compensating element and exhibits depth-dependent characteristics. When the sodiation reaction dominates over dehydration, the reaction intermediates undergo gradual electronic structure changes with increasing temperature, as indicated by the spectral features of TM3d-O2p hybrid states. Calcination duration is also a critical parameter governing the microstructure, surface reactivity, phase fraction distribution and electrochemical performance of the material. The optimal calcination duration was determined to be 18 hours at 850 °C under the conditions evaluated here. Calcination beyond this duration was found to be detrimental to electrochemical performance due to Na and O loss and heterogeneous sodium distribution throughout the particles. Our work sheds light on the complex crystallographic-chemical-microstructural evolution of sodium ion layered oxide cathodes and provides insight into precisely tuning material properties which are intimately linked to battery performances.

25 ENERGY STORAGE↗
NASA NTRS2002

Dynamics of Sheared Granular Materials

This work focuses on the properties of sheared granular materials near the jamming transition. The project currently involves two aspects. The first of these is an experiment that is a prototype for a planned ISS (International Space Station) flight. The second is discrete element simulations (DES) that can give insight into the behavior one might expect in a reduced-g environment. The experimental arrangement consists of an annular channel that contains the granular material. One surface, say the upper surface, rotates so as to shear the material contained in the annulus. The lower surface controls the mean density/mean stress on the sample through an actuator or other control system. A novel feature under development is the ability to 'thermalize' the layer, i.e. create a larger amount of random motion in the material, by using the actuating system to provide vibrations as well control the mean volume of the annulus. The stress states of the system are determined by transducers on the non-rotating wall. These measure both shear and normal components of the stress on different size scales. Here, the idea is to characterize the system as the density varies through values spanning dense almost solid to relatively mobile granular states. This transition regime encompasses the regime usually thought of as the glass transition, and/or the jamming transition. Motivation for this experiment springs from ideas of a granular glass transition, a related jamming transition, and from recent experiments. In particular, we note recent experiments carried out by our group to characterize this type of transition and also to demonstrate/ characterize fluctuations in slowly sheared systems. These experiments give key insights into what one might expect in near-zero g. In particular, they show that the compressibility of granular systems diverges at a transition or critical point. It is this divergence, coupled to gravity, that makes it extremely difficult if not impossible to characterize the transition region in an earth-bound experiment. In the DE modeling, we analyze dynamics of a sheared granular system in Couette geometry in two (2D) and three (3D) space dimensions. Here, the idea is to both better understand what we might encounter in a reduced-g environment, and at a deeper level to deduce the physics of sheared systems in a density regime that has not been addressed by past experiments or simulations. One aspect of the simulations addresses sheared 2D system in zero-g environment. For low volume fractions, the expected dynamics of this type of system is relatively well understood. However, as the volume fraction is increased, the system undergoes a phase transition, as explained above. The DES concentrate on the evolution of the system as the solid volume fraction is slowly increased, and in particular on the behavior of very dense systems. For these configurations, the simulations show that polydispersity of the sheared particles is a crucial factor that determines the system response. Figures 1 and 2 below, that present the total force on each grain, show that even relatively small (10 %) nonuniformity of the size of the grains (expected in typical experiments) may lead to significant modifications of the system properties, such as velocity profiles, temperature, force propagation, and formation shear bands. The simulations are extended in a few other directions, in order to provide additional insight to the experimental system analyzed above. In one direction, both gravity, and driving due to vibrations are included. These simulations allow for predictions on the driving regime that is required in the experiments in order to analyze the jamming transition. Furthermore, direct comparison of experiments and DES will allow for verification of the modeling assumptions. We have also extended our modeling efforts to 3D. The (preliminary) results of these simulations of an annular system in zero-g environment will conclude the presentation.

Kondic, Lou↗
NASA NTRS2003

Numerical Simulations of Non-Newtonian Convection in Ice: Application to Europa

Numerical simulations of solid state convection in Europa's ice shell have so far been limited to consideration of Newtonian flow laws, where the viscosity of ice is strongly dependent upon temperature, predicting that a stagnant lid should form at the top (10-40%) of a convecting ice shell. Such large thicknesses seem to contradict estimates of the effective elastic thickness of Europa s ice shell during its geologically active period. Recent laboratory experiments characterize the rheology of ice as the sum of contributions from several temperature and strain rate-dependent creep mechanisms. We present the results of numerical simulations of convection within Europa s ice shell using the finite-element model Citcom, applying the non-Newtonian rheology of grain boundry sliding. Our calculations suggest a shallower brittle/ductile transition and larger interior convective velocities compared to Newtonian rheology. The flow field is time-dependent, with small, localized upwellings and downwellings at the thermal boundary layers that have minimal topographic expression at the surface.

Barr, A. C.↗
NASA NTRS2010

Comparisons of Supergranule Properties from SDO/HMI with Other Datasets

While supergranules, a component of solar convection, have been well studied through the use of Dopplergrams, other datasets also exhibit these features. Quiet Sun magnetograms show local magnetic field elements distributed around the boundaries of supergranule cells, notably clustering at the common apex points of adjacent cells, while more solid cellular features are seen near active regions. Ca II K images are notable for exhibiting the chromospheric network representing a cellular distribution of local magnetic field lines across the solar disk that coincides with supergranulation boundaries. Measurements at 304 A further above the solar surface also show a similar pattern to the chromospheric network, but the boundaries are more nebulous in nature. While previous observations of these different solar features were obtained with a variety of instruments, SDO provides a single platform, from which the relevant data products at a high cadence and high-definition image quality are delivered. The images may also be cross-referenced due to their coincidental time of observation. We present images of these different solar features from HMI & AIA and use them to make composite images of supergranules at different atmospheric layers in which they manifest. We also compare each data product to equivalent data from previous observations, for example HMI magnetograms with those from MDI.

Pesnell, William Dean↗