Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability properties”

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 181 records · Page 10

Kinetic-ballooning-bifurcation in tokamak pedestals across shaping and aspect-ratio

We use a new gyrokinetic threshold model to predict a bifurcation in tokamak pedestal width-height scalings that depends strongly on plasma shaping and aspect-ratio. The bifurcation arises from the first and second stability properties of kinetic-ballooning-modes that yields wide and narrow pedestal branches, expanding the space of accessible pedestal widths and heights. The wide branch offers potential for edge-localized-mode-free pedestals with high core pressure. For negative triangularity, low-aspect-ratio configurations are predicted to give steeper pedestals than conventional-aspect-ratio. Both wide and narrow branches have been attained in tokamak experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A stable corridor for toroidal plasma compression

In this study, a toroidal plasma compressed by a collapsing flux conserver is analyzed to reveal stable scenarios of operation to high compression ratios. The resistive and ideal MHD stability is calculated in full toroidal geometry, using the asymptotic matching method in realistic conditions, and comparing with nonlinear simulations. The near edge current profile, controlled by toroidal field ramping during compression, is shown to be critical to stability due to coupling between poloidal components of the least stable mode. The extension of a length of shaft on axis is also found to be critical at high compression, as the resulting good curvature region in magnetic field stabilizes pressure driven modes that would otherwise be unstable. This work extends from previous studies, which initially showed the existence of a stable scenario, to include findings of more extensive stable zones, detailed effects of geometry, and nonlinear simulations of the instabilities. The nonlinear simulations of the compression are consistent with the linear analyses, confirming both the conservation and stability properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Finite magnetic well effects on resistive and drift-resistive ballooning modes in a shaped tokamak

The impact of plasma shaping through magnetic well modifications on the stability of resistive ballooning modes (RBMs) in tokamaks is analysed, also including finite diamagnetic flows. Furthermore, various limiting cases of the dispersion relation, obtained by matching the ballooning equation across the ideal and resistive layers, are analysed. It is found that stability is generally improved by the combination of vertical elongation and positive triangularity, although, in some cases, the growth rate of the unstable mode can be enhanced by these effects. Usually, vertically elongated plasmas with no triangularity are prone to exhibit worse stability properties. A value for the critical β above which RBMs are driven unstable is identified, and a connection with type-III ELM activity is established.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stochastic fluctuations in relativistic fluids: Causality, stability, and the information current

We develop a general formalism for introducing stochastic fluctuations around thermodynamic equilibrium which takes into account, for the first time, recent developments in the causality and stability properties of relativistic hydrodynamic theories. The method is valid for any covariantly stable theory of relativistic viscous fluid dynamics derived from a covariant maximum entropy principle. We illustrate the formalism with some applications, showing how it could be used to consistently introduce fluctuations in a model of relativistic heat diffusion and in conformally invariant Israel-Stewart theory in a general hydrodynamic frame. Furthermore, the latter example is used to study the hydrodynamic frame dependence of the symmetric two-point function of fluctuations of the energy-momentum tensor.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Composite Fermi Liquid at Zero Magnetic Field in Twisted MoTe 2

Here, the pursuit of exotic phases of matter outside of the extreme conditions of a quantizing magnetic field is a long-standing quest of solid state physics. Recent experiments have observed spontaneous valley polarization and fractional Chern insulators in zero magnetic field in twisted bilayers of MoTe 2 , at partial filling of the topological valence band (ν =-2/3 and -3/5). We study the topological valence band at half filling, using exact diagonalization and density matrix renormalization group calculations. We discover a composite Fermi liquid (CFL) phase even at zero magnetic field that covers a large portion of the phase diagram near twist angle ~3.6°. The CFL is a non-Fermi liquid phase with metallic behavior despite the absence of Landau quasiparticles. We discuss experimental implications including the competition between the CFL and a Fermi liquid, which can be tuned with a displacement field. The topological valence band has excellent quantum geometry over a wide range of twist angles and a small bandwidth that is, remarkably, reduced by interactions. These key properties stabilize the exotic zero field quantum Hall phases. Finally, we present an optical signature involving “extinguished” optical responses that detects Chern bands with ideal quantum geometry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part II: Nonlinearities and DAEs

Fully implicit Runge--Kutta (IRK) methods have many desirable accuracy and stability properties as time integration schemes, but high-order IRK methods are not commonly used in practice with large-scale numerical PDEs because of the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic framework for solving the nonlinear equations that arise from IRK methods (and discontinuous Galerkin discretizations in time) applied to nonlinear numerical PDEs, including PDEs with algebraic constraints. Several new linearizations of the nonlinear IRK equations are developed, offering faster and more robust convergence than the often-considered simplified Newton, as well as an effective preconditioner for the true Jacobian if exact Newton iterations are desired. Inverting these linearizations requires solving a set of block 2 x 2 systems. Under quite general assumptions, it is proven that the preconditioned 2 x 2 operator's condition number is bounded by a small constant close to one, independent of the spatial discretization, spatial mesh, and time step, and with only weak dependence on the number of stages or integration accuracy. Moreover, the new method is built using the same preconditioners needed for backward Euler-type time stepping schemes, so can be readily added to existing codes. The new methods are applied to several challenging fluid flow problems, including the compressible Euler and Navier--Stokes equations, and the vorticity-streamfunction formulation of the incompressible Euler and Navier--Stokes equations. Up to 10th-order accuracy is demonstrated using Gauss IRK, while in all cases fourth-order Gauss IRK requires roughly half the number of preconditioner applications as required by standard Singly diagonally implicit Runge--Kutta methods.

97 MATHEMATICS AND COMPUTING↗

FAR3d

The FAR3d model calculates the linear and nonlinear stability properties of energetic particle driven Alfven instabilities for both tokamak and stellarator plasma confinement devices using gyro-landau closure methods. This is an important fundamental physics problem for existing fusion energy experiments and for future fusion reactors.

Varela, Jacobo↗

A Virtual Element Method for the Full MHD system

In this manuscript we present a novel discretization for the incompressible MHD system. Our approach follows the framework of the Virtual Element Method and offers two main advantages. The method can be implemented in unstructured meshes making it highly versatile and capable of handling a wide array of problems involving interfaces, free-boundaries or adaptive refinements on the mesh. The second advantage involves the divergence of the magnetic field, our approach guarantees that it remains solenoidal. We include a theoretical proof of the condition on the magnetic field as well as energy estimates and a well-posedness study. The latter sheds light as to the stability properties of the method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling of Macroscopic Dynamics in Three-Dimensional Plasma Configurations (Final Technical Report)

A numerical simulation-based study to investigate the stability properties of inherently three-dimensional plasma magnetic configurations is presented. These configurations, called stellarators and heliotrons, use external helical coils to provide the twisting magnetic field that is required to confine ionized plasma particles. In contrast to the tokamak configuration, which nominally has the property of axisymmetry, electrical current does not need to be induced in stellarator and heliotron plasma. Experiments have shown that these three-dimensional configurations are remarkably robust to macroscopic instability, even when driven to relatively high levels of plasma pressure. Our understanding of how the plasma is able to withstand macroscopic instability will be enhanced by the work proposed here. The effort will generalize existing numerical modeling capability, and it will apply the new capability to model macroscopic dynamics in these three-dimensional configurations. The development work has been applied to the NIMROD simulation code (https://nimrodteam.org) that was originally based on axisymmetric problem domains. It has been generalized to accommodate the three-dimensional shaping of the plasma volume in stellarators and heliotrons. Its use of equilibrium data has also been generalized, and the implementation allows different numerical representation of the magnetic field. Analysis has been performed to inform the selection of numerical methods, and a method to accelerate the iterative solution of the necessary algebraic systems has been developed and implemented. Verification of the new NIMSTELL code on benchmark problems has also been an important part of the effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NSTX-U National Research Program: White Paper in Response to Call from FESAC Sub-Committee

Both scientific and technical innovation is needed for the realization of an attractive engineering solution for a timely and cost-effective Pilot Plant, the design and construction of which is the overarching recommendation of the FESAC Long Range Plan, and the 2021 NASEM Pilot Plant reports, which underpin the Bold Decadal Vision. The two most significant plasma physics gaps to close for a Compact Pilot Plant (CPP) are core confinement improvement and heat flux mitigation, neither of which have been closed in an integrated fashion for any planned fusion power production device. High core confinement and stability are essential for producing majority self-driven plasmas in CPPs with reduced size and auxiliary heating power requirements, with an improvement in confinement being the major driver for cost reduction of a CPP. The National Spherical Tokamak Experiment - Upgrade (NSTX-U) is a unique low aspect ratio research facility that will address the fundamental challenge of developing the science and technology basis for a CPP design that integrates high core and edge confinement with the ability to mitigate very high incident heat fluxes. NSTX-U capabilities will enable the high performance, already achieved on NSTX, to extend into physics regimes much closer to those anticipated in Spherical Tokamak (ST)-based CPPs. These confinement and stability properties will be assessed by a full complement of diagnostics and analysis tools, which will also aid in the development of the underlying theory and predictive models needed for further optimization. Both conventional and transformative heat flux mitigation methods, such as liquid lithium plasma-facing components, will be developed and tested in-situ in NSTX-U at incident heat fluxes of ~100 MW/m 2 , and will inform plans and reduce risk for a subsequent major upgrade to the device to fully heated, high-Z wall and full liquid lithium divertor capability, a technology that potentially could then be implemented on any magnetic confinement device at any aspect ratio. NSTX-U research is fully complementary to programs performed on other STs, nationally and internationally. Furthermore, NSTX-U research has a direct connection to the private sector by informing design choices for future power production facilities being developed by these companies. The NSTX-U program will operate as a national User Facility, with collaborating researchers, engineers, and graduate students from 19 outside institutions, and open to participation and experiments led by researchers from both public and private entities. The research program will advance workforce development through training of young scientists, engineers, and technicians, and it will also serve for further diagnostic innovation, especially for high heat flux and high-Z wall environments, and implementation of advanced artificial intelligence (AI) for plasma and heat flux control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Functionalization of Cathode–Electrolyte Interface with Ionic Liquids for High-Performance Quasi-Solid-State Lithium–Sulfur Batteries: A Low-Sulfur Loading Study

We introduce a quasi-solid-state electrolyte lithium-sulfur (Li–S) battery (QSSEB) based on a novel Li-argyrodite solid-state electrolyte (SSE), Super P–Sulfur cathode, and Li-anode. The cathode was prepared using a water-based carboxymethyl cellulose (CMC) solution and styrene butadiene rubber (SBR) as the binder while Li6PS5F0.5Cl0.5 SSE was synthesized using a solvent-based process, via the introduction of LiF into the argyrodite crystal structure, which enhances both the ionic conductivity and interface-stabilizing properties of the SSE. Ionic liquids (IL) were prepared using lithium bis(trifluoromethyl sulfonyl)imide (LiTFSI) as the salt, with pre-mixed pyrrolidinium bis(trifluoromethyl sulfonyl)imide (PYR) as solvent and 1,3-dioxolane (DOL) as diluent, and they were used to wet the SSE–electrode interfaces. The effect of IL dilution, the co-solvent amount, the LiTFSI concentration, the C rate at which the batteries are tested and the effect of the introduction of SSE in the cathode, were systematically studied and optimized to develop a QSSEB with higher capacity retention and cyclability. Interfacial reactions occurring at the cathode–SSE interface during cycling were also investigated using electrochemical impedance spectroscopy, cyclic voltammetry, and X-ray photoelectron spectroscopy supported by ab initio molecular dynamics simulations. This work offers a new insight into the intimate interfacial contacts between the SSE and carbon–sulfur cathodes, which are critical for improving the electrochemical performance of quasi-solid-state lithium–sulfur batteries.

25 ENERGY STORAGE↗

Viscosity Analysis of Battery Electrode Slurry

We report the effects of component ratios and mixing time on electrode slurry viscosity. Three component quantities were varied: active material (graphite), conductive material (carbon black), and polymer binder (carboxymethyl cellulose, CMC). The slurries demonstrated shear-thinning behavior, and suspension properties stabilized after a relatively short mixing duration. However, micrographs of the slurries suggested their internal structures did not stabilize after the same mixing time. Increasing the content of polymer binder CMC caused the greatest viscosity increase compared to that of carbon black and graphite.

25 ENERGY STORAGE↗

The Kelvin–Helmholtz Instability at the Boundary of Relativistic Magnetized Jets

We study the linear stability of a planar interface separating two fluids in relative motion, focusing on conditions appropriate for the boundaries of relativistic jets. The jet is magnetically dominated, whereas the ambient wind is gas-pressure-dominated. We derive the most general form of the dispersion relation and provide an analytical approximation of its solution for an ambient sound speed much smaller than the jet Alfvén speed v A , as appropriate for realistic systems. The stability properties are chiefly determined by the angle ψ between the wavevector and the jet magnetic field. For ψ = π/2, magnetic tension plays no role, and our solution resembles the one of a gas-pressure-dominated jet. Here, only sub-Alfvénic jets are unstable ($0\lt {M}_{e}\equiv (v/{v}_{{\rm{A}}})\cos \theta \lt 1$ , where v is the shear velocity and θ the angle between the velocity and the wavevector). For ψ = 0, the free energy in the velocity shear needs to overcome the magnetic tension, and only super-Alfvénic jets are unstable ($1\lt {M}_{e}\lt \sqrt{{(1+{{\rm{\Gamma }}}_{w}^{2})/[1+({v}_{{\rm{A}}}/c)}^{2}{{\rm{\Gamma }}}_{w}^{2}]}$, with Γ w the wind adiabatic index). Our results have important implications for the propagation and emission of relativistic magnetized jets.

79 ASTRONOMY AND ASTROPHYSICS↗

Research in the Restricted Problems of Three and Four Bodies Final Scientific Report

Seven studies have been conducted on research in the existence and nature of solutions of the restricted problems of three and four bodies. The details and results of five of these research investigations have already been published, and the latest two studies will be published shortly. A complete bibliography of publications is included in this report. This research has been primarily qualitative and has yielded new information on the behavior of trajectories near the libration points in the Earth-Moon-Sun and Sun-Jupiter-Saturn systems, and on the existence of periodic trajectories about the libration points of the circular and elliptical restricted four-body models. We have also implemented Birkhoff's normalization process for conservative and nonconservative Hamiltonian systems with equilibrium points. This makes available a technique for analyzing stability properties of certain nonlinear dynamical systems, and we have applied this technique to the circular and elliptical restricted three-body models. A related study was also conducted to determine the feasibility of using cislunar periodic trajectories for various space missions. Preliminary results suggest that this concept is attractive for space flight safety operations in cislunar space. Results of this research will be of interest to mathematicians, particularly those working in ordinary differential equations, dynamical systems and celestial mechanics; to astronomers; and to space guidance and mission analysts.

Richards, Paul B.↗

Error compensation for hybrid-computer solution of linear differential equations

Z-transform technique compensates for digital transport delay and digital-to-analog hold. Method determines best values for compensation constants in multi-step and Taylor series projections. Technique also provides hybrid-calculation error compared to continuous exact solution, plus system stability properties.

Kemp, N. H.↗

Research in plasma physics

Three aspects of barium ion cloud dynamics are discussed. First, the effect of the ratio of ion cloud conductivity to background ionospheric conductivity on the motion of barium ion clouds is investigated and compared with observations of barium ion clouds. This study led to the suggestion that the conjugate ionosphere participates in the dynamics of barium ion clouds. Second, analytic work on the deformation of ion clouds is presented. Third, a linearized stability theory was extended to include the effect of the finite extent of an ion cloud, as well as the effect of the ratio of ion cloud to ionospheric conductivities. The stability properties of a plasma with contra-streaming ion beams parallel to a magnetic field are investigated. The results are interpreted in terms of parameters appropriate for collisionless shock waves. It is found that this particular instability can be operative only if the up-stream Alfven Mach number exceeds 5.5.

Source record↗