Search NASA⌕ Search

SEARCH · Search NASA

Results for “self configuration”

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

Small tensor product distributed active space (STP-DAS) framework for relativistic and non-relativistic multiconfiguration calculations: Scaling from 10 9 on a laptop to 10 12 determinants on a supercomputer

Despite the power and flexibility of configuration interaction (CI) based methods in computational chemistry, their broader application is limited by an exponential increase in both computational and storage requirements, particularly due to the substantial memory needed for excitation lists that are crucial for scalable parallel computing. Here, the objective of this work is to develop a new CI framework, namely, the small tensor product distributed active space (STP-DAS) framework, aimed at drastically reducing memory demands for extensive CI calculations on individual workstations or laptops, while simultaneously enhancing scalability for extensive parallel computing. Moreover, the STP-DAS framework can support various CI-based techniques, such as complete active space (CAS), restricted active space, generalized active space, multireference CI, and multireference perturbation theory, applicable to both relativistic (two- and four-component) and non-relativistic theories, thus extending the utility of CI methods in computational research. We conducted benchmark studies on a supercomputer to evaluate the storage needs, parallel scalability, and communication downtime using a realistic exact-two-component CASCI (X2C-CASCI) approach, covering a range of determinants from 10 9 to 10 12 . Additionally, we performed large X2C-CASCI calculations on a single laptop and examined how the STP-DAS partitioning affects performance.

Complete-active space self-consistent field↗
DOE OSTI2025

Methods to Observe Tribological Failures in Self-Mated Steel Contacts

Scuffing, a type of wear found in highly stressed or poorly lubricated contacts, is characterized by a rapid increase in friction and severe plastic deformation of the near-surface material. Scuffing has proven difficult to study because it initiates unpredictably, progresses rapidly, and typically develops within an inaccessible contact interface. Although there have been successful in-situ studies of scuffing in real-time, the transparent counter body needed for these studies changes the interactions between the surfaces and the lubricant, which affects the scuffing process in unknown ways. This paper describes the development of X-ray-compatible tribometry to study the scuffing of self-mated steels in-situ and in real-time. The method uses a crossed cylinders configuration with a thin (500 μm thick) stationary component and a small (≈200 μm) contact width to maximize X-ray interactions with atoms within the stress field generated by the contact. The resulting instrument and method are used to benchmark the scuffing response of self-mated 52,100 steel under tribologically challenging ‘oil-off’ lubrication conditions. The results demonstrate reliable scuffing in this configuration despite the relatively small contact areas and loads used. Following scuffing, gross plastic deformation was observed on both surfaces along with significant subsurface grain refinement and flow only on the stationary surface, which experienced constant contact. Interestingly, high friction initiated at specific locations of the migratory surface, which experienced intermittent contact, and then propagated across the track over time, suggesting that local conditions of the migratory surface dominated friction leading into the failure event.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development of a second COMB dipole magnet with STAR wires

Rare-Earth Barium Copper Oxide (REBCO) coated conductors are an attractive option for application in high-field accelerator magnets due to their high critical field and the convenience of fabrication without heat treatment compared to some other superconductors. Fermilab plays a vital role in superconducting accelerator magnet R&D under the framework of the U.S. Magnet Development Program. An integral part of that program is the accelerator magnet development based on REBCO conductors to demonstrate self-fields of 5 T or greater compatible with operation in hybrid configurations to generate fields beyond 16 T for future High Energy Physics (HEP) applications. A small, two-layer REBCO dipole magnet based on the Conductor on Molded Barrel (COMB) magnet technology, developed at Fermilab with Symmetric Tape Round (STAR®) wires from AMPeers LLC was recently fabricated and tested in liquid helium, demonstrating the design feasibility. A second dipole magnet based on the same technology is currently under development at Fermilab. The coil has a 60-mm clear bore and a 120-mm OD and consists of six layers of STAR® conductor, which allows to significantly increase the magnetic field compared to the previous design, and potentially probe the 5 T self-field target. This presentation reports the progress in the magnet development.

Kashikhin, Vadim [Fermilab]↗
DOE OSTI2024

Comment on ‘Resonant inelastic x-ray scattering from U 3 O 8 and UN’

In a recent manuscript, Lawrence Bright et al (2023 J. Phys.: Condens. Matter 35 175501) reported the resonant inelastic x-ray scattering spectra of U 3 O 8 , as well as UN. Their goal was to identify electronic multiplets associated with a 5f 1 configuration with ground state 2 F 5/2 . Complete active space self-consistent field with spin–orbit coupling (CASSCF-SOC) predicted that 2 F 5/2 transitions should be observable at 190 and 328 meV. However, these energies were not accessible in their experiment. Here they suggested that the recent inelastic neutron scattering results of Miskowiec et al (2021 Phys. Rev. B 103 205101) could have been sensitive to these transitions. Here we show that transitions of this possible origin appear in that dataset near 198, 262, 362, and potentially 448 meV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Strategies to Obtain Reliable Energy Landscapes from Embedded Multireference Correlated Wavefunction Methods for Surface Reactions

Embedded correlated wavefunction (ECW) theory is a powerful tool for studying ground- and excited-state reaction mechanisms and associated energetics in heterogeneous catalysis. Several factors are important to obtaining reliable ECW energies, critically the construction of consistent active spaces (ASs) along reaction pathways when using a multireference correlated wavefunction (CW) method that relies on a subset of orbital spaces in the configuration interaction expansion to account for static electron correlation, e.g., complete AS self-consistent field theory, in addition to the adequate partitioning of the system into a cluster and environment, as well as the choice of a suitable basis set and number of states included in excited-state simulations. Here, in this work, we conducted a series of systematic studies to develop best-practice guidelines for ground- and excited-state ECW theory simulations, utilizing the decomposition of NH 3 on Pd(111) as an example. We determine that ECW theory results are relatively insensitive to cluster size, the aug-cc-pVDZ basis set provides an adequate compromise between computational complexity and accuracy, and that a fixed-clean-surface approximation holds well for the derivation of the embedding potential. Additionally, we demonstrate that a merging approach, which involves generating ASs from the molecular fragments at each configuration, is preferable to a creeping approach, which utilizes ASs from adjacent structures as an initial guess, for the generation of consistent potential energy curves involving open-d-shell metal surfaces, and, finally, we show that it is essential to include bands of excited states in their entirety when simulating excited-state reaction pathways.

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

Real-Time GPU-Accelerated OFDR With an Integrated Auxiliary Interferometer

A GPU-accelerated optical frequency domain reflectometry (OFDR) system with an improved integrated auxiliary interferometer is proposed. Unlike conventional approaches that require separate auxiliary interferometers and multiple detection channels, the proposed OFDR system embeds this functionality directly into the signal via an intentional beat component. This enables self-calibration of laser nonlinearity while maintaining a cost-effective hardware configuration. Building on this simplified configuration, the system leverages GPU acceleration with an NVIDIA RTX 4070 Ti to achieve real-time performance, delivering high-throughput signal processing for continuous OFDR interrogation. The signal processing pipeline comprises signal capture, resampling for nonlinearity compensation, and frequency shift computation, all optimized for parallel execution. Hardware benchmarking demonstrates substantial acceleration over CPU implementations, achieving up to a 45× speedup for resampling and frequency shift computations and enabling processing latencies below 30 ms. Thermal response validation is conducted under two complementary scenarios: localized heating using a water bath and cryogenic-temperature conditions using liquid nitrogen. Under localized heating, the system achieves an accuracy of 0.249 °C with a thermal sensitivity of 5.971 GHz/°C, while cryogenic-temperature validation demonstrates a frequency shift response with a sensitivity of 2.383 GHz/°C and an accuracy of 2.04 °C. The high acceleration of the proposed GPU-accelerated OFDR system and its accuracy are achieved by exploiting CUDA-based stride indexing, enabling efficient parallel segmentation and processing of large datasets without additional memory copies. The benchmarking results confirm the robustness, accuracy, and deployability of the proposed OFDR system across a wide temperature range, establishing it as a practical platform for real-time distributed fiber sensing in structurally dynamic environments.

Harb, Salah [Lawrence Berkeley National Laboratory↗
DOE OSTI2026

AutoLabs: cognitive multi-agent systems with self-correction for autonomous chemical experimentation

The automation of chemical research through self-driving laboratories (SDLs) promises to accelerate scientific discovery, yet the reliability and granular performance of the underlying AI agents remain critical, under-examined challenges. In this work, we introduce AutoLabs, a self-correcting, multi-agent architecture designed to autonomously translate natural-language instructions into executable protocols for a high-throughput liquid handler. The system engages users in dialogue, decomposes experimental goals into discrete tasks for specialized agents, performs tool-assisted stoichiometric calculations, and iteratively self-corrects its output before generating a hardware-ready file. We present a comprehensive evaluation framework featuring five benchmark experiments of increasing complexity, from simple sample preparation to multi-plate timed syntheses. Through a systematic ablation study of 20 agent configurations, we assess the impact of reasoning capacity, architectural design (single- vs. multi-agent), tool use, and self-correction mechanisms. Our results demonstrate that agent reasoning capacity is the most critical factor for success, reducing quantitative errors in chemical amounts (nRMSE) by over 85% in complex tasks. When combined with a multi-agent architecture and iterative self-correction, AutoLabs approaches expert-authored reference procedures on the benchmark (F1-score > 0.89) on challenging multi-plate syntheses. These findings establish a clear blueprint for developing robust and trustworthy AI partners for autonomous laboratories, highlighting the synergistic effects of modular design, advanced reasoning, and self-correction to ensure both performance and reliability in high-stakes scientific applications. Code: https://github.com/pnnl/autolabs

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

Modeling Strong Light-Matter Coupling in Correlated Systems: State-Averaged Cavity Quantum Electrodynamics Complete Active Space Self-Consistent Field Theory

The description of strongly correlated systems interacting with quantized cavity modes poses significant theoretical challenges due to the combinatorial scaling of electronic and photonic degrees of freedom. Recent advances addressing this complexity include cavity quantum electrodynamics (QED) generalizations of complete active space configuration interaction and density matrix renormalization group methods. In this work, we introduce a QED extension of state-averaged complete active space self-consistent field theory, which incorporates cavity-induced correlations through a second-order orbital optimization framework with robust convergence properties. The method is implemented using both photon number state and coherent state representations, with the latter showing robust origin invariance in the energies regardless of the completeness of the photonic Fock space. The implementation enables symmetry-free orbital relaxations to account for photon-mediated symmetry breaking in polaritonic systems. Numerical validation on lithium hydride, hydroxide anion, and magnesium hydride cation demonstrates that this method achieves significantly improved accuracy in modeling ground-state and polariton potential energy surfaces compared to QED-CASCI in a fixed orbital basis. In these studies, we reach sub-kcal/mol accuracy in potential energy surface in much smaller active spaces than are required for QED-CASCI. This advancement provides a more robust approach for studying cavity-altered chemical landscapes for ground and exited strongly coupled systems.

CASSCF↗
DOE OSTI2024

Fabrication and test of COMB dipole magnet with STAR® wires

Fermilab is working on superconducting accelerator magnet R&D under the framework of the U.S. MagnetDevelopment Program (USMDP). An integral part of that program is High Temperature Superconducting (HTS)accelerator magnet development to- demonstrate HTS self‐fields of 5 T or greater- be compatible with operation in hybrid LTS/HTS configuration• to probe fields beyond 16 T for future High Energy Physics (HEP) experiments

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Self-consistent solution of the Frank–Bilby equation for interfaces containing disconnections

The quantized Frank–Bilby equation can be used to identify interfacial line defect array configurations which relax the misorientation and/or misfit of a coherent crystalline interface. These line defect arrays may be comprised of dislocations and/or disconnections, which are interfacial steps with dislocation character. When an interface contains disconnections, solution of the quantized Frank–Bilby equation is complicated by the fact that the habit plane orientation is not known in advance because it depends on the unknown spacing of the disconnection array. We present a root-finding-based method for addressing this issue, enabling a self-consistent solution for arbitrary defect content. Our method has been implemented in an open-source code which enumerates all possible solutions given a list of candidate line defects. Two cases are presented employing the code: a misoriented FCC twin boundary and an FCC/BCC phase boundary with the Nishiyama-Wasserman orientation relationship. Both cases exhibit more than 10,000 solutions to the Frank–Bilby equation, with several hundred solutions categorized as ‘‘low energy’’ and thus plausible configurations for the actual interface. The resulting set of solutions can be utilized to predict and understand the properties of a given interface.

42 ENGINEERING↗
DOE OSTI2025

SICAS, a new integrated divertor-SOL-pedestal-core framework for self-consistent modeling of ion and impurity transport in fusion devices

We present a new framework for core-edge integration studies named SICAS (SOLPS-ITER coupled to ASTRA-STRAHL) which enables high fidelity simulations of the core, edge, and divertor regions encompassing the transport of ions as well as the impurities through the entire plasma domain. SICAS handles the exchanging of the particle and power fluxes as well as transport coefficients to ensure consistency through the codes. An overlapping region is defined from the inside of the separatrix of the plasma to the SOLPS-ITER core boundary allowing for matching profiles and fluxes between the two codes for a self-consistent approach. The results presented here demonstrated the flexibility of SICAS to simulate different configurations, scenarios, divertor geometries, and plasma species with good agreement with DIII-D experimental data. This tool opens new possibilities in integrated modeling of fusion devices integrating all relevant phenomena in the core and the divertor plasmas. These capabilities are required for the interpretation of current experiments as well as the design of new devices.

core-edge integration↗
DOE OSTI2025

Interactive Gas Chemistry for Enhanced Science Capabilities of the Energy Exascale Earth System Model Version 3

Atmospheric chemistry plays a crucial role in Earth system models (ESMs), controlling atmospheric composition and radiative balance; it is highly interactive with the physical climate, biogeochemical cycles, and human systems. However, it often imposes computational challenges in an ESM. Here we develop a full troposphere‐stratosphere interactive chemistry module for the US Department of Energy's Energy Exascale Earth System Model (E3SM). We intentionally build a streamlined module based on E3SM version 2 that interacts with other components and maintains all of major chemical and chemistry‐climate feedbacks. The module incorporates a new, highly efficient tracer advection scheme; linearization of stratospheric chemistry; and abridged tropospheric chemical mechanism with 28 reactive tracers. This new model, E3SM‐chem, can readily perform century‐long climate simulations of ozone, methane, and nitrous oxide based on emission scenarios as well as provide hourly budgets for the gas‐phase radicals that drive aerosol chemistry. We evaluate E3SM‐chem with an atmosphere‐only simulation as in the recent climate model intercomparison project (CMIP6) finding results similar to the other CMIP6 models. For the present‐day, E3SM‐chem matches the standard measurement metrics for stratospheric and tropospheric ozone, surface air quality, other key reactive gases like carbon monoxide, and the methane lifetime. Overall, E3SM‐chem maintains the climate fidelity of the baseline model while adding at most 20% to the computational cost of the atmosphere model. Hence, interactive chemistry can be a default configuration for long climate simulations at resolutions of 1° or finer, which is crucial for producing self‐consistent chemistry‐climate feedbacks that alter the climate system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Interpreting Mass and Radius Measurements of Neutron Stars with Dark Matter Halos

Abstract The high densities of neutron stars (NSs) could provide astrophysical locations for dark matter (DM) to accumulate. Depending on the DM model, these DM admixed NSs (DANSs) could have significantly different properties than pure baryonic NSs, accessible through X-ray observations of rotation-powered pulsars. We adopt the two-fluid formalism in general relativity to numerically simulate stable configurations of DANSs, assuming a fermionic equation of state (EOS) for the DM with repulsive self-interaction. The distribution of DM in the DANS as a halo affects the path of X-rays emitted from hot spots on the visible baryonic surface, causing notable changes in the pulse profile observed by telescopes such as NICER, compared to pure baryonic NSs. We explore how various DM models affect the DM mass distribution, leading to different types of dark halos. We quantify the deviation in observed X-ray flux from stars with each of these halos. We identify the pitfalls in interpreting mass and radius measurements of NSs inferred from electromagnetic radiation and constraining the baryonic matter EOS if these dark halos exist.

Shawqi, Shafayat (ORCID:0000000210956183)↗
DOE OSTI2026

Active space selection with self-healing diffusion Monte Carlo algorithms for periodic solids

Multideterminant Diffusion Monte Carlo (DMC) displays improved accuracy over single determinant DMC. Self-Healing Diffusion Monte Carlo (SHDMC) is a DMC based method that iteratively improves a multideterminant trial wavefunction. Although configuration interaction or complete active space (CAS) methods are very accurate and computationally feasible for many systems, they are not optimal for application to solids. SHDMC is accurate and designed for application to solids, so developing SHDMC based active space selection algorithms is a worthy endeavor. Here, we present and compare active space selection algorithms that are designed for use in conjunction with SHDMC, without relying on external approaches. For benchmarking, we calculated the ground state energy of a small unit cell of graphene and compared the results with a complete basis set extrapolated selected CI and a reference SHDMC trajectory. We found that systematically expanding the active space using an “auto-branching” algorithm optimally balances accuracy with computational practicality. To the best of our knowledge, this is the first work that demonstrates completely self-contained DMC-based active space selection algorithms that do not depend on external methods for determinant selection.

Spanedda, Nicole [ORNL]↗
DOE OSTI2026

Assessing the Impacts of Conformational Fluxionality on Copper(II/I) Electron Transfer Self-Exchange

Typical Cu(II/I) complexes exhibit hallmark structural changes during their electron transfer (ET) reactions that result from the (pseudo) Jahn-Teller distortions and changes in polarizability inherent to their d 10 / d 9 configurations. Given that such structural changes incur large reorganization energy penalties, slow the rates of ET are characteristic of these compounds. Notwithstanding, we recently reported a set of Cu(II/I) complexes that undergo significant and well-defined structural changes during their redox reactions yet exhibit rapid (> 10 5 M ‒1 s ‒1 ) ET self-exchange rate constants ( k 11 ). To explain these results we proposed a pre-equilibrium model in which inherent conformational fluxionality in one of the two oxidation states provides access to pathways involving lower reorganization energies during the ET event. Herein we report our results testing this hypothesis through the preparation and study of a homologous series of compounds exhibiting varying extents of con-formational fluxionality in the Cu(I) state. We characterize these compounds electrochemically, structurally, and by varia-ble temperature NMR spectroscopy to provide experimental evidence for increase fluxionality across the series. We then correlate the trend with increasing k 11 through NMR linewidth broadening experiments, further taking care to define the impacts of solvent impurities therein. Finally, the nature of the conformational rearrangements, and their relation to increased k 11 are explored computationally to reveal differential Boltzmann populations of conformers across the series. The cumulative results of these studies support a previously underappreciated strategy for overcoming barriers to slow ET kinetics: via the incorporation of conformational fluxionality.

Singh, Aditi↗
DOE OSTI2024

A molecular ground electronic state with an occupied 5g spinor—The superheavy (E125)F molecule

Fully relativistic calculations, primarily at the 4-component coupled-cluster singles and doubles with perturbative triples [CCSD(T)] level of theory with the Dirac–Coulomb (DC) Hamiltonian, have been carried out for the superheavy (E125)F molecule using large Gaussian basis sets. The electronic ground state is determined to have an [Og]8s 2 5g 1 6f 3 configuration on E125 with an Ω = 6 ground state and an 8p electron largely donated to F. A Mulliken population analysis indicates that the ground state is mainly ionic with a partial charge of +0.79 on E125 and a single sigma bond involving the F 2p and E125 8p spinors. The occupied g spinor is not involved in the bonding. With the largest basis set used in this work, the (0 K) dissociation energy was calculated at the DC-CCSD(T) level of theory to be 7.02 eV. Analogous calculations were also carried out for the E125 atom, both the neutral and its cation. The lowest energy electron configuration of E125 + , [Og]$8s$$^{2}_{1/2}$$5g$$^{1}_{7/2}$$6f$$^{3}_{5/2}$ with a J = 6 ground state, was found to be similar to that in (E125)F, while the neutral E125 atom has an [Og]$8s$$^{2}_{1/2}$$5g$$^{1}_{7/2}$$6f$$^{2}_{5/2}$$7d$$^{1}_{3/2}$$8p$$^{1}_{1/2}$ ground state electron configuration with a J = 17/2 ground state. The ionization energy (IE) of E125 is reported for the first time and is calculated to be 4.70 eV at the DC-CCSD(T) level of theory. Non-relativistic calculations were also carried out on the E125 atom and the (E125)F molecule. Here, the non-relativistic ground state of the E125 atom was calculated to have a 5g 5 ground state with an IE of just 3.4 eV. The net effect of relativity on (E125)F is to stabilize its bonding.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Characterizing impurity sourcing and transport in the high temperature boundary of DIII-D wide pedestal QH-mode plasmas

Wide Pedestal QH-mode (WPQH) plasmas in the DIII-D tokamak show a sheath limited SOL, where electron density and temperature remain nearly constant along field lines from the midplane to the divertor. Consequently, parallel gradients are weak. The first Langmuir probe measurements in this regime point to a high sheath temperature with target electron temperature T e up to ~150 eV, leading to high carbon self-sputtering. Midplane carbon densities from SOLPS-ITER modeling of these plasmas in a double-null configuration fall far below experimental measurements unless full drifts are activated. In the drift-dependent SOLPS-ITER modeling, anomalous poloidally uniform, radially varying transport is adjusted to match measured radial electron temperature and density profiles in the pedestal and SOL region. The resulting carbon density (C6+) matches measured carbon densities just inside the separatrix. The C2+ density near the outer strike point is also consistent with spectroscopic imaging. With the ion B×∇B drift towards the X-point, carbon in the lower divertors is redistributed from the private flux region to the high field side (HFS) and pushed upstream in the SOL by poloidal E×B drifts. The B × ∇B drift dominates the radial flow of carbon as it moves upstream. Simulations with reversed toroidal field (ion B × ∇B drift away from the X-point) show a radically different behavior, where carbon accumulates on the low field side in lower divertors, and is pushed towards the HFS in the upper divertors, indicating a strong effect of particle drifts on impurity distribution. In double-null configurations which are usually used in WPQH plasmas, these drift effects tend to counterbalance, so that the carbon density and Z eff are reduced at the outer mid-plane only modestly, by around 7% in the modeling. Further predictive modeling indicates that the carbon density can be significantly reduced by an order of magnitude by implementing a single-null shape with reversed toroidal field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Towards fully predictive gyrokinetic full- f simulations: validation and triangularity studies in TCV

Designing economical magnetic confinement fusion power plants motivates computational tools that can estimate plasma behavior from engineering parameters without direct reliance on experimental measurement of the plasma profiles. In this work, we present full-f global long-wavelength gyrokinetic simulations of edge and scrape-off layer turbulence in tokamaks that use only magnetic geometry, heating power, and particle inventory as inputs. Unlike many modeling approaches that employ free parameters fitted to experimental data, raising uncertainties when extrapolating to reactor scales. This approach directly simulates turbulence and resulting profiles through gyrokinetics without such empirical adjustments. This is achieved via an adaptive sourcing algorithm in Gkeyll that strictly controls energy injection and emulates particle sourcing due to neutral recycling. We show that the simulated kinetic profiles compare reasonably well with Thomson scattering and Langmuir probe data for Tokamak á Configuration Variable (TCV) discharge #65125, and that the simulations reproduce characteristic features such as blob transport and self-organized electric fields. Applying the same framework to study triangularity effects suggests mechanisms contributing to the improved confinement reported for negative triangularity (NT). Simulations of TCV discharges #65125 and #65130 indicate that NT increases the E x B flow shear (by about 20% in these cases), which correlates with reduced turbulent losses and a modest change in the distribution of power exhaust to the vessel wall. While the physical models contain approximations that can be refined in future work, the predictive capability demonstrated here, evolving multiple profile relaxation times with kinetic electron and ion models in hundreds of GPU hours, indicates the feasibility of using Gkeyll to support design studies of fusion devices.

Hoffmann, Antoine Cyril David [Princeton Plasma Ph↗