Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sharp interface”

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 181 records · Page 10

The Implementation of Framework for Improvement by Vertical Enhancement Into Energy Exascale Earth System Model

Abstract The low cloud bias in global climate models (GCMs) remains an unsolved problem. Coarse vertical resolution in GCMs has been suggested to be a significant cause of low cloud bias because planetary boundary layer parameterizations cannot resolve sharp temperature and moisture gradients often found at the top of subtropical stratocumulus layers. This work aims to ameliorate the low cloud problem by implementing a new computational method, the Framework for Improvement by Vertical Enhancement (FIVE), into the Energy Exascale Earth System Model (E3SM). Three physics schemes representing microphysics, radiation, and turbulence as well as vertical advection are interfaced to vertically enhanced physics (VEP), which allows for these processes to be computed on a higher vertical resolution grid compared to the rest of the E3SM model. We demonstrate the better representation of subtropical boundary layer clouds with FIVE while limiting additional computational cost from the increased number of levels. When the vertical resolution approaches the large eddy simulation‐like vertical resolution in VEP, the climatological low cloud amount shows a significant increase of more than 30% in the southeastern Pacific Ocean. Using FIVE to improve the representation of low‐level clouds does not come with any negative side effects associated with the simulation of mid‐ and high‐level cloud and precipitation, that can occur when running the full model at higher vertical resolution.

54 ENVIRONMENTAL SCIENCES↗

Alcove formation in dissolving cliffs driven by density inversion instability

We demonstrate conditions that give rise to cave-like features commonly found in dissolving cliffsides with a minimal two-phase physical model. Alcoves that are wider at the top and tapered at the bottom, with sharp-edged ceilings and sloping floors, are shown to develop on vertical solid surfaces dissolving in aqueous solvents. As evident from descending plumes, sufficiently large indentations evolve into alcoves as a result of the faster dissolution of the ceiling due to a solutal Rayleigh–Bénard density inversion instability. In contrast, defects of size below the boundary layer thickness set by the critical Rayleigh number smooth out, leading to stable planar interfaces. Furthermore, the ceiling recession rate and the alcove opening area evolution are shown to be given to first-order by the critical Rayleigh number. By tracking passive tracers in the fluid phase, we show that the alcoves are shaped by the detachment of the boundary layer flow and the appearance of a pinned vortex at the leading edge of the indentations. The attached boundary layer past the developing alcove is then found to lead to rounding of the other sides and the gradual sloping of the floor.

58 GEOSCIENCES↗

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↗

Winding θ and destructive interference of instantons

While the θ dependence of field theories is 2π periodic, the ground-state wavefunctions at θ and θ + 2π often belong to different classes of symmetry-protected topological states. When this is the case, a continuous change of the θ parameter can introduce an interface that supports a nontrivial field theory localized on the wall. We consider the 2d CP N-1 sigma model as an example and construct a weak-coupling setup of this interface theory by considering the small S1 compactification with nonzero winding θ parameter and a suitable symmetry-twisted boundary condition. This system has N classical vacua connected by fractional instantons, but the anomaly constraint tells us that the fractional-instanton amplitudes should vanish completely to have N-fold degeneracy at the quantum level. We show how this happens in this purely bosonic system, uncovering that the integration over the zero modes annihilates the fractional instanton amplitudes, in sharp contrast to what happens when the θ angle is constant. Moreover, we provide another explanation of this selection rule by showing that the N perturbative vacua acquire different charges under the global symmetry with the activation of the winding θ angle. We also demonstrate a similar destructive interference between instanton effects in the CP N-1 quantum mechanics with the Berry phase.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Ir(hkl) Surface Electrochemistry in a Nonadsorbing Acidic Medium

The fundamental properties of electrochemical materials depend on the multiple and often complex interactions between electrode surface sites and electrolyte species at the electrochemical interface. Despite Iridium use in electrolyzer systems, much of its surface electrochemistry remains underexplored. This study investigates the surface electrochemistry of Ir(111), Ir(100), and Ir(110) surfaces in acidic media. Using cyclic voltammetry and CO charge displacement experiments, we establish the charge states and adsorbate coverages as a function of the electrode potential, revealing the presence of hydrogen and hydroxyl co-adsorption at low potentials on (111), and almost no coverage of H ad on (110) facet. In situ Shell Isolated Nanoparticle Enhanced Raman Spectroscopy experiments provide direct evidence of the formation of key adsorbate species, such as hydrogen, hydroxyl, and oxygen, but most importantly, their interactions with interfacial water, confirmed by Density Functional Theory calculations. Our findings highlight the role of co-adsorption and interspecies interactions, with microkinetic adsorption voltammetry simulations corroborating the influence of lateral interactions on adsorption dynamics, particularly for Ir(100) where the OHad formation occurs as a sharp adsorption/desorption current. Our results underscores the importance of interfacial water and hydrogen bonding networks in shaping the electrochemical behavior on Ir surfaces, refining our baseline understanding of the Ir surface electrochemistry necessary for the development of advanced Ir-based electrochemical materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Regolith, rock and fluid distributions at the Upper Colorado River Basin via a multicomponent seismic imaging approach

In October 2018, we conducted an exploratory seismic survey of the East River Watershed Function Science Focus Area, near Crested Butte, CO. Our survey goals were to map seismic bedrock, regolith and sediment properties and relate them to geomorphic and hydrologic processes. We acquired ~12 km of seismic data along roads with a 72-channel, 1.25m spaced streamer, an average shot spacing of 2.5m and an average acquisition rate of 3 km/day. We also obtained ~3km of planted geophone data with geophone and shot spacing between 2.5 and 5 meters and an acquisition rate of 0.5 km/day. Vp results show a sharp transition between regolith and bedrock along hillslopes suggesting a similarly abrupt change in porosity. Vp anisotropy in the Mancos Shale is observable in the lower reaches of the watershed, but is absent in regions near crystalline intrusions where contact metamorphism has likely modified pore structure within the shale. Alluvium thickness varies along the river axis from 0 to > 20 meters reflecting variations in depositional processes as function of stream gradient. We show that when strong discontinuities at the regolith base are present, inhomogeneous waves that travel along subsurface interfaces are generated and can be measured on vertical component geophones. We present a summary manuscript in the report appendix.

58 GEOSCIENCES↗

M-Edge Spectroscopy of Transition Metals: Principles, Advances, and Applications

M-edge X-ray absorption spectroscopy (XAS), which probes 3p→3d transitions in first-row transition metals, provides detailed insights into oxidation states, spin-states, and local electronic structure with high element and orbital specificity. Operating in the extreme ultraviolet (XUV) region, this technique provides sharp multiplet-resolved features with high sensitivity to ligand field and covalency effects. Compared to K- and L-edge XAS, M-edge spectra exhibit significantly narrower full widths at half maximum (typically 0.3–0.5 eV versus >1 eV at the L-edge and >1.5–2 eV at the K-edge), owing to longer 3p core-hole lifetimes. M-edge measurements are also more surface-sensitive due to the lower photon energy range, making them particularly well-suited for probing thin films, interfaces, and surface-bound species. The advent of tabletop high-harmonic generation (HHG) sources has enabled femtosecond time-resolved M-edge measurements, allowing direct observation of ultrafast photoinduced processes such as charge transfer and spin crossover dynamics. This review presents an overview of the fundamental principles, experimental advances, and current theoretical approaches for interpreting M-edge spectra. We further discuss a range of applications in catalysis, materials science, and coordination chemistry, highlighting the technique’s growing impact and potential for future studies.

M-edge XANES↗

Moiréless correlations in ABCA graphene

Atomically thin van der Waals materials stacked with an interlayer twist have proven to be an excellent platform toward achieving gate-tunable correlated phenomena linked to the formation of flat electronic bands. In this work we demonstrate the formation of emergent correlated phases in multilayer rhombohedral graphene––a simple material that also exhibits a flat electronic band edge but without the need of having a moiré superlattice induced by twisted van der Waals layers. We show that two layers of bilayer graphene that are twisted by an arbitrary tiny angle host large (micrometer-scale) regions of uniform rhombohedral four-layer (ABCA) graphene that can be independently studied. Scanning tunneling spectroscopy reveals that ABCA graphene hosts an unprecedentedly sharp van Hove singularity of 3–5-meV half-width. We demonstrate that when this van Hove singularity straddles the Fermi level, a correlated many-body gap emerges with peak-to-peak value of 9.5 meV at charge neutrality. Mean-field theoretical calculations for model with short-ranged interactions indicate that two primary candidates for the appearance of this broken symmetry state are a charge-transfer excitonic insulator and a ferrimagnet. Finally, we show that ABCA graphene hosts surface topological helical edge states at natural interfaces with ABAB graphene which can be turned on and off with gate voltage, implying that small-angle twisted double-bilayer graphene is an ideal programmable topological quantum material.

42 ENGINEERING↗

Off-state magnetoresistance in long-channel germanium Schottky-barrier MOSFETs

An increasing magnetic field perpendicular to an undoped semiconductor surface at low temperature is known to strengthen the binding of localized electrons to stationary ions, as the wavefunction's tails evolve from exponential to Gaussian. It is also known that application of a high bias voltage to a depleted semiconductor can liberate bound charge and induce a large drop in electrical resistance. We connect these established results to experimental electrical transport measurements on off-state germanium Schottky-barrier metal–oxide–semiconductor field-effect transistor (MOSFETs) with an aluminum oxide insulating dielectric and platinum germanide contacts. We make measurements at the three distinct orientations of the magnetic field with respect to the substrate and the current. At 6 K, we observe sharp attenuation of current by more than 2 orders of magnitude, within 60 mT, at a crossover magnetic field perpendicular to the substrate. A 1 T magnetic field attenuates the current by more than 4 orders of magnitude. The strength of the attenuation and the value of the crossover field are controlled by both the gate–source and drain–source voltages. The attenuation is much weaker when the magnetic field is parallel to the current. Finally, we orient the magnetic field parallel to the substrate, but perpendicular to the current, allowing us to distinguish charge hopping at the oxide interface from charge hopping in the bulk. In conclusion, this large off-state magnetoresistance can be exploited for cryogenic magnetic- and photo-detection, and for high-bias, low-leakage MOSFETs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Microscale mechanical modeling of deformable geomaterials with dynamic contacts based on the numerical manifold method

Abstract Micromechanical modeling of geomaterials is challenging because of the complex geometry of discontinuities and potentially large number of deformable material bodies that contact each other dynamically. In this study, we have developed a numerical approach for micromechanical analysis of deformable geomaterials with dynamic contacts. In our approach, we detect contacts among multiple blocks with arbitrary shapes, enforce different contact constraints for three different contact states of separated, bonded, and sliding, and iterate within each time step to ensure convergence of contact states. With these features, we are able to simulate the dynamic contact evolution at the microscale for realistic geomaterials having arbitrary shapes of grains and interfaces. We demonstrate the capability with several examples, including a rough fracture with different geometric surface asperity characteristics, settling of clay aggregates, compaction of a loosely packed sand, and failure of an intact marble sample. With our model, we are able to accurately analyze (1) large displacements and/or deformation, (2) the process of high stress accumulated at contact areas, (3) the failure of a mineral cemented rock samples under high stress, and (4) post-failure fragmentation. The analysis highlights the importance of accurately capturing (1) the sequential evolution of geomaterials responding to stress as motion, deformation, and high stress; (2) large geometric features outside the norms (such as large asperities and sharp corners) as such features can dominate the micromechanical behavior; and (3) different mechanical behavior between loosely packed and tightly packed granular systems.

58 GEOSCIENCES↗

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗

SNL-NJOY-2016

This is a wrapper around the LANL NJOY-2016 code that interfaces with extended code capabilities and minor modification of the LANL NJOY-2016 code that supports modeling radiation damage to materials. The NJOY Nuclear Data Processing System is a modular computer code designed to read evaluated data in ENDF format, transform the data in various ways, and output the results as libraries designed to be used in various applications. The wrapper provided here permits Sandia-specific control parameters to be used in the input data file. The modifications incorporated here enhance the ability of NJOY to address material damage response functions for many materials, e.g. to include the NRT and arc-dpa forms of the damage energy in addition to the default NJOY-2016 implementation of the sharp-threshold Kinchin-Pease threshold energy treatment. All of the modifications provided here are being made available to the GitHub-based NJOY-2016 code. As useful enhancements found here are incorporated into the baseline NJOY-2016 code, they will be eliminated from this version so as to maintain our compatibility with the baseline NJOY-2016 code. SAND2020-13060 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA00035

Griffin, Patrick↗

Neural Network-Enhanced Reproducing Kernel Particle Method for Image-Based Multiphysics Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant stresses during charge/discharge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), image-based modeling techniques are used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electrochemical-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in capturing crack propagation due to mesh dependency. Additionally, commonly used damage models, such as the continuous damage model and the cohesive zone model, often have steep tradeoffs between discontinuous field accuracy and computational expense. In this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure. Reference: [1] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, pp 4422-4454, https://doi.org/10.1002/nme.7040, 2022.

damage modeling↗

Selectivity in Yttrium Manganese Oxide Synthesis via Local Chemical Potentials in Hyperdimensional Phase Space

In sharp contrast to molecular synthesis, materials synthesis is generally presumed to lack selectivity. The few known methods of designing selectivity in solid-state reactions have limited scope, such as topotactic reactions or strain stabilization. This contribution describes a general approach for searching large chemical spaces to identify selective reactions. This novel approach explains the ability of a nominally “innocent” Na 2 CO 3 precursor to enable the metathesis synthesis of single-phase Y 2 Mn 2 O 7 : an outcome that was previously only accomplished at extreme pressures and which cannot be achieved with closely related precursors of Li 2 CO 3 and K 2 CO 3 under identical conditions. By calculating the required change in chemical potential across all possible reactant-product interfaces in an expanded chemical space including Y, Mn, O, alkali metals, and halogens, using thermodynamic parameters obtained from density functional theory calculations, we identify reactions that minimize the thermodynamic competition from intermediates. In this manner, only the Na-based intermediates minimize the distance in the hyperdimensional chemical potential space to Y 2 Mn 2 O 7 , thus providing selective access to a phase which was previously thought to be metastable. Experimental evidence validating this mechanism for pathway-dependent selectivity is provided by intermediates identified from in situ synchrotron-based crystallographic analysis. Finally, this approach of calculating chemical potential distances in hyperdimensional compositional spaces provides a general method for designing selective solid-state syntheses that will be useful for gaining access to metastable phases and for identifying reaction pathways that can reduce the synthesis temperature, and cost, of technological materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Neural Network Enhanced RKPM for Electrochemical-Mechanical Coupled Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant charge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. For Li-ion batteries, for example, lithium moving between electrodes during charging and discharging process causes expansion and contraction of grains, and the strongly anisotropic and nonlinearly [Li]-dependent grain material properties can cause grains to expand into and contract away from each other, leading to chemo-mechanical cracking. In the first part of this work, a RKPM based computational framework for solving the coupled solid-phase lithium conservation with Fickian diffusion and the lithium concentration dependent anisotropic mechanical problem subjected to a highly nonlinear Butler-Volmer boundary condition is introduced. The choice of RKPM completeness conditions for lithium concentration and mechanical deformation fields, and the variational consistency condition for the domain integration of the coupled problem is first determined. In the second part of this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material, by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that is necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure.

damage modeling↗

Magnetic excitations in strained infinite-layer nickelate PrNiO 2 films

Strongly correlated materials respond sensitively to external perturbations such as strain, pressure, and doping. In the recently discovered superconducting infinite-layer nickelates, the superconducting transition temperature can be enhanced via only ~ 1% compressive strain-tuning with the root of such enhancement still being elusive. Using resonant inelastic x-ray scattering (RIXS), we investigate the magnetic excitations in infinite-layer PrNiO 2 thin films grown on two different substrates, namely SrTiO 3 (STO) and (LaAlO 3 ) 0.3 (Sr 2 TaAlO 6 ) 0.7 (LSAT) enforcing different strain on the nickelates films. The magnon bandwidth of PrNiO 2 shows only marginal response to strain-tuning, in sharp contrast to the enhancement of the superconducting transition temperature T c in the doped superconducting samples. These results suggest the bandwidth of spin excitations of the parent compounds is similar under strain while T c in the doped ones is not, and thus provide important empirics for the understanding of superconductivity in infinite-layer nickelates.

36 MATERIALS SCIENCE↗

Engineering Edge States of Graphene Nanoribbons for Narrow-Band Photoluminescence

Solid-state narrow-band light emitters are on-demand for quantum optoelectronics. Current approaches based on defect engineering in low-dimensional materials usually introduce a broad range of emission centers. Here, we report narrow-band light emission from covalent heterostructures fused to the edges of graphene nanoribbons (GNRs) by controllable on-surface reactions from molecular precursors. Two types of heterojunction (HJ) states are realized by sequentially synthesizing GNRs and graphene nanodots (GNDs) and then coupling them together. HJs between armchair GNDs and armchair edges of the GNR are coherent and give rise to narrow-band photoluminescence. In contrast, HJs between the armchair GNDs and the zigzag ends of GNRs are defective and give rise to nonradiative states near the Fermi level. At low temperatures, sharp photoluminescence emissions with peak energy range from 2.03 to 2.08 eV and line widths of 2–5 meV are observed. Furthermore, the radiative HJ states are uniform, and the optical transition energy is controlled by the band gaps of GNRs and GNDs. As these HJs can be synthesized in a large quantity with atomic precision, this finding highlights a route to programmable and deterministic creation of quantum light emitters.

36 MATERIALS SCIENCE↗