Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “test matrix”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Mineralized sclerites in the gorgonian coral Leptogorgia chilensis as a natural jamming system

The soft corals (Cnidaria, Octocorallia), a diverse group of colonial marine invertebrates, can reversibly tune their body stiffness in response to external stimuli. This capability is attributed to their dynamic skeletal systems, which consist of thousands of mineralized skeletal elements, called sclerites, embedded within a gel-like matrix that swells/deswells and unjams/jams the sclerites, thus modulating skeletal stiffness. While sclerite morphology is widely used for species identification, its role in the mechanical performance of a soft coral’s skeletal system is largely unknown. Here, we investigated structure-jamming relationships in sclerite-based skeletal architectures using the red gorgonian octocoral Leptogorgia chilensis as a model system. The sclerites of L. chilensis exhibit a shaft-like geometry with two axial branches and two sets of triradiate side branches, which are aligned with the crystallographic symmetry of the constituent magnesium-containing calcite. By combining multiscale three-dimensional (3D) structural characterization, parametric geometrical modeling, 3D printing, mechanical testing, and discrete element simulations, we demonstrate how sclerite geometry achieves a balanced jamming performance in terms of stiffness, weight, strength, and fracture resistance in comparison to alternative geometries parametrically modified from the native sclerites (e.g., changes in the length and number of side branches). Here, we also found that these performance metrics are achieved through the effective interlocking among side and axial branches, which is further enhanced by the fractal-like microscopic spikes on the branch tips. The findings in this natural jamming system offer insights for designing synthetic mechanotunable material architectures for a wide range of applications, from soft robotics to mechanical dampeners.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Exploring Hilbert space on a budget: Novel benchmark set and performance metric for testing electronic structure methods in the regime of strong correlation

This work explores the ability of classical electronic structure methods to efficiently represent (compress) the information content of full configuration interaction (FCI) wave functions. We introduce a benchmark set of four hydrogen model systems of different dimensionalities and distinctive electronic structures: a 1D chain, a 1D ring, a 2D triangular lattice, and a 3D close-packed pyramid. To assess the ability of a computational method to produce accurate and compact wave functions, we introduce the accuracy volume, a metric that measures the number of variational parameters necessary to achieve a target energy error. Using this metric and the hydrogen models, we examine the performance of three classical deterministic methods: (i) selected configuration interaction (sCI) realized both via an a posteriori (ap-sCI) and variational selection of the most important determinants, (ii) an a posteriori singular value decomposition (SVD) of the FCI tensor (SVD-FCI), and (iii) the matrix product state representation obtained via the density matrix renormalization group (DMRG). We find that the DMRG generally gives the most efficient wave function representation for all systems, particularly in the 1D chain with a localized basis. For the 2D and 3D systems, all methods (except DMRG) perform best with a delocalized basis, and the efficiency of sCI and SVD-FCI is closer to that of DMRG. For larger analogs of the models, the DMRG consistently requires the fewest parameters but still scales exponentially in 2D and 3D systems, and the performance of SVD-FCI is essentially equivalent to that of ap-sCI.

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

Reactive Flows in Porous Media: Challenges in Theoretical and Numerical Methods

We review theoretical and computational research, primarily from the past 10 years, addressing the flow of reactive fluids in porous media. The focus is on systems where chemical reactions at the solid–fluid interface cause dissolution of the surrounding porous matrix, creating nonlinear feedback mechanisms that can often lead to greatly enhanced permeability. Here, we discuss insights into the evolution of geological forms that can be inferred from these feedback mechanisms, as well as some geotechnical applications such as enhanced oil recovery, hydraulic fracturing, and carbon sequestration. Until recently, most practical applications of reactive transport have been based on Darcy-scale modeling, where averaged equations for the flow and reactant transport are solved. We summarize the successes and limitations of volume averaging, which leads to Darcy-scale equations, as an introduction to pore-scale modeling. Pore-scale modeling is computationally intensive but offers new insights as well as tests of averaging theories and pore-network models. We include recent research devoted to validation of pore-scale simulations, particularly the use of visual observations from microfluidic experiments.

wormhole formation↗
DOE OSTI2025

Rapid Bayesian High Entropy Alloy Designs Fabricated via Wire Arc Additive Manufacturing

Purpose: This project seeks to demonstrate a new high-throughput (rapid) alloy design technique applied to creating new high entropy alloys (HEAs) for extreme environments. High entropy alloys shift the design paradigm from being focused on a single principal element (e.g. nickel-based alloys) to target alloys that include high atomic fractions (X >10%) of multiple elements. These HEA materials can exhibit sluggish diffusion and enhanced corrosion resistance, ideal for potential applications in advanced ultra supercritical (A-USC) steam cycles for power generation. Scope: The addition of multiple elements in high atomic fractions creates an enormous design space that cannot easily be investigated by traditional material design strategies such as designed of experiments (DOE). This project utilizes a Bayesian machine learning algorithm that has been modified to work with calculation of phase diagrams (CALPHAD) software. This Bayesian algorithm reduces manual inputs and increase the likelihood of achieving an optimal solution. Compositional inputs to this algorithm will be assessed using existing material property models for high temperature strength and corrosion resistance. The target for alloy performance will be a 15% (~100 ⁰C) increase in allowable service temperature beyond heat-resistant stainless steels while maintaining or improving alloy cost and corrosion resistance. Haynes 230 was selected as a baseline, which is 57 wt% Ni with 22 wt% Cr 14 wt% W, and 2 wt% Mo as solid solution strengtheners. In addition to rapid design via Bayesian machine learning, the alloys were rapidly fabricated using a multi-wire arc additive manufacturing (mWAAM) technique which allows for precise control of alloy composition and assessing of alloy design “windows” to study composition effects. Build speeds for wire-arc additive processes are among the highest for additive technologies enabling rapid and reliable sample fabrication when compared to conventional methods such as arc button melting. The mWAAM samples will be rapidly characterized via instrumented indentation for room temperature modulus and strength and for elevated temperature strength via hot hardness tests. After being screened with hardness testing, potential alloys will be further evaluated with conventional microscopy techniques including scanning electron microscopy (SEM) and transmission electron microscopy (TEM) to assess agreement with modeling results. The most promising compositions will also be evaluated by printing full sized tensile specimens for mechanical behavior tests at elevated temperatures. Results: Bayesian machine learning of a single performance function was initially used to optimize five performance metrics: 1) single phase stability, 2) yield strength, 3) creep resistance (low diffusion coefficient), 4) freezing range (weldability), and 5) material cost. The single performance function was suboptimal as assumptions had to be made about the results while formulating the optimization. A goal-oriented Bayesian optimization strategy (Hanaoka, 2021) was implemented with CALPHAD for use with the five metrics above. This multi-objective Bayesian optimization (MOBO) enabled the design of NiCrCoFe alloys with V and W additions. A base composition of NiCoCr was selected as Ni provides a stable FCC matrix, Cr aids corrosion/oxidation resistance, and Co is a solid-solutions strengthener that also improves creep by increasing the activation energy. Fe helps reduce diffusion coefficients and cost. Finally, V and W were selected for their reasonable solubility and high atomic misfit to aid in solid solution strengthening. Cracking of the mWAAM specimens was an early issue, and the Easton solidification cracking model (Easton et al., 2014a) was selected for addition to the MOBO function. High performing alloys fabricated by mWAAM included Ni 28 Cr 25 Co 26 Fe 15 V 8 and Ni 62 Cr 18 Co 1 Fe 3 W 15 . It was observed that even after adapting the mWAAM process for W, the W did not fully dissolve. To fully evaluate the Ni 62 Cr 18 Co 1 Fe 3 W 15 composition, a cored wire (80-20 NiCr sheath/powder core) was manufactured and printed via WAAM, and HIP’ing was utilized to homogenize and densify the printed alloy. The V and W alloys produced met metrics 1 (solid solution), 4 (solidification cracking), and 5 (cost). However, an unmodeled mechanism of thermal stress cracking was identified in the WAAM produced materials, perhaps exacerbated by the lack of grain boundary strengthening elements (B, C). Conclusions & Recommendations: A high-throughput (rapid) alloy design technique was applied to designing and manufacturing new high entropy alloys (HEAs) for extreme environments utilizing MOBO and mWAAM. The developed process was rapid and effective in addressing the mechanisms included in the model. The lack of grain boundary strengthening element additions (e.g., B, C) was a simplification that likely produced thermal stress cracking that turned into a large part of the investigation. Additions on the order of 0.005 wt% B and 0.05 wt% C likely would have minimized thermal stress grain boundary cracking. Overall, the high throughput design strategy is promising for rapid design of metrics-driven alloys for advanced ultra supercritical (A-USC) steam cycles for power generation. The MOBO and mWAAM process could be commercialized to accelerate metrics-driven alloy design. In addition, the cored-wire process utilized for scale-up is a promising high-volume process for WAAM alloy development and scale-up.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Long-term permeability evolution of shale seal rocks with argon and scCO 2

In this work, we conducted a suite of experiments to evaluate the long-term permeability evolution of shale rocks under constant effective stress, before, during, and after interaction with supercritical carbon dioxide (scCO 2 ). To do so, we measured the time-dependent evolution of argon permeabilities for the pre- and post-reacted samples (i.e., before and after long-term interaction with scCO 2 ). In addition, we obtained permeability evolution during long-term interaction with scCO 2 . The samples showed either relatively constant permeabilities or a moderate decrease during pre reaction long-term argon tests. The permeability evolution during long-term CO 2 tests showed continuous increase, continuous decrease, or cycles of increase/decrease in permeability. The long term response of the samples to CO 2 included phenomena such as 20 (i) salt precipitation, (ii) swelling-induced cracks, and (iii) carbonate dissolution. While it is obvious that salt precipitation and swelling-induced cracks decrease and increase the permeability, respectively, the sample response to carbonate dissolution proved to be more complex and may increase or decrease the permeability. The permeability evolution during post-reaction long-term argon injection is also affected by the contribution of each of these three phenomena, during long-term interaction with CO 2 . We observe increase, decrease, and constant permeability evolution during post-reaction argon tests. Our experiments reveal that the initial permeability of the samples plays a significant role on the long-term permeability response of shales in the presence argon and scCO 2 fluids. This study shows, when shales are hydraulically fractured with CO 2 , their initial permeability has a more significant role than their permeability evolution over time.

58 GEOSCIENCES↗
DOE OSTI2022

Transmission Electron Microscopy based Characterization of a U-20Pu-10Zr Fuel Irradiated in Experimental Breeder Reactor-II

U-Pu-Zr metallic fuels are important fuel candidate for future advanced and/or test reactors. To better understand fuel performance and guide future fuel development, a U-20Pu-10Zr (in weight) metallic fuel irradiated in the Experimental Breeder Reactor was revisited using advanced electron microscopies. This fuel was irradiated to a burnup of 7.6 % at a peak cladding temperature of 521 °C. In this research, several transmission electron microscopy samples were extracted from different but representative radial locations of the fuel cross section using focused ion beam technique. Phase identification and chemical analysis with sub-micron spatial resolution was carried out using scanning transmission electron microscopy. Multiple phenomena that are critical to fuel performance, such as Zr phase and its distribution in a different matrix, were revealed in an unpreceded and highly detailed manner. An improved understanding of fuel reconstruction in U-Pu-Zr metallic fuel under reactor irradiation is provided.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Orbit response measurements with oscillating trajectories

Optical properties of components in a linac can be reconstructed from analysis of the orbit response matrix, composed from responses of BPMs to changes in currents of dipole correctors. An efficient way to record such response, discussed in this note, is to excite the corrector current with a sine wave and record its Fourier component of the BPM spectrum. Because of a narrow band of the excitation signal, response to multiple correctors can be recorded at the same time without significantly affecting the signal-to-noise ratio. For that, the correctors are oscillated at different frequencies, and the corresponding BPM Fourier components are recorded. In an ion (non-relativistic) linac, similar measurements can be performed for the longitudinal plane as well. Examples of measurements at the PIP2IT test linac are shown.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

A Solubility Comparison of Cyclotetramethylenetetranitramine (HMX) and Potential Mock-HMX Candidates in Solvents Relevant to Slurry Formulating Polymer-Bonded Materials

Understanding solubility behavior of organic crystalline materials is critical in the development of slurry-formulated plastic bonded composites. In the field of high explosives (HE), a common method for obtaining high density, greater than 95 % theoretical maximum density, is to encapsulate HE in a polymer matrix resulting in the formation of plastic bonded explosive (PBX) prills. Prill formation is highly nuanced and requires understanding of particle size distributions, solvent, antisolvent, and binder system as well as interactions thereof. In this paper, in an effort to create mock-HE that mimics HE in density and mechanical strength, solubility of candidate molecules was compared to cyclotetramethylenetetranitramine (HMX) in binder solvents relevant to prill formation. The results indicate that the mock candidates chosen for this study have different solubilities compared to HMX in water and both organic solvents, 2-butanone and ethyl acetate, investigated. With these solubility differences known, prills to match PBXs based on HMX can be formulated, pressed and tested.

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

Intrinsic quality factor extraction of multi-port cavity with arbitrary coupling

In this work, we derived S-parameter based expressions for the intrinsic quality factor of an arbitrary coupled multiport microwave cavity and a non-ideal test fixture. Practical accuracy limitations of the obtained expressions specifically for superconducting accelerator cavities were evaluated both analytically and with the simulation software. The resulting formulas can be used to extract the intrinsic quality factor of normal-conducting and superconducting cavities directly from the calibrated S-parameter measurements.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Probability of Initiation in Neutron Transport

We discuss the numerical solution of the nonlinear integro-differential equation for the probability of a divergent neutron chain in a stationary system (i.e., the probability of initiation (POI)). We follow the development described in Bell’s classic paper on the stochastic theory of neutron transport. As noted by Bell, the linearized form of this equation resembles the linear adjoint neutron transport equation. A matrix formalism for the discretized steady state (or forward) neutron equation in slab geometry is first developed and is then used to derive the discrete adjoint equation. A main advantage of this discrete development is that the resulting discrete adjoint equation does not depend upon how the multigroup cross sections for the forward problem are obtained. That is, we derive the discrete adjoint directly from the discrete forward equations rather than discretizing directly the adjoint equation. This also guarantees that the discrete adjoint operator is consistent with the inner product used to define the adjoint operator. We discuss three approaches for the numerical solution of the POI equations, and present numerical results on several test problems. The three solution methods are a simple fixed-point iteration, a second approach that is akin to a nonlinear Power iteration, and a third approach which uses a Newton-Krylov nonlinear solver. We also give sufficient conditions to guarantee the existence and uniqueness of nontrivial solutions to our discrete POI equations when the discrete system is supercritical, and that only the trivial solution exists when the discrete system is subcritical. Our approach is modeled after the analysis presented for the continuous POI equations by Mokhtar-Kharroubi and Jarmouni-Idrissi, and by Pazy and Rabinowitz.

42 ENGINEERING↗
DOE OSTI2021

Distribution Grid Modeling Using Smart Meter Data

The knowledge of distribution grid models, including topologies and line impedances, is essential for grid monitoring, control and protection. However, such information is often unavailable, incomplete or outdated. The increasing deployment of smart meters (SMs) provides a unique opportunity to tackle this issue. This paper proposes a two-stage framework for distribution grid modeling using SM data. In the first stage, the network topology is identified by reconstructing a weighted Laplacian matrix of distribution networks. In the second stage, a least absolute deviations (LAD) regression model is developed for estimating line impedance of a single branch based on the nonlinear (inverse) power flow model, wherein a conductor library is leveraged to narrow down the solution space. The LAD regression model is originally a mixed-integer nonlinear program whose continuous relaxation is still non-convex. Furthermore, we specially address its convex relaxation and discuss the exactness. The modified regression model is then embedded within a bottom-up sweep algorithm to achieve the identification across the network in a branch-wise manner. Numerical results on the IEEE 13-bus, 37-bus and 69-bus test feeders validate the effectiveness of the proposed methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Efficient ptychographic phase retrieval via a matrix-free Levenberg-Marquardt algorithm

The phase retrieval problem, where one aims to recover a complex-valued image from far-field intensity measurements, is a classic problem encountered in a range of imaging applications. Modern phase retrieval approaches usually rely on gradient descent methods in a nonlinear minimization framework. Calculating closed-form gradients for use in these methods is tedious work, and formulating second order derivatives is even more laborious. Additionally, second order techniques often require the storage and inversion of large matrices of partial derivatives, with memory requirements that can be prohibitive for data-rich imaging modalities. We use a reverse-mode automatic differentiation (AD) framework to implement an efficient matrix-free version of the Levenberg-Marquardt (LM) algorithm, a longstanding method that finds popular use in nonlinear least-square minimization problems but which has seen little use in phase retrieval. Furthermore, we extend the basic LM algorithm so that it can be applied for more general constrained optimization problems (including phase retrieval problems) beyond just the least-square applications. Since we use AD, we only need to specify the physics-based forward model for a specific imaging application; the first and second-order derivative terms are calculated automatically through matrix-vector products, without explicitly forming the large Jacobian or Gauss-Newton matrices typically required for the LM method. We demonstrate that this algorithm can be used to solve both the unconstrained ptychographic object retrieval problem and the constrained “blind” ptychographic object and probe retrieval problems, under the popular Gaussian noise model as well as the Poisson noise model. We compare this algorithm to state-of-the-art first order ptychographic reconstruction methods to demonstrate empirically that this method outperforms best-in-class first-order methods: it provides excellent convergence guarantees with (in many cases) a superlinear rate of convergence, all with a computational cost comparable to, or lower than, the tested first-order algorithms.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Ground-state-based model reduction with unitary circuits

Here, we present a method to numerically obtain low-energy effective models based on a unitary transformation of the ground state. The algorithm finds a unitary circuit that transforms the ground state of the original model to a projected wavefunction with only the low-energy degrees of freedom. The effective model can then be derived using the unitary transformation encoded in the circuit. We test our method on the one-dimensional and two-dimensional square-lattice Hubbard model at half-filling, and obtain more accurate effective spin models than the standard perturbative approach.

Hubbard model↗
DOE OSTI2026

Neural posterior unfolding

Differential cross section measurements are the currency of scientific exchange in particle and nuclear physics. A key challenge for these analyses is the correction for detector distortions, known as deconvolution or unfolding. Binned unfolding of cross section measurements traditionally rely on the regularized inversion of the response matrix that represents the detector response, mapping pre-detector (`particle level') observables to post-detector (`detector level') observables. In this paper we introduce Neural Posterior Unfolding, a modern, Bayesian approach that leverages normalizing flows for unfolding. By using normalizing flows for neural posterior estimation, NPU offers several key advantages including implicit regularization through the neural network architecture, fast amortized inference that eliminates the need for repeated retraining, and direct access to the full uncertainty in the unfolded result. In addition to introducing NPU, we implement a classical Bayesian unfolding method called Fully Bayesian Unfolding (FBU) in modern Python so it can also be studied. These tools are validated on simple Gaussian examples and then tested on simulated jet substructure examples from the Large Hadron Collider (LHC). We find that the Bayesian methods are effective and worth additional development to be analysis ready for cross section measurements at the LHC and beyond.

Analysis and statistical methods↗
DOE OSTI2021

Development of Multi-Filament Textured-Powder Bi-2212/Ag Wire With Enhanced Local Area Ratio

The performance and cost of Bi-2212/Ag wire is limited by the large fraction of Ag matrix (~3:1) that is required in the oxide-powder-in-tube fabrication process. An alternative fabrication process is being developed in which fine-powder Bi-2212 is uniaxially compressed to form bars with a thin Ag foil sheath. The fine powder naturally textures (aligns the a-b planes perpendicular to the direction of compaction) with texture >80% using 200 MPa compression. A billet is formed by stacking trapezoidal- cross-section bars in a symmetric 8-12-16 pattern around an Ag rod and enclosed in an Ag-wall extrusion can. Here, the billet is extruded and drawn to fine wire. Results are presented on present status of the development and testing.

17 WIND ENERGY↗
DOE OSTI2021

Sub-leading structures in superconformal indices: subdominant saddles and logarithmic contributions

We systematically study various sub-leading structures in the superconformal index of N = 4 supersymmetric Yang-Mills theory with SU(N) gauge group. We concentrate in the superconformal index description as a matrix model of elliptic gamma functions and in the Bethe-Ansatz presentation. Our saddle-point approximation goes beyond the Cardy-like limit and we uncover various saddles governed by a matrix model corresponding to SU(N) Chern-Simons theory. The dominant saddle, however, leads to perfect agreement with the Bethe-Ansatz approach. We also determine the logarithmic correction to the superconformal index to be log N, finding precise agreement between the saddle-point and Bethe-Ansatz approaches in their respective approximations. We generalize the two approaches to cover a large class of 4d N = 1 superconformal theories. We find that also in this case both approximations agree all the way down to a universal contribution of the form log N. The universality of this last result constitutes a robust signature of this ultraviolet description of asymptotically AdS 5 black holes and could be tested by low-energy IIB supergravity.

1/N expansion↗
DOE OSTI2025

Code Coverage Status of the ARC Code PERSENT

The Argonne Reactor Code (ARC) software system supports users in their fast reactor design goals by providing neutronic, thermal-hydraulic, and structural analysis capabilities. PERSENT fulfills the role of generating reactivity coefficients for a given time point of a REBUS calculation usable in a point kinetics based safety analysis capability. PERSENT also provides a sensitivity coefficient capability on eigenvalue, reactivity worth, and several other key coefficients that are used in the follow-on safety analysis. Given a co-variance matrix, PERSENT can carry out the uncertainty quantification to indicate the amount of error in the reactivity coefficients derived from the errors in the cross section measurements. With continued improvement of computational resources, many of the geometry modeling capabilities in DIF3D that were primarily used in low order schemes are not really needed anymore. Today, the diffusion and transport capabilities of DIF3D-VARIANT are primarily used in the reactor design process with some scattered usage of DIF3D-FD and DIF3D-Nodal. PERSENT is part of the ARC code system and is built around DIF3D-VARIANT and the flux solution it provides. The purpose of the present work is to identify a set of test problems for PERSENT and assess the code coverage of PERSENT for those test problems. PERSENT treats the DIF3D executable as an external executable and thus the code coverage considerations only need to focus on the PERSENT source code and only a fraction of the connected modules in the existing ARC software library. The goal is to document what parts of the existing PERSENT code are touched by the set of test problems and which are not. Because the verification work done on PERSENT was focused on the most common uses of PERSENT for fast reactor analysis, the code coverage assessment of those capabilities is the highest priority. This will ensure that nothing is being missed by the existing verification test problems that users of PERSENT rely upon. The code coverage analysis of PERSENT was performed with the Code Coverage Tool of the Intel Fortran compiler which requires modifications to the compilation of PERSENT. The detailed coverage tables are given for each submodule of PERSENT. Most of the uncovered parts/files could be easily ignored because they are either for error message and debugging output or not needed by PERSENT today. Only a few uncovered parts of PERSENT deserve extending the verification test suite.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Investigation of tearing mode stability near ideal stability boundaries via asymptotic matching techniques

A number of improvements to the TJ toroidal tearing mode code [Fitzpatrick, Phys. Plasmas 31, 102507 (2024)] are documented. The TJ code is also successfully benchmarked against the STRIDE toroidal tearing mode code [Glasser and Koleman, Phys. Plasmas 25, 082502 (2018)]. Finally, the new capabilities of the TJ code are used to investigate the stability of tearing modes in tokamak plasmas as an ideal stability boundary, associated with either an external-kink or an internal-kink mode, is approached. All elements of the tearing stability matrix are found to tend to infinity as an ideal stability boundary is approached. Furthermore, as the stability boundary is approached, the eigenfunctions of the various tearing modes in the plasma, which are decoupled by sheared plasma rotation, are all found to morph into that of the marginally stable ideal mode. However, the growth-rates and real frequencies of the various “ideal-tearing-modes” are different from one another. Moreover, the growth-rate of the ideal-tearing-mode that reconnects magnetic flux at the rational surface that lies closest to the edge of the plasma is the one that tends to a very large value as the stability boundary is approached. A relatively simple test for ideal stability that is capable of detecting stability boundaries for external-kink and internal-kink modes, even in the presence of a very close-fitting ideal wall, is described and verified.

FOS: Physical sciences↗