Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum modeling”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Solvated Nuclear–Electronic Orbital Structure and Dynamics

Nonadiabatic dynamical processes such as proton-coupled electron transfer and excited state intramolecular proton transfer have been the subject of much research. One of the promising theoretical methods to describe these processes is the nuclear–electronic orbital (NEO) approach. This approach inherently accounts for nuclear quantum effects within quantum chemistry calculations, and it has recently been extended to directly simulate nonadiabatic processes with the development of real-time NEO methods. These processes can also be significantly dependent on the surrounding chemical environment, however, and capturing the effects of the environment is often necessary for analyzing experimentally relevant systems. This work couples the NEO density functional theory and real-time time-dependent density functional theory approaches with solvation through the polarizable continuum model. Here, the effects of this coupling are investigated for ground state properties, solvent-dependent vibrational frequencies, and direct excited state intramolecular proton transfer dynamics.

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

Anion–Cation–Anion Ion Triplet Characterization by Computation and Photoelectron Spectroscopy

Ion triplets of the chloride salts of two commonly used weakly coordinating cations are reported (i.e., Cl – ·NMe 4 + Cl – (1 – ) and Cl – ·PPh 4 + Cl – (2 – )). Negative ion photoelectron spectra at 20 K afford vertical and adiabatic detachment energies of 5.18 and 5.00 eV (1 – ) and 5.03 and 4.70 eV (2 – ), respectively. These results are well reproduced by coupled cluster calculations with single, double, and perturbative triple excitations (CCSD(T)) whereas M06-2X is systematically too small by ~0.3 eV (i.e., 7 kcal mol –1 ). The structures of both 1 – and 2 – have five or six C–H···Cl – interactions that stabilize these cluster anions by 32 (1 – ) and 27 (2 – ) kcal mol –1 as given by their chloride dissociation enthalpies. In conclusion, these values drop to 7.4 and 3.8 kcal mol –1 in dichloromethane based up conductor-like polarizable continuum model calculations and suggest that X – ·M + X – ion triplets with a weakly coordinating cation maybe the reactive form of salts under some conditions.

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

On the Nature of Field-Enhanced Water Dissociation in Bipolar Membranes

Bipolar membranes (BPMs) possess the potential to optimize pH environments for electrochemical synthesis applications when employed in reverse bias. Unfortunately, the performance of BPMs in reverse bias has long been limited by the rate of water dissociation (WD) occurring at the interface of the BPM. Herein, we develop a continuum model of the BPM that agrees with experiment to understand and enhance WD catalyst performance by considering multiple kinetic pathways for WD in the BPM junction catalyst layer. Here, the model reveals that WD catalysts with a more highly alkaline or acidic pH at the point of zero charge (pH PZC ) exhibit accelerated WD kinetics because the more acidic or alkaline pH PZC catalysts possess greater surface charge, enhancing the local electric field and rate of WD. The model is then employed to explore the sensitivity of the BPM performance to various BPM physical parameters. Finally, the model is used to simulate the operation of bimetallic WD catalysts, demonstrating that an optimal bimetallic catalyst has an acidic pH PZC catalyst matched with the cation-exchange layer and an alkaline pH PZC catalyst matched with the anion-exchange layer. The study provides insight into the operation of BPM WD catalysts and gives direction toward the development of next-generation WD catalysts for optimal BPM performance under water-splitting and related conditions.

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

Solvent Induced Proton Polarization within the Nuclear–Electronic Orbital Framework

To explicitly account for nuclear quantum effects and solvent environments in simulations of chemical processes, the nuclear-electronic orbital approach is coupled with a polarizable continuum model (PCM). This NEO–PCM approach is used to explore the influence of solvation on nuclear polarization through applications to a water dimer and a set of protonated water tetramers. Nuclear polarization in these species is analyzed in terms of changes in proton density and oxygen–hydrogen bond length. Solvation is shown to enhance nuclear polarization with increasing dielectric constant. For the water dimer, the internal, hydrogen-bonded proton is shown to polarize more than the external, free proton. Moreover, proton quantization leads to greater solvent polarization through their mutual polarization. Furthermore, these calculations highlight the complex interplay among electronic, nuclear, and solvent polarization in chemical systems.

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

Stress-Dependent Chemo-Mechanical Performance of Amorphous Si Anodes for Li-Ion Batteries upon Lithiation

Alloying-type anodes are significantly governed by their chemo-mechanical performance during the electrochemical cycling. The reaction-induced huge volumetric change of these anodes may cause material degradation and failure under mechanical constraints. Here, we investigate the stress-dependent lithiation behavior of amorphous Si (a-Si) anodes using molecular dynamics simulations. It is indicated that a-Si anodes can sustain higher hydrostatic stress than biaxial/uniaxial ones without the occurrence of mechanical failure. Thermodynamic and electrochemical calculations demonstrate that although the lithiation procedure also affects the thermodynamic stability of a-Si anodes, it is mainly dominated by the external mean stresses. Compressive stress is confirmed to destabilize a-Si anodes and further trigger their capacity fading. Compared with our atomistic simulations, previous continuum models underestimate the open-cell potentials of a-Si anodes, due to their ignored large volumetric deformation at higher stresses and Li concentrations. Finally, this computational study provides the intensive atomic-level understanding of the stress-dependent lithiation behavior of a-Si anodes.

25 ENERGY STORAGE↗
DOE OSTI2020

Understanding Multi-Ion Transport Mechanisms in Bipolar Membranes

Bipolar membranes (BPMs) have the potential to become critical components in electrochemical devices for a variety of electrolysis and electrosynthesis applications. Because they can operate under large pH gradients, BPMs enable favorable environments for electrocatalysis at the individual electrodes. Critical to the implementation of BPMs in these devices is understanding the kinetics of water dissociation that occurs within the BPM as well as the co- and counter-ion crossover through the BPM, which both present significant obstacles to developing efficient and stable BPM-electrolyzers. In this study, a continuum model of multi-ion transport in a BPM is developed and fit to experimental data. Specifically, concentration profiles are determined for all ionic species, and the importance of a water-dissociation catalyst is demonstrated. The model describes internal concentration polarization and co- and counter-ion crossover in BPMs, determining the mode of transport for ions within the BPM and revealing the significance of salt-ion crossover when operated with pH gradients relevant to electrolysis and electrosynthesis. Lastly, a sensitivity analysis reveals that the performance and lifetime of BPMs can be improved substantially by using of thinner dissociation catalysts, managing water transport, modulating the thickness of the individual layers in the BPM to control salt-ion crossover, and increasing the ion-exchange capacity of the ion-exchange layers in order to amplify the water-dissociation kinetics at the interface.

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

Insights into the Li Diffusion Mechanism in Si/C Composite Anodes for Lithium-Ion Batteries

Recently, Si/C composite materials have attracted enormous research interest as the most promising candidates for the anodes of next-generation lithium-ion batteries, owing to their high energy density and mechanical buffering property. However, the fundamental mechanism of Li diffusion behavior in various Si/C composite materials remains unclear, with our understanding limited by experimental techniques and continuum modeling methodologies. Herein, the atomic behavior of Li diffusion in the Si/C composite material is studied within the framework of density functional theory. Two representative structural mixing formats, that is, simple mixture mode and core-shell mode, are modeled and compared. We discover that the carbon material increases Li diffusion in silicon from 7.75 x 10 -5 to 2.097 x 10 -4 cm 2 /s. The boost is about 50% more obvious in the mixture mode, while the core-shell structure shows more dependence on the atomic structures of the carbon layer. These results offer new insights into Li diffusion behavior in Si/C composites and unlock the enhancing mechanism for Li diffusion in Si/C. Furthermore, this understanding facilitates the modeling of batteries with composite anodes and will guide the corresponding structure designs for robust and high-energy-density batteries.

25 ENERGY STORAGE↗
DOE OSTI2023

Theory and Simulation of Metal–Insulator–Semiconductor (MIS) Photoelectrodes

A metal–insulator–semiconductor (MIS) structure is an attractive photoelectrode-catalyst architecture for promoting photoelectrochemical reactions, such as the formation of H 2 by proton reduction. The metal catalyzes the generation of H 2 using electrons generated by photon absorption and charge separation in the semiconductor. The insulator layer between the metal and the semiconductor protects the latter element from photo-corrosion and, also, significantly impacts the photovoltage at the metal surface. Understanding how the insulator layer determines the photovoltage and what properties lead to high photovoltages is critical to the development of MIS structures for solar-to-chemical energy conversion. Herein, we present a continuum model for charge-carrier transport from the semiconductor to the metal with an emphasis on mechanisms of charge transport across the insulator. The polarization curves and photovoltages predicted by this model for a Pt/HfO 2 /p-Si MIS structure at different HfO 2 thicknesses agree well with experimentally measured data. The simulations reveal how insulator properties (i.e., thickness and band structure) affect band bending near the semiconductor/insulator interface and how tuning them can lead to operation closer to the maximally attainable photovoltage, the flat-band potential. This phenomenon is understood by considering the change in tunneling resistance with insulator properties. The model shows that the best MIS performance is attained with highly symmetric semiconductor/insulator band offsets (e.g., BeO, MgO, SiO 2 , HfO 2 , or ZrO 2 deposited on Si) and a low to moderate insulator thickness (e.g., between 0.8 and 1.5 nm). Beyond 1.5 nm, the density of filled interfacial trap sites is high and significantly limits the photovoltage and the solar-to-chemical conversion rate. These conclusions are true for photocathodes and photoanodes. This understanding provides critical insight into the phenomena enhancing and limiting photoelectrode performance and how this phenomenon is influenced by insulator properties. The study gives guidance toward the development of next-generation insulators for MIS structures that achieve high performance.

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

Dynamic Boundary Layer Simulation of Pulsed CO 2 Electrolysis on a Copper Catalyst

Pulsed electrolysis has been demonstrated to improve the faradaic efficiency (FE) to C 2+ products during the electrochemical reduction of CO 2 over a Cu catalyst, but the nature of this enhancement is poorly understood. Herein, we developed a time-dependent continuum model of pulsed CO 2 electrolysis on Cu in 0.1 M CsHCO 3 that faithfully represents the experimentally observed effects of pulsed electrolysis. This work shows that pulsing results in dynamic changes in the pH and CO 2 concentration near the Cu surface, which lead to an enhanced C 2+ FE as a consequence of repeatedly accessing a transient state of heightened pH and CO 2 concentration at high cathodic overpotential. Using these insights, a variety of pulse shapes were explored to establish operating conditions that maximize the rate of C 2+ product formation and minimize the rates of H 2 and C 1 product formation.

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

Electrolytic Methane Production from Reactive Carbon Solutions

Here, we report an electrochemical reactor that converts 3.0 M KHCO 3 into methane at the cathode, and oxidizes water at the anode. The molar ratio of methane product to unreacted CO 2 gas (defined herein as "methane yield") was measured to be 34% at a partial current density of 120 mA cm -2 . The highest previously reported CO 2 -to-methane yield is 3%. Our reactor achieved this improvement in methane yield because it is fed with 3.0 M KHCO 3 , a type of reactive carbon solution, rather than gaseous CO 2 . The reactor uses H + delivered by a bipolar membrane to form CO 2 at the cathode. This CO 2 is subsequently reduced into methane. A cationic surfactant added to the catholyte suppressed hydrogen evolution and increased methane formation. A 1D continuum model confirmed that H + from the membrane promotes the formation of methane over multicarbon products at the cathode. These findings present design principles for electrochemical methane synthesis.

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

Asymmetric Bipolar Membrane for High Current Density Electrodialysis Operation with Exceptional Stability

Bipolar membranes (BPMs) enable isolated acidic/alkaline regions in electrochemical devices, facilitating optimized environments for electrochemical separations and catalysis. For economic viability, BPMs must attain stable, high current density operation with low overpotentials in a freestanding configuration. We report an asymmetric, graphene oxide (GrOx)-catalyzed BPM capable of freestanding electrodialysis operation at 1 A cm–2 with overpotentials <250 mV. Use of a thin anion-exchange layer improves water transport while maintaining near unity Faradaic efficiency for acid and base generation. Voltage stability exceeding 1100 h with an average drift of 70 μV/h at 80 mA cm–2 and 100 h with an average drift of −300 μV/h at 500 mA cm–2 and implementation in an electrodialysis stack demonstrate real-world applicability. Continuum modeling reveals that water dissociation in GrOx BPMs is both catalyzed and electric-field enhanced, where low pK a moieties on GrOx enhance local electric fields and high pK a moieties serve as active sites for surface-catalyzed water dissociation. These results establish commercially viable BPM electrodialysis and provide fundamental insight to advance design of next-generation devices.

Lucas, Éowyn↗
DOE OSTI2026

Spatially Mapping the CO 2 Alkaline Sorbent Diffuse Microenvironment Using Operando Raman Spectroscopy

When designing a chemical process, the local balance of transport and kinetics, collectively referred to as the diffuse microenvironment, plays a critical role in performance but is difficult to directly observe. This work demonstrates a method of two-dimensional spatial chemical mapping of the diffuse microenvironment in the context of alkali metal hydroxide direct air capture of carbon dioxide using a custom operando gas-absorption flow cell along with confocal Raman spectroscopy. Notably, we observe the concentration boundary layer near the gas–liquid interface and elucidate the interplay of carbonate and bicarbonate ions within it while inferring local hydroxide depletion through continuum modeling. These first of their kind observations provide a technique to compare the performance of direct air capture solvents based on diffuse microenvironment dynamics while also providing metrics important for air contactor design such as boundary layer thickness. Overall, this work showcases a new experimental platform to study interfacial diffuse microenvironments in and outside of the field of direct air capture of carbon dioxide.

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

Pulsed Electrolysis Promotes Catalyst Activity in Dilute CO 2 Streams

Industrial CO 2 streams vary widely in composition, from pure to as low as 3%, posing challenges for purification or direct conversion. Electrochemical reduction offers a route for converting dilute CO 2 streams but faces severe mass transport limitations. This study demonstrates that pulsed electrolysis effectively overcomes these limitations, enhancing CO 2 electroreduction across variable feed compositions and current densities, particularly at low CO 2 concentrations and high current densities. At 25% CO 2 and 400 mA cm −2 , pulsing improved selectivity from 25.6 to 78.6%, production rate from 13.7 to 21.0 mol m −2 h −1 , and energy productivity from 0.77 to 2.59 mol kWh −1 . A dynamic, multiphysics continuum model confirms a 64% increase in CO 2 concentration within the catalyst layer during pulsing, resolving the transient chemical microenvironment. These findings establish pulsed electrolysis as a viable strategy for converting dilute industrial CO 2 streams into valuable feedstocks, bypassing costly pre-separation.

Orfali, Dania Muhieddine [New York University (NYU↗
DOE OSTI2025

Diffusion and Deformation Mechanism Maps in Li Metal from Atomistic Simulations

The rate of Li transport in the Li metal anode is important for the operation of Li metal-solid state batteries. Here, transport rates due to diffusion and creep are predicted in Li using atomistic simulations. First, molecular dynamics is used to estimate the rate of Li diffusion along dislocations and in grain boundary triple junctions. By combining this data with that from a prior study of grain boundary diffusion the dominant mechanisms and rates of self-diffusion in Li polycrystals are predicted as a function of grain size, grain shape, dislocation density, and temperature. Second, the dominant creep mechanisms are predicted and used to estimate critical current densities and void annihilation times. Grain boundary sliding and coble creep are the dominant mechanisms for micron-sized grains. Lastly, a continuum model for interfacial contact loss reveals that high dislocation densities of ∼10 12 /cm 2 enable achieving battery performance targets for Li grain sizes of ∼10 μm.

Batteries↗
DOE OSTI2026

Electrostatic Gating of Ionic Conductance through Heterogeneous van der Waals Nanopores

Nanofluidic ionic transistors typically require gate voltages above 1 V and operate only at submillimolar ionic strengths, limiting their biocompatible applications. We demonstrate ionic transistors consisting of single sub-10 nm nanopores drilled in van der Waals (vdW) heterostructures with internal gate electrodes made of few-layer graphene. These devices deliver up to 10-fold current modulation at gate voltages as low as 0.3 V in 10 mM KCl, and ∼2-fold modulation at near-physiological 100 mM KCl. Baseline conductance with no gate shows surface-charge-dominated transport below 100 mM KCl, consistent with negatively charged hBN walls and ∼5 nm opening of the pores. The surface charge and the electrochemical asymmetry introduced by the three-electrode configuration govern the device’s behavior: negative gate voltage (V G ) enriches ionic concentrations and enhances current, whereas positive V G induces a local depletion zone that suppresses transport. The current modulation by V G is dependent on the polarity of the transmembrane potential and leads to ion current rectification. Molecular dynamics simulations of a nanopore in a hBN–graphene–hBN stack reveal confinement and surface charge-dependent suppression of the relative permittivity of interfacial water. Continuum modeling with radially varying interfacial water permittivity reproduces the asymmetric I–V characteristics and explains how the embedded gate sculpts local potential and ion concentrations. By enabling sub-0.5 V control of ionic transport at up to 100 mM salt concentrations, these devices address a key need in nanofluidics to create low-power ionic circuits and biosensing.

Materials science↗
DOE OSTI2020

Gating of Hydrophobic Nanopores with Large Anions

In this report, understanding ion transport in nanoporous materials is critical to a wide variety of energy and environmental technologies, ranging from ion-selective membranes, drug delivery, and biosensing, to ion batteries and supercapacitors. While nanoscale transport is often described by continuum models that rely on a point charge description for ions and a homogeneous dielectric medium for the solvent, here, we show that transport of aqueous solutions at a hydrophobic interface can be highly dependent on the size and hydration strength of the solvated ions. Specifically, measurements of ion current through single silicon nitride nanopores that contain a hydrophobic-hydrophilic junction show that transport properties are dependent not only on applied voltage but also the type of anion. We find that in Cl - -contained solutions the nanopores only conducted ionic current above a negative voltage threshold. On the other hand, introduction of large polarizable anions, such as Br - and I - , facilitated the pore wetting, making the pore conductive at all examined voltages. Molecular dynamics simulations revealed that the large anions, Br - and I - , have a weaker solvation shell compared to that of Cl - , and consequently were prone to migrate from the aqueous solution to the hydrophobic surface, leading to the anion accumulation responsible for pore wetting. The results are essential for designing nanoporous systems that are selective to ions of the same charge, for realization of ion-induced wetting in hydrophobic pores, as well as for fundamental understanding on the role of ion hydration shell on the properties of solid/liquid interfaces.

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

Modification at the C2'-O-Position with 2-Methylbenzothiophene Induces Unique Structural Changes and Thermal Transitions on Duplexes of RNA and DNA

Oligonucleotides can be chemically modified for a variety of applications that include their use as biomaterials, in therapeutics, or as tools to understand biochemical processes, among others. This work focuses on the functionalization of oligonucleotides of RNA and DNA (12- or 14-nucleotides long) with methylbenzothiophene (BT), at the C2'-O-position, which led to unique structural features. Circular dichroism (CD) analyses showed that positioning the BT units on one strand led to significant thermal destabilization, while duplexes where each strand contained 4-BT rings formed a distinct arrangement with cooperativity/interactions among the modifications (evidenced from the appearance of a band with positive ellipticity at 235 nm). Interestingly, the structural arrays displayed increased duplex stabilization (>10 °C higher than the canonical analogue) as a function of [Na + ] with an unexpected structural rearrangement at temperatures above 50 °C. Density functional theory–polarizable continuum model (DFT-PCM) calculations were carried out, and the analyses were in agreement with induced structural changes as a function of salt content. A model was proposed where the hydrophobic surface allows for an internal nucleobase rearrangement into a more thermodynamically stable structure, before undergoing full denaturation, with increased heat. While this behavior is not common, B- to Z-form duplex transitions can occur and are dependent on parameters that were probed in this work, i.e., temperature, nature of modification, or ionic content. To take advantage of this phenomenon, we probed the ability of the modified duplexes to be recognized by Zα (an RNA binding protein that targets Z-form RNA) via electrophoretic analysis and CD. Interestingly, the protein did not bind to canonical duplexes of DNA or RNA; however, it recognized the modified duplexes, in a [monovalent/divalent salt] dependent manner. Overall, the findings describe methodology to attain unique structural motifs of modified duplexes of DNA or RNA, and control their behavior as a function of salt concentration. While their affinity to RNA binding proteins, and the corresponding mechanism of action, requires further exploration, the tunable properties can be of potential use to study this, and other, types of modifications. The novel arrays that formed, under the conditions described herein, provide a useful way to explore the structure and behavior of modified oligonucleotides, in general.

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

Spectroscopy of a tunable moiré system with a correlated and topological flat band

Abstract Moiré superlattices created by the twisted stacking of two-dimensional crystals can host electronic bands with flat energy dispersion in which enhanced interactions promote correlated electron states. The twisted double bilayer graphene (TDBG), where two Bernal bilayer graphene are stacked with a twist angle, is such a moiré system with tunable flat bands. Here, we use gate-tuned scanning tunneling spectroscopy to directly demonstrate the tunability of the band structure of TDBG with an electric field and to show spectroscopic signatures of electronic correlations and topology for its flat band. Our spectroscopic experiments are in agreement with a continuum model of TDBG band structure and reveal signatures of a correlated insulator gap at partial filling of its isolated flat band. The topological properties of this flat band are probed with the application of a magnetic field, which leads to valley polarization and the splitting of Chern bands with a large effective g-factor.

36 MATERIALS SCIENCE↗