SEARCH · Search NASA

Results for “Linear stability theory”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Propane Dehydrogenation on Pt x Zn y Active Sites in Silicalite‐1

Abstract The improvement of Pt‐based catalysts for propane dehydrogenation (PDH) has progressed by recent investigations that have identified Zn as a promising promoter for Pt subnanometer catalysts. It is desirable to gain insights into the structure, stability, and activity of such active sites and the factors that influence them, such as Zn : Pt ratio, Pt coordination and nuclearity. Here, we employ density functional theory and microkinetic simulations to investigate the stability of Pt x Zn y ( x =1–3, y=0–3) active sites grafted on silanols of Silicalite‐1 and the PDH activity of Pt. We find that the coordination of a Pt atom to a nest of grafted Zn(II) atoms increases the stability of the Pt 1 Zn y sites, whose activity is similar for y=0–2 and drops dramatically for y>2. We further demonstrate, via linear scaling relations and microkinetic simulations, that the turnover frequency obeys a volcano law as a function of propylene binding strength. The Pt 2 Zn 1 and Pt 3 Zn 1 sites are stable and exhibit activity similar to Pt 1 Zn 2 , but only Pt 1 Zn 2 manifests reaction kinetics consistent with experimental data, strongly suggesting the active site composition in the synthesized catalyst samples. The methodology presented here suggests a general strategy for deducing active site information such as composition through simple kinetic experiments.

Liu, Yilang↗

DOE OSTI2024

Propane Dehydrogenation on Pt x Zn y Active Sites in Silicalite‐1

Liu, Yilang↗

DOE OSTI2024

Neutrino fast flavor oscillations with moments: Linear stability analysis and application to neutron star mergers

Providing an accurate modeling of neutrino physics in dense astrophysical environments such as binary neutron star mergers presents a challenge for hydrodynamic simulations. Nevertheless, understanding how flavor transformation can occur and affect the dynamics, the mass ejection, and the nucleosynthesis will need to be achieved in the future. Computationally expensive, large-scale simulations frequently evolve the first classical angular moments of the neutrino distributions. By promoting these quantities to matrices in flavor space, we develop a linear stability analysis of fast flavor oscillations using only the first two “quantum” moments, which notably requires generalizing the classical closure relations that appropriately truncate the hierarchy of moment equations in order to treat quantum flavor coherence. After showing the efficiency of this method on a well-understood test situation, we perform a systematic search of the occurrence of fast flavor instabilities in a neutron star merger simulation. Here, we discuss the successes and shortcomings of moment linear stability analysis, as this framework provides a time-efficient way to design and study better closure prescriptions in the future.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

A fluctuating hydrodynamics model for nanoscale surfactant-laden interfaces

A multispecies diffuse interface model is formulated in a fluctuating hydrodynamics framework for the purpose of simulating surfactant interfaces at the nanoscale. The model generalizes previous work to ternary mixtures, employing a Cahn-Hilliard free energy density combined with incompressible, isothermal fluctuating hydrodynamics where dissipative fluxes include both deterministic and stochastic terms. The intermolecular parameters in the free energy are chosen such that one species acts as a partially miscible surfactant. From Laplace pressure measurements, we show that in this model the surface tension decreases linearly with surfactant concentration, leading to Marangoni convection for interfaces with concentration gradients. In the capillary wave spectrum for interfaces with and without surfactant, we find that for the former, the spectrum deviates significantly from classical capillary wave theory, presumably due to Gibbs elasticity. In non-equilibrium simulations of the Rayleigh-Plateau instability, deterministic simulations showed that the surfactant delays pinching of a fluid cylinder into droplets. However, stochastic simulations indicate that thermal fluctuations disrupt the surfactant's stabilizing effect. Similarly, the spreading of a patch of surfactant, driven by Marangoni convection, was found to be partially suppressed by thermal fluctuations.

Capillary waves↗

DOE OSTI2021

High-Entropy Ceramic Coatings: Transformative New Materials for Environmentally-Compatible Thin-Film Insulators against High-T Molten Salts

An emerging new class of promising materials, High-Entropy Ceramics (HECs) or Compositionally-Complex Ceramics (CCCs), have been successfully fabricated in the last three years. This project evaluates HECs/CCCs in the concentrated solar power (CSP) context. In this project, we have identified a couple of simple, and key descriptors governing thermal conductivity in HECs. Size disorder among the cations was found to have a strong influence on the thermal conductivity across multiple systems and throughout a vast compositional space. Since identifying this, we have optimized previous compositions to achieve low thermal conductivities of $k$~$1-1.4Wm{-1}K{-1}$. We have further reduced our selection of promising compositions based on project guidelines and applied a theory that can potentially predict reactivity with molten salts to identify a small set of HECs that are suitable for further molten salt heat transfer media (HTM) stability testing. Many compositions reveal poor resistance to molten chloride salt after being exposed for 24 h compared to the reference structural alloy, Haynes 230. However, after 100 h, a few HECs/CCCs show strong parabolic rate constant effects and provide greater resistance to chloride salt compared to Haynes 230 alloy, which shows more linear kinetics. Postmortem characterization of the microstructure reveals the parabolic rate constant originates from a dense reaction layer form at the interface. This feature is absent in Haynes 230, which only shows grain boundary infiltration. The HECs/CCCs were also put in contact with the metal alloy under Ar for 800 h at 800°C and show negligible interaction. Lastly, parabolic rate constant HECs/CCCs also have low thermal conductivity $(<2 Wm^{-1} K^{-1})$ at 800°C. These specimens satisfy all project requirements and demonstrate their promising application in next-generation CSP.

14 SOLAR ENERGY↗

DOE OSTI2021

Artificial intelligence “sees” split electrons

Chemical bonds between atoms are stabilized by the exchange-correlation (xc) energy, a quantum-mechanical effect in which “social distancing” by electrons lowers their electrostatic repulsion energy. Kohn-Sham density functional theory (DFT) states that the electron density determines this xc energy, but the density functional must be approximated. Furthermore, this is usually done by satisfying exact constraints of the exact functional (making the approximation predictive), by fitting to data (making it interpolative), or both. Two exact constraints—the ensemble-based piecewise linear variation of the total energy with respect to fractional electron number and fractional electron z-component of spin —require hard-to-control nonlocality. On page 1385 of this issue, Kirkpatrick et al. have taken a big step toward more accurate predictions for chemistry through the machine learning of molecular data plus the fractional charge and spin constraints, expressed as data that a machine can learn.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Boundary Corrections for Kernel Approximation to Differential Operators

The kernel-based approach to operator approximation for partial differential equations has been shown to be unconditionally stable for linear PDEs and numerically exhibit unconditional stability for non-linear PDEs. These methods have the same computational cost as an explicit finite difference scheme but can exhibit order reduction at boundaries. In previous work on periodic domains, order reduction was addressed, yielding high-order accuracy. The issue addressed in this work is the elimination of order reduction of the kernel-based approach for a more general set of boundary conditions. Further, we consider the case of both first and second order operators. To demonstrate the theory, we provide not only the mathematical proofs but also experimental results by applying various boundary conditions to different types of equations. The results agree with the theory, demonstrating a systematic path to high order for kernel-based methods on bounded domains.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Fragmentation of faceted crystalline wires

The reduction of surface energy is a well-known driving force for the fragmentation of nanowires. Fragmentation can limit the service life of those small low-dimensional building blocks or can be controlled to induce beneficial shape changes. For isotropic surface energy, fragmentation is analogous to the classic Rayleigh-Plateau instability of liquid jets. However, commonly synthesized crystalline nanowires have strongly anisotropic surface energies and exhibit facets coinciding with cusps in the Wulff plot. Depending on growth orientation, different fragmentation behaviors have been seen in nanowires. Using phase-field simulations, we show that fragmentation of faceted nanowires with cubic crystal symmetry may occur by a finite-amplitude nonlinear instability, as opposed to a Rayleigh-Plateau-like linear instability, depending on nanowire growth orientation. Here, we carry out a weakly nonlinear analysis based on sharp interface theory to characterize the faceted nanowire shape corresponding to a surface-energy saddle point. The analysis predicts that the minimum amplitude of a periodic shape perturbation to trigger fragmentation increases with cusp strength but decreases inversely proportionally to the perturbation wavelength for long wavelengths, in good quantitative agreement with phase-field predictions. The results provide the theoretical foundation to predict nanowire stability as a function of length and surface energy in diverse applications.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Review—Development of Highly Active and Durable Hybrid Compressive Platinum Lattice Catalysts for Polymer Electrolyte Membrane Fuel Cells: Mathematical Modeling and Experimental Work

This review provides a comprehensive overview on the development of highly active and durable platinum catalysts with ultra-low Pt loadings for polymer electrolyte membrane fuel cells (PEMFCs) through a combined mathematical modeling and experimental work. First, simulation techniques were applied to evaluate the validity of the Tafel approximation for the calculation of the mass activity (MA) and specific activity (SA). A one-dimensional agglomeration model was developed and solved to understand the effects of exchange current density, porosity, agglomerate size, Nafion ® film thickness, and Pt loading on the MA and SA. High porosity (> 60%) and agglomerations at high Pt loadings cause the loss of the Tafel approximation and consequently the decrease in MA and SA. A new structure parameter was introduced to estimate the real porous structure using the fractal theory. The volumetric catalyst density was corrected by the fractal dimension (measured by Hg porosimetry), which gave a good agreement with the experimental values. The loading-dependent Tafel equation was then derived, which contains both the utilization and the non-linear scaling factor. Second, activated carbon composite support (ACCS) with optimized surface area, porosity, pore size, and pore size distribution was developed. The hydrophilic/hydrophobic ratio, structural properties (amorphous/crystalline ratio), and the number of active sites were optimized through metal-catalyzed pyrolysis. Stability of ACCS and Pt/ACCS were evaluated using an accelerated stress test (AST). The results indicated that Pt/ACCS showed no significant loss of MA and power density after 5,000 cycles at 1.0–1.5 V, while the commercial Pt/C catalysts showed drastic losses of MA and power density. Finally, monolayers of compressed Pt (core–shell-type Pt 3 Co 1 ) catalysts were structured by diffusing Co atoms (previously embedded in ACCS) into Pt. Compressive Pt lattice (Pt * ) catalysts were synthesized through an annealing procedure developed at the University of South Carolina (USC). The Pt * /ACCS catalyst showed high initial power density (rated) of 0.174 g Pt kW −1 and high stability (24 mV loss) at 0.8 A cm −2 after 30,000 cycles (0.6–1.0 V). The outstanding performance of Pt * /ACCS is due to the synergistic effect of ACCS and compressive Pt * lattice.

25 ENERGY STORAGE↗

DOE OSTI2025

Metastability of stratified magnetohydrostatic equilibria and their relaxation

Motivated by explosive releases of energy in fusion, space and astrophysical plasmas, we consider the nonlinear stability of stratified magnetohydrodynamic equilibria against two-dimensional interchanges of straight magnetic-flux tubes. We demonstrate that, even within this restricted class of dynamics, the linear stability of an equilibrium does not guarantee its nonlinear stability: equilibria can be metastable. We show that the minimum-energy state accessible to a metastable equilibrium under non-diffusive two-dimensional dynamics can be found by solving a combinatorial optimisation problem. These minimum-energy states are, to good approximation, the final states reached by our simulations of destabilised metastable equilibria for which turbulent mixing is suppressed by viscosity. To predict the result of fully turbulent relaxation, we construct a statistical mechanical theory based on the maximisation of Boltzmann's mixing entropy. This theory is analogous to the Lynden-Bell statistical mechanics of collisionless stellar systems and plasma, and to the Robert–Sommeria–Miller theory of two-dimensional vortex turbulence. Our theory reproduces well the results of our numerical simulations for sufficiently large perturbations to the metastable equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Hybrid simulations of sub-cyclotron compressional and global Alfven Eigenmode stability in spherical tokamaks

A comprehensive numerical study has been conducted in order to investigate the stability of beam-driven, sub-cyclotron frequency compressional (CAE) and global (GAE) Alfven Eigenmodes in low aspect ratio plasmas for a wide range of beam parameters. The presence of CAEs and GAEs has previously been linked to anomalous electron temperature profile flattening at high beam power in NSTX experiments, prompting further examination of the conditions for their excitation. Linear simulations are performed with the hybrid MHD-kinetic initial value code HYM in order to capture the general Doppler-shifted cyclotron resonance that drives the modes. Three distinct types of modes are found in simulations -- co-CAEs, cntr-GAEs, and co-GAEs -- with differing spectral and stability properties. The simulations reveal that unstable GAEs are more ubiquitous than unstable CAEs, consistent with experimental observations, as they are excited at lower beam energies and generally have larger growth rates. Local analytic theory is used to explain key features of the simulation results, including the preferential excitation of different modes based on beam injection geometry and the growth rate dependence on the beam injection velocity, critical velocity, and degree of velocity space anisotropy. The background damping rate is inferred from simulations and estimated analytically for relevant sources not present in the simulation model, indicating that co-CAEs are closer to marginal stability than modes driven by the cyclotron resonances.

Compressional Alfven Eigenmodes↗

DOE OSTI2021

Effect of polycyclosilane microstructure on thermal properties

Thermal characterization of polysilanes has focused on the influence of organic side chains, whereas little is understood about the influence of silane backbone microstructure on thermal stability, phase properties, and pyrolysis. To address this knowledge gap, we prepared three distinct polycyclosilanes: linear polymers synthesized from the cyclosilane building blocks 1,4Si 6 and 1,3Si 6 , as well as a cyclic polymer of 1,3Si 6 . Thermal properties across the temperature range 25 to 600 °C were investigated using differential scanning calromietry (DSC) and thermogravimetric analysis (TGA). We found differences between linear and cyclic materials, including a phase transition unique to the cyclic polymer and lower rates of mass loss during pyrolysis. Density functional theory (DFT) calculations provided insight into microstructure-dependent pyrolysis.

36 MATERIALS SCIENCE↗

DOE OSTI2021

The ultrahigh pressure stability of silver: An experimental and theoretical study

We measured the atomic volume of Ag in a toroidal diamond anvil cell to a maximum pressure of 416 GPa and calculated the atomic volume and elastic constants of Ag up to 750 and 460 GPa, respectively. Our density functional theory calculations at 0 K utilize an all-electron fully relativistic method and agree well with our volume measurements, particularly at pressures above ~75 GPa. In this work, we corrected our experimental results for non-hydrostaticity using a line shift analysis, and the resulting Vinet equation of state (EOS) parameters are reported. We find that the uniaxial stress sustained by Ag increases linearly up to 4.5 GPa at a pressure of 416 GPa. Our experimental results indicate that the fcc structure of Ag remains stable to at least 416 GPa at room temperature. Our theoretical results show that C44 increases as pressure increases and reaches a maximum at ~100 GPa above which it begins to decrease, a sign that the fcc structure of Ag is becoming unstable, and at V/V 0 = 0.30, the bcc structure is lower in energy than fcc.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Hybrid simulations of sub-cyclotron compressional and global Alfvén eigenmode stability in spherical tokamaks

A comprehensive numerical study has been conducted in order to investigate the stability of beam-driven, sub-cyclotron-frequency compressional Alfvén eigenmodes (CAEs) and global Alfvén eigenmodes (GAEs) in low-aspect-ratio plasmas for a wide range of beam parameters. The presence of CAEs and GAEs has previously been linked to anomalous electron temperature profile flattening at high beam powers in NSTX experiments, prompting a further examination of the conditions necessary for their excitation. Linear simulations have been performed with the hybrid MHD–kinetic initial value code HYM in order to capture the general Doppler-shifted cyclotron resonance that drives the modes. Additionally, three distinct types of modes were found in the simulations—co-CAEs, cntr-GAEs, and co-GAEs—with differing spectral and stability properties. The simulations revealed that unstable GAEs are more ubiquitous than unstable CAEs, which is consistent with experimental observations, as they are excited at lower beam energies and generally have larger growth rates. Local analytic theory is used to explain key features of the simulation results, including the preferential excitation of different modes based on beam injection geometry and the growth rate dependence on the beam injection velocity, critical velocity, and degree of velocity space anisotropy. The background damping rate is inferred from the simulations and analytically estimated for relevant sources absent from the simulation model, indicating that co-CAEs are closer to marginal stability than modes driven by the cyclotron resonances.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Electron– K -phonon interaction in twisted bilayer graphene

We develop an analytic theory to describe the interaction between electrons and K phonons and study its influence on superconductivity in the bare bands of twisted bilayer graphene (TBG). We find that, due to symmetry and the two-center approximation, only one optical K phonon (160meV) of graphene is responsible for the intervalley electron-phonon interaction. Here, this phonon has recently been found in angular-resolved photoemission spectroscopy to be responsible for replicas of the TBG flat bands. By projecting the interaction to the TBG flat bands, we perform the full symmetry analysis of the phonon-mediated attractive interaction and pairing channels in the Chern basis, and show that several channels are guaranteed to have gapless order parameters. From the linearized gap equations, we find that the highest T c pairing induced by this phonon is a singlet gapped s-wave inter-Chern-band order parameter, followed closely by a gapless nematic d-wave intra-Chern-band order parameter. We justify these results analytically, using the topological heavy-fermion mapping of TBG which has allowed us to obtain an analytic form of a phonon-mediated attractive interaction and to analytically solve the linearized and T = 0 gap equations. For the intra-Chern-band channel, the nematic state with nodes is shown to be stabilized in the chiral flat-band limit. While the flat-band Coulomb interaction can be screened sufficiently enough—around the Van Hove singularities—to allow for electron-phonon based superconductivity, it is unlikely that this effect can be maintained in the lower density of states excitation bands around the correlated insulator states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Heterotic instantons for monad and extension bundles

We consider non-perturbative superpotentials from world-sheet instantons wrapped on holomorphic genus zero curves in heterotic string theory. These superpotential contributions feature prominently in moduli stabilization and large field axion inflation, which makes their presence or absence, as well as their functional dependence on moduli, an important issue. We develop geometric methods to compute the instanton superpotentials for heterotic string theory with monad and extension bundles. Using our methods, we find a variety of examples with a non-vanishing superpotential. In view of standard vanishing theorems, we speculate that these results are likely to be attributed to the non-compactness of the instanton moduli space. We test this proposal, for the case of monad bundles, by considering gauged linear sigma models where compactness of the instanton moduli space can be explicitly checked. In all such cases, we find that the geometric results are consistent with the vanishing theorems. Surprisingly, linearly dependent Pfaffians even arise for cases with a non-compact instanton moduli space. This suggests some gauged linear sigma models with a non-compact instanton moduli space may still have a vanishing instanton superpotential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Importance of resistivity on edge-localized mode onset in spherical tokamaks

In this work, we explore the impact of realistic plasma resistivity on the linear stability of peeling-ballooning (PB) modes in tokamak discharges with low-aspect ratio. For this study we consider discharges that are subject to edge-localized modes (ELMs) in the National Spherical Torus Experiment (NSTX). Employing the state of the art extended-magnetohydrodynamic (MHD) code M3D-C1 it is demonstrated that non-ideal effects can significantly affect PB stability thresholds in NSTX discharges. In particular, robust resistive PB modes are found to exist well before the ideal PB stability threshold is met. These novel results can explain why ideal-MHD theory often does not accurately describe ELM onset in spherical torus configurations, and also present a valuable basis for the development of a predictive model for ELMs in low-aspect ratio tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Configurational stability of a crack propagating in a material with mode-dependent fracture energy – Part II: Drift of fracture facets in mixed-mode I+II+III

In earlier papers (Leblond et al., 2011; 2019), we presented linear stability analyses of the coplanar propagation of a crack loaded in mixed-mode I+III, based on a “double” propagation criterion combining (Griffith, 1920)’s energetic condition and (Goldstein and Salganik, 1974)’s principle of local symmetry. The difference between the two papers was that in the more recent one, the local value of the critical energy-release-rate was no longer considered as a constant, but heuristically allowed to depend upon the ratio of the local mode III to mode I stress intensity factors. This led to a much improved, qualitatively acceptable agreement of theory and experiments, for the “threshold” value of the ratio of the unperturbed mode III to mode I stress intensity factors, above which coplanar propagation becomes unstable. In this paper, the analysis is extended to the case where a small additional mode II loading component is present in the initially planar configuration of the crack, generating a small, general kink of this crack from the moment it is applied. The main new effect resulting from presence of such a loading component is that the instability modes present above the threshold must drift along the crack front during its propagation. This prediction may be useful for future theoretical interpretations of a number of experiments where such a drifting motion was indeed observed.

36 MATERIALS SCIENCE↗