Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface coupling”

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 631 records · Page 35

Tuning the interfacial spin-orbit coupling with ferroelectricity

Detection and manipulation of spin current lie in the core of spintronics. Here we report an active control of a net spin Hall angle, θ SHE (net), in Pt at an interface with a ferroelectric material PZT (PbZr 0.2 Ti 0.8 O 3 ), using its ferroelectric polarization. The spin Hall angle in the ultra-thin Pt layer is measured using the inverse spin Hall effect with a pulsed tunneling current from a ferromagnetic La 0.67 Sr 0.33 MnO 3 electrode. The effect of the ferroelectric polarization on θ SHE (net) is enhanced when the thickness of the Pt layer is reduced. When the Pt layer is thinner than 6 nm, switching the ferroelectric polarization even changes the sign of θ SHE (net). This is attributed to the reversed polarity of the spin Hall angle in the 1 st -layer Pt at the PZT/Pt interface when the ferroelectric polarization is inverted, as supported by the first-principles calculations. These findings suggest a route for designing future energy efficient spin-orbitronic devices using ferroelectric control.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Adsorption of Non-Ionic Surfactant and Monoclonal Antibody on Siliconized Surface Studied by Neutron Reflectometry

The adsorption of monoclonal antibodies (mAbs) on hydrophobic surfaces is known to cause protein aggregation and degradation. Therefore, surfactants, such as Poloxamer 188, are widely used in therapeutic formulations to stabilize mAbs and protect mAbs from interacting with liquid-solid interfaces. Here, the adsorption of Poloxamer 188, one mAb and their competitive adsorption on a model hydrophobic siliconized surface is investigated with neutron scattering coupled with contrast variation to determine the molecular structure of adsorbed layers for each case. Small angle neutron scattering measurements of the affinity of Poloxamer 188 to this mAb indicate that there is negligible binding at these solution conditions. Neutron reflectometry measurements of the mAb show irreversible adsorption on the siliconized surface, which cannot be washed off with neat buffer. Poloxamer 188 can be adsorbed on the surface already occupied by mAb, which enables partial removal of some adsorbed mAb by washing with buffer. The adsorption of the surfactant introduces significant conformational changes for mAb molecules that remain on the surface. In contrast, if the siliconized surface is first saturated with the surfactant, no adsorption of mAb is observed. Competitive adsorption of mAb and Poloxamer 188 from solution leads to a surface dominantly occupied with surfactant molecules, whereas only a minor amount of mAb absorbs. These findings clearly indicate that Poloxamer 188 can protect against mAb adsorption as well as modify the adsorbed conformation of previously adsorbed mAb.

36 MATERIALS SCIENCE↗

Snowflake relocated Cu 2 O electrocatalyst on an Ag backbone template for the production of liquid C 2+ chemicals from CO 2

In this study, we performed the CO 2 reduction reaction (CO 2 RR) using a structural composite catalyst of cuprous oxide (Cu 2 O) and silver (Ag) that was simultaneously electrodeposited. While the underneath Ag electrodeposits maintained their spiky backbone structures even after the CO 2 RR, the Cu 2 O deposits were reduced to Cu(111) and relocated on the backbone template. Here, the structural changes in Cu 2 O to Cu increase the active area of the Cu–Ag interface, resulting in a remarkable production rate of 125.01 μmol h -1 of liquid C 2+ chemicals via the stabilization of the C–C coupling of the key intermediate species of acetaldehyde. This study provides new insights into designing a bimetallic catalyst for producing sustainable C 2+ products from CO 2 without any selectivity towards the production of methane.

36 MATERIALS SCIENCE↗

Mapping causal patterns in crystalline solids

The evolution of the atomic structures of the combinatorial library of Sm-substituted thin film BiFeO 3 along the phase transition boundary from the ferroelectric rhombohedral phase to the non-ferroelectric orthorhombic phase is explored using scanning transmission electron microscopy. Localized properties, including polarization, lattice parameter, and chemical composition, are parameterized from atomic-scale imaging, and their causal relationships are reconstructed using a linear non-Gaussian acyclic model. This approach is further extended to explore the spatial variability of the causal coupling using the sliding window transform method, which revealed that new causal relationships emerged at both the expected locations, such as domain walls and interfaces, and at additional regions forming clusters in the vicinity of the walls or spatially distributed features. While the exact physical origins of these relationships are unclear, they likely represent nanophase-separated regions in the morphotropic phase boundaries. Overall, we posit that an in-depth understanding of complex disordered materials away from thermodynamic equilibrium necessitates understanding not only the generative processes that can lead to observed microscopic states but also the causal links between multiple interacting subsystems.

Causal inference↗

Observations and seismoacoustic simulations of earthquake-generated infrasound waves in non-epicentral regions

We analysed infrasound waves associated with the Gyeongju earthquake (ML 5.8) that occurred on 2016 September 12, in the southeastern Korean Peninsula. For infrasound wave detection, the Progressive Multichannel Correlation method was applied to the infrasound data set recorded at seven arrays operating in South Korea at epicentral distances ranging from 178 to 472 km. Based on the back-projection method constrained by array-dependent celerity and azimuth deviation models, the source regions were identified in both the epicentral and non-epicentral regions. Remarkably, the non-epicentral secondary sources of this earthquake were located in regions with shallow water depths: (i) the western coastal area in the Yellow Sea and (ii) the shallow ocean basin and bank in the East Sea. The location results obtained from the earthquake could be corroborated through its foreshock (ML 5.1), yielding location results consistent with those of the main shock. The generation of infrasound waves over shallow water depths was fortuitously validated by direct recordings of dominant single-frequency (~0.3 Hz) infrasound waves at close range via temporary sensors near the ocean basin and bank. We interpreted that low-frequency infrasound signals could be generated from interactions among the ocean floor, shallow seawater and atmosphere. We performed numerical simulations of seismoacoustic fields to predict ground motions on the seafloor and acoustic transmission efficiency between the water and air interface. The simulations quantified the energy transfer through different media and clarified our observational results. We found that because this solid Earth‒water‒atmosphere coupled air wave has a relatively low frequency (~0.3 Hz), it can survive propagation over long distances compared with high-frequency infrasound waves generated in inland and mountain regions. In this study, we extend our understanding of water‒atmosphere coupling and the monitoring framework for earthquake-associated non-epicentral infrasound waves, encompassing not only inland ground shaking but also shallow sea regions located far from the epicentre.

Earthquake monitoring and test-ban treaty verifica↗

Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes

Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multimessenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of fluid angular momentum, as well as requiring fewer numbers of computational cells. However, the use of spherical coordinates to numerically solve hyperbolic partial differential equations can result in severe Courant-Friedrichs-Lewy (CFL) stability condition time step limitations, which can make simulations prohibitively expensive. This paper addresses this issue for the numerical solution of coupled spacetime and general relativistic magnetohydrodynamics evolutions by introducing a double fast Fourier transform (FFT) filter and implementing it within the fully message passing interface (mpi)-parallelized sphericalnr framework in the einstein toolkit. In conclusion, we demonstrate the effectiveness and robustness of the filtering algorithm by applying it to a number of challenging code tests, and show that it passes these tests effectively, demonstrating convergence while also increasing the time step significantly compared to unfiltered simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluating Climate Models with the CLIVAR 2020 ENSO Metrics Package

El Niño–Southern Oscillation (ENSO) is the dominant mode of interannual climate variability on the planet, with far-reaching global impacts. It is therefore key to evaluate ENSO simulations in state-of-the-art numerical models used to study past, present, and future climate. Recently, the Pacific Region Panel of the International Climate and Ocean: Variability, Predictability and Change (CLIVAR) Project, as a part of the World Climate Research Programme (WCRP), led a community-wide effort to evaluate the simulation of ENSO variability, teleconnections, and processes in climate models. The new CLIVAR 2020 ENSO metrics package enables model diagnosis, comparison, and evaluation to 1) highlight aspects that need improvement; 2) monitor progress across model generations; 3) help in selecting models that are well suited for particular analyses; 4) reveal links between various model biases, illuminating the impacts of those biases on ENSO and its sensitivity to climate change; and to 5) advance ENSO literacy. By interfacing with existing model evaluation tools, the ENSO metrics package enables rapid analysis of multipetabyte databases of simulations, such as those generated by the Coupled Model Intercomparison Project phases 5 (CMIP5) and 6 (CMIP6). The CMIP6 models are found to significantly outperform those from CMIP5 for 8 out of 24 ENSO-relevant metrics, with most CMIP6 models showing improved tropical Pacific seasonality and ENSO teleconnections. Only one ENSO metric is significantly degraded in CMIP6, namely, the coupling between the ocean surface and subsurface temperature anomalies, while the majority of metrics remain unchanged.

54 ENVIRONMENTAL SCIENCES↗

Real-Time Drilling Optimization System for Improved Overall Rate of Penetration and Reduced Cost Per Foot in Geothermal Drilling

The key to success in geothermal drilling is economic feasibility, and a major cost in the development of geothermal resources is the actual drilling of the wells. In this project, a real-time drilling optimization system for geothermal drilling was developed. The system couples three individual components while drilling. The first component is a drill stem vibration analysis model, the second is Mechanical Specific Energy (MSE) analyses, and the third is a detailed PDC Rate of Penetration (ROP) drill bit model for optimum RPM and WOB combinations. The benefit of the coupled system is that the range of WOB and RPM could be selected to avoid drill stem vibrations. Secondly, MSE is used as an efficiency measure and the detailed PDC drill bit model ensures the drill bit does not endure temperatures that exceed the temperature at which the PDC cutters experience accelerated wear. The new detailed PDC bit model is based on rock/bit interaction that physically tracks the PDC cutter wear flats as the bit drills ahead giving the capability to calculate the temperature being generated underneath the worn cutters to better advise on operational parameters to avoid accelerated cutter wear and failure and to ensure that operational parameters are applied so that overall ROP is maximized. By combining the drill stem vibrations and the detailed PDC bit cutter wear and “safe” non-accelerated cutter wear temperature and optimum ranges of operating parameters, it results in higher ROP and lower cost drilling. Single cutter PDC testing performed in different lithologies at Sandia was utilized to verify the PDC cutter forces and depth of cut for new and worn cutters. Based on single cutter PDC temperature modeling, verification using single cutter data from the testing done by National Oilwell Varco (NOV) was performed. Sandia’s Hard-Rock Drilling Facility (HRDF) was utilized to test different drill bit configurations with different cutter designs and wear status with different induced modes of vibration to obtain the critical bit RPM/WOB ranges resulting in ineffective drilling and low ROP. The collected test data were further used to verify and calibrate the full hole PDC ROP model that was developed based on single cutter interaction data. A full coupled drill stem vibration model was formulated and verified with geothermal field data from the Chocolate Mountain Aerial Gunnery Range (CMAGR). A graphical user interface (GUI) was developed using Tkinter library in the computer programming language Python, which integrates all the developed models in one system. The developed system consists mainly of the PDC ROP model, PDC bit wear model, PDC cutter temperature model, Mechanical Specific Energy (MSE) model, and drillstring vibration model integrated into one system. The developed system can be used for both, post well analysis and real-time optimization using different criteria such as ROP maximization or MSE minimization. The software uses Differential Evolution Algorithm (DEA) to find optimum values for operational parameters based on last foot drilled while avoiding the drillstring vibration and cutter temperature critical operating parameters.

15 GEOTHERMAL ENERGY↗

Improving the Capabilities and Computational Efficiency of the RTE+RRTMGP Radiation Code (Final Report)

This report details progress on the RTE+RRTMGP radiation codes made during the period of performance. RTE+RRTMGP is a set of codes for computing radiative fluxes in planetary atmospheres. RRTMGP uses a k-distribution to provide an optical description (absorption and possibly Rayleigh optical depth) of the gaseous atmosphere, along with the relevant source functions, on a pre-determined spectral grid given temperatures, pressures, and gas concentration. RTE computes fluxes given spectrally-resolved optical descriptions and source functions. Spectrally-resolved fluxes are summarized (“reduced”) via a user extensible class. The initial release of the code and the design choices are described in Pincus et al. 2019; the codes are available on Github. Although RRTMGP was based on current (at the time) empirical spectroscopic data, RTE and RRTMGP were developed in large part to modernize software practices. The design focused on flexibility broadly interpreted: by separating code from data and allowing data to drive computation; in coupling to the host model (e.g. the coupling of clouds to radiative fluxes is user-controlled); with respect to programming languages (computational tasks are accessed via widely-compatible C interfaces); and with respect to hardware (the codes run on a range of CPU and GPU architectures). The code also puts an emphasis on modularity and clarity. RTE+RRTMGP v1.0 was released in September 20219. This award supported the evolution of the RTE+RRTMGP code base to support greater flexibility, accuracy, and efficiency.

54 ENVIRONMENTAL SCIENCES↗

FY25 Cellular Silicone Foam Friction Testing Summary Report

This effort was a collaboration between LANL groups MPA-CINT, E-13, and V-14 to characterize various properties of foam materials for programmatic work. Specifically, MPA-CINT staff members’ contribution to this effort was to perform friction testing of cellular silicone foams of two different porosities with respect to two engineering materials with predefined surface roughnesses. The selected foams (henceforth serving as the sled material interface exposed to motion) had porosities of 67% and 70% with respect to fabricated values. The other materials used in the friction testing measurement couple (static track material) were 304 Stainless Steel (304 SS) and Very High-Temperature Glass-Mica Ceramic, known as Macor®.

36 MATERIALS SCIENCE↗

Compact, high-efficiency accelerators driven by low-voltage solid-state amplifiers

A compact particle accelerator can include two or more cavities disposed along an axis of the particle accelerator, each of which is coupled to two or more drivers. The accelerator can also include a power supply coupled to the two or more drivers such that a particle beam traveling along the axis is accelerated. The power supply can be an interface with a commercial power outlet, battery power, or a combination thereof depending upon the use case. Example configurations of the accelerator include hand held or mobile devices that are capable of delivering up to and greater than a 1 MeV electron beam.

42 ENGINEERING↗

Compact, high-efficiency accelerators driven by low-voltage solid-state amplifiers

A compact particle accelerator can include two or more cavities disposed along an axis of the particle accelerator, each of which is coupled to two or more drivers. The accelerator can also include a power supply coupled to the two or more drivers such that a particle beam traveling along the axis is accelerated. The power supply can be an interface with a commercial power outlet, battery power, or a combination thereof depending upon the use case. Example configurations of the accelerator include hand held or mobile devices that are capable of delivering up to and greater than a 1 MeV electron beam.

Nguyen, Dinh Cong↗

Structural metamaterials comprising interpenetrating lattices

Multi-body interpenetrating lattices comprise two or more lattices that interlace or interpenetrate through the same volume without any direct physical connection to each other, wherein energy transfer is controlled by surface interactions. As a result, multifunctional or composite-like responses can be achieved by additive manufacturing of the interpenetrating lattices, even with only a single print material, with programmable interface-dominated properties. As a result, the interpenetrating lattices can have unique mechanical properties, including improved toughness, multi-stable/negative stiffness, and electromechanical coupling.

Boyce, Brad↗

Characterization of self-magnetic pinch radiographic diode performance on RITS-6 at Sandia National Laboratories. II. Coupling between the inductive voltage adder and the SMP load

The self-magnetic pinch (SMP) diode is a type of radiographic diode used to generate an intense electron beam for radiographic applications. At Sandia National Laboratories, SMP was the diode load for the six-cavity radiographic integrated test stand inductive voltage adder (IVA) driver operated in a magnetically insulated transmission line (MITL). The MITL contributes a flow current in addition to the current generated within the diode itself. Extensive experiments with a MITL of 40Ω load impedance [T. J. Renk et al., Phys. Plasmas 29, 023105 (2022)] indicate that the additional flow current leads to results similar to what might be expected from a conventional high-voltage interface driver, where flow current is not present. However, when the MITL flow impedance was increased to 80 Ω, qualitatively different diode behavior was observed. This includes large retrapping waves suggestive of an initial coupling to low impedance as well as diode current decreasing with time even as the total current does not. A key observation is that the driver generates total current (flow + diode) consistent with the flow impedance of the MITL used. The case is made in this paper that the 80Ω MITL experiments detailed here can only be understood when the IVA-MITL-SMP diode is considered as a total system. Additionally, the constraint of fixed total current plus the relatively high flow impedance limits the ability of the diode (whether SMP or other type) to act as an independent load. An unexpected new result is that in tracking the behavior of the electron strike angle on the converter as a function of time, we observed that the conventional cIV x “Radiographic” radiation scaling (where x~2.2) begins to break down for voltages above 8 MV, and cubic scaling is required to recover accurate angle tracking.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A multiphysics coupling framework for exascale simulation of fracture evolution in subsurface energy applications

Predicting the evolution of fractured media is challenging due to coupled thermal, hydrological, chemical and mechanical processes that occur over a broad range of spatial scales, from the microscopic pore scale to field scale. We present a software framework and scientific workflow that couples the pore scale flow and reactive transport simulator Chombo-Crunch with the field scale geomechanics solver in GEOS to simulate fracture evolution in subsurface fluid-rock systems. This new multiphysics coupling capability comprises several novel features. An HDF5 data schema for coupling fracture positions between the two codes is employed and leverages the coarse resolution of the GEOS mechanics solver which limits the size of data coupled, and is, thus, not taxed by data resulting from the high resolution pore scale Chombo-Crunch solver. The coupling framework requires tracking of both before and after coarse nodal positions in GEOS as well as the resolved embedded boundary in Chombo-Crunch. We accomplished this by developing an approach to geometry generation that tracks the fracture interface between the two different methodologies. The GEOS quadrilateral mesh is converted to triangles which are organized into bins and an accessible tree structure; the nodes are then mapped to the Chombo representation using a continuous signed distance function that determines locations inside, on and outside of the fracture boundary. The GEOS positions are retained in memory on the Chombo-Crunch side of the coupling. The time stepping cadence for coupled multiphysics processes of flow, transport, reactions and mechanics is stable and demonstrates temporal reach to experimental time scales. The approach is validated by demonstration of 9 days of simulated time of a core flood experiment with fracture aperture evolution due to invasion of carbonated brine in wellbore-cement and sandstone. We also demonstrate usage of exascale computing resources by simulating a high resolution version of the validation problem on OLCF Frontier.

97 MATHEMATICS AND COMPUTING↗

Tunable n- and p-Type Behavior of Ruddlesden–Popper and Dion–Jacobson 2D Perovskites

Interfacial band alignment in two-dimensional (2D)/three-dimensional (3D) perovskite heterostructures is important for maximizing separation, extraction, and collection of charge carriers, which, in turn, stabilize the perovskite layer during solar cell operation. Despite the wide use of 2D spacer cations for stabilizing perovskite solar cells, spacer cation exchange across the 2D/3D interface induces structural transformation and degradation in performance and stability. We have now examined the electrochemical and photoelectrochemical behavior of 2D Ruddlesden–Popper (RP) and Dion–Jacobson (DJ) perovskites to assess the impact of binding configuration on the interfacial charge separation and their operational stability. In-situ open circuit potential (OCP) responses of 2D perovskites under light illumination, paired with redox couples (ferrocene/ferrocenium), reveal different semiconducting polarizabilities of BA-RP (n-type) and BDA-DJ (p-type). The tunability of the Fermi energy level in 2D perovskites discussed in this study offers insight into the design of 2D/3D interfaces for optimizing charge transfer and enhancing charge neutrality and stability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Redox Couples Control Band Bending, Photovoltage, and Quasi-Fermi Levels in Tungsten Oxide (WO 3 ) Photoanodes

Tungsten oxide (WO 3 ) is a well-known photoanode and photocatalyst for photoelectrochemical (PEC) water oxidation. Because the compound has a deep valence band, it can facilitate the oxygen evolution reaction without added cocatalysts, and it can drive the oxidation of species with much higher electrochemical potentials, including the conversion of water to hydrogen peroxide, sulfate to persulfate, and iodate to meta-periodate. Here, we use the liquid vibrating Kelvin probe surface photovoltage (liquid VK-SPV) technique in combination with open circuit potential (OCP) and photoelectrochemical (PEC) scans to assess the possibility of reaching such oxidizing potentials in aqueous electrolytes and at open circuit. Here, this is done by mapping the quasi-Fermi levels of electrons and holes at the interfaces as a function of the light intensity. Nanostructured WO 3 photoelectrodes for this purpose were fabricated by thermal annealing of a tungstic acid solution on fluorine-doped tin oxide. Electrochemical measurements are conducted at open circuit and 400 nm LED light illumination in electrolytes containing fast (O 2 /H 2 O 2 ), slow (O 2 /H 2 O), and very oxidizing (NaIO 4 /NaIO 3 ) redox couples. Photovoltage values scale with the light intensity and with the built-in potential for each redox couple and reach values up to 0.61 V under 20 mW cm –2 illumination for the NaIO 4 electrolyte. This shows that the photoelectrodes behave like Schottky-type diodes whose maximum possible energy output is determined mainly by the built-in voltage of each junction. For slow redox couples, the quasi-Fermi level of the holes increases with light intensity due to hole accumulation at the WO 3 –liquid interface. For example, for the O 2 /H 2 O electrolyte, interfacial hole accumulation and removal occur on the 90–300 s time scale. For the fast hole acceptor H 2 O 2 , on the other hand, the quasi-Fermi level of the photoholes is pinned to the electrochemical potential of the O 2 /H 2 O 2 couple. This limits the energy conversion efficiency of the electrode. Overall, these results reveal the influence of charge transfer thermodynamics and kinetics on the photovoltage of WO 3 . Furthermore, the work further establishes VK-SPV as a contactless method to observe the photovoltage, carrier dynamics, and quasi-Fermi levels of semiconductor-liquid junctions.

Electrodes↗

Macroscale and Nanoscale Photoelectrochemical Behavior of p-Type Si(111) Covered by a Single Layer of Graphene or Hexagonal Boron Nitride

Two-dimensional (2D) materials may enable a general approach to the introduction of a dipole at a semiconductor surface as well as control over other properties of the double layer at a semiconductor/liquid interface. Vastly different properties can be found in the 2D materials currently studied due in part to the range of the distribution of density-of-states. In this work, the open-circuit voltage (V oc ) of p-Si–H, p-Si/Gr (graphene), and p-Si/h-BN (hexagonal boron nitride) in contact with a series of one-electron outer-sphere redox couples was investigated by macroscale measurements as well as by scanning electrochemical cell microscopy (SECCM). The band gaps of Gr and h-BN (0–5.97 eV) encompass the wide range of band gaps for 2D materials, so these interfaces (p-Si/Gr and p-Si/h-BN) serve as useful references to understand the behavior of 2D materials more generally. The value of Voc shifted with respect to the effective potential of the contacting solution, with slopes (ΔV oc /ΔE Eff ) of -0.27 and -0.38 for p-Si/Gr and p-Si/h-BN, respectively, indicating that band bending at the p-Si/h-BN and p-Si/Gr interfaces responds at least partially to changes in the electrochemical potential of the contacting liquid electrolyte. Additionally, SECCM is shown to be an effective method to interrogate the nanoscale photoelectrochemical behavior of an interface, showing little spatial variance over scales exceeding the grain size of the CVD-grown 2D materials in this work. The measurements demonstrated that the polycrystalline nature of the 2D materials had little effect on the results and confirmed that the macroscale measurements reflected the junction behavior at the nanoscale.

2D material↗