Search NASA⌕ Search

SEARCH · Search NASA

Results for “minimum”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Directional dependence of equilibrated TLD-400 chips in multiple radiation fields

Fielding of thermoluminescent dosimeters (TLDs) for measurement of photon radiation dose in experiments is the standard practice for γ-irradiation facilities, pulsed power x-ray facilities, and reactor facilities at Sandia National Laboratories (SNL). Due to the high-dose experimental conditions and the mixed ( 1 n, γ) fields in these facilities, SNL radiation metrologists have historically used CaF 2 :Mn TLDs (also known as TLD-400). Recent inquiries to the radiation metrology staff have raised concerns that the aluminum-equilibrated TLD-400 chips may exhibit a directionally dependent response. The metrologists were asked whether the dose measured by the chip may be impacted by the angle of incidence on the equilibrated TLD. To provide a thorough answer to this query, a set of adjoint Monte Carlo radiation transport calculations was performed for three different equilibrated TLDs as well as bare TLD chips using the Integrated Tiger Series (ITS) code. The key feature of the adjoint calculations performed in this study is that each photon escaping the modeled geometry was tallied into angular bins to provide information on both the energy and angular dependence of the equilibrated chip. After the adjoint Monte Carlo calculations were completed, the resulting energy-dependent response function for each angular bin was convolved with multiple photon energy spectra representing various radiation facilities at SNL. The radiation facilities selected for analysis span a range of photon energies from approximately 1 keV up to approximately 20 MeV. Thus, the presented results are applicable to a wide variety of radiation facilities around the world. Although the bare TLD-400 chip was expected to display the largest variation due to the photon angle of impact on the dosimeter, the dosimeter with the thinnest aluminum equilibrator (SNL thin equilibrated TLD-400) was determined to have the biggest differential between the impact angle with the maximum dose ( D max ) and the impact angle with the minimum dose ( D min ). However, the SNL normal equilibrated TLD-400 chip demonstrated a dramatic reduction in that differential between maximum and minimum dose angles. The reduction in this differential is one of the dominant factors in experimenters’ choice to field these dosimeters at SNL radiation facilities. The results from the PNNL (Hanford) energy-flattening field capsule are consistent with the previous publications. The differential dose responses due to the photon impact angle indicate that experimenters should strive to field their TLD-400 dosimeters in a consistent manner to avoid additional uncertainty in the measurements based on dosimeter orientation.

Adjoint radiation transport↗
DOE OSTI2026

Thermomechanical analysis and modeling of a high-temperature light trapping planar cavity solar receiver

The development of durable particle-based high-temperature solar receivers is critical for advancing concentrating solar-thermal (CST) technologies to enable high-efficiency power generation and industrial process heat. Here, this study presents a computational framework to evaluate the thermomechanical performance of a proposed enclosed light-trapping planar cavity receiver designed for particle-based thermal energy systems. The receiver incorporates absorptive cavities and fluidized particle-bed channels to enhance heat capture and reduce thermal losses. Finite element analysis (FEA) is employed to assess stress, strain, and creep-fatigue behavior under concentrated solar flux using realistic thermal boundary conditions derived from coupled system models and experimental assembly parameters. The analysis investigates the influence of particle-to-wall heat transfer coefficients (HTC) ranging from 800 to 1800 W/m 2 .K on the thermomechanical response of six candidate high-temperature alloys: Alloy 740H, Alloy 282, Alloy 617, 316H, Alloy 230, and 800H. Results show that increasing HTC reduces thermal gradients, leading to lower stresses and strains and extended minimum predicted creep life. While all materials satisfy fatigue life requirements under the investigated conditions, significant differences in creep resistance are observed. Alloy 740H consistently exhibits the longest minimum predicted creep life and the most favorable durability margins, followed by Alloy 282, with the remaining materials showing reduced creep resistance under identical loading. The reported creep lives are conservative lower-bound estimates intended for comparative material evaluation. This framework highlights the critical roles of material selection and geometry optimization in improving mechanical durability and reliability of solar-thermal receivers, forming a foundation for future experimental validation and design optimization.

14 SOLAR ENERGY↗
DOE OSTI2025

Demonstration and technoeconomic analysis of dodecanol production from acetate using metabolically engineered Escherichia coli

In a circular bioeconomy, the one-way conversion of petroleum to chemicals and CO 2 is replaced with processes that reduce CO 2 to energy carriers and useful materials that are returned to CO 2 upon combustion. A circular bioeconomy that relies on photosynthesis to generate sugars as the chief energy carrier and precursor to chemical building blocks has yet to overcome many recalcitrant aspects of plant-based photosynthesis, namely, high feedstock costs, arable land scarcity, food competition, and fertilizer overuse. Acetate is a potential sustainable energy carrier because it can be produced from CO 2 either electrocatalytically or by acetogens via the Wood-Ljungdahl pathway. Here, in this work, we conducted a metabolic engineering study of Escherichia coli 's ability to convert acetate into dodecanol as a model oleochemical product. We performed techno-economic and life cycle analyses to determine break-even points with alternative fossil fuel-based strategies and identified critical process performance parameters for supporting an industrial acetate-based bioprocess. These analyses showed that oleochemical yield is the primary driver of minimum oleochemical selling price and carbon intensity. Therefore, to increase yield on acetate, we deleted the aceBAK operon, which facilitates funneling of acetate into biomass instead of product. We performed additional strain engineering to increase flux towards dodecanol and increase acetate uptake. Finally, we demonstrated increased yield in controlled bioreactors, improving from 13% of the maximum theoretical yield to 37%. Rigorous uncertainty analyses assuming a range of market conditions and future technological performances resulted in 88% and 37% of simulated scenarios having lower carbon intensities than fossil fuel-based routes and lower minimum selling prices than the market price.

Acetate↗
DOE OSTI2025

Metastability of stratified magnetohydrostatic equilibria and their relaxation

Motivated by explosive releases of energy in fusion, space and astrophysical plasmas, we consider the nonlinear stability of stratified magnetohydrodynamic equilibria against two-dimensional interchanges of straight magnetic-flux tubes. We demonstrate that, even within this restricted class of dynamics, the linear stability of an equilibrium does not guarantee its nonlinear stability: equilibria can be metastable. We show that the minimum-energy state accessible to a metastable equilibrium under non-diffusive two-dimensional dynamics can be found by solving a combinatorial optimisation problem. These minimum-energy states are, to good approximation, the final states reached by our simulations of destabilised metastable equilibria for which turbulent mixing is suppressed by viscosity. To predict the result of fully turbulent relaxation, we construct a statistical mechanical theory based on the maximisation of Boltzmann's mixing entropy. This theory is analogous to the Lynden-Bell statistical mechanics of collisionless stellar systems and plasma, and to the Robert–Sommeria–Miller theory of two-dimensional vortex turbulence. Our theory reproduces well the results of our numerical simulations for sufficiently large perturbations to the metastable equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

(NH 3 (CH 2 ) 7 NH 3 ) 2 Sn 3 I 10 , a Vacancy-Ordered Three-Dimensional Tin(II) Perovskite-Derived Semiconductor

Ordering vacancies in hybrid Sn(II) halide semiconductors provides a strategy for preventing uncontrolled oxidation and formation of mobile holes. In this study, we report the structure and optical and electronic properties of (NH 3 (CH 2 ) 7 NH 3 ) 2 Sn 3 I 10 , a vacancy-ordered perovskite derivative with three-dimensional inorganic connectivity. The crystal structure resembles that of a Dion–Jacobson layered perovskite derivative, but with [SnI 5 ] square pyramids bridging the layers. UV–vis diffuse reflectance spectroscopy reveals a sharp onset of light absorption at 1.86(1) eV with the photoluminescence emission maximum at 1.90(1) eV. However, the maximum excitation occurs from 3.42 to 3.81 eV (325 to 370 nm), revealing a significant Stokes shift of 1.3 eV. The electronic properties determined from dark and time-resolved microwave conductivity measurements reveal a minimum carrier mobility of 4.3 × 10 –2 cm 2 V –1 s –1 and a maximum carrier density of 5.96 × 10 16 cm –3 , a uniquely low value for a hybrid Sn(II) halide semiconductor. The transport behavior in combination with first-principles calculations of the electronic band structure and dielectric permittivity suggest polaron-mediated electronic transport, yet the photogenerated carriers have a fast and fluence-dependent nonradiative recombination rate, suggestive of localized “defect-like” states at the band edge. The observed photoluminescence is most consistent with single-ion-like behavior of an asymmetric Sn(II) environment. Together, these results suggest that defect ordering presents a strategy for the reduction of mobile charge carriers at equilibrium.

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

Navigating Large Chemical Spaces Using Graph Theory and Integer Programming

Navigating and analyzing large chemical spaces are necessary to accelerate the design and discovery of new molecules and chemical processes. In this work, we introduce a computational framework that integrates graph theory and integer programming to enable the efficient navigation of large chemical spaces. Our framework represents the chemical space as a graph, wherein nodes represent molecules and edges represent the degree of similarity or connectivity based on domain-specific information. Using the graph representation, we identify representative molecules by computing the so-called minimum dominating set (MDS), which in our context is the minimum set of molecules that is connected to all other molecules. We present a suite of solution strategies for the MDS problem including heuristic and rigorous integer programming (IP) approaches. We show that these approaches allow us to capture physicochemical properties and domain-specific logic and constraints, facilitating the identification of molecules with the target properties. We demonstrate the effectiveness of the proposed approach by navigating the chemical space of per- and polyfluoroalkyl substances (PFAS); this comprises approximately 15,000 molecular structures. We compare our framework against traditional dimensionality reduction and clustering methods such as t-SNE and K-means clustering.

Chemical structure↗
DOE OSTI2026

Photoelectron Spectroscopic Determination of the Interfacial Energetics of Metal Oxide Protection Layers on p-InP Photocathodes

The interfacial energetics between p-type InP and a series of metal oxides, including TiO 2 , Nb 2 O 5 , Ta 2 O 5 , and HfO 2 , were evaluated using X-ray photoelectron spectroscopy (XPS), ultraviolet photoelectron spectroscopy (UPS), and optical absorption spectroscopy. The energy of the conduction-band minimum ( E cb ) of TiO 2 and Nb 2 O 5 was more negative (i.e., further from the vacuum level) than the conduction-band minimum at the surface of InP ( E cb,s ), whereas E cb for Ta 2 O 5 and HfO 2 was more positive than E cb,s for InP. The data are consistent with the electrochemical behavior of p-InP coated with various metal oxide candidate protection layers, with TiO 2 and Nb 2 O 5 facilitating interfacial transfer of photogenerated minority-carrier electrons in p-InP photocathodes, and Ta 2 O 5 and HfO 2 blocking photogenerated electrons in p-InP from readily transferring across the oxide-coated photocathodes. The energy of the valence-band maximum ( E vb ) for all of the oxides was much more negative than E vb,s for InP, consistent with observations that these protection layers effectively block hole transport and consequently suppress oxidative degradation of the underlying p-InP photocathodes.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2025

Characterizing the Potential for Sustainable Azelaic Acid Production from High-Oleic Vegetable Oil Using Two-Step Oxidative Cleavage

Azelaic acid is a renewable monomer conventionally produced via the energy-intensive ozonolysis of oleic acid. Recent advancements have enabled the use of high-oleic vegetable oils (rather than tallow-derived oleic acid) and replaced ozonolysis with two-step oxidative cleavage using hydrogen and oxygen. Although this shift would improve process safety, the financial viability and environmental implications remain uncertain. In this study, we characterized the sustainability of azelaic acid production from high-oleic vegetable oil using two-step oxidative cleavage. Process design, simulation, technoeconomic analysis (TEA), and life cycle assessment (LCA) were executed under uncertainty using BioSTEAM. The modeled system produces azelaic acid at a market-competitive minimum selling price (MSP) of 8.32 [4.93−13.34] $\$kg$ −1 (median 5th−95th percentiles), below the minimum estimated market price of 9.93 $\$kg$ −1 . Further, it has the potential to approach carbon neutrality (0.0 [−5.5 to 5.6] kg of CO 2 -eq kg −1 ) under displacement allocation. Improvements to dihydroxylation (86 to 99%) and oxidative cleavage conversions (93 to 99%) would reduce MSP to $\$5.24$ kg −1 and carbon intensity to −1.90 kg of CO 2 -eq kg −1 (displacement). Additionally, increasing the feedstock triolein content (75 to 85%) lowers MSP by $\$0.82$ kg −1 . Overall, this research demonstrates the potential for financially viable production of azelaic acid from vegetable oils and the utility of agile TEA/LCA.

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

Climate, Hydrology, and Nutrients Control the Seasonality of Si Concentrations in Rivers

Abstract The seasonal behavior of fluvial dissolved silica (DSi) concentrations, termedDSi regime, mediates the timing of DSi delivery to downstream waters and thus governs river biogeochemical function and aquatic community condition. Previous work identified five distinct DSi regimes across rivers spanning the Northern Hemisphere, with many rivers exhibiting multiple DSi regimes over time. Several potential drivers of DSi regime behavior have been identified at small scales, including climate, land cover, and lithology, and yet the large‐scale spatiotemporal controls on DSi regimes have not been identified. We evaluate the role of environmental variables on the behavior of DSi regimes in nearly 200 rivers across the Northern Hemisphere using random forest models. Our models aim to elucidate the controls that give rise to (a) average DSi regime behavior, (b) interannual variability in DSi regime behavior (i.e., Annual DSi regime), and (c) controls on DSi regime shape (i.e., minimum and maximum DSi concentrations). Average DSi regime behavior across the period of record was classified accurately 59% of the time, whereas Annual DSi regime behavior was classified accurately 80% of the time. Climate and primary productivity variables were important in predicting Average DSi regime behavior, whereas climate and hydrologic variables were important in predicting Annual DSi regime behavior. Median nitrogen and phosphorus concentrations were important drivers of minimum and maximum DSi concentrations, indicating that these macronutrients may be important for seasonal DSi drawdown and rebound. Our findings demonstrate that fluctuations in climate, hydrology, and nutrient availability of rivers shape the temporal availability of fluvial DSi.

Environmental Sciences & Ecology↗
DOE OSTI2025

Global heliospheric termination shock strength in the solar–interstellar interaction

A heliospheric termination shock (HTS) surrounds our Solar System at approximately 100 astronomical units from the Sun, where the expanding solar wind (SW) is compressed and heated before encountering the interstellar medium. HTS-accelerated particles govern the pressure balance with the interstellar medium, but little is known about the global properties of the HTS beyond in situ measurements from Voyager in only two directions of the sky. Here we fill this gap by extracting the HTS strength using particle-in-cell, test particle and magnetohydrodynamic simulations, constrained by Interstellar Boundary Explorer observations of energetic neutral atoms produced from HTS-accelerated particles. Our results reveal there is a higher compression near the poles during solar minimum compared with solar maximum due to the higher Mach number flow. North–south asymmetries arise from the disparate evolution of the polar coronal holes, while minimum compression near the flanks is probably due to SW slowing from mass loading over a greater distance to the HTS. The results imply a strong connection between the HTS strength and the SW and interstellar medium dynamics.

Astrophysical plasmas↗
DOE OSTI2025

Dislocation-Grain Boundary Interaction Dataset for FCC Cu

Interactions between dislocations and grain boundaries play a major role in controlling the strength and ductility of structural materials. Experimentally, assessing and probing geometric and stress-based criteria at the local level for dislocation transmission through grain boundaries remains challenging. Therefore, there have been many efforts to systematically generate datasets of dislocation-grain boundary interactions (DGI) via computational models such as molecular dynamics simulations. So far, most DGI datasets have focused only on the subset of nominal minimum-energy grain boundary structures, which limits their applicability, especially to materials processed far from equilibrium. We present a comprehensive database of dislocation-grain boundary interactions for edge, screw, and 60° mixed dislocation with 330 <110> and 257 <112> symmetric tilt grain boundaries (total of 587) in FCC Cu consisting of 73 minimum-energy grain boundary structures and 514 metastable structures. The dataset contains the outcomes for 5234 unique interactions for various dislocation types, grain boundary structures, and applied shear stresses.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Enhanced energy gain through higher-order resonances during direct laser acceleration with superluminal phase velocity

Ultra-high intensity laser–plasma interactions can produce ultra-relativistic electrons via direct laser acceleration, assisted by quasi-static plasma magnetic and electric fields. These fields transversely confine electron motion and induce betatron oscillations. The net energy gain is strongly influenced by the interplay between two frequencies: the betatron frequency and the frequency of laser field oscillations experienced by the electron. Prior work has shown that energy gain is enabled by a resonance between the betatron oscillations and the oscillations of the laser field. In particular, higher-order resonances occur when the laser field completes multiple cycles during one betatron oscillation, allowing additional regimes of energy transfer beyond the fundamental (betatron) resonance. In this work, we demonstrate that such resonances become particularly effective when the laser's phase velocity is superluminal. Although the two frequencies generally evolve differently with increasing electron energy—leading to detuning—a superluminal phase velocity introduces a non-monotonic frequency ratio with a global minimum. This minimum allows sustained frequency matching over a broad energy range, thereby enabling enhanced energy gain. As the phase velocity increases, the betatron resonance becomes ineffective due to premature frequency detuning. At the same time, higher-order resonances become increasingly effective, emerging as the dominant mechanisms for enhanced energy gain in this regime of direct laser acceleration.

Laser plasma interactions↗
DOE OSTI2025

On resistive interchange, double tearing, and resonant and non-resonant infernal modes in spherical tokamak negative central shear discharges

We examine the linear and nonlinear stability of a sequence of reversed shear NSTX equilibria with the same toroidal current density and pressure profile but with different toroidal field strengths. All equilibria have two q = 2 surfaces and have q > 1 everywhere. For this sequence, which all have β P ∼ 0.5, as the minimum value of q decreases below about 1.4, the localized resistive interchange criterion is strongly violated, and the unstable mode with toroidal mode number n=1 changes its dominant poloidal mode number m from a (m,n)=(2,1) double tearing mode to a non-resonant (1,1) infernal mode. This (1,1) mode nonlinearly flattens both the current density and the pressure near the magnetic axis. The higher-n unstable modes are all resistive infernal modes localized around the surface where q has a minimum and the magnetic shear vanishes. They generally saturate nonlinearly at low amplitude but can cause magnetic surface breakup near the axis.

Spitzer resistivity↗
DOE OSTI2025

Electron Cyclotron Emission–Based Separatrix Identification in ITER with OMFIT Synthetic Modeling

Accurate determination of the separatrix location is essential for understanding edge plasma behavior and optimizing confinement in tokamaks, especially in next-generation devices such as ITER. In this study, a synthetic microwave diagnostics module was developed and implemented in the OMFIT framework to assess the feasibility of an electron cyclotron emission–based separatrix detection method in ITER plasmas. Simulations were carried out using ITER H-mode equilibrium scenarios with different plasma density profiles and different pedestal widths. Here, the results show that the electron emission temperature profiles consistently exhibit an inversion pattern near the edge, with a well-defined minimum point that could serve as a proxy for the separatrix location. However, unlike in the DIII-D, the minimum point in ITER is systematically offset by approximately 2 cm into the scrape-off layer, independent of density or pedestal width, which is within the radial resolution range (2 to 5 cm) determined by the 500-MHz channel spacing. While the method does not provide the exact separatrix location, it offers a reliable indicator of the boundary region and has potential applications for real-time boundary monitoring in ITER and other future fusion devices.

Electron cyclotron emission↗
DOE OSTI2025

Robustness of quasi-symmetry along parametric boundary variation

Quasi-symmetry is a guiding principle to modern stellarator optimization for improved plasma confinement. However, the robustness of optimized configurations, which can be crucial for maintaining performance under diverse engineering constraints and practical limitations, has received relatively little attention. Here, in this work, we present various case studies on this robustness, by investigating variations in 1/ v neoclassical transport when plasma configurations are smoothly altered across distinct optimized targets. These targets, optimized from different families—quasi-axisymmetric, quasi-helical, and quasi-isodynamic—are approximately matched in major radius as part of a flexible stellarator design. Our study shows that an optimized target does not always represent a local minimum in transport and that the robustness of a local minimum when present can vary significantly. Furthermore, there are configurations which belong to no established families but have transport levels as low as those of optimized targets. These results highlight the importance of conducting extended searches with key parametric variations around optimized configurations, to ensure its robustness as well as flexibility if desired.

neoclassical transport↗
DOE OSTI2024

Warm inflation with a heavy QCD axion

We propose the first model of warm inflation in which the particle production emerges directly from coupling the inflaton to Standard Model particles. Warm inflation, an early epoch of sustained accelerated expansion at finite temperature, is a compelling alternative to cold inflation, with distinct predictions for inflationary observables such as the amplitude of fluctuations, the spectral tilt, the tensor-to-scalar ratio, and non-gaussianities. In our model a heavy QCD axion acts as the warm inflaton whose coupling to Standard Model gluons sources the thermal bath during warm inflation. Axion-like couplings to non-Abelian gauge bosons have been considered before as a successful microphysical theory with emerging thermal friction that can maintain finite temperature during inflation via sphaleron heating. However, the presence of light fermions charged under the non-Abelian group suppresses particle production, hindering a realization of warm inflation by coupling to QCD. We point out that the Standard Model quarks can be heavy during warm inflation if the Higgs field resides in a high-energy second minimum which restores efficient sphaleron heating. A subsequent large reheating temperature is required to allow the Higgs field to relax to its electroweak minimum. Exploring a scenario in which hybrid warm inflation provides the large reheating temperature, we show that future collider and beam dump experiments have discovery potential for a heavy QCD axion taking the role of the warm inflaton.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Non-linear dependence of ion heat flux on plasma density at the L–H transition of JET NBI-heated deuterium–tritium plasmas

Abstract Recent JET D–T campaigns opened the possibility of unique isotope studies to investigate the L–H transition physics in view of reactor plasmas and to study the origin of the observed power threshold minimum. In the present paper, we characterise L–H transitions in the low and high-density branches of JET NBI-heated D–T plasmas. As discussed in the paper, L–H transition has been hypothesised to be determined by the transport power losses of plasma ions, i.e. the so-called ion heat flux ( Q i ). We present the first power balance analysis of JET NBI-heated D–T plasmas to evaluate the ion heat flux at the transition. Due to the experimental setting being similar to previous JET D experiments, we also directly compare the results, discussing the isotope effect and similarities between datasets. First, we find an isotope effect between D and D–T Q i , with a lower Q i in D–T plasmas. We confirm that the ion heat flux deviates from density linearity compared to the linear trend observed in wave-heated D plasmas of other tokamaks. The deviation we observe in NBI-heated L–H transitions happens at an isotope-dependent density. Plasma edge rotation correlates with Q i deviation from density linearity in the low-density branch. However, further investigations would be required to assess the role of rotation on Q i and the power threshold minimum at JET. At low plasma density, NBI power dominates Q i , while increasing the density makes the equipartition power dominant. We finally compare our results with hypotheses proposed from evidence in other tokamaks to present a complete overview of ion heat flux analyses in D and D–T NBI-heated plasmas at JET.

Vincenzi, P. (ORCID:0000000251564354)↗
DOE OSTI2025

Global stellarator coil optimization with quadratic constraints and objectives

Most present stellarator designs are produced by costly two-stage optimization: the first for an optimized equilibrium, and the second for a coil design reproducing its magnetic configuration. Few proxies for coil complexity and forces exist at the equilibrium stage. Rapid initial state finding for both stages is a topic of active research. Most present convex coil optimization codes use the least square winding surface method by Merkel (NESCOIL), with recent improvements in conditioning, regularization, sparsity, and physics objectives. While elegant, the method is limited to modeling the norms of linear functions in coil current. We present QUADCOIL, a global coil optimization method that targets combinations of linear and quadratic functions of the current. It can directly constrain and/or minimize a wide range of physics objectives unavailable in NESCOIL and REGCOIL, including the Lorentz force, magnetic energy, curvature, field-current alignment, and the maximum density of a dipole array. QUADCOIL requires no initial guess and runs nearly $10$ 2 x faster than filament optimization. Integrating it in the equilibrium optimization stage can potentially exclude equilibria with difficult-to-design coils, without significantly increasing the computation time per iteration. QUADCOIL finds the exact, global minimum in a large parameter space when possible, and otherwise finds a well-performing approximate global minimum. It supports most regularization techniques developed for NESCOIL and REGCOIL. We demonstrate QUADCOIL’s effectiveness in coil topology control, minimizing non-convex penalties, and predicting filament coil complexity with three numerical examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗