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 271 records · Page 15

Improving ideal MHD equilibrium accuracy with physics-informed neural networks

We present a novel approach to compute three-dimensional magnetohydrodynamic equilibria with isotropic pressure profiles and nested surfaces by parametrizing Fourier modes with artificial neural networks (NNs). The full nonlinear global force residual of single equilibria across the volume in real space is then minimized with first order optimizers and compared to equilibria computed by conventional solvers. Already, we observe competitive computational cost to arrive at the same minimum residuals computable with existing codes. With increased computational cost, lower minima of the residual are computable with the NNs than with any other tested solver, establishing a new lower bound for the force residual. We use minimally complex NNs, and we expect significant improvements for solving not only single equilibria with NNs, but also for creating NN models valid over continuous distributions of equilibria.

ideal magnetohydrodynamics↗

Topological constraints and the existence of force-free fields

A fundamental problem in plasma theory is the question of the existence of MHD equilibria. The issue of topological constraints is of crucial importance for the problem of the existence of equilibria. Heuristic methods are used to discuss the coronal wrapping pattern. It is concluded that for a given set of footpoint positions the wrapping pattern in the corona is completely fixed. The topological constraints are included in the boundary conditions on the Euler potentials and impost no additional restrictions on possible equilibria. Although this does not prove that equilibria always exist, it does show that the force-free problem is not overdetermined and that existence of equilibria is still an open question.

Antiochos, S. K.↗

Driven resonant current amplification in self-organized plasma configurations with uniform λ and plasma pressure confinement

Driven resonant current amplification in uniform λ=$^{μoJ}_{B}$, self-organized plasma configurations with plasma pressure confinement is studied. This work extends previous analyses of force-free states (J ∥ B) to equilibria with plasma pressure gradients (∇P = J × B ≠ 0). An analytical treatment shows that in equilibria with uniform λ and a plasma pressure gradient that varies linearly with poloidal magnetic flux, driven resonances are either downshifted or upshifted from their force-free values. Additionally, a modification of off resonance current amplification is observed in pressure confining equilibria. Computational results add support to analytical predictions, and also include driven resonance behavior with equilibrium profiles not considered analytically. Furthermore, these results suggest not only persistence, but also modification of driven resonant current amplification in sustained spheromak configurations with plasma pressure confinement when compared to force-free states. By providing a basis for the existence of driven, high current amplification, pressure confining equilibria, this work provides insight into the accessibility of sustained spheromak states of most interest for potential fusion energy applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Feedforward equilibrium trajectory optimization with GSPulse

One of the common tasks required for designing new plasma scenarios or evaluating capabilities of a tokamak is to design the desired equilibria using a Grad-Shafranov (GS) equilibrium solver. However, most standard equilibrium solvers are time-independent and do not include dynamic effects such as plasma current flux consumption, induced vessel currents, or voltage constraints. Another class of tools, plasma equilibrium evolution simulators, do include time-dependent effects. These are generally structured to solve the forward problem of evolving the plasma equilibrium given feedback-controlled voltages. In this work, we introduce GSPulse, a novel algorithm for equilibrium trajectory optimization, that is more akin to a pulse planner than a pulse simulator. GSPulse includes time-dependent effects and solves the inverse problem: given a user-specified set of target equilibrium shapes, as well as limits on the coil currents and voltages, the optimizer returns trajectories of the voltages, currents, and achievable equilibria. This task is useful for scoping performance of a tokamak and exploring the space of achievable pulses. The computed equilibria satisfy both Grad-Shafranov force balance and axisymmetric circuit dynamics. The optimization is performed by restructuring the free-boundary equilibrium evolution equations into a form where it is computationally efficient to optimize the entire dynamic sequence. GSPulse can solve for hundreds of equilibria simultaneously within a few minutes. GSPulse has been validated against NSTX-U and MAST-U experiments and against SPARC feedback control simulations, and is being used to perform scenario design for SPARC. The computed trajectories can be used as feedforward inputs that are connected to the feedback controller to inform and improve feedback performance. The code for GSPulse is available open-source at github.com/jwai-cfs/GSPulse_public.

equilibrium↗

Validation of an integrated modeling framework for investigating 3D plasma responses in tokamak plasmas

As contemporary experimental tokamaks are pushed toward reactor-relevant operation, they provide essential testbeds for demonstrating and exploring ELM control strategies for deployment in future fusion pilot plants. Accurate predictions of full plasma responses are essential to guide and optimize these demonstrations. This paper introduces and validates an integrated modeling framework over a historical range of DIII-D operational space. The integrated modeling framework uses only scalar plasma parameters and optional reference boundary to self-consistently and flexibly scan through tokamak operational space and estimate the corresponding ELM-suppression relevant plasma response over 3D coil phase space. The framework generates tightly converged equilibria that satisfy a target set of plasma parameters (I p , β N , l i ), with kinetic profiles constrained by an EPED(NN)-computed pedestal and empirical core model. The plasma response of these modeled equilibria is calculated with GPEC. Validated against 55 distinct DIII-D equilibria, the framework consistently reproduces experimentally constrained equilibrium pressure, q, and other representative profiles using the time-varying information from only evolving scalar plasma parameters; plasma response validation is performed for the n=3 perturbation. In particular, the plasma response of a DIII-D discharge scanning q 95 in search of ELM suppression windows is quantitatively reproduced. By accurately reproducing experimental equilibria and plasma responses across wide parameter variation, the framework supports potential for synthetic parameter scans in key operational, stability, and plasma response dimensions to investigate RMP ELM-suppression experiments and inform predictive RMP scenario optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetohydrodynamic instabilities in coronal arcades

The MHD stability of coronal arcades is studied with and without a detached flux using a two-dimensional linear MHD stability numerical model. Two-dimensional magnetohydrostatic equilibria with and without gravity are computed. A coronal arcade without a detached flux tube is found to be stable for any magnetic shear and for any longitudinal mode. On the other hand, an arcade with a detached flux tube is unstable for perturbations with nonzero longitudinal wave number n and the instability mode structure and growth rate vary depending on the field twist and n. All the equilibria studied are stable to n = 0 perturbations. Gravity has a stabilizing effect on the equilibria. As the ratio lambda of the standard arcade width to the gravitational scale height increases from zero to 0.2, the m = 0 transverse mode growth rate decreases, but high m modes are stabilized. The equilibria studied here are completely stabilized for lambda larger than 0.33.

An, Chang-Hyuk↗

Multi-Anvil Experimentation Applied to Planetary Differentiation

Planets undergo differentiation that includes segregation of metal from silicate at high temperatures and pressures ranging from deep planetary core pressures (>300 GPa for Earth)to very shallow conditions of asteroids (<100 MPa).The multi-anvil solid media apparatus accesses the middle part of this range from 3 to 30 GPa –pressure relevant to the interior of Mercury, Venus, Earth, Earth’s Moon, and Mars. Early planets are thought to have experienced high temperatures from a combination of heat sources including radioactive decay, gravitational and accretional heating, and impact processes. These heating events led to melting of mantles and cores, thus requiring an understanding of solid-liquid equilibria in metal-silicate systems. In 2006 we established a multi-anvil facility at NASA-JSC combining an 880 ton press and a Kawai/Walker type module from Rockland Research. Our high PT work has been greatly facilitated by use of the COMPRES multi-anvil assemblies (1). Our recent work has included studies of element partitioning between liquid metal and liquid silicate(e.g., 2), as well as between minerals and melts(e.g., 3), both of which have led to better constraints on the timing and conditions of planetary differentiation(e.g. 4). Several examples involving sustained efforts will be summarized below and for the presentation. The distribution of siderophile (iron-loving) elements between core and mantle is controlled by metal-silicate equilibrium across a wide range of pressures. Therefore, experimentation across this pressure range helps to calibrate elemental partitioning models that can be applied to planets, and used to predict mantle chemistry and composition during planetary differentiation. Our studies have focused on a wide range of siderophile elements (refractory Ni, Co, W, Mo; volatile P, Ga, Cu, Sn, Sb; highly siderophile Au, Pd) that have constrained partitioning, valence, and isotopic fractionation, and applied to Earth, Moon, and Mars. When molten mantles (magma oceans) cool enough to initiate crystallization, the solids precipitate at depth and in large planets this involves high pressure phases like garnet, majorite, akimotoite, and ringwoodite. As these solids precipitate they can segregate from liquid by density contrasts, thus causing elemental fractionation which can be used to decipher timing of differentiation. Mineral/melt and metal/silicate equilibria in our lab have helped to better understand high pressure fractionation of isotopic parent/daughter pairs Hf/W, Mn/Cr, Pd/Ag, Pt/Os, Re/Os, and U/Pb, and their application to Earth, Moon and Mars. There remains great potential for multi-anvil experimentation to shed light on many pressure–dependent aspects of planetary evolution such as core formation, high pressure phase equilibria, redox equilibria, and volatile evolution and storage. References1. Leinenweber, K., et al.(2012) American Mineralogist,97, 353–368. 2.Righter, K., et al.(2020) Geochem. Persp. Lett.15, 1-6.3. Righter, K., et al.(2020) Met. Planet. Sci 55, 2741-2757.4. Righter, K., et al.(2020)Earth and Planetary Science Letters,552, 116590.

pressure↗

Thermodynamic modeling of adsorption at the liquid-solid interface

Adsorptive separation techniques are significantly energy efficient in comparison to conventional thermal separation techniques such as distillation. Despite extensive research and development activities undertaken for mixed-gas adsorption, the use of adsorption techniques for the separation of multicomponent liquid mixtures is still limited. A major barrier is the lack of accurate adsorption thermodynamic models, which form the scientific foundation of process simulation of such systems, making the translation to industrial scale challenging. In this work, we have rigorously computed the surface excess of adsorption for six binary liquid mixtures on silica gel at 303 K using the frameworks of the generalized Langmuir isotherm (gL) and the adsorbed solution theory (AST). The six binary liquid mixtures were formed by the pair-wise combinations of four components: benzene, 1,2-dichloroethane, cyclohexane, and n-heptane. We have based our calculations by considering simultaneous equilibria of three phases: saturated vapor phase, bulk liquid phase, and adsorbed phase. Further, the composition of the corresponding saturated vapor phase was estimated by the Nonrandom Two-Liquid activity coefficient model and experimental vapor-liquid equilibria data. The activity coefficients of the adsorbed phase, the central issue of multicomponent adsorption thermodynamics, were calculated using the adsorption Nonrandom Two-Liquid activity coefficient model. Devoid of simplifying assumptions, gL and AST provide rigorous thermodynamic frameworks for adsorption equilibria of multicomponent liquid mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solving the problem of overdetermination of quasisymmetric equilibrium solutions by near-axis expansions. I. Generalized force balance

It is well known that the process of construction of quasisymmetric magnetic fields in magnetostatic equilibrium with isotropic pressure suffers from the problem of overdetermination. This has led to the widespread belief that global quasisymmetric solutions are likely not to exist. Here we develop a general near-axis expansion procedure that does not rely on the assumption of magnetostatic equilibria with isotropic pressure. We then demonstrate that in equilibria with anisotropic pressure, it is possible to circumvent the problem of overdetermination and carry out the power-series solutions to higher order. This suggests, contrary to current belief, that the existence of globally quasisymmetric fields is likely if one relaxes the assumption of magnetostatic equilibria with isotropic pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stepped pressure equilibrium with relaxed flow and applications in reversed-field pinch plasmas

The multi-region relaxed magnetohydrodynamics (MRxMHD) has been successful in the construction of equilibria in three-dimensional (3D) configurations. In MRxMHD, the plasma is sliced into sub-volumes separated by ideal interfaces, each undergoing relaxation, allowing the formation of islands and chaos. The resulting equilibrium has a stepped pressure profile across sub-volumes. The stepped pressure equilibrium code (SPEC) was developed to calculate MRxMHD equilibria numerically. In this work, we have extended the SPEC code to compute MRxMHD equilibria with field-aligned flow and rotation, following the theoretical development to incorporate cross-helicity and angular momentum constraints. The code has been verified for convergence and compared to a Grad–Shafranov solver in 2D. We apply our new tool to study the flow profile change before and after the sawtooth crash of a reversed-field pinch discharge, in which data of the parallel flow is available. We find the promising result that under the constraints of cross-helicity and angular momentum, the parallel flow profile in post-crash SPEC equilibrium is flat in the plasma core and the amplitude of the flow matches experimental observations. Lastly, we provide an example equilibrium with a 3D helical field structure as the favored lower energy state. This will be the first 3D numerical equilibrium in which the flow effects are self-consistently calculated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Projected global stability of high beta MAST-U spherical tokamak plasmas

Assessment of the limits of stability of tokamak plasmas is essential to operation in high fusion performance ranges without disruption of the plasma current. Projected equilibria have been generated for the MAST-U spherical tokamak experiment, an upgrade of the previous MAST device, in order to prepare for operation. These equilibria are scanned in pressure and current profiles, and assessed with the DCON and MARS-F stability codes to find the so-called 'no-wall' beta limit, above which resistive wall mode instabilities can be expected in the absence of other stabilising effects. The no-wall limit was generally found to increase as plasma internal inductance increased. The equilibria are also assessed for the 'with-wall' limit, theoretically the highest achievable performance point, again with the DCON and MARS-F codes, including different approximate axisymmetric walls, and with the VALEN code which includes a 3D model of the surrounding conducting structure. Similar limits were found, despite the difference between the 2D and 3D codes in the treatment of the wall. Conducting passive stabilisation plates, which were newly installed in MAST-U, are in a region of significant mode perturbation when the plasma β N is sufficiently high and eddy currents are driven in these structures. Due to the increased stabilising effect of the wall in MAST-U vs. MAST, a significant gap exists between the approximate no-wall limits of β N / l i = 6.71 and 7.13, found from DCON and MARS-F respectively, and the with-wall limits of β N / l i = 8.23 and 8.53 for the equilibrium profiles analysed in this study. This opens a region of high beta operating space in MAST-U for potentially stable operation if non-ideal effects or active control can stabilise the resistive wall mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

FES 2023 Theory & Simulation Performance Target

The Report for the 2023 Theory & Simulation Performance Target (TSPT) is presented herein. The 1st, 2nd and 3rd Quarter Reports are given in Sec. I, Sec. II and Sec. III. These are unchanged from those previously submitted, except that in this final document the equations, figures and citations have been automatically renumbered. The 4th Quarter Report is given in Sec. IV. The following material presents verification calculations of several widely used three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium codes in strongly shaped configurations; for both fixed and free-boundary equilibria; for the so-called vacuum case and for the case of high plasma pressure; and, for magnetic fields that have nested magnetic surfaces and for fields that have a magnetic island. The configurations studied includes the recently discovered “precise-QA” case [Phys. Rev. Lett., 128, 035001 (2002)] and other highly optimized configurations. Because of the strong (extreme) shaping, the latter of these pose challenges to the 3D MHD codes. A new method is introduced for constructing coordinates in arbitrarily shaped toroidal domains, and this can improve the robustness of 3D equilibrium calculations. Beyond the equilibrium calculation, this report demonstrates the capability of performing non-linear extended-MHD modelling of optimized stellarator equilibria. This capability allows the assessment not only of linear stability thresholds but also of the nonlinear evolution of unstable equilibria to new configurations that might be far from the original configuration, with significantly different transport properties and magnetic topology. They also enable the self-consistent determination of pressure profiles in non-integrable magnetic topology through the use of highly anisotropic thermal-transport equations. The calculations presented here represent the first demonstration of extended-MHD calculations in strongly-shaped, optimized quasi-symmetric configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Influence of transient baroclinic eddies on planetary-scale waves

Low-order models with a mean zonal flow forced by zonally symmetric heating and one planetary-scale wave forced by topography have multiple flow equilibria. Two stable equilibria are characterized by strong zonal flow with low wave amplitude (high index circulation) and weak zonal flow with a high wave amplitude (low index circulation) fixed in phase with the topography. These two states presumably represent normal zonal circulation and a blocking configuration, respectively. When a shorter, baroclinically unstable wave is dded to this low-order model, the planetary-scale wave no longer stays in stable equilibrium states. Instead, the long wave remains in weather regimes or preferred regions in phase space. These regimes are, in general, different from the equilibria of the model with the planetary-scale wave only. Thus, the short unstable wave adds some randomness to planetary-scale circulation and changes its position with respect to the planetary-scale topography. The Reinhold and Pierrehumbert model is further explored here by adding a long-wave in the thermal forcing and exploring wider parameter space. When the symmetric thermal forcing is weak so that the short wave is stable, the planetary scale wave is fixed in amplitude and phase by the asymmetric forcing. However, when the symmetric forcing is increased to an unstable level, the planetary-scale wave becomes less organized. With asymmetric thermal forcing only (i.e., no topography) the long-wave is randomly distributed through phase space. Thermal forcing seems to be less effective in organizing weather regimes than topographic forcing.

Branscome, L. E.↗

Unified description of accretion flows around black holes

We provide a unified description of thermal equilibria of black hole accretion disks, including the newly discovered advection-dominated solutions. We classify the solutions on the basis of optical depth and impotance of advection cooling. We demonstrate that only four physically distinct topological types of equilibria exist. Two of the types correspond to optically thin and optically thick equilibria, while the other two types are distinguished by whether advection is negligible or dominant. A stable Shakura-Sunyaev disk exists only for accretion rates dot-M below a certain maximum. However, there is a critical viscosity parameter alpha(sub crit), which is a function of radius, such that for alpha greater than alpha(sub crit) advection-dominated solutions exist for all dot-M. Even when alpha less than alpha(sub crit), the advection-dominated solutions are available for a wide range of dot-M except for a gap around the Eddington rate. We therefore suggest that advection-dominated flows may be more common than standard thin disks in black hole systems. For certain ranges of radii and dot-M, no stable steady state solution is possible. In these cases, we suggest that limit cycle behavior may occur, leading to variability.

Chen, Xingming↗

The Origin of Monsoon Onset: Rotational ITCZ Attractors - Part 2

Through various specially designed numerical experiments with an aqua-planet general circulation model and theoretical arguments. Chao showed the existence of multiple quasi-equilibria of the intertropical convergence zone (ITCZ). He also showed that monsoon onset could be interpreted as an abrupt transition between the quasi-equilibria of the ITCZ. He further showed that the origin of these quasi-equilibria is related to two different types of attraction pulling the ITCZ in opposite directions. One type of attraction on the ITCZ is due to earth's rotation, which pulls the ITCZ toward the equator or two equatorial latitudes symmetric with respect to the equator depending on the choice of convection scheme, and the other due to the peak of the sea surface temperature (SST, which is given in the experiments a Gaussian profile in latitude and is uniform in longitude), which pulls the ITCZ toward a latitude just poleward of the SST peak. The strength of the attraction due to the earth's rotation has a highly nonlinear dependence on the latitude and that due to the SST peak has a linear (at least in a relative sense) dependence on the latitude.

Chao, Winston C.↗

Band Energy Dependence of Defect Formation in the Topological Semimetal Cd3As2

Cadmium Arsenide (Cd3As2) is a prototypical Dirac semimetal that manifests topological properties in a 3D bulk material. In defect-free Cd3As2, the Fermi level lies at a minimum in the density of states at the Dirac point, but experimentally it forms with excess electron carriers and an elevated EF, thereby masking the topological features. To computationally study the self-doping of Cd3As2, we combine density functional theory (DFT) calculations for defect formation energies with quasi-particle self-consistent GW (QSGW) electronic structure calculations. We demonstrate an innate dependence of the point defect formation energies on carrier concentrations and use the QSGW calculated density of states to extrapolate formation energies to arbitrary electron concentrations. This approach allows the quantitative modeling of thermodynamic defect equilibria in topological semimetals and is used to predict how Cd3As2 growth conditions affect the position of EF relative to the Dirac point.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗

Activity coefficient acquisition with thermodynamics–informed active learning for phase diagram construction

This work explores the use of thermodynamics-informed Gaussian processes (GPs) and active learning (AL) to model activity coefficients and construct phase diagrams. Relying on synthetic data generated from an excess Gibbs energy model, GPs were found to accurately describe the activity coefficients of several binary mixtures across large composition and temperature ranges. Moreover, GPs could estimate their own uncertainty and identify composition/temperature regions where activity coefficient data provide the most information to the models. This was leveraged to build AL algorithms targeted at modeling phase equilibria. In many cases, a single active-learning-acquired data point was sufficient to describe the phase diagrams studied. Lastly, the ability of AL to greatly reduce the amount of data needed to obtain accurate models was further verified on experimental case studies, namely individual ion activity coefficients, the solid–liquid and vapor–liquid equilibrium of deep eutectic solvents, and phase equilibria in ternary mixtures.

25 ENERGY STORAGE↗