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 307 records · Page 17

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Magic in twisted transition metal dichalcogenide bilayers

The long-wavelength moiré superlattices in twisted 2D structures have emerged as a highly tunable platform for strongly correlated electron physics. We study the moiré bands in twisted transition metal dichalcogenide homobilayers, focusing on WSe 2 , at small twist angles using a combination of first principles density functional theory, continuum modeling, and Hartree-Fock approximation. We reveal the rich physics at small twist angles θ < 4°, and identify a particular magic angle at which the top valence moiré band achieves almost perfect flatness. In the vicinity of this magic angle, we predict the realization of a generalized Kane-Mele model with a topological flat band, interaction-driven Haldane insulator, and Mott insulators at the filling of one hole per moiré unit cell. The combination of flat dispersion and uniformity of Berry curvature near the magic angle holds promise for realizing fractional quantum anomalous Hall effect at fractional filling. We also identify twist angles favorable for quantum spin Hall insulators and interaction-induced quantum anomalous Hall insulators at other integer fillings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ca-dimers, solvent layering, and dominant electrochemically active species in Ca(BH4)2 in THF

Abstract Divalent ions (Mg, Ca, and Zn) are being considered as competitive, safe, and earth-abundant alternatives to Li-ion electrochemistry, but present challenges for stable cycling due to undesirable interfacial phenomena. We explore the formation of electroactive species in the electrolyte Ca(BH 4 ) 2 ∣THF using molecular dynamics coupled with a continuum model of bulk and interfacial speciation. Free-energy analysis and unsupervised learning indicate a majority population of neutral Ca dimers and monomers with diverse molecular conformations and an order of magnitude lower concentration of the primary electroactive charged species – the monocation, $${\rm{CaBH}}_{4}^{+}$$ CaBH 4 + – produced via disproportionation of neutral complexes. Dense layering of THF molecules within ~1 nm of the electrode surface strongly modulates local electrolyte species populations. A dramatic increase in monocation population in this interfacial zone is induced at negative bias. We see no evidence for electrochemical activity of fully-solvated Ca 2+ . The consequences for performance are discussed in light of this molecular-scale insight.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Emergent actin flows explain distinct modes of gliding motility

During host infection, Toxoplasma gondii and related unicellular parasites move using gliding, which differs fundamentally from other known mechanisms of eukaryotic cell motility. Gliding is thought to be powered by a thin layer of flowing filamentous (F)-actin sandwiched between the plasma membrane and a myosin-covered inner membrane complex. How this surface actin layer drives the various gliding modes observed in experiments—helical, circular, twirling and patch, pendulum or rolling—is unclear. Here we suggest that F-actin flows arise through self-organization and develop a continuum model of emergent F-actin flow within the confines provided by Toxoplasma geometry. In the presence of F-actin turnover, our model predicts the emergence of a steady-state mode in which actin transport is largely directed rearward. Removing F-actin turnover leads to actin patches that recirculate up and down the cell, which we observe experimentally for drug-stabilized actin bundles in live Toxoplasma gondii parasites. These distinct self-organized actin states can account for observed gliding modes, illustrating how different forms of gliding motility can emerge as an intrinsic consequence of the self-organizing properties of F-actin flow in a confined geometry.

59 BASIC BIOLOGICAL SCIENCES↗

Experimental signature of layer skyrmions and implications for band topology in twisted WSe 2 bilayers

Twisted homobilayers of transition-metal dichalcogenides have been established as an ideal platform for studying strong correlation phenomena, as exemplified by the recent discovery of fractional Chern insulator states in twisted MoTe 2 , as well as Chern insulators and unconventional superconductivity in twisted WSe 2 . In these systems, a non-trivial topology in the strongly layer-hybridized regime can arise from a spatial patterning of interlayer tunnelling amplitudes and layer-dependent potentials that yields a lattice of layer skyrmions. Here, in this work, we report experimental signatures of skyrmion textures in the layer degree of freedom of rhombohedral-stacked twisted WSe 2 homobilayers. Using scanning tunnelling spectroscopy that separately resolves the Γ-valley and K-valley moiré electronic states, we reveal opposite layer polarizations of the K valley at two different lattice sites with opposite stacking order within the moiré unit cell. These findings are consistent with the theoretically predicted layer-skyrmion texture. We also use our experimental results to parameterize a common continuum model of moiré bands in twisted bilayers, further establishing a direct correlation between the shape of the local density of states in real space and the topology of the topmost moiré band.

electronic properties and materials↗

Flow of long chain hydrocarbons through carbon nanotubes (CNTs)

The pressure-driven flow of long-chain hydrocarbons in nanosized pores is important in energy, environmental, biological, and pharmaceutical applications. This paper examines the flow of hexane, heptane, and decane in carbon nanotubes (CNTs) of pore diameters 1–8 nm using molecular dynamic simulations. Enhancement of water flow in CNTs in comparison to rates predicted by continuum models has been well established in the literature. Our work was intended to observe if molecular dynamic simulations of hydrocarbon flow in CNTs produced similar enhancements. We used the OPLS-AA force field to simulate the hydrocarbons and the CNTs. Our simulations predicted the bulk densities of the hydrocarbons to be within 3% of the literature values. Molecular sizes and shapes of the hydrocarbon molecules compared to the pore size create interesting density patterns for smaller sized CNTs. We observed moderate flow enhancements for all the hydrocarbons (1–100) flowing through small-sized CNTs. For very small CNTs the larger hydrocarbons were forced to flow in a cork-screw fashion. As a result of this flow orientation, the larger molecules flowed as effectively (similar enhancements) as the smaller hydrocarbons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Emergence of Lanes and Turbulent-like Motion in Active Spinner Fluid

Assemblies of self-rotating particles are gaining interest as a novel realization of active matter with unique collective behaviors such as edge currents and non-trivial dynamic states. Here, we develop a continuum model for a system of fluid-embedded spinners by coarse-graining the equations of motion of the discrete particles. We apply the model to explore mixtures of clockwise and counterclockwise rotating spinners. We find that the dynamics is sensitive to fluid inertia: in the inertialess system, after transient turbulent-like motion the spinners segregate and form steady traffic lanes. At small but finite Reynolds number instead, the turbulent-like motion persists and the system exhibits a chirality breaking transition leading to a single rotation sense state. Our results shed light on the dynamic behavior of non-equilibrium materials exemplified by active spinners.

36 MATERIALS SCIENCE↗

Concentration and velocity profiles in a polymeric lithium-ion battery electrolyte

Predictive knowledge of ion transport in electrolytes which bridges microscopic and macroscopic length scales is imperative to design new ion conductors and to simulate device performance. Here, we employed a novel approach combining operando X-ray photon correlation spectroscopy, X-ray absorption microscopy, continuum modelling, and molecular dynamics simulations to probe the ion transport in a baseline polymeric lithium-ion battery electrolyte. In a Li/PEO–LiTFSI/Li symmetric cell under polarization, we determined and rationalized microscopic properties including local electrolyte velocities and ion correlations and connected this insight to measured and simulated macroscopic ion concentration gradients. By relating our results across length scales, we suggest a fairly concentration-independent transference number of about 0.2. Our study shows the broad applicability of operando X-ray photon correlation spectroscopy to the understanding of dynamic phenomena.

25 ENERGY STORAGE↗

An interplay between a hydrogen atmosphere and dislocation characteristics in BCC Fe from time-averaged molecular dynamics

The interplay between hydrogen and dislocations (e.g., core and elastic energies, and dislocation–dislocation interactions) has implications on hydrogen embrittlement but is poorly understood. Continuum models of hydrogen enhanced local plasticity have not considered the effect of hydrogen on dislocation core energies. Energy minimization atomistic simulations can only resolve dislocation core energies in hydrogen-free systems because hydrogen motion is omitted so hydrogen atmosphere formation can’t occur. Additionally, previous studies focused more on face-centered-cubic than body-centered-cubic metals. Discrete dislocation dynamics studies of hydrogen–dislocation interactions assume isotropic elasticity, but the validity of this assumption isn’t understood. Here, we perform time-averaged molecular dynamics simulations to study the effect of hydrogen on dislocation energies in body-centered-cubic iron for several dislocation character angles. We see atmosphere formation and highly converged dislocation energies. We find that hydrogen reduces dislocation core energies but can increase or decrease elastic energies of isolated dislocations and dislocation–dislocation interaction energies depending on character angle. We also find that isotropic elasticity can be well fitted to dislocation energies obtained from simulations if the isotropic elastic constants are not constrained to their anisotropic counterparts. These results are relevant to ongoing efforts in understanding hydrogen embrittlement and provide a foundation for future work in this field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Anion-exchange membranes with internal microchannels for water control in CO 2 electrolysis

Electrochemical reduction of carbon dioxide (CO 2 R) poses substantial promise to convert abundant feedstocks (water and CO 2 ) to value-added chemicals and fuels using solely renewable energy. However, recent membrane-electrode assembly (MEA) devices that have been demonstrated to achieve high rates of CO 2 R are limited by water management within the cell, due to both consumption of water by the CO 2 R reaction and electro-osmotic fluxes that transport water from the cathode to the anode. Additionally, crossover of potassium (K + ) ions poses concern at high current densities where saturation and precipitation of the salt ions can degrade cell performance. Herein, a device architecture incorporating an anion-exchange membrane (AEM) with internal water channels to mitigate MEA dehydration is proposed and demonstrated. A macroscale, two-dimensional continuum model is used to assess water fluxes and local water content within the modified MEA, as well as to determine the optimal channel geometry and composition. The modified AEMs are then fabricated and tested experimentally, demonstrating that the internal channels can both reduce K + cation crossover as well as improve AEM conductivity and therefore overall cell performance. This work demonstrates the promise of these materials, and operando water-management strategies in general, in handling some of the major hurdles in the development of MEA devices for CO 2 R.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cryptate binding energies towards high throughput chelator design: metadynamics ensembles with cluster–continuum solvation

A tiered forcefield/semiempirical/meta-GGA pipeline together with a thermodynamic scheme designed with error cancellation in mind was developed to calculate binding energies of [2.2.2] cryptate complexes of mono- and divalent cations. Stable complexes of Na, K, Rb, Ca, Zn and Pb were generated, revealing consistent cation–N lengths but highly variable cation–O lengths and an amine stacking mechanism potentially augmenting the cation size selectivity. Metadynamics, used for searching the high-dimensional potential energy surface, together with a cluster–continuum model for affordable – yet accurate – solvation modeling, enabled the discovery of more stable geometries than those previously reported. Similar solvation energy curve shapes for lone vs. coordinated ions enabled rapid solvation convergence via the cancellation of errors stemming from finite cluster sizes. In conclusion, an R 2 of 0.850 vs. experimental aqueous binding energies was obtained, validating this scheme as the backbone of a high-throughput workflow for chelator design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗