Search NASA⌕ Search

SEARCH · Search NASA

Results for “stability analysis”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Spinodal electronic phase separation during insulator-metal transitions

Electronic phase transitions such as insulator-metal transitions are common in strongly correlated systems. Here, using a combination of thermodynamic linear-stability analysis and phase-field simulations and employing VO 2 as a prototypical example, we predict that an insulator-metal transition driven by photoexcitation may involve an intermediate, modulated charge density state with a temperature-dependent characteristic wavelength. It is shown that such an intermediate two-phase electronic state is formed through a spinodal mechanism and that its formation can be generic for insulator-metal transitions driven by fast stimuli. This transient electronic phase separation is expected to stimulate future experimental and computational efforts.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Elastically mediated interactions between grain boundaries and precipitates in two-phase coherent solids

In this work, we investigate analytically and numerically the interaction between grain boundaries and second phase precipitates in two-phase coherent solids in the presence of misfit strain. Our numerical study uses amplitude equations that describe the interaction of composition and stress [R. Spatschek and A. Karma, Phys. Rev. B 81, 214201 (2010)] and free-energies corresponding to two-dimensional hexagonal and three-dimensional BCC crystal structures that exhibit isotropic and anisotropic elastic properties, respectively. We consider two experimentally motivated geometries where (i) a lamellar precipitate nucleates along a planar grain boundary that is centered inside the precipitate, and (ii) a circular precipitate nucleates inside a grain at a finite distance to an initially planar grain boundary. For the first geometry, we find that the grain boundary becomes morphologically unstable due to the combination of long-range elastic interaction between the grain boundary and compositional domain boundaries, and shear-coupled grain boundary motion. We characterize this instability analytically by extending the linear stability analysis carried out recently [P.-A. Geslin, Y.-C. Xu, and A. Karma, Phys. Rev. Lett. 114, 105501 (2015)] to the more general case of elastic anisotropy. The analysis predicts that elastic anisotropy hinders but does not suppress the instability. Simulations also reveal that, in a well-developed nonlinear regime, this instability can lead to the break-up of low-angle grain boundaries when the misfit strain exceeds a threshold that depends on the grain boundary misorientation. For the second geometry, simulations show that the elastic interaction between an initially planar grain boundary and an adjacent circular precipitate causes the precipitate to migrate to and anchor at the grain boundary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Comparison of momentum transport models for numerical relativity

The main problems of nonvacuum numerical relativity, compact binary mergers and stellar collapse, involve hydromagnetic instabilities and turbulent flows, so that kinetic energy at small scales leads to mean effects at large scale that drive the secular evolution. Notable among these effects is momentum transport. We investigate two models of this transport effect, a relativistic Navier-Stokes system and a turbulent mean stress model, that are similar to all of the prescriptions that have been attempted to date for treating subgrid effects on binary neutron star mergers and their aftermath. Our investigation involves both stability analysis and numerical experimentation on star and disk systems. We also begin the investigation of the effects of particle and heat transport on postmerger simulations. Here, we find that correct handling of turbulent heating is crucial for avoiding unphysical instabilities. Given such appropriate handling, the evolution of a differentially rotating star and the accretion rate of a disk are reassuringly insensitive to the choice of prescription. However, disk outflows can be sensitive to the choice of method, even for the same effective viscous strength. We also consider the effects of eddy diffusion in the evolution of an accretion disk and show that it can interestingly affect the composition of outflows.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Three-flavor collective neutrino conversions with multi-azimuthal-angle instability in an electron-capture supernova model

We investigate the multi-azimuthal angle (MAA) effect on collective neutrino oscillation by considering the three-dimensional neutrino momentum distribution in a realistic electron-capture supernova model with an 8.8 M ⊙ progenitor. We find that the MAA effect induces collective flavor conversions at epochs when it is completely suppressed under the axial-symmetric approximation. This novel activity is switched on/off by the growth of the MAA instability and imprints additional time evolution in the expected neutrino event rate. We validate our results by extending the linear stability analysis into the three-flavor scheme including mixing angles, and confirm that the onset of collective neutrino oscillation matches the steep growth of flavor instability. Furthermore, we discuss how the MAA effect alters neutrino detection at Super-Kamiokande and DUNE.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Entanglement and correlations in fast collective neutrino flavor oscillations

Collective neutrino oscillations play a crucial role in transporting lepton flavor in astrophysical settings like supernovae and neutron star binary merger remnants, which are characterized by large neutrino densities. In these settings, simulations in the mean-field approximation show that neutrino-neutrino interactions can overtake vacuum oscillations and give rise to fast collective flavor evolution on time-scales t ∝ μ –1 , with μ proportional to the local neutrino density. In this work, we study the full out-of-equilibrium flavor dynamics in simple multi-angle geometries displaying fast oscillations in the mean field linear stability analysis. Focusing on simple initial conditions, we analyze the production of pair correlations and entanglement in the complete many-body-dynamics as a function of the number N of neutrinos in the system, for up to thousands of neutrinos. Similarly to simpler geometries with only two neutrino beams, we identify three regimes: stable configurations with vanishing flavor oscillations, marginally unstable configurations with evolution occurring on long time scales τ ≈ µ –1 $\sqrt{N}$, and unstable configurations showing flavor evolution on short time scales τ ≈ μ –1 log(N). Here, we present evidence that these fast collective modes are generated by the same dynamical phase transition which leads to the slow bipolar oscillations, establishing a connection between these two phenomena and explaining the difference in their time scales. We conclude by discussing a semi-classical approximation which reproduces the entanglement entropy at short to medium time scales and can be potentially useful in situations with more complicated geometries where classical simulation methods starts to become inefficient.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Resolution requirements for numerical modeling of neutrino quantum kinetics

Neutrino quantum kinetics is a rapidly evolving field in computational astrophysics, with a primary focus on collective neutrino oscillations in core-collapse supernovae and postmerger phases of binary neutron star mergers. In recent years, there has been considerable debate concerning resolution dependence in numerical simulations. In this paper, we conduct a comprehensive resolution study in both angular- and spatial directions by using two independent schemes of quantum kinetic neutrino transport: finite volume and pseudospectral methods. We complement our discussion by linear stability analysis including inhomogeneous modes. Our result suggests that decreasing spatial resolutions underestimates the growth of flavor instability, and then leads to wrong asymptotic states of flavor conversions, which potentially has a critical impact on astrophysical consequences. We further delve into numerical results of low resolution simulations, that reveals the underlying mechanism responsible for numerical artifacts caused by insufficient resolutions. Finally, this study settles the debate on requirements of resolutions and serves as a guideline for numerical modeling of quantum kinetic neutrino transport.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Quantum maximum entropy closure for small flavor coherence

Quantum angular moment transport schemes are an important avenue toward describing neutrino flavor mixing phenomena in dense astrophysical environments such as supernovae and merging neutron stars. Successful implementation will require new closure relations that go beyond those used in classical transport. In this paper, we derive the first analytic expression for a quantum M1 closure, valid in the limit of small flavor coherence, based on the maximum entropy principle. We verify that the resulting closure relation has the appropriate limits and characteristic speeds in the diffusive and free-streaming regimes. We then use this new closure in a moment linear stability analysis to search for fast flavor instabilities in a binary neutron star merger simulation and find better results as compared with previously designed, ad hoc , semiclassical closures.

astrophysical & cosmological simulations↗
DOE OSTI2017

Failures of homogeneous and isotropic cosmologies in extended quasidilaton massive gravity

We analyze the extended quasidilaton massive gravity model around a Friedmann-Lemaître-Robertson-Walker cosmological background. We present a careful stability analysis of asymptotic fixed points. We find that the traditional fixed point cannot be approached dynamically, except from a perfectly fine-tuned initial condition involving both the quasidilaton and the Hubble parameter. A less-well examined fixed-point solution, where the time derivative of the zeroth Stückelberg field vanishes ${\stackrel{{}}{{\phi}}}^{0}=0$, encounters no such difficulty, and the fixed point is an attractor in some finite region of initial conditions. We examine the question of the presence of a Boulware-Deser ghost in the theory. We show that the additional constraint that generically allows for the elimination of the Boulware-Deser mode is only present under special initial conditions. We find that the only possibility corresponds to the traditional fixed point and the initial conditions are the same fine-tuned conditions that allow the fixed point to be approached dynamically.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Eulerian formulation of the interacting particle representation model of homogeneous turbulence

The Interacting Particle Representation Model (IPRM) of homogeneous turbulence incorporates information about the morphology of turbulent structures within the con nes of a one-point model. In the original formulation [Kassinos & Reynolds, Center for Turbulence Research: Annual Research Briefs, 31{51, (1996)], the IPRM was developed in a Lagrangian setting by evolving second moments of velocity conditional on a given gradient vector. In the present work, the IPRM is re-formulated in an Eulerian framework and evolution equations are developed for the marginal PDFs. Eulerian methods avoid the issues associated with statistical estimators used by Lagrangian approaches, such as slow convergence. A specific emphasis of this work is to use the IPRM to examine the long time evolution of homogeneous turbulence. We first describe the derivation of the marginal PDF in spherical coordinates, which reduces the number of independent variables and the cost associated with Eulerian simulations of PDF models. Next, a numerical method based on radial basis functions over a spherical domain is adapted to the IPRM. Finally, results obtained with the new Eulerian solution method are thoroughly analyzed. The sensitivity of the Eulerian simulations to parameters of the numerical scheme, such as the size of the time step and the shape parameter of the radial basis functions, is examined. A comparison between Eulerian and Lagrangian simulations is performed to discern the capabilities of each of the methods. Finally, a linear stability analysis based on the eigenvalues of the discrete differential operators is carried out for both the new Eulerian solution method and the original Lagrangian approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Morphological Instability of Grain Boundaries in Two-Phase Coherent Solids

We show both computationally and analytically that grain boundaries that exhibit shear-coupled motion become morphologically unstable in solid alloys that phase separate into coherent domains of distinct chemical compositions. Here, we carry out simulations of continuum models demonstrating that this instability is mediated by long-range elastic interaction between compositional domains and grain boundaries. In addition, we perform a linear stability analysis that predicts the range of unstable wavelengths in good quantitative agreement with simulations. In nonlinear stages, this pattern-forming instability leads to the breakup of low-angle grain boundaries, thereby strongly impacting microstructural evolution in a wide range of phase-separating materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Controlling Viscous Fingering Using Time-Dependent Strategies

Control and stabilization of viscous fingering of immiscible fluids impacts a wide variety of pressure-driven multiphase flows. Here, we report theoretical and experimental results on time-dependent control strategy by manipulating the gap thickness b(t) in a lifting Hele-Shaw cell in the power-law form b(t) = b 1 t 1/7 . Experimental results show good quantitative agreement with the predictions of linear stability analysis. Furthermore, by choosing the value of a single time-independent control parameter we can either totally suppress the viscous fingering instability or maintain a series of non-splitting viscous fingers during the fluid displacement process. Besides the gap thickness of a Hele-Shaw cell, in principle, time-dependent control strategies can also be placed on the injection rate, viscosity of the displaced fluid, and interfacial tensions between the two fluids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Crack Front Segmentation and Facet Coarsening in Mixed-Mode Fracture

A planar crack generically segments into an array of “daughter cracks” shaped as tilted facets when loaded with both a tensile stress normal to the crack plane (mode I) and a shear stress parallel to the crack front (mode III). We investigate facet propagation and coarsening using in situ microscopy observations of fracture surfaces at different stages of quasistatic mixed-mode crack propagation and phase-field simulations. The results demonstrate that the bifurcation from propagating a planar to segmented crack front is strongly subcritical, reconciling previous theoretical predictions of linear stability analysis with experimental observations. As a result, they further show that facet coarsening is a self-similar process driven by a spatial period-doubling instability of facet arrays.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Invisible Anchors Trap Particles in Branching Junctions

Here, we combine numerical simulations and an analytic approach to show that the capture of finite, inertial particles during flow in branching junctions is due to invisible, anchor-shaped three-dimensional flow structures. These Reynolds-number-dependent anchors define trapping regions that confine particles to the junction. For a wide range of Stokes numbers, these structures occupy a large part of the flow domain. For flow in a V-shaped junction, at a critical Stokes number, we observe a topological transition due to the merger of two anchors into one. Lastly, from a stability analysis, we identify the parameter region of particle sizes and densities where capture due to anchors occurs.

42 ENGINEERING↗
DOE OSTI2022

Intrinsic ferromagnetism and restrictive thermodynamic stability in MA$_{2}$N$_{4}$ and Janus VSiGeN$_{4}$ monolayers

The seminal experimental discovery of the remarkably stable MoSi 2 N 4 monolayer has led to a handful of predicted magnetic two-dimensional (2D) materials in the MA 2 Z 4 family (M = transition metals, A = Si, Ge, and Z = N, P, As). These magnetic monolayers were predicted to be dynamically stable, but none of them has been synthesized to date. In this Research Letter, from first-principles thermodynamic stability analysis, we demonstrate that only the nitrides are thermodynamically stable and this occurs under N-rich conditions. Based on this finding, we propose two ferromagnetic, semiconducting Janus monolayers in the family: VSiGeN 4 and VSiSnN 4 . They are both dynamically and thermally stable, but only the former is thermodynamically stable. Intriguingly, Janus VSiGeN 4 and VSiSnN 4 monolayers show weak in-plane anisotropy compared with the VSi 2 N 4 monolayer. Furthermore, these two emerging Janus magnetic semiconductors offer opportunities for studying 2D magnetism and spin control for spintronics applications.

2-dimensional systems↗
DOE OSTI2023

Two-dimensional Dirac semimetal based on the alkaline earth metal CaP 3

Using an evolutionary algorithm in combination with first-principles density-functional theory calculations, we identify a two-dimensional (2D) CaP 3 monolayer as a new Dirac semimetal due to inversion and nonsymmorphic spatial symmetries of the structure. This new topological material, composed of light elements, exhibits high structural stability (higher than the phase known in the literature), which is confirmed by thermodynamic and kinetic stability analysis. Moreover, it satisfies the electron filling criteria, so that its Dirac state is located near the Fermi level. The existence of the Dirac state predicted by the theoretical symmetry analysis is also confirmed by first-principles electronic band structure calculations. We find that the energy position of the Dirac state can be tuned by strain, while the Dirac state is unstable against an external electric field since it breaks the spatial inversion symmetry. In conclusion, our findings should be instrumental in the development of 2D Dirac fermions based on light elements for their application in nanoelectronic devices and topological electronics.

2-dimensional systems↗
DOE OSTI2025

Competing mechanisms at vibrated interfaces of density-contrast fluids

Fluid-fluid interfacial instability and subsequent fluid mixing are ubiquitous in nature and engineering. The hydrodynamic instability of fluid interfaces has long centered on the pressure gradient-driven long-wavelength Rayleigh-Taylor instability and the resonance-induced short-wavelength Faraday instability. However, neither instability alone can explain the dynamics when both mechanisms are present. We identify a previously unseen multi-modal instability emerging from their coexistence. When the denser fluid is polydimethylsiloxane, the mixed region at a high density contrast (Atwood number = 0.9) spans a vibration amplitude range approximately twice the gravitational acceleration. Using Floquet stability analysis, we show how vibrations govern transitions between the RT and Faraday instabilities, leading to contention between these instabilities rather than resonant enhancement. Here, the initial transient growth is represented by the exponential modal growth of the most unstable Floquet exponent, along with its accompanying periodic behavior. Direct numerical simulations validate these findings and track interface breakup into the multiscale and nonlinear regimes. Specifically, we show that growing RT modes nonlinearly suppresses Faraday responses even when the initial growth rate of the Faraday instability is 3.63 times that of RT, so a bidirectional competition hinders their sustained coexistence.

Direct numerical simulations↗
DOE OSTI2025

Rotated-Droop Control for Enhanced Stability and Power Decoupling in Microgrids With Complex Line Impedances

Classical droop control in microgrids predominantly assumes inductive line impedance, which simplifies implementation but causes power coupling and steady-state errors in systems with resistive and inductive lines. This paper proposes a rotated-droop control strategy for grid-forming inverters that reformulates the power equations by incorporating the magnitude and angle of the line impedance within a rotated reference frame. This method enhances the decoupling of active and reactive power dynamics without increasing complexity or requiring communication links. A small-signal state-space model was created to capture the dynamic behavior of the system under varying impedance parameters, preserving the original droop gains by rotating the power control structure. This allows eigenvalue-based stability analysis and enhances damping and transient performance. Simulation and experimental results validated the improved power-sharing performance, faster response, and robustness of the proposed method under different impedance conditions. This approach maintains the decentralized structure of the conventional droop control while enabling greater adaptability and scalability, making it suitable for modern inverter-based microgrids with dynamic topologies.

Campo-Ossa, Daniel Dario [Univ. of Puerto Rico, Ag↗