SEARCH · Search NASA

Results for “phase stability 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.

SourceResearch connectionFrom yearThrough year

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

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Nonlinear stability and control study of highly maneuverable high performance aircraft

This project is intended to research and develop new nonlinear methodologies for the control and stability analysis of high-performance, high angle-of-attack aircraft such as HARV (F18). Past research (reported in our Phase 1, 2, and 3 progress reports) is summarized and more details of final Phase 3 research is provided. While research emphasis is on nonlinear control, other tasks such as associated model development, system identification, stability analysis, and simulation are performed in some detail as well. An overview of various models that were investigated for different purposes such as an approximate model reference for control adaptation, as well as another model for accurate rigid-body longitudinal motion is provided. Only a very cursory analysis was made relative to type 8 (flexible body dynamics). Standard nonlinear longitudinal airframe dynamics (type 7) with the available modified F18 stability derivatives, thrust vectoring, actuator dynamics, and control constraints are utilized for simulated flight evaluation of derived controller performance in all cases studied.

Mohler, R. R.↗

NASA NTRS2004

Two-Fluid Models and Interfacial Area Transport in Microgravity Condition

The objective of the present study is to develop a two-fluid model formulation with interfacial area transport equation applicable for microgravity conditions. The new model is expected to make a leapfrog improvement by furnishing the constitutive relations for the interfacial interaction terms with the interfacial area transport equation, which can dynamically model the changes of the interfacial structures. In the first year of this three-year project supported by the U.S. NASA, Office of Biological and Physics Research, the primary focus is to design and construct a ground-based, microgravity two-phase flow simulation facility, in which two immiscible fluids with close density will be used. In predicting the two-phase flow behaviors in any two-phase flow system, the interfacial transfer terms are among the most essential factors in the modeling. These interfacial transfer terms in a two-fluid model specify the rate of phase change, momentum exchange, and energy transfer at the interface between the two phases. For the two-phase flow under the microgravity condition, the stability of the fluid particle interface and the interfacial structures are quite different from those under normal gravity condition. The flow structure may not reach an equilibrium condition and the two fluids may be loosely coupled such that the inertia terms of each fluid should be considered separately by use of the two-fluid model. Previous studies indicated that, unless phase-interaction terms are accurately modeled in the two-fluid model, the complex modeling does not necessarily warrant an accurate solution.

Ishii, Mamoru↗

NASA NTRS1986

Effect of surface ionization on wetting layers

A surface ionization model due to Langmuir is generalized to liquid mixtures of polar and nonpolar components in contact with ionizable substrates. When a predominantly nonpolar mixture is near a miscibility gap, thick wetting layers of the conjugate polar phase form on the substrate. Such charged layers can be much thicker than similar wetting layers stabilized by dispersion forces. This model may explain the 0.4- to 0.6-micron-thick wetting layers formed in stirred mixtures of nitromethane and carbon disulfide in contact with glass.

Kayser, R. F.↗

DOE OSTI2020

Fluorescence Enhancement in the Solid State by Isolating Perylene Fluorophores in Metal–Organic Frameworks

Polycyclic aromatic hydrocarbons such as perylene, pyrene, and their derivatives are highly emissive fluorophores in solution. Furthermore, the practical applications of these materials in the field of molecular electronic and light emitting devices are often hindered by self-quenching effects due to the formation of non-fluorescent aggregates in concentrated solutions or in the solid state. Herein, we demonstrate that aggregation caused quenching of perylenes can be minimalized by molecular incorporation into metalorganic frameworks (MOFs). This study utilized a stable Zr 6 cluster-based MOF, UiO-67, as a matrix. Linear linkers containing the photo-responsive moieties were designed and incorporated into the parent UiO-67 scaffold through the partially replacement of the non-fluorescent linkers of a similar length, forming mixed-linker MOFs. Here, the average distance between perylene moieties was tuned by changing the linker ratios, thus controlling the fluorescence intensity, emission wavelength, and quantum yield. Molecular modeling was further adopted to correlate the number of isolated perylene linkers within the framework to the ratio between the two linkers, therefore rationalizing the change in the observed fluorescent properties. Taking advantage of the tunable fluorescence, inherent porosity, and high chemical stability of this MOF platform, it was applied as a fluorescent sensor for oxygen detection in the gas phase, a model reaction, showing fast response and good recyclability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Variational study of the Kitaev-Heisenberg-Gamma model

Here we compute the low-energy excitation spectrum and the dynamical spin structure factor of the Kitaev-Heisenberg-Gamma model through a variational approach based on the exact fractionalized excitations of the pure Kitaev honeycomb model. This novel approach reveals the physical reason for the asymmetric stability of the Kitaev spin liquid phases around the ferromagnetic and antiferromagnetic Kitaev limits. Moreover, we demonstrate that the fractionalized excitations form bound states in specific regions of each Kitaev spin liquid phase and that certain phase transitions induced by Heisenberg and Gamma interactions are driven by the condensation of such a bound state. Remarkably, this bound state appears as a sharp mode in the dynamical spin structure factor, while its condensation patterns at the appropriate phase transitions provide a simple explanation for the magnetically ordered phases surrounding each Kitaev spin liquid phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Fictitious phase separation in Li layered oxides driven by electro-autocatalysis

Layered oxides widely used as lithium-ion battery electrodes are designed to be cycled under conditions that avoid phase transitions. Although the desired single-phase composition ranges are well established near equilibrium, operando diffraction studies on many-particle porous electrodes have suggested phase separation during delithiation. Notably, the separation is not always observed, and never during lithiation. These anomalies have been attributed to irreversible processes during the first delithiation or reversible concentration-dependent diffusion. However, these explanations are not consistent with all experimental observations such as rate and path dependencies and particle-by-particle lithium concentration changes. Here, we show that the apparent phase separation is a dynamical artefact occurring in a many-particle system driven by autocatalytic electrochemical reactions, that is, an interfacial exchange current that increases with the extent of delithiation. We experimentally validate this population-dynamics model using the single-phase material L ix (Ni 1/3 Mn 1/3 Co 1/3 )O 2 (0.5 < x < 1) and demonstrate generality with other transition-metal compositions. Operando diffraction and nanoscale oxidation-state mapping unambiguously prove that this fictitious phase separation is a repeatable non-equilibrium effect. We quantitatively confirm the theory with multiple-datastream-driven model extraction. More generally, our study experimentally demonstrates the control of ensemble stability by electro-autocatalysis, highlighting the importance of population dynamics in battery electrodes (even non-phase-separating ones).

36 MATERIALS SCIENCE↗

DOE OSTI2023

Role of Fe in long-range ordered Ni 2 Cr precipitates in Ni-Cr-Fe model alloys during isothermal aging

The precipitation of new phases during long-term service at elevated temperatures is a concern for the thermal stability of engineering alloys. In Ni-Cr-based alloys, e.g., Alloy 690, the formation of long-range ordered Ni 2 Cr causes embrittlement and may impact the lifetime of nuclear power plant components. In this work, we quantify the formation and evolution of Ni 2 Cr precipitation in eleven Ni-Cr-Fe model alloys with 0, 5, 7, and 10wt % Fe contents, and with Ni/Cr atomic ratios of 1.8, 2.0, 2.2, 2.4. These alloys were isothermally aged up to 10,000h at temperatures between 330 and 475°C. The alloys were characterized by synchrotron-based x-ray diffraction and Vickers hardness testing to quantify Ni 2 Cr precipitate size, and the impact of precipitate size on the mechanical properties as a function of Fe content. After 10,000 h of aging at 475°C and 418°C, the formation of Ni 2 Cr was observed in all alloys with 0 and 5wt % Fe. After 10,000h of aging at 418°C, Ni 2 Cr precipitates were also observed in the 7wt % Fe containing Ni/Cr=2.0 sample. Further, no clear evidence of Ni 2 Cr was observed in any of the 10wt % Fe samples at any time and temperature combination. We find that the face-centered cubic matrix lattice contraction and Vickers hardness are correlated with the Ni 2 Cr formation. The greatest change in hardness and lattice contraction occurs in stoichiometric alloys (Ni/Cr=2.0) with 0wt % Fe at 475°C. The rate of change in the material properties for the 5wt % Fe alloys is reduced, however the magnitude of changes is similar to 0wt % Fe alloys. A precipitation hardening model developed for Ni-Cr alloys based on critical resolved shear stress with weakly coupled dislocations shows a clear link between Ni 2 Cr precipitate size and hardness. This trend held across all alloys with Ni 2 Cr formation regardless of Fe concentration. This important structure-property relationship can potentially help define Ni-Cr-Fe-based component lifetimes directly through an understanding of how Ni 2 Cr formation impacts mechanical properties as a function of Fe content.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Unraveling the transformation pathway of the 𝛽 to 𝛾 phase transition in Ga 2 ⁢O 3 from atomistic simulations

Defect spinel 𝛾−Ga 2 ⁢O 3 is the least stable polymorph of Ga 2 ⁢O 3 , so its frequent appearance as a structural defect within or on the surface of monoclinic 𝛽−Ga 2 ⁢O 3 remains a mystery. Through first-principles calculations, we explore potential pathways for the phase transition from 𝛽−Ga 2⁢ O 3 to 𝛾−Ga 2 ⁢O 3 , and examine two key driving forces: tensile strain and Ga deficiency. When configurational entropy contributions to phase energies are included, the 𝛾 phase becomes energetically competitive with the 𝛽 phase, with the free energy difference between these phases diminishing even further under Ga-deficient conditions. Notably, a stability crossover occurs at room temperature at high vacancy concentrations ([V$^{3−}_{Ga}$]>3%) . A simple model 𝛽 → 𝛾 transformation pathway is identified, comprising two primary reactions, that enables the formation of the 𝛾 phase via simultaneous migration of Ga atoms from tetrahedral lattice sites to octahedral interstitial positions. The transformation barriers are prohibitively large in pristine Ga 2 ⁢O 3 , but can be substantially reduced by: (1) the presence of Ga vacancies, (2) elongational strains along the crystallographic 𝑎-axis, and (3) when volumetric relaxations are possible during transformation. These results elucidate prior experimental observations, where 𝛾−Ga 2⁢ O 3 is seen on damaged surfaces or in highly 𝑛-type 𝛽−Ga 2⁢ O 3 environments, which support Ga deficiency and mechanical strain. The insights into the driving forces and mechanisms of 𝛾−Ga 2⁢ O 3 formation enhance understanding of how localized strain and nonequilibrium defect concentrations may facilitate its formation from the 𝛽 phase.

Defects↗

DOE OSTI2024

First Principles Evaluation of Phase Stability in the In-Sn Binary System

Abstract The In-Sn binary alloy system exhibits several unusual features that challenge crystallographic and thermodynamic expectations. We combine first principles total energy calculation with simple thermodynamic modeling to address two key points. First, we evaluate energies along the Bain path to interpret the discontinuous transition between the phases α-In (Pearson type tI2) and β-In 3 Sn (also Pearson type tI2) that are identical in symmetry. Second, we demonstrate that the solid solution phases β-In 3 Sn and γ-InSn 4 (Pearson type hP1) exist at high temperatures only, and they exhibit eutectoid decompositions at low temperatures.

Chemistry↗

NASA NTRS1970

Space shuttle: MMC phase A space shuttle modified orbiter: Investigation of stability characteristics and control effectiveness

Trisonic wind tunnel stability tests on space shuttle modified orbiter model

Sarver, D.↗

NASA NTRS2015

Robust Stability Analysis of the Space Launch System Control Design: A Singular Value Approach

Classical stability analysis consists of breaking the feedback loops one at a time and determining separately how much gain or phase variations would destabilize the stable nominal feedback system. For typical launch vehicle control design, classical control techniques are generally employed. In addition to stability margins, frequency domain Monte Carlo methods are used to evaluate the robustness of the design. However, such techniques were developed for Single-Input-Single-Output (SISO) systems and do not take into consideration the off-diagonal terms in the transfer function matrix of Multi-Input-Multi-Output (MIMO) systems. Robust stability analysis techniques such as H(sub infinity) and mu are applicable to MIMO systems but have not been adopted as standard practices within the launch vehicle controls community. This paper took advantage of a simple singular-value-based MIMO stability margin evaluation method based on work done by Mukhopadhyay and Newsom and applied it to the SLS high-fidelity dynamics model. The method computes a simultaneous multi-loop gain and phase margin that could be related back to classical margins. The results presented in this paper suggest that for the SLS system, traditional SISO stability margins are similar to the MIMO margins. This additional level of verification provides confidence in the robustness of the control design.

Pei, Jing↗

DOE OSTI2021

Three-Dimensional Free Gas Flow Focuses Basin-Wide Microbial Methane to Concentrated Methane Hydrate Reservoirs in Geological System

Here, we present a systematic model that links the generation, migration, phase partitioning, and accumulation of methane into a closed loop as the sediment is deposited from the seafloor and buried through the base of hydrate stability zone (BHSZ). In our model, methane is generated by biodegradation of organic carbon in muds. Hydrate does not form and methane is not trapped until a coarse-grained layer is deposited, because the small pores prevent hydrate formation in muds. Instead, methane diffuses into sands/silts where methane solidifies into hydrate. As hydrate-bearing sands/silts pass through the BHSZ during sediment burial, methane hydrate dissociates, and releases free gas. The released and the newly generated free gas below the BHSZ concentrates into a vertical/dipping zone with low capillary entry pressure and high permeability and flows upward driven its buoyancy. When free gas reaches the hydrate stability zone (HSZ), capillary forces drive free gas to flow laterally, preferentially enter sands/silts, feed hydrate growth, and elevate hydrate saturation. With three-dimensional focused free gas flow, microbial methane that is generated from a much larger fetch area of the entire basin, both above and below the BHSZ, is concentrated into coarse-grained layers at structural closures for hydrate formation. Our model illustrates how geological evolution, microbial methane generation, and gas flow by buoyancy couple to generate concentrated hydrate deposits in geological system. These insights can be used to explore for high-concentration methane hydrate and are important for understanding the methane budget and carbon cycle under the seafloor.

58 GEOSCIENCES↗

NASA NTRS1994

Parameter dimension of turbulence-induced phase errors and its effects on estimation in phase diversity

The method of phase diversity has been used in the context of incoherent imaging to estimate jointly an object that is being imaged and phase aberrations induced by atmospheric turbulence. The method requires a parametric model for the phase-aberration function. Typically, the parameters are coefficients to a finite set of basis functions. Care must be taken in selecting a parameterization that properly balances accuracy in the representation of the phase-aberration function with stability in the estimates. It is well known that over parameterization can result in unstable estimates. Thus a certain amount of model mismatch is often desirable. We derive expressions that quantify the bias and variance in object and aberration estimates as a function of parameter dimension.

Thelen, Brian J.↗

DOE OSTI2025

Enhancing stability, magnetic anisotropy, and coercivity of manganese aluminum: Machine learning, ab initio , and micromagnetic modeling

The binary manganese aluminum (MnAl) alloy with L⁢1 0 crystal structure is a promising rare earth (RE) element-free permanent magnetic material because of its exceptional magnetic properties. However, experimentally synthesizing it in a stable bulk form is extremely challenging. Here, in this study, an alternative method of stabilizing the material, a pathway for experimental synthesis and validation, is proposed and theoretically verified. This is done by partially substituting Mn and Al sites with Fe and Ni and identifying its enhanced phase stability, saturation magnetization density, magnetic anisotropy, and coercivity from density functional theory (DFT), machine learning (ML) crystal graph convolution neural network (CGCNN), and micro-magnetic modeling. When considering a fixed 50% Ni, the magnetic anisotropy increases with the increasing Fe content but decreases the formation energy. The calculated formation energies, elastic constants, and phonon frequencies demonstrate that the binary and quaternary compositions are stable. Most importantly, in 50% Fe and Ni-substituted-equiatomic phase, magnetic anisotropy constants and saturation magnetization density increase by 56% and 23% as compared to the MnAl. Further, the coercivity of the equiatomic phase predicted with micro-magnetic modeling is higher by 17% than the parent compound.

Bhandari, Churna [Ames National Laboratory, and Io↗

DOE OSTI2022

Hubbard-corrected oxide formation enthalpies without adjustable parameters

A density functional theory (DFT) approach to computing transition metal oxide heat of formation without adjustable parameters is presented. Different degrees of d-electron localization in oxides are treated within the DFT+U approach with site-dependent, first-principles Hubbard U-parameters obtained from linear response theory, and delocalized states in the metallic phases are treated without Hubbard corrections. Comparison of relative stabilities of these differently treated phases is enabled by a local d-electron density matrix-dependent model, which was found by genetic programming against experimental reference formation enthalpies. This mathematically simple model does not explicitly depend on the Hubbard-corrected ionic species and is shown to reproduce the heats of formation of the Mott insulators Ca 2 RuO 4 and Y 2 Ru 2 O 7 within ~3% of experimental results, where the experimental training data did not contain Ru oxides. This newly developed method thus absolves from the need for element-specific corrections fitted to experiments in existing Hubbard-corrected approaches to the prediction of reaction energies of transition metal oxides and metals. The absence of fitting parameters opens up here the possibility to predict relative thermodynamic stabilities and reaction energies involving d-states of varying degree of localization at transition metal oxide interfaces and defects, where site-dependent U-parameters will be particularly important and devising a fitting scheme against experimental data with predictive power would be exceedingly difficult.

transition metal oxides↗

NASA NTRS1974

On the steady flow of gas from the nuclei of Seyfert galaxies

Heating of gas energy released during explosive events produces hot winds in Seyfert nuclei. Steady-state solutions indicate that supersonic velocities up to about 800 km per sec are obtainable. Lyman continuum radiation emitted by the wind ionizes and maintains temperature of cooler gas clouds seen in emission, and dominates energy losses during early stages of expansion. A study of stability of the thermal equilibrium state leads to a two-component model: a hot stable phase (the wind), in pressure equilibrium with a cooler stable phase (clouds). Extreme variations in cloud electron densities are a consequence of the steep decrease with radius of external pressure. Cloud dynamics are related to the radius of formation and internal parameters. Permitted lines of H are emitted by high-density clouds formed at the base of the wind, while forbidden lines arise in radially moving clouds, produced throughout the wind.

Wolfe, A. M.↗

NASA NTRS1976

Compensated control loops for a 30-cm ion thruster

The vaporizer dynamic control characteristics of a 30-cm diameter mercury ion thruster were determined by operating the thruster in an open loop steady state mode and then introducing a small sinusoidal signal on the main, cathode, or neutralizer vaporizer current and observing the response of the beam current, discharge voltage, and neutralizer keeper voltage, respectively. This was done over a range of frequencies and operating conditions. From these data, Bode plots for gain and phase were made and mathematical models were obtained. The Bode plots and mathematical models were analyzed for stability and appropriate compensation networks determined. The compensated control loops were incorporated into a power processor and operated with a thruster. The time responses of the compensated loops to changes in set points and recovery from arc conditions are presented.

Robson, R. R.↗

NASA NTRS1976

Compensated control loops for a 30-cm ion thruster

Robson, R. R.↗