Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linearity preservation”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Anti-symmetric and positivity preserving formulation of a spectral method for Vlasov-Poisson equations

We analyze the anti-symmetric properties of a spectral discretization for the one-dimensional Vlasov-Poisson equations. The discretization is based on a spectral expansion in velocity with the symmetrically weighted Hermite basis functions, central finite differencing in space, and an implicit Runge Kutta integrator in time. The proposed discretization preserves the anti-symmetric structure of the advection operator in the Vlasov equation, resulting in a stable numerical method. We apply such discretization to two formulations: the canonical Vlasov-Poisson equations and their continuously transformed square-root representation. The latter preserves the positivity of the particle distribution function. We derive analytically the conservation properties of both formulations, including particle number, momentum, and energy, which are verified numerically on the following benchmark problems: manufactured solution, linear and nonlinear Landau damping, two-stream instability, bump-on-tail instability, and ion-acoustic wave.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A Trajectory-Driven Algorithm for Differentiating SRB Measures on Unstable Manifolds

Sinai-Ruelle-Bowen (SRB) measures are limiting stationary distributions describing the statistical behavior of chaotic dynamical systems. Directional derivatives of SRB measure densities conditioned on unstable manifolds are critical in the sensitivity analysis of hyperbolic chaos. These derivatives, known as the SRB density gradients, are by-products of the regularization of Lebesgue integrals appearing in the original linear response expression. In this paper, we propose a novel trajectory- driven algorithm for computing the SRB density gradient defined for systems with high-dimensional unstable manifolds. We apply the concept of measure preservation together with the chain rule on smooth manifolds. Due to the recursive one-step nature of our derivations, the proposed procedure is memory-efficient and can be naturally integrated with existing Monte Carlo schemes widely used in computational chaotic dynamics. Here, we numerically show the exponential convergence of our scheme, analyze the computational cost, and present its use in the context of Monte Carlo integration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Structure and Fc-Effector Function of Rhesusized Variants of Human Anti-HIV-1 IgG1s

Passive transfer of monoclonal antibodies (mAbs) of human origin into Non-Human Primates (NHPs), especially those which function predominantly by a Fc-effector mechanism, requires an a priori preparation step, in which the human mAb is reengineered to an equivalent NHP IgG subclass. This can be achieved by changing both the Fc and Fab sequence while simultaneously maintaining the epitope specificity of the parent antibody. This Ab reengineering process, referred to as rhesusization, can be challenging because the simple grafting of the complementarity determining regions (CDRs) into an NHP IgG subclass may impact the functionality of the mAb. Here we describe the successful rhesusization of a set of human mAbs targeting HIV-1 envelope (Env) epitopes involved in potent Fc-effector function against the virus. This set includes a mAb targeting a linear gp120 V1V2 epitope isolated from a RV144 vaccinee, a gp120 conformational epitope within the Cluster A region isolated from a RV305 vaccinated individual, and a linear gp41 epitope within the immunodominant Cys-loop region commonly targeted by most HIV-1 infected individuals. Structural analyses confirm that the rhesusized variants bind their respective Env antigens with almost identical specificity preserving epitope footprints and most antigen-Fab atomic contacts with constant regions folded as in control RM IgG1s. In addition, functional analyses confirm preservation of the Fc effector function of the rhesusized mAbs including the ability to mediate Antibody Dependent Cell-mediated Cytotoxicity (ADCC) and antibody dependent cellular phagocytosis by monocytes (ADCP) and neutrophils (ADNP) with potencies comparable to native macaque antibodies of similar specificity. While the antibodies chosen here are relevant for the examination of the correlates of protection in HIV-1 vaccine trials, the methods used are generally applicable to antibodies for other purposes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

hi_class background evolution, initial conditions and approximation schemes

Cosmological datasets have great potential to elucidate the nature of dark energy and test gravity on the largest scales available to observation. Theoretical predictions can be computed with hi_class (www.hiclass-code.net), an accurate, fast and flexible code for linear cosmology, incorporating a wide range of dark energy theories and modifications to general relativity. We introduce three new functionalities into hi_class: (1) Support for models based on covariant Lagrangians, including a constraint-preserving integration scheme for the background evolution and a series of worked-out examples: Galileon, nKGB, quintessence (monomial, tracker) and Brans-Dicke. (2) Consistent initial conditions for the scalar-field perturbations in the deep radiation era, identifying the conditions under which modified-gravity isocurvature perturbations may grow faster than adiabatic modes leading to a loss of predictivity. (3) An automated quasistatic approximation scheme allowing order-of-magnitude improvement in computing performance without sacrificing accuracy for wide classes of models. These enhancements bring the treatment of dark energy and modified gravity models to the level of detail comparable to software tools restricted to standard $\Lambda$CDM cosmologies. Finally, the hi_class code is publicly available (https://github.com/miguelzuma/hi_class_public), ready to explore current data and prepare for next-generation experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Note on the Probability of Initiation Problem

The purpose of this note is to draw attention to Pál’s derivation of the Pál Bell equation which is used by LLNL and LANL to model Probability of Initiation (POI) problems. Both laboratories credit the equation to Bell and Lee. Although the equation formulated by Bell and Lee and the equation derived by Pál appear to be the same, the derivation by Bell and Lee is flawed. Scattering, which was absent in Bell’s original formulation of a POI problem, was later introduced by Bell and Lee by simply replacing a fission term with a scattering term. Such a replacement, however, leads to an inaccurate average number v̅ of neutrons that is produced by a fission event. By approaching a POI problem from probabilistic point of view, Pál formulates the problem as a branching process that determines v̅ correctly. There are six reasons for writing this report. The first is to clarify the failing in Bell and Lee’s derivation. The second is to establish the origin of the equation used by the laboratories. The third is to elucidate the fixed point technique formulated by Mokhtar-Kharroubi for solving the Pál-Bell equation. The fourth is to derive the α eigenvalue that is associated with a time dependent Pál-Bell equation. The fifth is to show in Appendix A that the non-linear fission operator of the Pál-Bell equation is descended from the operator of the Galton-Watson process of stochastic theory. The sixth is to describe in Appendix B a second order positivity preserving technique that accelerates the convergence of the first order fixed point method developed by Mokhtar-Kharroubi.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Linear stability analysis via simulated annealing and accelerated relaxation

Simulated annealing (SA) is a kind of relaxation method for finding equilibria of Hamiltonian systems. A set of evolution equations is solved with SA, which is derived from the original Hamiltonian system so that the energy of the system changes monotonically while preserving Casimir invariants inherent to noncanonical Hamiltonian systems. The energy extremum reached by SA is an equilibrium. Since SA searches for an energy extremum, it can also be used for stability analysis when initiated from a state where a perturbation is added to an equilibrium. The procedure of the stability analysis is explained, and some examples are shown. Because the time evolution is computationally time consuming, efficient relaxation is necessary for SA to be practically useful. An acceleration method is developed by introducing time dependence in the symmetric kernel used in the double bracket, which is part of the SA formulation described here. An explicit formulation for low-beta reduced magnetohydrodynamics (MHD) in cylindrical geometry is presented. In conclusion, since SA for low-beta reduced MHD has two advection fields that relax, it is important to balance the orders of magnitude of these advection fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Dynamic, symmetry-preserving, and hardware-adaptable circuits for quantum computing many-body states and correlators of the Anderson impurity model

We present a hardware-reconfigurable ansatz on N q -qubits for the variational preparation of many-body states of the Anderson impurity model (AIM) with N imp + N bath = N q /2 sites, which conserves total charge and spin z component within each variational search subspace. The many-body ground state of the AIM is determined as the minimum over all minima of O(N$^2_ q$) distinct charge-spin sectors. Hamiltonian expectation values are shown to require ω(N q ) < N meas. $\leqslant$ O(N imp N bath ) symmetry-preserving, parallelizable measurement circuits, each amenable to postselection. To obtain the one-particle impurity Green’s function we show how initial Krylov vectors can be computed via midcircuit measurement and how Lanczos iterations can be computed using the symmetry-preserving ansatz. For a single-impurity Anderson model with a number of bath sites increasing from one to seven, we show using numerical emulation that the ease of variational ground-state preparation is suggestive of linear scaling in circuit depth and subquartic scaling in optimizer complexity. We therefore expect that, combined with time-dependent methods for Green’s function computation, our ansatz provides a useful tool to account for electronic correlations on early fault-tolerant processors. Finally, with a view towards computing real materials properties of interest like magnetic susceptibilities and electron-hole propagators, we provide a straightforward method to compute many-body, time-dependent correlation functions using a combination of time evolution, midcircuit measurement-conditioned operations, and the Hadamard test.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Lorentz symmetry in ghost-free massive gravity

The role of Lorentz symmetry in ghost-free massive gravity is studied, emphasizing features emerging in approximately Minkowski spacetime. The static extrema and saddle points of the potential are determined and their Lorentz properties identified. Solutions preserving Lorentz invariance and ones breaking four of the six Lorentz generators are constructed. Locally, globally, and absolutely stable Lorentz-invariant extrema are found to exist for certain parameter ranges of the potential. Gravitational waves in the linearized theory are investigated. Deviations of the fiducial metric from the Minkowski metric are shown to lead to pentarefringence of the five wave polarizations, which can include superluminal modes and subluminal modes with negative energies in certain observer frames. The Newton limit of ghost-free massive gravity is explored. The propagator is constructed and used to obtain the gravitational potential energy between two point masses. Here, the result extends the Fierz-Pauli limit to include corrections generically breaking both rotation and boost invariance.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Simulating large one-dimensional neutral-atom quantum systems

While abstract models of quantum computation assume a closed system of two-level states, practical quantum devices inevitably couple to the environment in some way, creating sources of noise. Understanding the tolerance to noise of specific quantum algorithms run on specific devices is important for determining the feasibility of quantum computing in the current noisy intermediate-scale quantum era. Of particular interest is understanding the noise sensitivity of these devices as more qubits are added to the system. Classical simulations are a useful tool to understand the effects of this noise, but direct classical simulations of open quantum systems are burdened by an exponentially growing cost in the number of qubits and a large local Hilbert space dimension. For onedimensional, shallow circuits, using tensor networks can replace this exponential cost with a linear one and simulate far wider systems than what would normally be available. In this paper, we describe a tensor network simulation of a neutral atom quantum system under the presence of noise, while introducing a purity-preserving truncation technique that compromises between the simplicity of the matrix product state and the positivity of the matrix product density operator. We apply this simulation to a near-optimized iteration of the quantum approximate optimization algorithm on a transverse field Ising model in order to investigate the influence of large system sizes on the performance of the algorithm. We find that while circuits with a large number of qubits fail more often under noise that depletes the qubit population, their outputs on a successful measurement are just as robust under Rydberg atom dissipation or qubit dephasing as smaller systems. However, such circuits might not perform as well under coherent multiqubit errors such as Rydberg atom crosstalk. We also find that the optimized parameters are especially robust to noise, suggesting that a noisier quantum system can be used to find the optimal parameters before switching to a cleaner system for measurements of observables.

Allen, James↗
DOE OSTI2022

Nonnegative canonical tensor decomposition with linear constraints: nnCANDELINC

Abstract There is an emerging interest for tensor factorization applications in big‐data analytics and machine learning. To speed up the factorization of extra‐large datasets, organized in multidimensional arrays (also known as tensors), easy to compute compression‐based tensor representations, such as, Tucker and tensor train formats, are used to approximate the initial large‐tensor. Further, tensor factorization is used to extract latent features that can facilitate discoveries of new mechanisms and signatures hidden in the data, where the explainability of the latent features is of principal importance. Nonnegative tensor factorization extracts latent features that are naturally sparse and parts of the data, which makes them easily interpretable. However, to take into account available domain knowledge and subject matter expertise, often additional constraints need to be imposed, which lead us to canonical decomposition with linear constraints (CANDELINC), a canonical polyadic decomposition with rank deficient factors. In CANDELINC, Tucker compression is used as a preprocessing step, which lead to a larger residual error but to more explainable latent features. Here, we propose a nonnegative CANDELINC (nnCANDELINC) accomplished via a specific nonnegative Tucker decomposition; we refer to as minimal or canonical nonnegative Tucker. We derive several results required to understand the specificity of nnCANDELINC, focusing on the difficulties of preserving the nonnegative rank of a tensor to its Tucker core and comparing the real valued to nonnegative case. Finally, we demonstrate nnCANDELINC performance on synthetic and real‐world examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Experimental Limit on Nonlinear State-Dependent Terms in Quantum Theory

Linear time evolution is one of the fundamental postulates of quantum theory. Past theoretical attempts to introduce nonlinearity into quantum evolution have violated causality. However, a recent theory has introduced nonlinear state-dependent terms in quantum field theory, preserving causality [D. E. Kaplan and S. Rajendran, Phys. Rev. D 105, 055002 (2022)]. We report the results of an experiment that searches for such terms. Our approach, inspired by the Everett many-worlds interpretation of quantum theory, correlates a binary macroscopic classical voltage with the outcome of a projective measurement of a quantum bit, prepared in a coherent superposition state. Measurement results are recorded in a bit string, which is used to control a voltage switch. Presence of a nonzero voltage reading in cases of no applied voltage is the experimental signature of a nonlinear state-dependent shift of the electromagnetic field operator. We implement blinded measurement and data analysis with three control bit strings. Control of systematic effects is realized by producing one of the control bit strings with a classical random-bit generator. The other two bit strings are generated by measurements performed on a superconducting qubit in an IBM Quantum processor and on a N 15 nuclear spin in a nitrogen-vacancy center in diamond. Our measurements find no evidence for electromagnetic quantum state-dependent nonlinearity. We set a bound on the parameter that quantifies this nonlinearity | ε γ | < 4.7 × 10 - 11 , at 90% confidence level.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Investigation of modified uni-traveling carrier photodiode for cryogenic microwave photonic links

Quantum devices present the potential for unparalleled computing and communications capabilities; however, the cryogenic temperatures required to successfully control and read out many qubit platforms can prove to be very challenging to scale. Recently, there has emerged an interest in using microwave photonics to deliver control signals down to ultracold stages via optical fiber, thereby reducing thermal load and facilitating dense wavelength multiplexing. Photodetectors can then convert this optical energy to electrical signals for qubit control. The fidelity of the quantum operations of interest therefore depend heavily upon the characteristics of the photodiode, yet experimental demonstrations of fiber-coupled photodetection systems at low temperatures are relatively few in number, leaving important open questions regarding how specific detectors may perform in real-world cryogenic settings. In this work, we examine a highly linear modified uni-traveling carrier photodiode (MUTC-PD) under C-band illumination (1530–1565 nm) at three temperature regimes (300 K, 80 K, and ∼4 K) and multiple bias conditions. Our findings of reduced responsivity but preserved bandwidth are consistent with previous studies, while our saturation tests suggest a variety of potential applications for MUTC-PDs in cryogenic microwave photonics with and without electrical bias. Overall, our results should provide a valuable foundation for the continued and expanding use of this detector technology in quantum information processing.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Tuning dispersity of linear polymers and polymeric brushes grown from nanoparticles by atom transfer radical polymerization

Molecular weight distribution imposes considerable influence on the properties of polymers, making it an important parameter, impacting morphology and structural behavior of polymeric materials. Atom transfer radical polymerization (ATRP) has established itself as a powerful tool to prepare polymers with predetermined molecular weight, preserved chain-end functionality, and low dispersity. More recently, ATRP has also been shown to provide a means to deliberately broaden molecular weight distributions, and, via retaining living chain-ends, to enable the formation of block copolymers with designed block dispersity, featuring new microstructures and potentially attractive properties. Here, similar methodologies have been developed to facilitate tuning of the dispersity of polymeric brushes grown from nanoparticles thus resulting in hybrid materials with enhanced fracture toughness and high inorganic content. Recent advances have given access to brush architectures comprised of uni- and bimodal block copolymers with unique morphologies along with interesting mechanical, thermal, and optical properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Fully Homomorphic Encryption

This code implements a Fully Homomorphic Encryption (FHE) system, enabling secure computation on encrypted data without requiring decryption. It supports encryption, decryption, and homomorphic operations like matrix multiplication and addition. This code is adaptable for integrating FHE into linear-time invariant (LTI) systems, including digital control and filtering. With proper configuration from subject matter expertise, encrypted system parameters and signals can be manipulated to perform tasks like state updates, output calculations, and convolution in the encrypted domain. By preserving the structure of LTI systems while ensuring privacy, the framework facilitates secure applications in areas such as autonomous systems, signal processing, and industrial automation. The code initializes the encryption system using parameters provided in the env dictionary. These parameters include the ciphertext modulus, key dimension, plaintext fixed-point scaling factor, and noise bound. During initialization, a secret key is generated, which is essential for encrypting and decrypting data securely. The modular design allows users to tailor these parameters to specific use cases or security requirements. The code implements multiple cryptographic schemes. The learning with errors (LWE) encryption method encodes cleartext message to their plaintext fixed-point representation then encrypted into ciphertext space with additive noise. This noise ensures the security of the scheme, relying on the computational hardness of the LWE problem. The code also includes the Gentry-Sahai-Waters (GSW) scheme based off the LWE problem. Homomorphic matrix multiplication is performed between the LWE and GSW to encrypted data. This is achieved using a decomposition function on the LWE ciphertext during the multiplication operation. For higher-dimensional data, the code includes a method to encrypt entire matrices (GSWMat) using GSW encryption. These encrypted matrices can then be used for homomorphic matrix multiplications (MatMult). The decryption function uses the secret key to recover the original plaintext, removing the added noise and scaling that was originally applied during encryption.

Lois, Roberts [Idaho National Laboratory (INL), Id↗
DOE OSTI2025

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

We propose a second-order temporally implicit, fourth-order-accurate spatial discretization scheme for the strongly anisotropic heat transport equation characteristic of hot, fusion-grade plasmas. Following Du Toit et al. (2018), the scheme transforms mixed-derivative diffusion fluxes (which are responsible for the lack of a discrete maximum principle) into nonlinear advective fluxes, amenable to nonlinear-solver-friendly monotonicity-preserving limiters. The scheme enables accurate multi-dimensional heat transport simulations with up to seven orders of magnitude of heat-transport-coefficient anisotropies with low cross-field numerical error pollution and excellent algorithmic performance, with the number of linear iterations scaling very weakly with grid resolution and grid anisotropy, and scaling with the square-root of the implicit timestep. We propose a multigrid preconditioning strategy based on a lower-order approximation that renders the scheme efficient and scalable under grid refinement. Several numerical tests are presented that display the expected spatial convergence rates and strong algorithmic performance, including fully nonlinear magnetohydrodynamics simulations of kink instabilities in a Bennett pinch in 2D helical geometry and of ITER in 3D toroidal geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Dielectric-Engineered Monolayer MoS 2 Memtransistors for Brain-Inspired Computing with High Recognition Accuracy

Two-dimensional transition metal dichalcogenides (2D-TMDs)-based memtransistors have emerged as promising candidates for neuromorphic hardware due to their exceptional ability to emulate synaptic behavior. However, many existing 2D-TMDs memtransistors rely on polycrystalline channels with grain boundaries or defects introduced through postgrowth treatments, raising concerns about material integrity and the preservation of intrinsic properties. Here, in this work, we demonstrate a monocrystalline monolayer MoS 2 memtransistor fabricated on a silicon nitride (SiN X ) substrate, achieving a large resistive switching ratio of 10 4 , a dynamic range exceeding 90, along with highly linear and symmetric weight updates, minimal cycle-to-cycle variability, and low device-to-device variability. These attributes are critical for enabling high-performance neuromorphic hardware. Based on experimental data, we further show that these artificial synapses enable a recognition accuracy of more than 97% on the MNIST handwritten digits data set. Our findings present a straightforward approach to realizing 2D-TMDs memtransistors through dielectric engineering, offering a promising platform for next-generation neuromorphic computing systems.

2D TMDs↗
DOE OSTI2025

IM3 Projected U.S. Western Interconnection Grid Stress Dataset

This dataset provides projected grid stress and reliability results (including all model inputs and outputs from GO WEST and TEP) for Integrated Multisector, Multiscale Modeling (IM3) Phase 2 simulations across eight different scenarios for the U.S. Western Interconnection through 2055. The scenarios include combinations of two Shared Socioeconomic Pathways (SSP3 and SSP5) with four high-resolution climate projections specific to the United States from a set of Thermodynamic Global Warming (TGW) simulations. These climate projections include "hotter" and "cooler" variants for two Representative Concentration Pathways (RCP4.5 and RCP8.5). The resulting eight simulations are: rcp45cooler_ssp3 rcp45cooler_ssp5 rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85cooler_ssp3 rcp85cooler_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 GO WEST is an open-source power grid modeling framework for the U.S. Western Interconnection, which allows users to tailor the model depending on their research study and science questions. It covers 28 balancing authorities (BAs) and 12 states in U.S. Western Interconnection. GO WEST allows users to select different number of nodes and come up with a simplified network by utilizing 10,000 nodal topology of the U.S. Western Interconnection (ACTIVSg10k). Users can select different number of nodes, mathematical formulations (linear programming vs. mixed-integer linear programming), transmission line limit scaling factors, and hurdle rate scaling factors. GO WEST offers a unit commitment and economic dispatch (UC/ED) module to simulate grid operations on an hourly scale. In this sense, users can calibrate and validate their model versions by comparing model outputs to historical datasets. TEP is an open-source transmission capacity expansion model, built on the GO WEST framework. It utilizes linear programming to optimize transmission capacity addition investment on existing lines within the GO WEST framework. The TEP model only increases the thermal capacity of existing transmission lines and does not add new lines to the system, which leaves the topology preserved. In order to use TEP model, users need to create scenarios with the GO WEST framework. Please refer to README file for a detailed description of the dataset including individual files and references.

Capacity Expansion Model↗
DOE OSTI2023

On the relationship between the multi-region relaxed variational principle and resistive inner-layer theory

We show that the variational energy principle of the multi-region relaxed magnetohydrodynamic (MRxMHD) model can be used to predict finite-pressure linear tearing instabilities. In this model, the plasma volume is sliced into sub-volumes separated by 'ideal interfaces', and in each volume the magnetic field relaxes to a Taylor state, where the pressure gradient $\nabla p = 0$. The MRxMHD model is implemented in the Stepped-Pressure Equilibrium Code (SPEC) so that the equilibrium solution in each region is computed while preserving the force balance across the interfaces. As SPEC computes the Hessian matrix (a discretized stability matrix), the stability of an MRxMHD equilibrium can also be computed with SPEC. In this article, using SPEC, we investigate the effect of local pressure gradients and the $\nabla p = 0$ in the vicinity of the resonant surface of a tearing mode. For low-beta plasma, we have been able to illustrate a relationship between the resistive singular-layer theory (Coppi et al 1966 Nucl. Fusion6 101; Glasser et al 1975 Phys. Fluids18 875–88) and the MRxMHD model. Within the singular layer, the volume-averaged magnetic helicity and the flux-averaged toroidal flux are shown to be the invariants for the linear tearing modes in SPEC simulations. Our technique to compute MRxMHD stability is first tested numerically in a cylindrical tokamak and its application in toroidal geometry is demonstrated. We demonstrate an agreement between the stability boundary obtained with SPEC simulation and the resistive inner-layer theories.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗