Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shape analysis”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

End-exclusion zones in strongly stretched, molten polymer brushes of arbitrary shape

Theories of strongly stretched polymer brushes, particularly the parabolic brush theory, are valuable for providing analytically tractable predictions for the thermodynamic behavior of surface-grafted polymers in a wide range of settings. Furthermore, the parabolic brush limit fails to describe polymers grafted to convex curved substrates, such as the surfaces of spherical nanoparticles or the interfaces of strongly segregated block copolymers. It has previously been shown that strongly stretched curved brushes require a boundary layer devoid of free chain ends, requiring modifications of the theoretical analysis. While this “end-exclusion zone” has been successfully incorporated into the descriptions of brushes grafted onto the outer surfaces of cylinders and spheres, the behavior of brushes on surfaces of arbitrary curvature has not yet been studied. We present a formulation of the strong-stretching theory for molten brushes on the surfaces of arbitrary curvature and identify four distinct regimes of interest for which brushes are predicted to possess end-exclusion zones, notably including regimes of positive mean curvature but negative Gaussian curvature. Through numerical solutions of the strong-stretching brush equations, we report predicted scaling of the size of the end-exclusion zone, the chain end distribution, the chain polarization, and the free energy of stretching with mean and Gaussian surface curvatures. Through these results, we present a comprehensive picture of how the brush geometry influences the end-exclusion zones and exact strong-stretching free energies, which can be applied, for example, to model the full spectrum of brush geometries encountered in block copolymer melt assembly.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Carrot rhamnogalacturonan I structure and composition changed during 2017 in California

Detailed composition, fine structure, and function were determined for carrot pomace fiber remaining after commercial juice extraction from carrots harvested at four different times during 2017, in different California, USA locations. This pomace was treated with commercial enzymes and the changes in carbohydrate composition were analyzed. The enzyme-treated carrot pomace contained pectic oligosaccharides that were branched, unsaturated, partially acetylated and methylated rhamnogalacturonan I (RG-I) with arabinogalacto-oligosaccharide side chains as well as free galacto-oligosaccharides and galacturono-oligosaccharides detected by MALDI-TOF MS and NMR analysis. The molecular weight and degree of galacto-oligosaccharide polymerization decreased while water solubility increased following enzymatic treatment. While carrot soluble dietary fiber remained relatively constant, insoluble dietary fiber and total dietary fiber increased during the year. The carrot pomace fiber structure changed between the June and September harvests so that it was more susceptible to enzymatic degradation, and its shape became more linear. The diversity of carrot RG-I oligosaccharide structures was the greatest in the enzyme-treated pomace harvested in September and the RG-I became almost completely debranched by the enzyme treatment for the pomace harvested in December. The structural properties of these carrot pomace fiber fractions suggest that they can function as a low viscosity beverage food ingredient, and the rhamnogalacturonan I structure indicated the potential for anti-viral, immuno-modulatory and prebiotic properties. Furthermore, the rhamnogalacturonan structural composition changes during the growing season highlights the need for bioactivity monitoring of carrot pomace obtained from different sources, and harvested at different times during the year.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Chemolithoautotroph distributions across the subsurface of a convergent margin

Subducting oceanic crusts release fluids rich in biologically relevant compounds into the overriding plate, fueling subsurface chemolithoautotrophic ecosystems. To understand the impact of subsurface geochemistry on microbial communities, we collected fluid and sediments from 14 natural springs across a ~200 km transect across the Costa Rican convergent margin and performed shotgun metagenomics. The resulting 404 metagenome-assembled genomes (MAGs) cluster into geologically distinct regions based on MAG abundance patterns: outer forearc-only (25% of total relative abundance), forearc/arc-only (38% of total relative abundance), and delocalized (37% of total relative abundance) clusters. In the outer forearc, Thermodesulfovibrionia, Candidatus Bipolaricaulia, and Firmicutes have hydrogenotrophic sulfate reduction and Wood-Ljungdahl (WL) carbon fixation pathways. In the forearc/arc, Anaerolineae, Ca. Bipolaricaulia, and Thermodesulfovibrionia have sulfur oxidation, nitrogen cycling, microaerophilic respiration, and WL, while Aquificae have aerobic sulfur oxidation and reverse tricarboxylic acid carbon fixation pathway. Transformation-based canonical correspondence analysis shows that MAG distribution corresponds to concentrations of aluminum, iron, nickel, dissolved inorganic carbon, and phosphate. While delocalized MAGs appear surface-derived, the subsurface chemolithoautotrophic, metabolic, and taxonomic landscape varies by the availability of minerals/metals and volcanically derived inorganic carbon. Furthermore, the WL pathway persists across all samples, suggesting that this versatile, energy-efficient carbon fixation pathway helps shape convergent margin subsurface ecosystems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Large-eddy simulations of turbulent flows in arrays of helical- and straight-bladed vertical-axis wind turbines

Effects of helical-shaped blades on the flow characteristics and power production of finite-length wind farms composed of vertical-axis wind turbines (VAWTs) are studied numerically using large-eddy simulation (LES). Two helical-bladed VAWTs (with opposite blade twist angles) are studied against one straight-bladed VAWT in different array configurations with coarse, intermediate, and tight spacings. Statistical analysis of the LES data shows that the helical-bladed VAWTs can improve the mean power production in the fully developed region of the array by about 4.94%–7.33% compared with the corresponding straight-bladed VAWT cases. The helical-bladed VAWTs also cover the azimuth angle more smoothly during the rotation, resulting in about 47.6%–60.1% reduction in the temporal fluctuation of the VAWT power output. Using the helical-bladed VAWTs also reduces the fatigue load on the structure by significantly reducing the spanwise bending moment (relative to the bottom base), which may improve the longevity of the VAWT system to reduce the long-term maintenance cost.

17 WIND ENERGY↗
DOE OSTI2026

The Role of Nuclear Data Sensitivities in Prompt α-Eigenvalue Predictions of Delayed Critical Benchmarks

Alpha (α) eigenvalues, which describe the logarithmic time derivative of the neutron population in a multiplying system, are integral to time-dependent behavior and diagnostic applications. However, uncertainties in the evaluated nuclear data can significantly impact the accuracy of transport simulations for such quantities. This work explores the use of machine learning models to predict two key outputs, α-eigenvalues and keff bias, using input features derived from α-eigenvalue sensitivities to nuclear data. The criticality safety benchmark models used in this study come from the International Handbook of Evaluated Criticality Safety Benchmark Experiments. Three models, random forest, XGBoost, and NGBoost, are trained on both energy-resolved and energy-summed α sensitivities. For the α-eigenvalue bias prediction, NGBoost achieved the highest R 2 (0.9476) using energy-resolved features, while XGBoost performed best using summed sensitivities. In contrast, when predicting the keff bias, all the models showed moderate predictive capability (best R 2 ≈ 0.72), as the mapping from the static α-sensitivities to the static keff bias was less direct. SHAP (SHapley Additive exPlanations) analysis was used to interpret the model predictions. Across both prediction tasks, the features associated with neutron capture [H-1 (n, γ)], uranium scattering reactions (such as 235 U elastic/inelastic), and actinide capture/fission reactions (such as 239 Pu and 234 U) were consistently identified as the most impactful. This highlights the key role of specific nuclear reactions and energy ranges in shaping both time-dependent and steady-state criticality behavior. These results demonstrated that α-sensitivities, despite being computed for time-dependent metrics, can provide valuable insights for predicting both α-eigenvalues and the keff bias. Moreover, machine learning models offer a promising pathway for uncovering important nuclear data dependencies and guiding future data evaluation efforts.

Nuclear data↗
DOE OSTI2023

EV charging site day-ahead load prediction in a synthetic environment for RL based grid-informed charging

Ensuring grid health in the face of increasing demand for power is an emerging challenge especially due to transportation electrification. A free market approach to influencing electric vehicle (EV) load through grid-informed hourly dynamic pricing is introduced in this work. The setting of charging price is done by a reinforcement learning (RL) agent that learns the complicated dynamics by interacting with a synthetic environment. This synthetic environment is a combination of distribution feeder simulation, EV charger user behavior dynamics, and EV charging simulation. A key module in this synthetic environment involves obtaining the day-ahead charging profile of EV charging stations based on real-world past data. The day-ahead prediction is also useful in other traditional optimizations related to EV charge scheduling. The proposed approach involves using EV charging data from two different past time horizons – one to determine the shape of the daily profile and the other to determine a scaling value to capture actual energy consumption. Real-world charging data over many years from the ACN charging network has been used to demonstrate the ability to predict the day-ahead profile with only charge session data. Both Python and MATLAB have been used for data cleaning, processing, analysis, and prediction.

Suryanarayana, Harish↗
DOE OSTI2025

Novel CHI3L1 ‐Associated Angiogenic Phenotypes Define Glioma Microenvironments: Insights From Multi‐Omics Integration

ABSTRACT The CHI3L1 signaling pathway significantly influences glioma angiogenesis, but its role in the tumor microenvironment (TME) remains elusive. We propose a novelCHI3L1‐associated vascular phenotype classification for glioma through integrative analyses of multiple datasets with bulk and single‐cell transcriptome, genomics, digital pathology, and clinical data. We investigated the biological characteristics, genomic alterations, therapeutic vulnerabilities, and immune profiles within these phenotypes through a comprehensive multi‐omics approach. We constructed the vascular‐related risk (VR) score based onCHI3L1‐associated vascular signatures (CAVS) identified by machine learning algorithms. Utilizing unsupervised consensus clustering, gliomas were stratified into three distinct vascular phenotypes: Cluster A, marked by high vascularization and stromal activation with a relatively low levels of tumor‐infiltrating lymphocytes (TILs); Cluster B, characterized by moderate vascularization and stromal activity, coupled with a high density of TILs; and Cluster C, defined by low vascularization and sparse immune cell infiltration. We observed that the CAVS effectively indicated glioma‐associated angiogenesis and immune suppression by single‐cell RNA‐seq analysis. Moreover, the high‐VR‐score group exhibited enhanced angiogenic activity, reduced immune response, resistance to immunotherapy, and poorer clinical outcomes. The VR score independently predicted glioma prognosis and, combined with a nomogram, provided a robust clinical decision‐making tool. Potential drug prediction based on transcription factors for high‐risk patients was also performed. Our study reveals thatCHI3L1‐associated vascular phenotypes shape distinct immune landscapes in gliomas, offering insights for optimizing therapeutic strategies to improve patient outcomes.

Oncology↗
DOE OSTI2024

Optimization of a High Pressure Turbine Blade and Sector-Based Annular Rig Design for Supercritical CO2 Power Cycle Representative Testing

Abstract As part of the ongoing research into the design of hardware for zero emission cycles, a first-stage high-pressure turbine (HPT) blade is optimized for a 300 MWe supercritical CO2 (sCO2) power cycle using the surrogate-assisted genetic algorithm optimizer in Numeca FINE/Design three-dimensional with objectives of increasing efficiency and decreasing heat load to the blade. Supercritical CO2 property tables are constructed from NIST REFPROP data for the condensable gas simulation in FINE/Turbo. A detailed mesh sensitivity study is performed for a baseline design to identify the proper-grid refinement and efficiently allocate resources for the optimization. Seventy design variables are selected for the initial population generation. Self-organizing maps are then used to focus the design variables on the most important ones affecting the objective functions. The optimization results in approximately 3000 three-dimensional Reynolds Averaged Navier Stokes simulations of different blade shapes with increases in efficiency of up to 0.85% and decreases in heat load of 14%. Families of blade shapes are identified for experimental testing in an annular rig at the Purdue Experimental Turbine Aerothermal Laboratory. A design to adapt the annular cascade for testing optimized geometries is introduced, which features eccentric radius sectors allowing for scaled-up geometries of sCO2 optimized blade profiles to be tested at design cycle representative conditions at high Reynolds numbers in dry air. Analysis into the effects of Reynolds number, working fluid, and geometric relations are presented to prove the efficacy of the test method.

Engineering↗
DOE OSTI2022

Functional Traits Resolve Mechanisms Governing the Assembly and Distribution of Nitrogen-Cycling Microbial Communities in the Global Ocean

Microorganisms drive much of the marine nitrogen (N) cycle, which jointly controls the primary production in the global ocean. However, our understanding of the microbial communities driving the global ocean N cycle remains fragmented. Focusing on “who is doing what, where, and how?”, this study draws a clear picture describing the global biogeography of marine N-cycling microbial communities by utilizing the Tara Oceans shotgun metagenomes. The marine N-cycling communities are highly variable taxonomically but relatively even at the functional trait level, showing clear functional redundancy properties. The functional traits and taxonomic groups are shaped by the same set of geo-environmental factors, among which, depth is the major factor impacting marine N-cycling communities, differentiating mesopelagic from epipelagic communities. Latitudinal diversity gradients and distance-decay relationships are observed for taxonomic groups, but rarely or weakly for functional traits. The composition of functional traits is strongly deterministic as revealed by null model analysis, while a higher degree of stochasticity is observed for taxonomic composition. Integrating multiple lines of evidence, in addition to drawing a biogeographic picture of marine N-cycling communities, this study also demonstrated an essential microbial ecological theory—determinism governs the assembly of microbial communities performing essential biogeochemical processes; the environment selects functional traits rather than taxonomic groups; functional redundancy underlies stochastic taxonomic community assembly.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

River sinuosity describes a continuum between randomness and ordered growth

River channels are among the most common landscape features on Earth. An essential characteristic of channels is sinuosity: their tendency to take a circuitous path, which is quantified as along-stream length divided by straight-line length. River sinuosity is interpreted as a characteristic that either forms randomly at channel inception or develops over time as meander bends migrate. Studies tend to assume the latter and thus have used river sinuosity as a proxy for both modern and ancient environmental factors including climate, tectonics, vegetation, and geologic structure. But no quantitative criterion for planform expression has distinguished between random, initial sinuosity and that developed by ordered growth through channel migration. This ambiguity calls into question the utility of river sinuosity for understanding Earth’s history. We propose a quantitative framework to reconcile these competing explanations for river sinuosity. Using a coupled analysis of modeled and natural channels, we show that while a majority of observed sinuosity is consistent with randomness and limited channel migration, rivers with sinuosity ≥1.5 likely formed their geometry through sustained, ordered growth due to channel migration. This criterion frames a null hypothesis for river sinuosity that can be applied to evaluate the significance of environmental interpretations in landscapes shaped by rivers. The quantitative link between sinuosity and channel migration further informs strategies for preservation and restoration of riparian habitat and guides predictions of fluvial deposits in the rock record and in remotely sensed environments from the seafloor to planetary surfaces.

58 GEOSCIENCES↗
DOE OSTI2024

Evaluation of Additively Manufactured Monolithic SiC and SiC Ceramic Matrix Composites for Concentrating Solar Receiver Applications

Supported by an award from the Solar Technology Office, US Department of Energy, GE Aerospace Research in collaboration with Heliogen Holdings Inc, is engaged in the development of ultra-High Operating Temperature SiC-matrix Solar Thermal Air Receiver (HOTSSTAR). We report the results of our study of the effective heat transfer characteristic of several candidate structure motifs, or feature geometries, made of SiC via additive manufacturing in a simulative laboratory test. The SiC structure motifs studied include different permutations of three-dimensional periodic lattices and defined shapes. The solar-thermal simulating laboratory test setup is constructed using a 4kW CO2 laser system with beam shaping optics to apply radiative heating power on one face of 2”-diameter cylindrical feature specimens representing the structure motifs of interest for receiver element design. Using the test setup, simulative test conditions representative of a concentrated solar flux of up to ~2000 suns could be achieved in the lab tests under varying air flow through the test structure. A numerical analysis scheme is developed to extract an effective or compound heat transfer coefficient representative of the test structure under steady-state heat flow conditions.

14 SOLAR ENERGY↗
DOE OSTI2024

Tuning Perovskites’ Hydration-Induced Chemical Expansion with Octahedral Tilt Angles

Hydration-induced strains in proton-conducting oxides compromise chemo-mechanical stability when these materials are applied in protonic ceramic electrochemical cells. Here, to develop design principles for zero-strain materials, we systematically studied the hydration coefficients of chemical expansion (CCEs) in perovskite (Sr, Ba)(Ce, Zr, Y)O 3–x solid solutions with in situ dilatometry and thermogravimetric analysis in the range of 430–630 °C. By including and decoupling a wide range of tolerance factors and lattice parameters, we were able to identify a minimum in hydration CCEs (0–0.02) at intermediate tolerance factor values (t ≈ 0.95). Conversely, despite expectations of lower CCEs in larger unit cells, no general trend in CCE versus lattice parameter was found, and opposite trends could be seen for Sr(Ce, Zr, Y)O 3–x versus Ba(Ce, Zr, Y)O 3–x separately. In situ neutron diffraction (ND) enabled atomistic insight. Upon decreasing t, chemical strain anisotropy increased, but this trend did not match the U-shaped dependence of macroscopic CCEs on t. Instead, perovskites with intermediate t, hosting intermediate octahedral tilt angles in the nominally dry state, underwent the largest change in the B–O–B angles during hydration. Accommodating hydration through decreasing B–O–B angles is beneficial because it does not result in large lattice parameter changes. We propose an intermediate tolerance factor as a simple structural descriptor to enable near-zero hydration strains in proton-conducting perovskites.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Arctic Mixed-Phase Cloud Base Ice Precipitation Properties Over the NSA Site

Cloud-climate feedbacks are still the greatest source of uncertainty in current climate projections. Arctic clouds, which are predominantly stratiform and supercooled, often long-lived, and nearly continuously precipitate ice particles, contribute roughly 10% of the uncertainty attributed to the global cloud feedback. This arctic cloud uncertainty is driven by incomplete observational and theoretical knowledge required to estimate and explain the state and active processes occurring in those clouds. A focus on ice precipitation properties at arctic cloud base rather than the surface deconfounds the product of cloud condensate sink processes from the influence of the atmospheric thermodynamic state below cloud base, rendering cloud-base properties a more appealing target for inference and evaluation of model simulations. This data set provides more than 1800 samples of cloud-base ice precipitation properties over Utqiagvik, North Slope of Alaska, all of which were retrieved using the synthesis of ARM radar and lidar measurements. The retrieved ice precipitation variables in this data set include, among others, the ice number concentration, water content, PSD parameters, precipitation rate, mass-weighted fall velocity, vertical air motion, and effective radius, all of which are highly valuable for model evaluation and a general understanding of polar cloud sink processes. Each variable sample includes its mean value and associated uncertainty. Additional variables based on ARM measurements (liquid layer statistics, etc.) are included in this data set. The retrieval algorithm and analysis of this data set are described in Silber (JGR, 2023, https://doi.org/10.1029/2022JD038202).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Granular Flow in Novel Octet Shape–Based Lattice Frame Material

Moving packed-bed heat exchangers in concentrated solar power (CSP) plants involves heat transfer between heated falling particles and supercritical carbon dioxide. The overall effective thermal conductivity of the moving packed bed and particle-side channel contact resistances are still the bottlenecks in achieving the desirable thermal transport levels. To this end, a novel moving packed bed heat exchanger consisting of an Octet lattice packed between the walls of the particle-side channel is proposed in this study. Granular flow analysis in Octet lattice moving packed bed heat exchanger (OLHX) was conducted through experiments and discrete element method (DEM)-based numerical simulations. The experimental images clearly demonstrated stagnation regions upstream of lattice fibers, void regions downstream of the fiber junctions, and wavy-type unobstructed flow on the lateral sides of the fibers. DEM simulations were successful in capturing all these critical flow phenomena. Larger flow velocities were observed on the lateral sides of the fibers in the simulations. Also, when the particles in the silo were emptied, the final images showed an accumulation of particles on the inter-fiber as well as fiber–channel wall junctions. Moreover, the fiber connections resulted in some regions devoid of particle contact on the channel endwall, which means that these regions would suffer from poor thermal exchange. Finally, the overall mass flowrate increased with increasing porosity for a fixed particle diameter.

14 SOLAR ENERGY↗
DOE OSTI2020

Analysis of the effect of spherical aberrations on a solenoid scan

Solenoid scans are often used to determine the beam emittance of an intense relativistic electron beam (IREB). In a solenoid scan, the beam size is measured at a location downstream of the solenoid as a function of the strength of the solenoid magnetic field. The emittance is then inferred using a beam optics model of the transport of the beam through the solenoid magnet to the imaging location. This requires precise knowledge of the field strength and shape as a function of the current supplied to the solenoid. The beam optics model for an IREB needs to include an accurate model of the magnetic field. Unfortunately, most beam optics models do not include the effect of spherical aberrations which are inherent in a solenoid magnet. Spherical aberrations are manifested as an essentially quadratic increase in the solenoid magnetic field as a function of radius. One consequence of this is that the focal length will be dependent on the beam radius. This note examines the sensitivity of a typical solenoid scan performed on DARHT axis 2 to spherical aberrations. The sensitivity to spherical aberrations is examined as a function of the beam envelope size and the drift distance from the solenoid magnet to the imaging station.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Can classical DEM simultaneously capture compressibility and flowability of milled biomass?

Accurate prediction of the rheological behavior of biomass is essential for the design and operation of hoppers, feeders, and storage systems in biorefineries. This study examines whether the classical, coarse-grained discrete element method (DEM) formulation can simultaneously reproduce the compressibility and flowability of milled herbaceous biomass, using Miscanthus × giganteus as a representative material. The model represents particles as rigid spheres interacting through Hertz-Mindlin elastic-frictional contacts augmented with an area-dependent cohesion term. Laboratory cyclic compression and wedge-shaped hopper discharge experiments were used as calibration benchmarks. Although the model can independently reproduce each behavior by appropriately tuning particle Young's modulus E and cohesion energy density k, an extensive parametric investigation comprising more than 600 simulations reveals that the optimal parameter regions for compression and hopper flow are distinct and non-overlapping in (E, k) space. Surrogate surface analysis further shows that the corresponding objective-function valleys exhibit similar trends but are approximately parallel and spatially offset, precluding a unified calibration within the explored domain. Sensitivity analysis indicates that compressibility is governed predominantly by stiffness and cohesion, whereas the slope of the mass flow rate-opening relation in hopper discharge is primarily controlled by tangential friction. Extensions incorporating particle size distribution and clumped-sphere representations do not eliminate the incompatibility. These results systematically reveal, for the first time, the structural limitation of simplified DEM formulations in representing biomass rheological behavior, underscoring the necessity for models incorporating additional physical mechanisms, such as particle deformability or enhanced interlocking, to achieve unified predictive capability for biomass handling behavior.

09 BIOMASS FUELS↗
DOE OSTI2021

DNA origami single crystals with Wulff shapes

DNA origami technology has proven to be an excellent tool for precisely manipulating molecules and colloidal elements in a three-dimensional manner. However, fabrication of single crystals with well-defined facets from highly programmable, complex DNA origami units is a great challenge. Here, we report the successful fabrication of DNA origami single crystals with Wulff shapes and high yield. By regulating the symmetries and binding modes of the DNA origami building blocks, the crystalline shapes can be designed and well-controlled. The single crystals are then used to induce precise growth of an ultrathin layer of silica on the edges, resulting in mechanically reinforced silica-DNA hybrid structures that preserve the details of the single crystals without distortion. The silica-infused microcrystals can be directly observed in the dry state, which allows meticulous analysis of the crystal facets and tomographic 3D reconstruction of the single crystals by high-resolution electron microscopy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Real-time plasma equilibrium reconstruction and shape control for the MAST Upgrade tokamak

Real-time magnetic control has been developed to deliver precise control of multiple plasma shape parameters for advanced divertor configurations, including double-null, Super-X, X-point target and X-divertor for the first time on the MAST Upgrade (MAST-U) spherical tokamak. Successful real-time magnetic equilibrium control of different plasma shape variables has been accomplished in the 2022–2023 MAST-U experimental campaign for the advanced MAST-U divertor configurations. Application of the MAST-U boundary reconstruction algorithm, LEMUR, is described and compared with off-line equilibrium reconstruction and diagnostic measurements. The process of design and verification of the axisymmetric magnetic control schemes using a suite of control analysis tools (known collectively as TokSys) is also described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗