Search NASA⌕ Search

SEARCH · Search NASA

Results for “Equilibria”

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 73 records · Page 4

Machine Learning Applications to Metal-Silicate Equilibria and their Insights into Core Formation

An extensive number of studies have experimentally investigated how elements distribute between metal and silicate phases, to better constrain core-mantle chemical equilibrium. Here, we present a new database compiling all (to our knowledge) experimental data on liquid metal-silicate partitioning from 118 peer-reviewed publications. We applied various machine learning techniques to gain further insights into these partitioning equilibria and their dependencies. We performed a network analysis to investigate the relationship between experiments and partition coefficients, which enables visualizing gaps in the experimental dataset and biases related to varying experimental conditions and analytical setup. In addition, semi-empirical thermodynamic models are commonly used to extrapolate these chemical reactions to the wide range of pressure, temperature and compositional conditions of planetary differentiation. These models are based on linear regressions that assume continuous relationship between partition coefficients and experimental variables. Here, we considered random forest regressions, which are algorithms based on ensembles of decision trees and does not consider continuous effects of each variable. The application of this regression significantly improves the prediction of metal-silicate partitioning for several elements including Ni, Si and Cr. We will show how this new approach improves our understanding of elemental exchange between metal and silicate and their implications for the Earth’s core formation.

siderophile element↗

Defect Equilibria from First Principles: From Widegap Oxides to Topological Semimetals

Materials functionality and performance is rarely determined by the ideal crystal alone but is usually affected by formation of imperfections and the solution of impurities. In some applications, such as solar thermochemical hydrogen generation, defect formation is the fundamentally enabling mechanism of the desired functionality. In other cases, such as Cd3As2 topological semimetals, unintentional self-doping presents an obstacle to the access to the unique electronic properties. In either case, a quantitative understanding of the relevant defect mechanism is essential for developing design strategies. This presentation will touch upon numerous aspects in the computational simulation of defect equilibria, including non-equilibrium design strategies, the coupling of solid state and gas-phase reactions, dopant-defect and defect-defect interactions, both attractive and repulsive, the accuracy of total energy functionals and electronic structure methods, and the role of the shape of the density of states for the charge balance condition and Fermi level position, as well as machine-learning prediction of defect energies (1). Specific materials systems include Ga2O3 (2), Cd3As2 (3), and (Sr,Ce)MnO3 (4). (1) M.D. Witman, A. Goyal, T. Ogitsu, A.H. McDaniel, S. Lany, Nat. Comput. Sci. 3, 675 (2023). (2) A. Goyal, A. Zakutayev, V. Stevanovic, S. Lany, J. Appl. Phys. 129, 245704 (2021). (3) C. Brooks, M. van Schilfgaarde, D. Pashov, J.N. Nelson, K. Alberi, D.S. Dessau, S. Lany, Phys. Rev. B 107, 224110 (2023). (4) A. Goyal, M.D. Sanders, R.P. O'Hayre, S. Lany, PRX Energy 3, 013008 (2024).

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,M↗

The Ce-Fe-Ti System: Phase Equilibria in the Fe-rich Corner at 1000 °C

The rare earth–iron–transition metal systems are of increasing interest in the search for novel permanent magnet phases. This work reports on the solidification behavior, phase equilibria and stability range of solid phases in the ternary Ce-Fe-Ti system, focusing on the iron-rich region (> 65 at.% Fe) of the isothermal section at 1000 °C. Further, two ternary phases Ce 1.03 Fe 12-x Ti x (x = 0.87–1.02) with ThMn 12 structure type and Ce 3.06 Fe 27.6 Ti 1.4 with Nd 3 Fe 29 structure type were observed. Magnetic measurements of Ce 1.03 Fe 12-x Ti x and Ce 3.06 Fe 27.6 Ti 1.4 revealed ferromagnetic ordering with Curie temperatures of 550 K and 327 K, respectively.

36 MATERIALS SCIENCE↗

Re-investigation of high-temperature phase equilibria in Fe-rich Sm–Fe–Ti alloys

Recently renewed attempts to develop high-performance rare-earth-lean permanent magnets based on the Sm(Fe,Ti) 12 compound have drawn attention to the limited knowledge about the high-temperature phase equilibria in the Sm–Fe–Ti system. Experimental investigation of equilibrated alloys with electron probe microanalysis, X-ray diffraction and thermomagnetic analysis revealed several inaccuracies in the currently accepted phase relations at 1000°C and allowed for a revision of the Fe-rich corner of the Sm–Fe–Ti phase diagram. The Sm(Fe,Ti) 12 and Sm 3 (Fe,Ti) 29 phases were found to have more extended Ti ranges of 5.1–9.7 at% and 2.8–6.9 at%, respectively. With increasing of the Ti content, the Curie temperature of the Sm(Fe,Ti) 12 remains nearly constant at 306–312°C, whereas that of the Sm 3 (Fe,Ti) 29 increases from 188°C to 207°C. The low-titanium Sm 3 (Fe,Ti) 29 phase equilibrates not only with the Sm(Fe,Ti) 12 and Sm2(Fe,Ti) 17 phases, but also with (α-Fe) solid solution. Newly demonstrated equilibrium between Sm 2 (Fe,Ti) 17 and TiFe 2 phases makes impossible the earlier reported equilibrium between the Sm 3 (Fe,Ti) 29 and Sm(Fe,Ti) 11 phases. Because of an invariant reaction at 1000 °C, the revised phase diagram also features a class II four-phase equilibrium Sm 3 (Fe,Ti) 29 + TiFe 2 + Sm(Fe,Ti) 12 + Sm 2 (Fe,Ti) 17 . Peritectic decomposition of the Sm(Fe,Ti) 11 phase, which occurs either at 1075°C or at 1087°C, was found to have among its products the Sm(Fe,Ti) 12 phase. Although no such equilibration was attempted, it must be possible to obtain above 1087°C a two-phase state composed of the Sm(Fe,Ti) 12 phase and a liquid – which is important for manufacturing of the Sm(Fe,Ti) 12 - based permanent magnets via the liquid-phase sintering.

36 MATERIALS SCIENCE↗

Phase equilibria in iron-rich Sm–Fe–Ti and Sm–(Fe,Co)–Ti alloys at 1100–1200 °C

Iron-rich corners of ternary Sm–Fe–Ti phase diagram at 1100 °C and quasi-ternary Sm–Fe 0.8 Co 0.2 –Ti phase diagrams at 1100 and 1200 °C are constructed based on experimental investigation of equilibrated alloys with electron probe microanalysis, X-ray diffraction and thermomagnetic analysis. In addition, the upper boundaries of the temperature ranges of Sm(Fe,[Co,]Ti) 12 and Sm 3 (Fe,[Co,]Ti) 29 phases are determined with differential thermal analysis to update earlier rough estimates. The existence of a high-temperature phase of the Th 2 Ni 17 type, originally reported by Ivanova et al. [J. Alloys Compd. 224 (1995) 29], is confirmed. In the Sm–Fe 0.8 Co 0.2 –Ti system, the composition and equilibria of this hexagonal phase are established for 1200 °C; it is Sm-depleted (≈9.8 at.% Sm) compared to the 2:17 stoichiometry and it coexists with the rhombohedral 2:17 phase. The magnetic anisotropy of the cobalt-substituted Th 2 Ni 17 -type phase is planar, with the easy magnetization direction parallel to [100]. Equilibrium between a Sm-rich liquid phase and the 1:12 phase, which is important for the development of new high-performance permanent magnets, is absent up to 1000 °C, but does exist at 1100 °C (for the 1:12 phase with at least 8.7–8.9 at.% Ti) and at 1200 °C (for the 1:12 phase with as little as 7.4 at.% Ti). The development of magnets may be complicated, however, by an observed tendency of the high-temperature liquid to solidify into ferromagnetic phases including the Th 2 Ni 17 -type phase. The Curie temperatures of the α-(Fe,Ti), Sm(Fe,Ti) 12 , Sm 3 (Fe,Ti) 29 and rhombohedral Sm 2 (Fe,Ti) 17 phases are not only increased by the partial Co substitution for Fe, but their dependence on the Ti concentration is changed by this Co substitution from positive (or, for the 1:12 phase, zero) to negative values.

36 MATERIALS SCIENCE↗

Simple, general, realistic, robust, analytic tokamak equilibria. Part 2. Pedestals and flow

Part 1 described a wide range of analytic tokamak equilibria modelling smooth limiter surfaces, double- and single-null divertor surfaces, arbitrary aspect ratio, elongation, triangularity and beta. Part 2 generalizes the analysis to further include edge pedestals and toroidal flow. Specifically, edge pedestals are allowed in the pressure, pressure gradient and toroidal current density. Also, an edge-localized contribution to the bootstrap current is treated. In terms of flow, analytic solutions are obtained for two cases: a $\gamma = 2$ adiabatic and a $\gamma = \infty $ incompressible energy conservation relation.

Physics↗

Kinetic-ballooning-mode turbulence in low-average-magnetic-shear equilibria

In this work, kinetic-ballooning-mode (KBM) turbulence is studied via gyrokinetic flux-tube simulations in three magnetic equilibria that exhibit small average magnetic shear: the Helically Symmetric eXperiment (HSX), the helical-axis Heliotron-J and a circular tokamak geometry. For HSX, the onset of KBM being the dominant instability at low wavenumber occurs at a critical value of normalized plasma pressure $\beta ^{\rm KBM}_{\rm crit}$ that is an order of magnitude smaller than the magnetohydrodynamic (MHD) ballooning limit $\beta ^{\rm MHD}_{\rm crit}$ when a strong ion temperature gradient (ITG) is present. However, $\beta ^{\rm KBM}_{\rm crit}$ increases and approaches the MHD ballooning limit as the ITG tends to zero. For these configurations, $\beta ^{\rm KBM}_{\rm crit}$ also increases as the magnitude of the average magnetic shear increases, regardless of the sign of the normalized magnetic shear. Simulations of Heliotron-J and a circular axisymmetric geometry display behaviour similar to HSX with respect to $\beta ^{\rm KBM}_{\rm crit}$ . Despite large KBM growth rates at long wavelengths in HSX, saturation of KBM turbulence with $\beta > \beta _{\rm crit}^{\rm KBM}$ is achievable in HSX and results in lower heat transport relative to the electrostatic limit by a factor of roughly five. Nonlinear simulations also show that KBM transport dominates the dynamics when KBMs are destabilized linearly, even if KBM growth rates are subdominant to ITG growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Asymptotic quasisymmetric high-beta three-dimensional magnetohydrodynamic equilibria near axisymmetry

Quasisymmetry (QS), a hidden symmetry of the magnetic field strength, is known to support nested flux surfaces and provide superior particle confinement in stellarators. In this work, we study the ideal magnetohydrodynamic (MHD) equilibrium and stability of high-beta plasma in a large-aspect-ratio stellarator. In particular, we show that the lowest-order description of a near-axisymmetric equilibrium vastly simplifies the problem of three-dimensional quasisymmetric MHD equilibria, which can be reduced to a standard elliptic Grad–Shafranov equation for the flux function. We show that any large-aspect-ratio tokamak, deformed periodically in the vertical direction, is a stellarator with approximate volumetric QS. We discuss exact analytical solutions and numerical benchmarks. Finally, we discuss the ideal ballooning and interchange stability of some of our equilibrium configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Narrow operator models of stellarator equilibria in Fourier Zernike basis

Numerical computation of the ideal magnetohydrodynamic (MHD) equilibrium magnetic field is at the base of stellarator optimisation and provides the starting point for solving more sophisticated partial differential equations like transport or turbulence models. Conventional approaches solve for a single stationary point of the ideal MHD equations, which is fully defined by three invariants and the numerical scheme employed by the solver. We present the first numerical approach that can solve for a continuous distribution of equilibria with fixed boundary and rotational transform, varying only the pressure invariant. This approach minimises the force residual by optimising parameters of multilayer perceptrons that map from a scalar pressure multiplier to the Fourier Zernike basis as implemented in the modern stellarator equilibrium solver DESC.

fusion plasma↗

Vortex fluidic regulated phospholipid equilibria involving liposomes down to sub-micelle size assemblies

Conventional channel-based microfluidic platforms have gained prominence in controlling the bottom-up formation of phospholipid based nanostructures including liposomes. However, there are challenges in the production of liposomes from rapidly scalable processes. These have been overcome using a vortex fluidic device (VFD), which is a thin film microfluidic platform rather than channel-based, affording ~110 nm diameter liposomes. The high yielding and high throughput continuous flow process has a 45° tilted rapidly rotating glass tube with an inner hydrophobic surface. Processing is also possible in the confined mode of operation which is effective for labelling pre-VFD-prepared liposomes with fluorophore tags for subsequent mechanistic studies on the fate of liposomes under shear stress in the VFD. In situ small-angle neutron scattering (SANS) established the co-existence of liposomes ~110 nm with small rafts, micelles, distorted micelles, or sub-micelle size assemblies of phospholipid, for increasing rotation speeds. The equilibria between these smaller entities and ~110 nm liposomes for a specific rotational speed of the tube is consistent with the spatial arrangement and dimensionality of topological fluid flow regimes in the VFD. The prevalence for the formation of ~110 nm diameter liposomes establishes that this is typically the most stable structure from the bottom-up self-assembly of the phospholipid and is in accord with dimensions of exosomes.

60 APPLIED LIFE SCIENCES↗

Shape equilibria of vesicles with rigid planar inclusions

Motivated by recent studies of two-phase lipid vesicles possessing 2D solid domains integrated within a fluid bilayer phase, we study the shape equilibria of closed vesicles possessing a single planar, circular inclusion.

Chemistry↗

All-in-one probe for exploring self-organized two-fluid equilibria in toroidal plasmas

This paper presents the development of an all-in-one probe to simultaneously measure all components of the generalized Ohm’s law in reversed-field pinch plasmas and tokamaks. The polyhedral configuration of the Mach probe is achieved through the specific arrangement, angle, and depth of the collimator channel apertures drilled into the surface of a hollow boron nitride cylinder encasing it. This probe includes a central Mach probe to assess the ion velocity field in three dimensions. Initial tests at the RELAX and Madison Symmetric Torus machines have confirmed the probe’s effectiveness, revealing an octahedron form similar to a tetrahedron. The probe seems to function correctly and is expected to facilitate the empirical validation of two-fluid equilibria at the periphery of toroidal plasmas.

Instruments & Instrumentation↗

Rigidity of MHD equilibria to smooth incompressible ideal motion near resonant surfaces

In ideal MHD, the magnetic flux is advected by the plasma motion, freezing flux-surfaces into the flow. An MHD equilibrium is reached when the flow relaxes and force balance is achieved. We ask what classes of MHD equilibria can be accessed from a given initial state via smooth incompressible ideal motion. It is found that certain boundary displacements are formally not supported. This follows from yet another investigation of the Hahm–Kulsrud–Taylor (HKT) problem, which highlights the resonant behaviour near a rational layer formed by a set of degenerate critical points in the flux-function. When trying to retain the mirror symmetry of the flux-function with respect to the resonant layer, the vector field that generates the volume-preserving diffeomorphism vanishes at the identity to all order in the time-like path parameter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the non-existence of stepped-pressure equilibria far from symmetry

The Stepped Pressure Equilibrium Code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 112502) has been successful in the construction of equilibria in 3D configurations that contain a mixture of flux surfaces, islands and chaotic magnetic field lines. In this model, the plasma is sliced into sub-volumes separated by ideal interfaces, and in each volume the magnetic field is a Beltrami field. In the cases where the system is far from possessing a continuous symmetry, such as in stellarators, the existence of solutions to a stepped-pressure equilibrium with given constraints, such as a multi-region relaxed MHD minimum energy state, is not guaranteed but is often taken for granted. Using SPEC, we have studied two different scenarios in which a solution fails to exist in a slab with analytic boundary perturbations. Here we found that with a large boundary perturbation, a certain interface becomes fractal, corresponding to the break up of a Kolmogorov–Arnold–Moser (KAM) surface. Moreover, an interface can only support a maximum pressure jump while a solution of the magnetic field consistent with the force balance condition can be found. An interface closer to break-up can support a smaller pressure jump. We discovered that the pressure jump can push the interface closer to being non-smooth through force balance, thus significantly decreasing the maximum pressure it can support. Our work shows that a convergence study must be performed on a SPEC equilibrium with interfaces close to break-up. These results may also provide insights into the choice of interfaces and have applications in finding out the maximum pressure a machine can support.

3D equilibrium↗

Single-stage stellarator optimization: combining coils with fixed boundary equilibria

We introduce a novel approach for the simultaneous optimization of plasma physics and coil engineering objectives using fixed-boundary equilibria that is computationally efficient and applicable to a broad range of vacuum and finite plasma pressure scenarios. Our approach treats the plasma boundary and coil shapes as independently optimized variables, penalizing the mismatch between the two using a quadratic flux term in the objective function. Four use cases are presented to demonstrate the effectiveness of the approach, including simple and complex stellarator geometries. As shown here, this method outperforms previous two-stage approaches, achieving smaller plasma objective function values when coils are taken into account.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Large-Signal Stability of Phase-Balanced Equilibria in Single-Phase Grid-Forming Inverter Systems

This article explores the setup where large numbers of single-phase grid-forming inverters with droop control across distribution networks self-organize into a stable and balanced system with 120° phase offsets across aggregates in the absence of balanced three-phase generating resources or external communication. A suite of circuit- and system-theoretic notions are leveraged to derive a dynamical model for phase-angle differences across aggregates of inverters connected in the three phases. Focusing on this model, large-signal stability is established and the region of attraction of the phase-balanced equilibria is determined with the aid of a Lyapunov function. Experimental validation for a bench-top prototype network is included to support the analytical developments. Altogether, the effort supports the vision of facilitating balanced operation of distribution networks with grid-forming inverters during service disruptions at the bulk transmission network.

Droop control↗

Thermochemical Modeling of Radionuclide Vapor-Liquid Equilibria in Sodium Pools for SFR Mechanistic Source Term Analysis

Mechanistic source term (MST) analysis of sodium fast reactors (SFR) requires understanding of various radionuclide (RN) transport phenomena influencing potential releases from the fuel to the environment. One such phenomenon includes the retention or release of RNs from the sodium coolant pool, representing the step after possible fuel failures and influencing transport to the cover gas region. Thermodynamic vapor-liquid equilibria (VLE) calculations were performed on systems representing SFR sodium pools containing oxygen impurities and radionuclide (RN) inventories. First, an assessment and recreation of a previously developed thermodynamic database was completed, including updates to thermodynamic parameters. The RN inventories used in VLE calculations represented hypothetical source terms that might be released to the pool during previously analyzed fuel failure scenarios. The calculations were performed for all possible combinations of sodium pool size (i.e., total oxygen) and number of failed fuel pins (i.e., total RNs). In this way, multiple ratios of the RN relative to the oxygen impurity (RN:O) were compared for their impact on RN volatility, which is discussed in terms of the vapor fraction (VF), defined as the fraction of the RN that is calculated to exist in the vapor phase at equilibrium above condensed phases of that element. Similar trends in VLE behavior are seen in the equilibrium calculation results for elements of similar chemistry, and for some element types, it was found that the RN:O ratio can be important due to oxide formation, which typically exist in the condensed phase.

Shahbazi, Shayan↗