Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “extended MHD plasmas”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

2D Magnetohydrodynamic Simulations of the Electrothermal Instability in Metallic Liners

The Virginia Tech (VT) Plasma Dynamics Laboratory Computational (PDCL) and Lawrence Livermore National Laboratory (LLNL) are performing two dimensional (2D) simulations of the electrothermal instability (ETI) using the LLNL multi-physics code Ares. Ares is a multi-physics arbitrary–Lagrangian-Eulerian (ALE) code developed by LLNL and is of particular use in studying magnetohydrodynamic (MHD) instabilities like the ETI due to its resistive MHD, magnetic diffusion, and radiative-hydrodynamics packages. Among its capabilities, it has the ability to model material strength, perform adaptive mesh refinement (AMR), and incorporate a wide variety of equations of state models and conductivity models. The 2D Ares simulation model created by VT-LLNL for studying the development and growth of the electrothermal instability has been configured with initial conditions based on the Mykonos Electrothermal Instability II (METI-II) experiments described by this grant and conducted by team members at the University of Nevada (UNR), the University of New Mexico (UNM), and Sandia National Laboratories. Previously, preliminary 2D Ares simulations of the ETI had been run to approximately 80ns. The rods in these simulations were initiated with sinusoidal perturbations at a similar order of magnitude to those measured on the aluminum rods used for the Mykonos experiment. This model has been improved by increasing the spatial resolution of the simulations and running the simulations further in time. In addition to the simulation run-times extending, the preliminary sinusoidal perturbation has been replaced with a perturbation derived from amplitude measurements by the experimental team, thereby correlating the simulation inputs better to the experimental runs. These new runs are capable of reaching 120ns of simulated time for the uncoated cases and to 200ns the 41 μm coated cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Investigation of Anchor Link Magnetohydrodynamic Effects in the Toroidally Symmetric Lead-Lithium (TSLL) Blanket Concept

This paper presents a comprehensive study on the magnetohydrodynamic (MHD) flow in a slotted channel with a cylindrical anchor link, which is an imperative component of the toroidally symmetric lead-lithium (TSLL) blanket concept. Following the validation of the MHD solver implemented in COMSOL and a proposed subtraction approach to compute the anchor link pressure drop, the effects of computational domain size on the pressure drop and velocity distribution are examined. The results show that the pressure drop associated with the anchor link and the maximum flow velocity follow an asymptotic trend as the domain width increases, with wall-induced pressure drop being more dominant than that of the anchor link. The velocity distribution analysis revealed the formation of an internal boundary layer extending along the magnetic field direction, which is a unique feature of the investigated MHD flow. An estimation of the total MHD pressure drop associated with the array of anchor links under the TSLL blanket conditions suggests ~0.3 MPa, which is significantly lower than the recommended maximum allowable blanket pressure drop of 2 MPa. In conclusion, this work offers valuable insights into the anchor link–associated MHD phenomena and serves as a foundation for further development of the TSLL blanket concept for future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Overview of the TCV tokamak program: Scientific progress and facility upgrades

The TCV tokamak is augmenting its unique historical capabilities (strong shaping, strong electron heating) with ion heating, additional electron heating compatible with high densities, and variable divertor geometry, in a multifaceted upgrade program designed to broaden its operational range without sacrificing its fundamental flexibility. The TCV program is rooted in a three-pronged approach aimed at ITER support, explorations towards DEMO, and fundamental research. A 1 MW, tangential neutral beam injector (NBI) was recently installed and promptly extended the TCV parameter range, with record ion temperatures and toroidal rotation velocities and measurable neutral-beam current drive. ITER-relevant scenario development has received particular attention, with strategies aimed at maximizing performance through optimized discharge trajectories to avoid MHD instabilities, such as peeling-ballooning and neoclassical tearing modes. Experiments on exhaust physics have focused particularly on detachment, a necessary step to a DEMO reactor, in a comprehensive set of conventional and advanced divertor concepts. The specific theoretical prediction of an enhanced radiation region between the two X-points in the low-field-side snowflake-minus configuration was experimentally confirmed. Fundamental investigations of the power decay length in the scrape-off layer (SOL) are progressing rapidly, again in widely varying configurations and in both D and He plasmas; in particular, the double decay length in L-mode limited plasmas was found to be replaced by a single length at high SOL resistivity. Experiments on disruption mitigation by massive gas injection and electron-cyclotron resonance heating (ECRH) have begun in earnest, in parallel with studies of runaway electron generation and control, in both stable and disruptive conditions; a quiescent runaway beam carrying the entire electrical current appears to develop in some cases. Developments in plasma control have benefited from progress in individual controller design and have evolved steadily towards controller integration, mostly within an environment supervised by a tokamak profile control simulator. TCV has demonstrated effective wall conditioning with ECRH in He in support of the preparations for JT-60SA operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Near-ideal relaxed MHD in slab geometry

We investigate the solutions of the relaxed magnetohydrodynamic (MHD) model (RxMHD) of R. Dewar and Z. Qu. This model generalizes Taylor relaxation by including the ideal Ohm's law constraint using an augmented Lagrangian method, providing a pathway to extend the multi-region relaxed MHD (MRxMHD) model. We present the first numerical solution of the RxMHD model by Dewar and Qu, demonstrating that it is mathematically well-defined and computationally feasible for constructing MHD equilibria in slab geometry. We also show that a cross-field flow can exist without enforcing an arbitrary constraint on the angular momentum, as is done in the case of MRxMHD with flow. Our results also demonstrate the self-organization of fully relaxed regions during the optimization, which was an important motivation behind developing this model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2008

Recent Insights into the Nature of Turbulence in the Solar Wind

During the past several years, studies of solar wind turbulence using data from Cluster and other spacecraft, and results from new numerical simulations, have revealed new aspects of solar wind turbulence. I will try to highlight some of that research. At the shortest length scales and highest frequencies, there is renewed interest in determining how the turbulence dissipates, e.g., whether by kinetic Alfven waves or whistler turbulence. Finding observational evidence for exponential damping of solar wind fluctuations has proven challenging. New studies using a combination of flux gate and search coil magnetometer data from Cluster have extended this search (in the spacecraft frame of reference) to more than 10 Hertz. New models and simulations are also being used to study the dissipation. A detailed study of fluctuations in the magnetosheath suggests that turbulent dissipation could be occurring at very thin current sheets as had been suggested by two-dimensional MHD simulations more than 20 years ago. Data from the four Cluster spacecraft, now at their maximum separation of 10,000 km provide new opportunities to investigate the symmetry properties, scale lengths, and the relative proportion of magnetic energy in parallel and perpendicular wave numbers of solar wind turbulence. By utilizing well-calibrated electron data, it has been possible to take advantage of the tetrahedral separation of Cluster in the solar wind near apogee to measure directly the compressibility and vorticity of the solar wind plasma.

Goldstein, Melvun L.↗
DOE OSTI2024

Low latency optical-based mode tracking with machine learning deployed on FPGAs on a tokamak

Active feedback control in magnetic confinement fusion devices is desirable to mitigate plasma instabilities and enable robust operation. Optical high-speed cameras provide a powerful, non-invasive diagnostic and can be suitable for these applications. Here, in this study, we process high-speed camera data, at rates exceeding 100 kfps, on in situ field-programmable gate array (FPGA) hardware to track magnetohydrodynamic (MHD) mode evolution and generate control signals in real time. Our system utilizes a convolutional neural network (CNN) model, which predicts the n = 1 MHD mode amplitude and phase using camera images with better accuracy than other tested non-deep-learning-based methods. By implementing this model directly within the standard FPGA readout hardware of the high-speed camera diagnostic, our mode tracking system achieves a total trigger-to-output latency of 17.6 μs and a throughput of up to 120 kfps. This study at the High Beta Tokamak-Extended Pulse (HBT-EP) experiment demonstrates an FPGA-based high-speed camera data acquisition and processing system, enabling application in real-time machine-learning-based tokamak diagnostic and control as well as potential applications in other scientific domains.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2007

Role of Ionospheric Plasmas in Earth's Magnetotail

This tutorial will summarize observations and theories indicating a prominent role of ionospheric plasma in the Earth's magnetotail. At the Global scale, I will argue that it is ionospheric plasma momentum and dynamic pressure that are responsible for the production of plasmoids, through the action of a transient near-Earth neutral or X-line, which serves to release excessive plasma pressure from the magnetotail field. Ionospheric plasma gains the momentum and energy to produce plasmoids and their related effects through its interaction with the solar wind, beginning at the dayside reconnection region and extending across the polar caps through the magnetotail lobes. This distant neutral line can be depicted as a feature much like that found in cometary magnetospheres, where disconnection limits the amount of IMF hung up on the cometary coma. On the other hand, the near-Earth neutral one can be seen as a feature unique to planets with an intrinsic magnetic field and internal source of plasma, the heating of which produces pressures too large to be restrained. Ionospheric plasmas also have other more local roles to play in the magnetotail. The circulation influences the composition of the plasma sheet, and the resultant wave environment, giving rise to reduced wave propagation speeds. Important heavy ion cyclotron resonances, and enhanced finite gyro-radius effects including non-adiabatic particle acceleration. At minimum, the presence of ionospheric plasma must influence the rate of reconnection via its enhanced mass density. Other non-MHD effects of ionospheric plasma presence are likely to be important but need much more investigation to be well understood. The MMS mission is designed to penetrate the subtle diffusion region physics that is involved, and its ability to observe ionospheric plasma involvement in reconnection will contribute significantly toward that goal.

Moore, Thomas E.↗
NASA NTRS1981

Alternative formulations of magnetospheric plasma electrodynamics

The fundamental equations of magnetospheric plasma electrodynamics are considered from a theoretical standpoint that stresses the basic equivalence of various seemingly different formal representations. The mathematical properties of vector fields are reviewed, and their implications in electrodynamics are studied. The irrotational and solenoidal parts of the electric field are associated with two physically distinct types of sources. Relativistic covariance and gauge invariance in electromagnetic theory are reviewed and discussed in the context of an approach in which the mathematical properties of vector fields are taken as primary concepts. Special attention is given to the use and interpretation of the Coulomb gauge potential functions. This choice of gauge is sometimes regarded with undue suspicion, possibly because of a certain paradox concerning causality. The paradox is discussed and resolved. Useful properties of the Coulomb gauge are identified. These need not be limited to the case of slow time variations and can extend beyond the limits of validity of ideal MHD theory.

Cragin, B. L.↗
DOE OSTI2021

Pressure drop in a prototypical 3D magnetohydrodynamic flow across contraction of a fusion blanket manifold

Predicting 3D magnetohydrodynamic (MHD) pressure losses associated with liquid metal flows in complex geometry ducts is required to design breeding blankets for fusion reactors. Of such components, manifolds exhibit major pressure losses. Recent correlations [T. Rhodeset al. Magnetohydrodynamic pressure drop and flow balancing of liquid metal flow in a prototypic fusion blanket manifold. Phys. Fluids. 2018; 30: 057101.] demonstrate promise for predicting pressure drops in electrically insulated inlet manifolds. In the present work, the extend to which these correlations can be applied to outlet manifolds, which feature sudden contractions, is investigated in both viscous-electromagnetic (VE) and inertial-electromagnetic (IE) regimes for Reynolds numbers 50 < Re < 1500, Hartmann numbers 2500 < Ha < 5475, and expansion/contraction ratio 10. Numerical computations have shown that the MHD pressure drops in the eutectic lead-lithium (PbLi) flows in contractions are almost identical to those in the flows featuring expansion with slightly higher magnitudes in the contraction cases. The small discripancy in the MHD pressure drop between contractions and expansions (<8%) suggests that the earlier obtained correlations for the 3D MHD pressure drop in a duct flow with a sudden expansion can also be applied to flows with a sudden contraction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Solving differential equations using deep neural networks

Recent work on solving partial differential equations (PDEs) with deep neural networks (DNNs) is presented. The paper reviews and extends some of these methods while carefully analyzing a fundamental feature in numerical PDEs and nonlinear analysis: irregular solutions. First, the Sod shock tube solution to the compressible Euler equations is discussed and analyzed. This analysis includes a comparison of a DNN-based approach with conventional finite element and finite volume methods, and demonstrates that the DNN is competitive in terms of degrees of freedom required for a given accuracy. Further, the DNN-based approach is extended to consider performance improvements and simultaneous parameter space exploration. Next, a shock solution to compressible magnetohydrodynamics (MHD) is solved for, and used in a scenario where experimental data is utilized to enhance a PDE system that is a priori insufficient to validate against the observed/experimental data. This is accomplished by enriching the model PDE system with source terms that are then inferred via supervised training with synthetic experimental data. The resulting DNN framework for PDEs enables straightforward system prototyping and natural integration of large data sets (be they synthetic or experimental), all while simultaneously enabling single-pass exploration of an entire parameter space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Observations of hydromagnetic turbulence in the solar wind

MHD turbulence is studied by analyzing magnetic field and plasma observations from Helios-1 and -2 at minimum solar activity. The steady conditions in the plasma flows and the magnetic field sector structure in 1975/1976 facilitate an investigation of the radial evolution of the turbulence from 0.29 to 1AU. In high speed streams the fluctuations in the solar wind velocity v and the magnetic field b are highly correlated (the correction coefficient almost being one), which indicates that the turbulence is mainly Alfvenic in high speed plasma. While some general fluctuation properties remain essentially unchanged from 0.29 to 1AU, power spectral analysis reveals a different frequency composition of the Alfvenic turbulence at different heliocentric distances. At 0.3AU much more 'high' frequency fluctuations contribute to the total power in the magnetic field and velocity fluctuations than at 1AU. The contributions of field magnitude fluctuations are found to be distance and frequency dependent. Magnetic field spectra with an extended frequency range up to 470Hz show certain frequency bands, where the steepness of the spectra is independent of the helicocentric distance.

Denskat, K. U.↗
DOE OSTI2024

Simulation of DIII-D disruption with argon pellet injection and runaway electron beam

Abstract The next generation of large tokamaks, including ITER, will be equipped with a disruption mitigation system (DMS) that can be activated if a disruption is deemed to be imminent. Introducing impurities by pellet (large or shattered) or massive gas injection has been shown to be an effective mitigation mechanism on many tokamaks. The goal of the mitigation is to lessen the thermal and electromagnetic loads from the disruption without generating enough high-energy (runaway) electrons to damage the device. Variations of this mitigation process with impurity injection are presently being tested on many experiments. We have modeled one such impurity injection experiment on DIII-D using the M3D-C1 nonlinear 3D extended MHD code (Jardin et al 2012 Comput. Sci. Discovery 6 014002), The model includes an argon large pellet injection and ablation model, impurity ionization, recombination, and radiation, and runaway electron formation and subsequent evolution, including both Dreicer and avalanche sources. We obtain reasonable agreement with the experimental results for the timescale of the thermal and current quench and for the magnitude of the runaway electron plateau formed during the mitigation. This is the first 3D full MHD simulation with pellets and REs to simulate the disruption process and it also provides a partial validation of the M3D-C1 DMS model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

On FIRE mode in KSTAR

We report on the status of Fast Ion Regulated Enhancement (FIRE) mode experiments in the Korea Superconducting Tokamak Advanced Research. This regime is being developed for high-performance, steady-state operation which features a stationary ion internal transport barrier, enabling a central ion temperature approaching 10 keV to be sustained for up to 50 s, without the need for delicate profile control and with no significant impurity accumulation. As its key novelty lies in the significant contribution of fast ions that stabilize core turbulence, the regime has been named FIRE mode. To achieve this regime, neutral beam injection is applied at moderate power levels near the L–H power threshold, while maintaining low plasma density to avoid the L–H transition. The scenario is typically established in diverted magnetic configurations. The core features of FIRE mode were investigated through power balance analysis and fluctuation measurements, revealing a clear transport bifurcation in the ion channel. At the plasma edge, FIRE mode occasionally exhibits I-mode characteristics, particularly in unfavorable magnetic null configurations with q 95 ∼ 4, including the presence of weakly coherent modes. In terms of MHD activity, sawtooth oscillations are observed but appear to be stabilized during the high-performance phase. Fast ion-driven Alfvénic eigenmodes (AEs), indicated by strong frequency chirping near 200 kHz, are also observed. Additionally, lower-frequency MHD activities, distinct from the fast ion-driven AEs, are present and have some influence on plasma performance. The enhancement of core confinement is primarily attributed to fast ion effects, which were evaluated with respect to dilution, alpha stabilization, and resonant interactions with turbulence using gyrokinetic analyses. Among these, the dilution effect was found to be the most dominant. The characteristics of FIRE mode were compared with those of other hot ion plasma scenarios, such as supershot and hot ion mode. While they share many similarities, FIRE mode is distinguished by the accessibility to conditions with T i ≈ T e , presence of I-mode edge features, and its long-duration sustainment. Predictive simulations of FIRE mode were performed using integrated transport modeling with TRIASSIC, employing the TGLF anomalous transport model. These simulations confirmed the critical role of fast ions in achieving this regime. The future prospect of FIRE mode for application in fusion reactors is discussed, with an emphasis on possibly extending the regime to higher density operation.

FIRE mode↗
NASA NTRS1995

The structure of the distanct geomagnetic tail during long periods of northward IMF

We have used a newly developed, parallelized, global MHD magnetosphere - ionosphere simulation model with a 400 R(sub E) long tail to study the evolution, structure, and dynamics of the distant magnetotail during extended periods of northward interplanetary magnetic field (IMF). We find that the tail evolves to a nearly time stationary structure about one solar wind transit time after the IMF turns northward. Four regions of different magnetic topology can be distinguished which extend at least to the end of the simulation box at 400 R(sub E). Besides lobe field lines and open solar wind field lines tailward of an X-line, there is a broad boundary layer of closed field lines which we call the tail flank boundary layer (TFBL). Just inside the TFBL there is a region of closed field loops. Besides the X-line we find two O-lines which are enclosed by the closed field loops and are roughly aligned with the tail axis. Together they form a U shaped separator between the northward and the southward plasma sheet fields.

Raeder, J.↗
DOE OSTI2022

Faster ablative Kelvin–Helmholtz instability growth in a magnetic field

Shear flows along a plasma interface will quickly grow unstable due to the Kelvin–Helmholtz instability. If there is a concurrent temperature gradient across the interface, higher modes are stabilized by the thermal diffusion. These ablative effects must be considered in, for example, jet features in inertial confinement fusion hot-spots, or plasma plumes in young supernovae. We show that magnetization of the plasma can greatly affect the instability, even if magnetic pressure is small. This is because electrons are localized by their gyromotion, reducing the heat flux and material ablation. We use a two-dimensional numerical extended-magnetohydrodynamics approach to assess this effect for dense fusion conditions. In comparison with the unmagnetized case, self-generated Biermann fields make only a minor difference to growth rates. However, simulations with a large 50 kT external field found that the growth rate of the least stable mode increased by 40%. In conclusion, this has implications for mix processes in Z-pinches and magnetized inertial confinement fusion concepts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Effects beyond ideal MHD on stability of wide and enhanced pedestal regimes in NSTX

Stability of edge-localized modes (ELMs) in spherical tokamaks is explored using the extended MHD model. Linear NIMROD simulations have been performed for three NSTX discharges 132543, 132588, and 141133, to investigate the role of resistivity, diffusivity, and shear flows on the onset of ELMs. The first discharge represents the wide pedestal regime and the later two discharges represent ELM-free enhanced pedestal H-mode. We first present the effect of toroidal rotation shear and find a flow shear destabilizing effect in these NSTX discharges. Simulations are also extended to include the two-fluid and ion gyroviscosity effects. Simulations show that the flow shear can shift the mode spectrum and alter the critical condition of ELM onset. We also uncover that ELM onset prediction in spherical tokamaks requires effects beyond MHD, in particular gyroviscosity and diamagnetic terms could stabilize Peeling-Ballooning modes consistent with the experimental observation of ELM-free regimes in NSTX. The findings give new insight into the nature of the interplay between resistivity, flows, and diamagnetic stabilization in ELM suppression and have potential applications to ELM control schemes in NSTX-U and next-generation spherical tokamaks. This study identifies the essential physical effects that must be included in future predictive and validation simulations.

NSTX↗
DOE OSTI2020

Free-boundary MRxMHD equilibrium calculations using the stepped-pressure equilibrium code

Here the stepped-pressure equilibrium code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 , 112 502) is extended to enable free-boundary multi-region relaxed magnetohydrodynamic (MRxMHD) equilibrium calculations. The vacuum field surrounding the plasma inside an arbitrary 'computational boundary', , is computed, and the virtual-casing principle is used iteratively to compute the normal field on so that the equilibrium is consistent with an externally produced magnetic field. Recent modifications to SPEC are described, such as the use of Chebyshev polynomials to describe the radial dependence of the magnetic vector potential, and a variety of free-boundary verification calculations are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗