Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffuse interface methods”

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

Creating Functional Oxynitride–Silicon Interfaces and SrNbO 2 N Thin Films for Photoelectrochemical Applications

Photoelectrochemical performance dependence upon absorption length, carrier diffusion length, and surface area of an oxynitride photoabsorber is investigated. How best to fabricate optical-quality thin films of bandgap-tunable oxynitrides is also discussed. We targeted the stoichiometric compound SrNbO 2 N as an optimal wide-bandgap photoabsorber (1.9 eV) for use with silicon (1.1 eV) in a tandem structure photoelectrochemical cell. Preparation of perovskite oxynitrides at high-temperature as isolated powders is often straightforward, but it is difficult to integrate them as thin films in tandem junction devices with low-temperature materials. Here we develop the first method to prepare optical-quality SrNbO 2 N thin-films of tunable thickness and roughness on single-crystal silicon substrate. This achievement required an interfacial layer of ultra-thin TaN to be used as a barrier to reduce the inter-diffusion of silicon and oxygen during oxynitride syn-thesis. We produced a variety of SrNbO 2 N film thicknesses (20-440 nm) on n + -Si(100) surfaces. Roughness factor (0.14-21) scaled with thickness. The intrinsic photoelectrochemical activity of these devices was evaluated using a low-barrier sacrificial electron donor. Photocurrent density and photovoltage revealed a significant (and non-linear) dependence on film thickness and roughness. Furthermore, absorption length, carrier diffusion length, and surface area were each found to play key roles. Balancing these is required for optimally performing devices.

08 HYDROGEN↗
DOE OSTI2021

Influence of Ink Formulation and Drying Conditions on Ionomer Distribution in High-Performance Roll-to-Roll-Coated Gas-Diffusion Electrodes

To enable mass production of fuel cell membrane electrode assemblies (MEAs) catalyst layers production will require continuous roll-to-roll (R2R) coating processes. Gas diffusion electrodes (GDEs) are advantageous for mass production because the catalyst layer can be directely coated on the microporous layer of the gas diffusion media without the need for a decal-transfer process. It is known that the water-to-alcohol ratio in the catalyst ink influences the interactions of the ionomer with the catalyst leading to different distributions of ionomer in spray-coated catalyst layers. It is also known that during drying of colloidal mixtures, like fuel cell inks, factors such as drying rate, particle size, and agglomeration influence how the materials distribute themselves throughout the thickness of the dired film. Thus far there have only been limited studies to understand how process conditions such as ink formulation and drying temperature influence the distribution of ionomer and catalyst coated using scalable methods. This understanding is especially important for GDEs since it is known that having a sufficient amount of ionomer at the catalyst layer-membrane interface is critical for high performance. In this study we have focused on determining how the ratio of water to 1-propanol in the catalyst ink ink and drying temperature influence the distribution of ionomer throughout the thickness of the catalyst layer. Using a combination of Kelvin probe and x-ray photoelectron spectroscopy we show that an ionomer-rich surface is promoted by a higher drying rate and a water-rich catalyst ink. In contrast, a 1-propanol catalyst ink leads to a lower concentration of ionomer on the top surface. Using x-ray computed tomography, we are able to characterize the ionomer distribution throughout the thickness of the layer. We find that, in addition to promoting an ionomer-rich top surface, water-rich inks lead to a more homogenous distribution of ionomer, whereas a 1-propanol-rich ink leads to a more irregular distribution. It is found that MEA performance is improved by selecting conditions and ink formulations that promote ionomer enrichment at the top surface to facilitate a good interface with the membrane. MEAs prepared with a 75 wt% water catalyst ink with a 0.9 I/C have equivalent performance to spray-coated GDEs. Critically, these R2R-coated GDEs do not need an additional ionomer overlayer like the spray-coated GDEs do, reducing the number of processing steps in a manufacturing setting. This work shows that with the appropriate selection of materials, ink formulation, and processing conditions gas-diffusion electrodes are a viable pathway for fuel cell manufacturing.

ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATION↗
DOE OSTI2024

Multilevel well modeling in aggregation-based nonlinear multigrid for multiphase flow in porous media

A full approximation scheme (FAS) nonlinear multigrid solver for two-phase flow and transport problems driven by wells with multiple perforations is developed here. It is an extension to our previous work on FAS solvers for diffusion and transport problems. The solver is applicable to discrete problems defined on unstructured grids as the coarsening algorithm is aggregation-based and algebraic. To construct coarse basis that can better capture the radial flow near wells, coarse grids in which perforated well cells are not near the coarse-element interface are desired. This is achieved by an aggregation algorithm proposed in this paper that makes use of the location of well cells in the cell-connectivity graph. Numerical examples in which the FAS solver is compared against Newton's method on benchmark problems are given. In particular, for a refined version of the SAIGUP model, the FAS solver is at least 35% faster than Newton's method for time steps with a CFL number greater than 10.

58 GEOSCIENCES↗
DOE OSTI2022

Dynamic heterogeneity and nanophase separation in rubber-toughened amine-cured highly cross-linked polymer networks

Solid state nuclear magnetic resonance (NMR) spectroscopy and small-to wide-angle X-ray scattering (SWAXS) methods were used to characterize the heterogeneous dynamics and polymer domain structure in rubber modified thermoset materials containing the diglycidyl ether of bisphenol A (DGEBA) epoxy resin and a mixture of Jeffamine reactive rubber and 4,4-diaminodicyclohexylmethane (PACM) amine curing agent. The polymer chain dynamics and morphologies as a function of the PACM/Jeffamine ratio were determined. Using dipolar-filtered NMR experiments, the resulting networks are shown to be composed of mobile and rigid regions that are separated on nanometer length scales, along with a dynamically immobilized interface region. Proton NMR spin diffusion experiments measured the dimensions of the mobile phase to range between 9 and 66 nm and varied with the relative PACM concentration. Solid state 13 C magic angle spinning NMR experiments show that the highly mobile phase is composed entirely of the dynamically flexible polyether chains of the Jeffamine rubber, the immobilized interface region is a mixture of DGEBA, PACM, and the Jeffamine rubber, with the PACM cross-linked to DGEBA predominantly residing in the rigid phase. The SWAXS results showed compositional nanophase separation spanning the 11–77 nm range. These measurements of the nanoscale compositional and dynamic heterogeneity provide molecular level insight into the very broad and controllable glass transition temperature distributions observed for these highly cross-linked polymer networks.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2023

COVIDScholar: An automated COVID-19 research aggregation and analysis platform

The ongoing COVID-19 pandemic produced far-reaching effects throughout society, and science is no exception. The scale, speed, and breadth of the scientific community’s COVID-19 response lead to the emergence of new research at the remarkable rate of more than 250 papers published per day. This posed a challenge for the scientific community as traditional methods of engagement with the literature were strained by the volume of new research being produced. Meanwhile, the urgency of response lead to an increasingly prominent role for preprint servers and a diffusion of relevant research through many channels simultaneously. These factors created a need for new tools to change the way scientific literature is organized and found by researchers. With this challenge in mind, we present an overview of COVIDScholar https://covidscholar.org/, an automated knowledge portal which utilizes natural language processing (NLP) that was built to meet these urgent needs. The search interface for this corpus of more than 260,000 research articles, patents, and clinical trials served more than 33,000 users at an average of 2,000 monthly active users and a peak of more than 8,600 weekly active users in the summer of 2020. Additionally, we include an analysis of trends in COVID-19 research over the course of the pandemic with a particular focus on the first 10 months, which represents a unique period of rapid worldwide shift in scientific attention.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

2D-EFICACY: Control of Metastable 2D Carbide[1]Chalcogenide Heterolayers: Strain and Moire Engineering

The experimental isolation of graphene led to the discovery of an entirely new world of two-dimensional (2D) materials in which the 2D nature often leads to emergent behaviors not seen in bulk systems. 2D transition metal dichalcogenides (TMDs) exhibit physico-chemical properties that depend on the transition metal, polymorph, thickness, and presence and type of defects. Recently, a group of thin (10-100nm) transition metal carbides (TMCs), such as Mo2C, has been synthesized that exhibit a thickness-dependent superconducting critical temperature (Tc). These thin TMCs are different from MXenes, another class of 2D materials consisting of few layers of nitrides or carbides (<5nm) produced by chemical etching and delamination. The goal of this renewal proposal is to combine experiment and computation to synthesize and elucidate the guiding principles that control the growth, orientation and strain of heterostacks of thin TMCs and TMDs composed with Nb, Ti and W. We expect to stabilize metastable hybrid phases of TMCs sandwiched between TMDs (H-TMD/Cs) with unprecedented physico-chemical properties. As part of previous DOE-funded work by the Terrones/Sinnott groups, thin (10-100 nm thick) Mo2C flakes were successfully synthesized by chemical vapor deposition (CVD). By subsequently exposing Mo2C to H2S, partial chalcogenization was achieved, resulting in heterostacks of MoCx phases and MoS2. The formation of MoS2 led to a deficiency of Mo atoms in the underlying Mo2C, resulting in an inhomogeneous phase change from α-Mo2C to γ’-MoCx and then to γ-MoC. The γ’-MoCx is a strained metastable phase and the heterostack of all three phases demonstrated an increased Tc relative to that of α-Mo2C, from 4 to 6K; its interleaved layered structure consisting of superconducting and semiconducting phases is ideal for future studies of Josephson junction series arrays. Moiré patterns in these heterostacked systems could result in new phenomena, as moiré patterns in bilayer graphene showed unconventional superconductivity and moiré excitons have been observed in twisted TMD heterobilayers. The scientific hypothesis of the proposed synergistic computational and experimental research is that orientation and strain control within confined thin metastable TMCs, sandwiched by stable phases of TMCs and layered TMDs, will depend on kinetic and thermodynamic “knobs” that include fast temperature changes, chalcogen diffusion through preferred crystallographic planes, reaction times, pressure, reactive atmosphere, precursors, and surfactants, which will also tailor properties such as superconductivity, magnetism, ferroelectricity, piezoelectricity, and catalytic performance. We will develop the guiding principles for the synthesis and stabilization of metastable H-TMD/Cs based on Nb, Ti and W. In order to validate the hypothesis, four tasks are proposed: The first task will synthesize ultra-thin TMCs based on Nb, W and Ti, by: 1) adapting the CVD method used for Mo2C, 2) plasma assisted CVD, 3) defect-mediated CVD processes, and 4) cryo-milling of carbide powders. The second task will accomplish the synthesis and basic physico-chemical characterizations of H-TMD/Cs by chalcogenization of the materials synthesized in task one, and by carbonization of TMDs. H-TMD/Cs will also be investigated for their suitability in energy conversion applications such as supercapacitors, Li and multivalent ion batteries, and electrocatalysts, topics of interest to DoE. These tasks will be carried out in close conjunction with density functional theory (DFT) calculations with insights into energetics, lattice parameters, stability, phase diagrams, band structures, and density of states of H-TMD/Cs. The third task will characterize and evaluate strain and moiré patterns at the interfaces of different H-TMD/Cs by high-resolution scanning transmission electron microscopy (HR-STEM), scanning tunneling microscopy (STM), and conductive tip atomic force microscopy. Nudged elastic band calculations with DFT will be performed to understand the chalcogen diffusion process, which will provide insights into the interfaces between different phases of TMCs and TMDs. The fourth task aims at quantifying the stability and dynamics of H-TMD/Cs by in-situ TEM and Raman studies under heating, strain, and electrical biasing. Phonon calculations using DFT will provide a basis for interpreting Raman spectra. This coherent framework involving synthesis, characterization, and computation will result in a broad scientific impact for energy related applications. The ability to develop new H-TMD/Cs will enhance a range of applications that include batteries, catalysts, switches, sensors, quantum computing components and smart coatings.

2-Dimensional materials↗
DOE OSTI2024

2D-EFICACY: Control of Metastable 2D Carbide Chalcogenide Heterolayers: Strain and Moire Engineering

The experimental isolation of graphene led to the discovery of an entirely new world of two-dimensional (2D) materials in which the 2D nature often leads to emergent behaviors not seen in bulk systems. 2D transition metal dichalcogenides (TMDs) exhibit physico-chemical properties that depend on the transition metal, polymorph, thickness, and presence and type of defects. Recently, a group of thin (10-100nm) transition metal carbides (TMCs), such as Mo 2 C, has been synthesized that exhibit a thickness-dependent superconducting critical temperature (Tc). These thin TMCs are different from MXenes, another class of 2D materials consisting of few layers of nitrides or carbides (<5nm) produced by chemical etching and delamination. The goal of this renewal proposal is to combine experiment and computation to synthesize and elucidate the guiding principles that control the growth, orientation and strain of heterostacks of thin TMCs and TMDs composed with Nb, Ti and W. We expect to stabilize metastable hybrid phases of TMCs sandwiched between TMDs (H-TMD/Cs) with unprecedented physico-chemical properties. As part of previous DOE-funded work by the Terrones/Sinnott groups, thin (10-100 nm thick) Mo 2 C flakes were successfully synthesized by chemical vapor deposition (CVD). By subsequently exposing Mo 2 C to H 2 S, partial chalcogenization was achieved, resulting in heterostacks of MoCx phases and MoS 2 . The formation of MoS 2 led to a deficiency of Mo atoms in the underlying Mo 2 C, resulting in an inhomogeneous phase change from α-Mo 2 C to γ’-MoCx and then to γ-MoC. The γ’-MoCx is a strained metastable phase and the heterostack of all three phases demonstrated an increased Tc relative to that of α-Mo 2 C, from 4 to 6K; its interleaved layered structure consisting of superconducting and semiconducting phases is ideal for future studies of Josephson junction series arrays. Moiré patterns in these heterostacked systems could result in new phenomena, as moiré patterns in bilayer graphene showed unconventional superconductivity and moiré excitons have been observed in twisted TMD heterobilayers. The scientific hypothesis of the proposed synergistic computational and experimental research is that orientation and strain control within confined thin metastable TMCs, sandwiched by stable phases of TMCs and layered TMDs, will depend on kinetic and thermodynamic “knobs” that include fast temperature changes, chalcogen diffusion through preferred crystallographic planes, reaction times, pressure, reactive atmosphere, precursors, and surfactants, which will also tailor properties such as superconductivity, magnetism, ferroelectricity, piezoelectricity, and catalytic performance. We will develop the guiding principles for the synthesis and stabilization of metastable H-TMD/Cs based on Nb, Ti and W. In order to validate the hypothesis, four tasks are proposed: The first task will synthesize ultra-thin TMCs based on Nb, W and Ti, by: 1) adapting the CVD method used for Mo 2 C, 2) plasma assisted CVD, 3) defect-mediated CVD processes, and 4) cryo-milling of carbide powders. The second task will accomplish the synthesis and basic physico-chemical characterizations of H-TMD/Cs by chalcogenization of the materials synthesized in task one, and by carbonization of TMDs. H-TMD/Cs will also be investigated for their suitability in energy conversion applications such as supercapacitors, Li and multivalent ion batteries, and electrocatalysts, topics of interest to DoE. These tasks will be carried out in close conjunction with density functional theory (DFT) calculations with insights into energetics, lattice parameters, stability, phase diagrams, band structures, and density of states of H-TMD/Cs. The third task will characterize and evaluate strain and moiré patterns at the interfaces of different H-TMD/Cs by high-resolution scanning transmission electron microscopy (HR-STEM), scanning tunneling microscopy (STM), and conductive tip atomic force microscopy. Nudged elastic band calculations with DFT will be performed to understand the chalcogen diffusion process, which will provide insights into the interfaces between different phases of TMCs and TMDs. The fourth task aims at quantifying the stability and dynamics of H-TMD/Cs by in-situ TEM and Raman studies under heating, strain, and electrical biasing. Phonon calculations using DFT will provide a basis for interpreting Raman spectra. This coherent framework involving synthesis, characterization, and computation will result in a broad scientific impact for energy related applications. The ability to develop new H-TMD/Cs will enhance a range of applications that include batteries, catalysts, switches, sensors, quantum computing components and smart coatings.

2-Dimensional Materials↗
DOE OSTI2025

In Situ Neutron Reflectometry Reveals the Interfacial Microenvironment Driving Electrochemical Ammonia Synthesis

Electrified interfaces are critical to the performance of energy systems and often demonstrate substantial complexity under operating conditions. A nanoscale understanding of the interfacial microenvironment, i.e., the solid-electrolyte interphase (SEI), in lithium-mediated nitrogen reduction (Li–N 2 R) is key for realizing efficient ammonia (NH 3 ) production. Herein, we used time-resolved neutron reflectometry (NR) to observe SEI formation under Li–N 2 R conditions. We found that the LiBF 4 -based electrolyte provided a substantially more well-defined SEI layer than previous SEI NR interrogations that used LiClO 4 , highlighting the underlying chemistry that dictates electrolyte design and enabling new NR-based studies. Using in situ NR, we found that the LiBF 4 -derived SEI under Li–N 2 R conditions comprises a thick, diffuse outer layer and a thin, compact inner layer at low current cycling (<2 mA/cm 2 ), revealing a structure which ex situ studies have not been able to probe. Increased current cycling and sustained current cycling led to the merging of the layers into a single-layer SEI. Here, we used isotope contrast methods with d 6 -EtOH and d 8 -THF to drive time-resolved tracking of SEI growth at low current cycling, revealing that the proton donor modifies the inner layer, and the solvent modifies the outer layer. Li dendritic growth was observed in the absence of a proton donor. Neutron absorption also indicated the presence of boron in the SEI, underscoring the value of neutron-based interrogation. Our results inform Li-based systems and reaction microenvironments, and these methods can be applied broadly to interfacial energy technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

An Interactive, Design and Educational Tool for Supersonic External-Compression Inlets

A workstation-based interactive design tool called VU-INLET was developed for the inviscid flow in rectangular, supersonic, external-compression inlets. VU-INLET solves for the flow conditions from free stream, through the supersonic compression ramps, across the terminal normal shock region and the subsonic diffuser to the engine face. It calculates the shock locations, the capture streamtube, and the additive drag of the inlet. The inlet geometry can be modified using a graphical user interface and the new flow conditions recalculated interactively. Free stream conditions and engine airflow can also be interactively varied and off-design performance evaluated. Flow results from VU-INLET can be saved to a file for a permanent record, and a series of help screens make the simulator easy to learn and use. This paper will detail the underlying assumptions of the models and the numerical methods used in the simulator.

Benson, Thomas J.↗
DOE OSTI2026

Multiscale Characterization of Electrode-Induced Degradation in Perovskite Solar Cells

The stability of metal-halide-perovskite (MHP) solar cells must be understood and improved for the commercial viability of MHP technologies. Here, we apply multiscale characterization methods to study degradation modes, specifically electrode corrosion, for p-i-n MHP partial device stacks and full devices that are stored in the dark under an inert atmosphere. Our multiscale characterization approaches include full-device electro-optical performance using current-voltage (JV) curves and spatial imaging with electroluminescence (EL) and photoluminescence (PL). We further correlate interface properties using cross-sectional Kelvin probe force microscopy, which maps the nanoscale electric field properties, and electron microscopy, which demonstrates structural and chemical features. Devices stored as a full device stack degrade primarily by metal (Ag) electrode diffusion into the absorber, with formation of AgI byproducts and Ag accumulation near the indium tin oxide (ITO) contact. This causes decomposition of the perovskite absorber domains, loss of the potential drop at the electron transport layer (ETL)/perovskite interface near the metal contact, and increased equivalent resistance at the perovskite/hole transport layer (HTL) interface near the ITO contact. The devices stored without metal show a different degradation pathway dominated by corrosion of the ITO, creating voids at the ITO electrode surface with diffusion of In and Sn into the absorber. We conclude that metal electrode-induced degradation is the most severe degradation pathway under dark storage, but that ITO corrosion and absorber instability must also be mitigated. We further demonstrate mitigation of these degradation pathways by changes to the device stack, including a SnO x blocking layer at the ETL side and replacing ITO with FTO at the HTL side. These results provide a useful demonstration of specific dark degradation pathways at each electrode interface, as well as a unique multiscale example that links degradation of chemical, structural, and electrical interface properties to the full-device electro-optical characteristics.

14 SOLAR ENERGY↗
DOE OSTI2022

Direct and Sensitive Electrochemical Evaluation of Pramipexole Using Graphitic Carbon Nitride (gCN) Sensor

Pramipexole (PMXL) belongs to the benzothiazole class of aromatic compounds and is used in treating Parkinson’s disease; however, overdosage leads to some abnormal effects that could trigger severe side effects. Therefore, it demands a sensitive analytical tool for trace level detection. In this work, we successfully developed an electrochemical sensor for the trace level detection of PMXL, using the voltammetric method. For the analysis, graphitic carbon nitride (gCN) was opted and synthesized by using a high-temperature thermal condensation method. The synthesized nanoparticles were employed for surface characterization, using transmission electron microscopy (TEM), X-ray diffraction (XRD), and atomic force microscopy (AFM) techniques. The electrochemical characterization of the material was evaluated by using the electrochemical impedance spectroscopy (EIS) technique to evaluate the solution–electrode interface property. The cyclic voltammetry (CV) behavior of PMXL displayed an anodic peak in the forward scan, indicating that PMXL underwent electrooxidation, and an enhanced detection peak with lower detection potential was achieved for gCN-modified carbon paste electrode (gCN·CPE). The influence of different parameters on the electrochemical behavior was analyzed, revealing the diffusion governing the electrode process with an equal number of hydronium ions and electron involvement. For the fabricated gCN·CPE, good linearity range was noticed from 0.05 to 500 µM, and a lower detection limit (L D ) of 0.012 µM was achieved for the selected concentration range (0.5 to 30 µM). Selectivity of the electrode in PMXL detection was investigated by conducting an interference study, while the tablet sample analysis demonstrates the sensitive and real-time application of the electrode. The good recovery values for the analysis illustrate the efficiency of the electrode for PMXL analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Locally developed electronic conduction in a yttria stabilized zirconia (YSZ) electrolyte for durable solid oxide fuel cells

In the present work, anode supported solid oxide fuel cells (SOFCs) based on a YSZ single layer electrolyte with locally developed electronic conduction near the anode side is fabricated to improve the cell durability under negative voltage conditions. A ceria added Ni/YSZ anode functional layer is utilized and co-sintered with the YSZ electrolyte at 1460 °C. Cerium ions are diffused into a limited region of the YSZ layer (~6 μm), establishing a ceria doped YSZ region. We successfully demonstrate that this YSZ single layer cell can operate stably for ~100 h without anode/electrolyte interface delamination, even under a negative voltage. The embedded probe technique and theoretical calculation reveal that the internal oxygen partial pressure (near the anode) is substantially reduced (stabilized) by means of local electronic conduction of ceria doped YSZ. Thus, we propose an SOFC electrolyte design comprising a YSZ single layer with a ceria added anode functional layer (locally formed ceria doped YSZ in the electrolyte) as a protection method against fuel cell stack degradation caused by negative voltage operation.

25 ENERGY STORAGE↗
DOE OSTI2023

A one-pot synthesis of high-density biofuels through bifunctional mesoporous zeolite-encapsulated Pd catalysts

Developing a powerful bifunctional catalyst for tandem reactions is essential to the future carbon neutrality by reducing the energy consumption in chemical industries. Herein, mesoporous zeolite-encapsulated palladium (Pd) nanoparticles (Pd@meso-ZSM-5) synthesized via emulsification-demulsification followed by a dry-gel transformation method were demonstrated to have a remarkable catalytic performance for a one-pot multiple tandem reaction of cyclic ketones (cyclopentanone, cyclohexanone) to bicyclic alkanes (bicyclopentane, bicyclohexane). Compared with supported catalysts (Pd/meso-ZSM-5) and microporous zeolite-encapsulated catalysts (Pd@ZSM-5) with a primary product of monocyclic alkane (cyclopentane, cyclohexane), Pd@meso-ZSM-5 shows much higher catalytic efficiency for the bicyclic alkanes synthesis, which was previously unattainable in the conventional two-step synthesis route. Controlled experiments and detailed characterizations show that mesoporosity provides sufficient space for the generation and diffusion of large molecular intermediates, and the intimate acid-Pd interface promotes the conversion of intermediates. As a result, this work proposes a design strategy for a mesoporous zeolite-encapsulated metal catalyst, which efficiently provides cooperative acid-hydrogenation catalysis.

09 BIOMASS FUELS↗
NASA NTRS1983

Advanced methods for preparation and characterization of infrared detector materials

Crystals were prepared by the Bridgman-Stockbarger method with a wide range of crystal growth rates and temperature gradients adequate to prevent constitutional supercooling under diffusion-limited, steady-state, growth conditions. The longitudinal compositional gradients for different growth conditions and alloy compositions were calculated and compared with experimental data to develop a quantitative model of solute redistribution during the crystal growth of the alloys. Measurements were performed to ascertain the effect of growth conditions on radial compositional gradients. The pseudobinary HgTe-CdTe constitutional phase diagram was determined by precision differential-thermal-analysis measurements and used to calculate the segregation coefficient of Cd as a function of x and interface temperature. Experiments were conducted to determine the ternary phase equilibria in selected regions of the Hg-Cd-Te constitutional phase diagram. Electron and hole mobilities as functions of temperature were analyzed to establish charge-carrier scattering probabilities. Computer algorithms specific to Hg(1-x)CdxTe were developed for calculations of the charge-carrier concentration, charge-carrier mobilities, Hall coefficient, and Dermi Fermi energy as functions of x, temperature, ionized donor and acceptor concentrations, and neutral defect concentrations.

Broerman, J. G.↗
DOE OSTI2024

The Effective Fragment Molecular Orbital Method: Achieving High Scalability and Accuracy for Large Systems

The effective fragment molecular orbital (EFMO) method has been developed to predict the total energy of a very large molecular system accurately (with respect to the underlying quantum mechanical method) and efficiently by taking advantage of the locality of strong chemical interactions and employing a two-level hierarchical parallelism. The accuracy of the EFMO method is partly attributed to the accurate and robust intermolecular interaction prediction between distant fragments, in particular, the many-body polarization and dispersion effects, which require the generation of static and dynamic polarizability tensors by solving the coupled perturbed Hartree–Fock (CPHF) and time-dependent HF (TDHF) equations, respectively. Solving the CPHF and TDHF equations is the main EFMO computational bottleneck due to the inefficient (serial) and I/O-intensive implementation of the CPHF and TDHF solvers. In this work, the efficiency and scalability of the EFMO method are significantly improved with a new CPU memory-based implementation for solving the CPHF and TDHF equations that are parallelized by either message passing interface (MPI) or hybrid MPI/OpenMP. Here, the accuracy of the EFMO method is demonstrated for both covalently bonded systems and noncovalently bound molecular clusters by systematically examining the effects of basis sets and a key distance-related cutoff parameter, R cut . R cut determines whether a fragment pair (dimer) is treated by the chosen ab initio method or calculated using the effective fragment potential (EFP) method (separated dimers). Decreasing the value of Rcut increases the number of separated (EFP) dimers, thereby decreasing the computational effort. It is demonstrated that excellent accuracy (<1 kcal/mol error per fragment) can be achieved when using a sufficiently large basis set with diffuse functions coupled with a small R cut value. With the new parallel implementation, the total EFMO wall time is substantially reduced, especially with a high number of MPI ranks. Given a sufficient workload, nearly ideal strong scaling is achieved for the CPHF and TDHF parts of the calculation. For the first time, EFMO calculations with the inclusion of long-range polarization and dispersion interactions on a hydrated mesoporous silica nanoparticle with explicit water solvent molecules (more than 15k atoms) are achieved on a massively parallel supercomputer using nearly 1000 physical nodes. In addition, EFMO calculations on the carbinolamine formation step of an amine-catalyzed aldol reaction at the nanoscale with explicit solvent effects are presented.

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

Two Dimensional Topology Optimization of Heat Exchangers with the Density and Level-Set Methods

We design heat exchangers using two topology optimization approaches: the density, i.e. volume fraction and level set methods. Our goal is to maximize the heat exchange between two fluids in separate channels while constraining the pressure drop across each channel. The heat exchanger is modeled with a coupled thermal-flow formulation. The flow is governed by an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation with high Peclet number. We solve one set of Stokes-Brinkman equations per fluid. Each Brinkman term in the flow equation serves to model the other phase as a solid, thereby preventing mixing. We first represent the solid and fluid phases using a volume fraction variable and apply a SIMP-like penalization in the Brinkman term to drive the optimization to a discrete design. The cost and constraint function derivatives are automatically calculated with the library pyadjoint and the optimization is performed by the Method of Moving Asymptotes. In a second optimization formulation, we use the level set approach to define the interface that separates the two fluids. Pyadjoint calculates the shape derivatives of the cost and constraint functions and the Hamilton-Jacobi advects the interface, allowing for topological changes. We present results in two dimensions and discuss the advantages and disadvantages of each approach.

42 ENGINEERING↗
DOE OSTI2024

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↗