Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flux optimization”

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 OSTI2022

Quantifying Northern High Latitude Gross Primary Productivity (GPP) Using Carbonyl Sulfide (OCS)

The northern high latitude (NHL, 40°N to 90°N) is where the second peak region of gross primary productivity (GPP) other than the tropics. The summer NHL GPP is about 80% of the tropical peak, but both regions are still highly uncertain. Carbonyl sulfide (OCS) provides an important proxy for photosynthetic carbon uptake. Here we optimize the OCS plant uptake fluxes across the NHL by fitting atmospheric concentration simulation with the GEOS-CHEM global transport model to the aircraft profiles acquired over Alaska during NASA's Carbon in Arctic Reservoirs Vulnerability Experiment (2012–2015). We use the empirical biome-specific linear relationship between OCS plant uptake flux and GPP to derive the six plant uptake OCS fluxes from different GPP data. Such GPP-based fluxes are used to drive the concentration simulations. We evaluate the simulations against the independent observations at two ground sites of Alaska. The optimized OCS fluxes suggest the NHL plant uptake OCS flux of -247 Gg S year -1 , about 25% stronger than the ensemble mean of the six GPP-based OCS fluxes. GPP-based OCS fluxes systematically underestimate the peak growing season across the NHL, while a subset of models predict early start of season in Alaska, consistent with previous studies of net ecosystem exchange. The OCS optimized GPP of 34 PgC yr -1 for NHL is also about 25% more than the ensembles mean from six GPP data. Further work is needed to fully understand the environmental and biotic drivers and quantify their rate of photosynthetic carbon uptake in Arctic ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Functional analysis of the methylerythritol phosphate pathway terminal enzymes IspG and IspH from Zymomonas mobilis

ABSTRACT Isoprenoids are a diverse family of compounds that are synthesized from two isomeric compounds, isopentenyl diphosphate and dimethylallyl diphosphate. In most bacteria, isoprenoids are produced from the essential methylerythritol phosphate (MEP) pathway. The terminal enzymes of the MEP pathway IspG and IspH are [4Fe-4S] cluster proteins, and in Zymomonas mobilis, the substrates of IspG and IspH accumulate in cells in response to O 2 , suggesting possible lability of their [4Fe-4S] clusters. Here, we show using complementation assays in Escherichia coli that even under anaerobic conditions, Z. mobilis IspG and IspH are not as functional as their E. coli counterparts, requiring higher levels of expression to rescue viability. A deficit of the sulfur utilization factor (SUF) Fe-S cluster biogenesis pathway did not explain the reduced function of Z. mobilis IspG and IspH since no improvement in viability was observed in E. coli expressing the Z. mobilis SUF pathway or having increased expression of the E. coli SUF pathway. Complementation of single and double mutants with various combinations of Z. mobilis and E. coli IspG and IspH indicated that optimal growth required the pairing of IspG and IspH from the same species. Furthermore, Z. mobilis IspH conferred an O 2 -sensitive growth defect to E. coli that could be partially rescued by co-expression of Z. mobilis IspG. In vitro analysis showed O 2 sensitivity of the [4Fe-4S] cluster of both Z. mobilis IspG and IspH. Altogether, our data indicate an important role of the cognate protein IspG in Z. mobilis IspH function under both aerobic and anaerobic conditions. IMPORTANCE Isoprenoids are one of the largest classes of natural products, exhibiting diversity in structure and function. They also include compounds that are essential for cellular life across the biological world. In bacteria, isoprenoids are derived from two precursors, isopentenyl diphosphate and dimethylallyl diphosphate, synthesized primarily by the methylerythritol phosphate pathway. The aerotolerant Z. mobilis has the potential for methylerythritol phosphate pathway engineering by diverting some of the glucose that is typically efficiently converted into ethanol to produce isoprenoid precursors to make bioproducts and biofuels. Our data revealed the surprising finding that Z. mobilis IspG and IspH need to be co-optimized to improve flux via the methyl erythritol phosphate pathway in part to evade the oxygen sensitivity of IspH.

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

Magneto-Structural Combined Dimensional and Topology Optimization of Interior Permanent Magnet Synchronous Machine Rotors

A magneto-structural combined dimensional and topology optimization approach for interior permanent magnet synchronous machine (IPMSM) rotor design is proposed using the Solid Isotropic Material Penalization (SIMP) density-based topology optimization method. This method optimizes the location and dimensions of the permanent magnets (PMs) embedded in the rotor, while simultaneously optimizing the lamination flux barriers using topology optimization.To allow for the location and dimensions of the PMs to change while retaining the same number of mesh elements and connectivity, the mesh must be deformed. Mesh deformation is accomplished using a Laplace smoothing technique. To increase the range of PM displacement and size changes, the non-PM areas of the rotor design domain are divided into sub-regions along whose boundaries the mesh nodes are redistributed also using a Laplace smoothing algorithm. Three examples of IPMSM rotor optimizations are presented to demonstrate the proposed approach.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Progress of Divertor Heat and Particle Flux Control in EAST for Advanced Steady-State Operation in the Last 10 Years

Active control of the excessively high heat and particle fluxes on the divertor target plates is of fundamental importance to the steady state operation of tokamaks, especially for fusion reactors. A series of experiments have been carried out on this critical issue to relieve the plasma-wall interactions in the experimental advanced superconducting tokamak (EAST) in the last ten years, not only contributing to the long pulse operation of EAST itself, but also providing physical understandings and potential techniques to the next-generation devices like ITER. Here we have characterized the power deposition pattern and broadened the divertor footprint width effectively. The plasma-wetted area is actively handled using either 3-dimentional edge magnetic topology or advanced plasma equilibrium, thereby peak heat flux around the strike point is reduced. Active control of detachment or radiation compatible with core plasma performance has progressed significantly in very recent years, with a series of active feedback control modules developed and utilized successfully, based on the divertor physics advances with both experiments and simulation. The upper divertor of EAST was upgraded from graphite to active water-cooling ITER-like tungsten in 2014, exhibiting much enhanced heat removal capability. As for the particle exhaust including both fueling and impurity particles, in addition to wall conditioning and impurity source control, the efficiency of particle flux exhaust is optimized by making full use of the divertor closure and the plasma drifts in both scrape-off layer and divertor volume. These heat and particle exhaust advances contribute greatly to a series of EAST achievements like H-mode operation over 100 s. A brief near-term plan on the integrated control of divertor plasma-wall interactions in long-time scale will also be introduced, aiming to provide favorable divertor operation solution for ITER and CFETR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Gyrokinetic profile prediction and validation of a negative triangularity plasma in ASDEX Upgrade

In this work, gyrokinetic simulations are performed with the CGYRO code (Candy et al 2016 J. Comput. Phys. 324 73–93) for a negative triangularity H-mode plasma in ASDEX Upgrade, and compared with experimental measurements. The PORTALS framework (Rodriguez-Fernandez et al 2024 Nucl. Fusion 64 076034) is used to accelerate the prediction of kinetic profiles for this plasma, using surrogate modeling and Bayesian optimization. Ion heat flux, electron heat flux, and electron particle flux are simultaneously matched across the simulated radial regime of the plasma (normalized radius $r/a = 0.35-0.90$), and the resulting ion temperature, electron temperature, and electron density profiles match well with the experimental profile data within this radial range. A synthetic Correlation Electron Cyclotron Emission diagnostic is applied to find well-matched electron temperature fluctuation properties between simulation and experiment. The flux-matched profiles provide a basis for investigation of the turbulence nature across the plasma radius, revealing the dominance of Trapped Electron Mode turbulence at $r/a = 0.35$, the dominance of Ion Temperature Gradient turbulence at $r/a = 0.55$, 0.75, and 0.83, and an instability boundary at $r/a = 0.90$.

gyrokinetic simulation↗
DOE OSTI2021

Spatial distributions of X CO 2 seasonal cycle amplitude and phase over northern high-latitude regions

Satellite-based observations of atmospheric carbon dioxide (CO 2 ) provide measurements in remote regions, such as the biologically sensitive but undersampled northern high latitudes, and are progressing toward true global data coverage. Recent improvements in satellite retrievals of total column-averaged dry air mole fractions of CO 2 (X CO 2 ) from the NASA Orbiting Carbon Observatory 2 (OCO-2) have allowed for unprecedented data coverage of northern high-latitude regions, while maintaining acceptable accuracy and consistency relative to ground-based observations, and finally providing sufficient data in spring and autumn for analysis of satellite-observed X CO 2 seasonal cycles across a majority of terrestrial northern high-latitude regions. Here, we present an analysis of X CO 2 seasonal cycles calculated from OCO-2 data for temperate, boreal, and tundra regions, subdivided into 5° latitude by 20° longitude zones. We quantify the seasonal cycle amplitudes (SCAs) and the annual half drawdown day (HDD). OCO-2 SCAs are in good agreement with ground-based observations at five high-latitude sites, and OCO-2 SCAs show very close agreement with SCAs calculated for model estimates of X CO 2 from the Copernicus Atmosphere Monitoring Services (CAMS) global inversion-optimized greenhouse gas flux model v19r1 and the CarbonTracker2019 model (CT2019B). Model estimates of X CO 2 from the GEOS-Chem CO 2 simulation version 12.7.2 with underlying biospheric fluxes from CarbonTracker2019 (GC-CT2019) yield SCAs of larger magnitude and spread over a larger range than those from CAMS, CT2019B, or OCO-2; however, GC-CT2019 SCAs still exhibit a very similar spatial distribution across northern high-latitude regions to that from CAMS, CT2019B, and OCO-2. Zones in the Asian boreal forest were found to have exceptionally large SCA and early HDD, and both OCO-2 data and model estimates yield a distinct longitudinal gradient of increasing SCA from west to east across the Eurasian continent. In northern high-latitude regions, spanning latitudes from 47 to 72° N, longitudinal gradients in both SCA and HDD are at least as pronounced as latitudinal gradients, suggesting a role for global atmospheric transport patterns in defining spatial distributions of X CO2 seasonality across these regions. GEOS-Chem surface contact tracers show that the largest X CO 2 SCAs occur in areas with the greatest contact with land surfaces, integrated over 15–30d. The correlation of X CO 2 SCA with these land surface contact tracers is stronger than the correlation of X CO 2 SCA with the SCA of CO 2 fluxes or the total annual CO 2 flux within each 5° latitude by 20° longitude zone. This indicates that accumulation of terrestrial CO 2 flux during atmospheric transport is a major driver of regional variations in X CO 2 SCA.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Fast physics-based launcher optimization for electron cyclotron current drive

With the increased urgency to design fusion pilot plants, fast optimization of electron cyclotron current drive (ECCD) launchers is paramount. Traditionally, this is done by coarsely sampling the 4D parameter space of possible launch conditions consisting of (1) the launch location (constrained to lie along the reactor vessel), (2) the launch frequency, (3) the toroidal launch angle, and (4) the poloidal launch angle. For each initial condition, a ray-tracing simulation is performed to evaluate the ECCD efficiency. Unfortunately, this approach often requires a large number of simulations (sometimes millions in extreme cases) to build up a dataset that adequately covers the plasma volume, which must then be repeated every time the design point changes. Here we adopt a different approach. Rather than launching rays from the plasma periphery and hoping for the best, we instead directly reconstruct the optimal ray for driving current at a given flux surface using a reduced physics model coupled with a commercial ray-tracing code. Repeating this throughout the plasma volume requires only hundreds of simulations, constituting a significant speedup. The new method is validated on two separate example tokamak profiles, and is shown to reliably drive localized current at the specified flux surface with the same optimal efficiency as obtained from the traditional approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Optimization of radiochromic film stacks to diagnose high-flux laser-accelerated proton beams

Here, we extend flatbed scanner calibrations of GafChromic EBT3, MD-V3, and HD-V2 radiochromic films using high-precision x-ray irradiation and monoenergetic proton bombardment. By computing a visibility parameter based on fractional errors, optimal dose ranges and transitions between film types are identified. The visibility analysis is used to design an ideal radiochromic film stack for the proton energy spectrum expected from the interaction of a petawatt laser with a cryogenic hydrogen jet target.

Curry, C. B. (ORCID:000000018756181X)↗
DOE OSTI2024

Optimization of nonlinear turbulence in stellarators

We present new stellarator equilibria that have been optimized for reduced turbulent transport using nonlinear gyrokinetic simulations within the optimization loop. The optimization routine involves coupling the pseudo-spectral GPU-native gyrokinetic code GX with the stellarator equilibrium and optimization code DESC . Since using GX allows for fast nonlinear simulations, we directly optimize for reduced nonlinear heat fluxes. To handle the noisy heat flux traces returned by these simulations, we employ the simultaneous perturbation stochastic approximation (SPSA) method that only uses two objective function evaluations for a simple estimate of the gradient. We show several examples that optimize for both reduced heat fluxes and good quasi-symmetry as a proxy for low neoclassical transport. Finally, we run full transport simulations using the T3D stellarator transport code to evaluate the changes in the macroscopic profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Efficient GPU Implementation of Automatic Differentiation for Computational Fluid Dynamics

Many scientific and engineering applications require repeated calculations of derivatives of output functions with respect to input parameters. Automatic Differentiation (AD) is a method that automates derivative calculations and can significantly speed up code development. In Computational Fluid Dynamics (CFD), derivatives of flux functions with respect to state variables (Jacobian) are needed for efficient solutions of the nonlinear governing equations. AD of flux functions on graphics processing units (GPUs) is challenging as flux computations involve many intermediate variables that create high register pressure and require significant memory traffic because of the need to store the derivatives. This paper presents a forward-mode AD method based on multivariate dual numbers that addresses these challenges and simultaneously reduces the floating-point operation count. The dimension of the multivariate dual numbers is optimized for performance. The flux computations are restructured to minimize the number of temporary variables and reduce register pressure. For effective utilization of memory bandwidth, shared memory is used to store the local flux Jacobian. This AD implementation is compared with several other Jacobian implementations on an NVIDIA V100 GPU (V100). For three-dimensional perfect-gas compressible-flow equations implemented in a practical CFD code, the AD implementation of a flux Jacobian based on multivariate dual numbers of dimension 5 outperforms all other GPU AD implementations on V100. Its performance is comparable with the optimized hand-differentiated version. Finally, the implementation achieves 75% of the peak floating-point throughput and 61 % of the peak global device memory bandwidth usage.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Development of an immersion fiber optic Raman probe for real-time analysis of molten materials

This study presents an advancement in high-temperature Raman spectroscopy, specifically for analyzing molten materials. It introduces an approach by integrating a fiber-optic Raman probe with a copper block protection system designed to endure extreme thermal conditions. The copper block features an open port designed to accommodate an external telescope with a 3 cm focal length, enabling Raman spectra collection in challenging high-temperature environments. A built-in gas channel ensures a continuous flow of argon gas to prevent flux intrusion. The robust copper block acts as a reliable shield, safeguarding the fiber-optic Raman probe within molten materials. This enhancement maintains the probe's integrity and significantly improves its resilience, making it ideal for rigorous investigations of molten substances. This advancement is particularly relevant in metallurgy, where flux materials impact production quality and efficiency. The ability to acquire Raman signals under elevated thermal conditions offers opportunities for studying molecular dynamics, compositional changes, and chemical interactions within molten substances. This introduced direct immersion probing technique has implications, benefiting both scientific and industrial fields. It holds promise for advancing research and exploration in various contexts, from fundamental scientific inquiries to practical applications in metallurgical processes, where flux materials are critical for optimizing production quality and efficiency. Furthermore, this approach enhances the capabilities of high-temperature Raman spectroscopy, making it a valuable tool for investigating molten materials and their properties in diverse settings.

Argon↗
DOE OSTI2025

Operando Neutron Imaging of Lithium Flux and Gradient Cathode Design for Enhanced Kinetics in High‐Loading All‐Solid‐State Li─S Batteries

All-solid-state Li–sulfur batteries (ASSLSBs) are considered promising candidates for next-generation energy storage owing to their inherent safety, high energy density, and abundant sulfur resources. However, slow redox kinetics greatly limit sulfur utilization during solid-solid sulfur reactions, leading to significant challenges to achieve efficient performance in high-mass-loading ASSLSBs. Here, operando neutron image is employed to directly visualize, for the first time, that sluggish Li + transport kinetics and the uneven distribution of Li + during cathodic reactions are critical factors limiting sulfur conversion. To address this issue, gradient cathode architectures comprising three and five layers are designed, in which catholyte concentrations are strategically varied to optimize Li-ion flux and enhance ionic conductivity of the whole composite cathode electrode. Operando neutron imaging distinctly visualizes and confirms that three-layer gradient approach significantly enhances Li-ion mobility, resulting in more uniform redox reactions and greatly improved sulfur utilization compared to traditional non-gradient structures. Consequently, the three-layer gradient cathode achieves superior rate performance and reduced electrode polarization at high sulfur mass loadings of 4.5 and 6.0 mg cm −2 . Furthermore, the applicability and scalability of this design are demonstrated in a five-layer gradient cathode architecture, achieving an impressive discharge specific capacity increase from 656 mAh g −1 (three-layer gradient) to 1232 mAh g −1 at 1/20 C for ultra-high sulfur loading of 7.5 mg cm −2 . In conclusion, this innovative gradient cathode design offers substantial advancements in understanding and overcoming Li-ion transport limitations, paving the way toward practical, high-energy-density ASSLSBs.

25 ENERGY STORAGE↗
DOE OSTI2024

Molecular beam epitaxy growth and characterization of GePb alloys

Pb based group-IV alloys such as GePb have been gaining interest as a potential alternative for infrared detectors, quantum materials, and high-speed electronic devices. Challenges remain in their growth due to the extremely low solid solubility of Pb in the Ge–Pb system. This paper reports molecular beam epitaxy growth of GePb alloy thin films on Ge(100) substrates. Effusion cells of Ge and Pb are used to control the flux ratio independently. The optimal substrate temperature is found to be near the thermocouple temperature of 300 °C based on the characterization of the grown films using high-resolution x-ray diffraction. A large change in the Ge:Pb beam equivalent pressure ratio from 10:1 to 1:1 results in only a minimal increase of the Pb composition from 0.74% to 2.84% as estimated from Raman spectroscopy and Rutherford backscattering spectrometry. Furthermore, scanning electron microscopy images show a large volume of Pb islands on the surface that form into either long trapezoidal rods or uniform droplets, with increasing Pb flux and growth time the density of Pb islands increased.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Neutrino oscillation prospects with a dual-baseline beam from BNL to SNOLAB and SURF

The Electron-Ion Collider (EIC) is a next-generation accelerator primarily designed to study the internal structure of nucleons through high-precision electron-hadron collisions. In this work, we explore the feasibility of employing a 1 MW fraction of the EIC proton beam to generate a high-intensity GeV-scale neutrino beam for long-baseline oscillation studies. We have simulated proton-target interactions and optimize the resulting neutrino fluxes for water-based liquid scintillator (WbLS) detectors located at distinct baselines of 900 km and at 2900 km. Oscillation analyses performed with GLoBES show that extended baselines allow access to multiple oscillation maxima, significantly enhancing sensitivity to leptonic CP violation. The study also examines the interplay between matter effects and the intrinsic CP violating phase in shaping observable asymmetries. We note that simplified systematics and no backgrounds are used in this analysis to establish the baseline physics potential. These results suggest that the EIC proton beam could provide a novel and complementary source for precision neutrino physics, extending the scientific reach of the EIC program.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Unveiling the complex configurational landscape of the intralayer cavities in a crystalline carbon nitride

The in-depth understanding of the reported photoelectrochemical properties of the layered carbon nitride, poly(triazine imide)/LiCl (PTI/LiCl), has been limited by the apparent disorder of the Li/H atoms within its framework. To understand and resolve the current structural ambiguities, an optimized one-step flux synthesis (470 °C, 36 h, LiCl/KCl flux) was used to prepare PTI/LiCl and deuterated-PTI/LiCl in high purity. Its structure was characterized by a combination of neutron/X-ray diffraction and transmission electron microscopy. The range of possible Li/H atomic configurations was enumerated for the first time and, combined with total energy calculations, reveals a more complex energetic landscape than previously considered. Experimental data were fitted against all possible structural models, exhibiting the most consistency with a new orthorhombic model (Sp. Grp. Ama2) that also has the lowest total energy. In addition, a new Cu(I)-containing PTI (PTI/CuCl) was prepared with the more strongly scattering Cu(I) cations in place of Li, and most closely matching with the partially-disorder structure in Cmc2 1 . Thus, a complex configurational landscape of PTI is revealed to consist of a number of ordered crystalline structures that are new potential synthetic targets, such as with the use of metal-exchange reactions.

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

Indium as a surfactant: Effects on growth morphology and background impurity in GaN films grown by ammonia-assisted molecular beam epitaxy

We report on the improvement of the surface morphology of c-plane GaN films grown at high growth rates (~1 µm/h) using ammonia molecular beam epitaxy through a series of growth optimizations as well as the introduction of indium as a surfactant. The indium surfactant was expected to help with the adatom mobility and, thus, provide smoother growth surfaces. Through a combination of varying V/III ratios, In flux, and growth temperatures, an optimal condition for surface morphology, characterized by atomic force microscopy, was achieved. At higher Ga fluxes for fast growth rates (~1 µm/h and beam equivalent pressures of ~5 × 10 –7 Torr), higher ammonia flows were necessary to preserve the surface morphology. In addition, indium was an effective surfactant—reducing the roughness and improving the overall surface morphology. However, excessive indium causes the surface morphology to degrade, potentially due to the enhancement of the Ga desorption from the surface as a result of the reaction of indium with ammonia for high indium fluxes. The indium surfactant also resulted in a reduction of background Si impurity concentrations in the film. These effects allow for the growth of thick drift layers with low background dopant concentrations for vertical GaN power devices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Optimal Heliostat Assignment Strategy for Multiple-Receiver Systems

Next-generation concentrating solar power concepts currently being developed under the U.S. Department of Energy's Gen3 program include several novel heliostat field and receiver arrangements. Among these concepts are solar fields utilizing multiple receivers or several distinct receiver flux control zones, and each of which has particular power or peak flux constraints. In this paper, we present a linear program that identifies optimal heliostat-to-receiver assignment strategies for multi-receiver or multi-zone solar fields based on detailed calculations made by an optics modeling package. Two distinct problems are addressed: namely, (1) a design problem that chooses the optimal set of heliostats to include in a final layout given receiver constraints and performance, and (2) an operations problem that assigns existing heliostats to a receiver, maximizing field power output while maintaining receiver power constraints. The optimization methodology is integrated into the concentrating solar power tower modeling package SolarPILOTTM, and results are presented showing that heliostat assignment strategies can impact both initial heliostat field layout characteristics and the cost of thermal energy produced by the field. The methodology succeeds in maintaining power requirements but at the expense of producing different flux profiles on each receiver surface. Consequently, the paper also discusses additional measures that are applied to ensure that desirable flux distribution properties are maintained, reducing local flux mean-absolute-deviation values from a baseline of 17% or greater to less than 5% compared to a desired reference flux profile.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2023

Direct stellarator coil optimization for nested magnetic surfaces with precise quasi-symmetry

We present a robust optimization algorithm for the design of electromagnetic coils that generate vacuum magnetic fields with nested flux surfaces and precise quasi-symmetry. The method is based on a bilevel optimization problem, where the outer coil optimization is constrained by a set of inner least squares optimization problems whose solutions describe magnetic surfaces. The outer optimization objective targets coils that generate a field with nested magnetic surfaces and good quasi-symmetry. The inner optimization problems identify magnetic surfaces when they exist, and approximate surfaces in the presence of magnetic islands or chaos. We show that this formulation can be used to heal islands and chaos, thus producing coils that result in magnetic fields with precise quasi-symmetry. We show that the method can be initialized with coils from the traditional two-stage coil design process, as well as coils from a near-axis expansion optimization. We present a numerical example where island chains are healed and quasi-symmetry is optimized up to surfaces with aspect ratio 6. Another numerical example illustrates that the aspect ratio of nested flux surfaces with optimized quasi-symmetry can be decreased from 6 to approximately 4. In conclusion, the last example shows that our approach is robust and a cold-start using coils from a near-axis expansion optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗