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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

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↗

DOE OSTI2025

Electromagnetic coil optimization for reduced Lorentz forces

Abstract The reduction of magnetic forces on electromagnetic coils is an important consideration in the design of high-field devices such as the stellarator or tokamak. Unfortunately, these forces may be too time-consuming to evaluate by conventional finite element modeling within an optimization loop. Although mutual forces can be computed rapidly by approximating large-bore coils as infinitely thin, this approximation does not hold for self-forces as it leads to an unphysical divergence. Recently, a novel reduced model for the self-field, self-force, and self-inductance of electromagnetic coils based on filamentary models was rigorously derived and demonstrated to be highly accurate and numerically efficient to evaluate (Hurwitz et al 2024 IEEE Trans. Magn. 60 7001614). In this paper, we present an implementation of the reduced self-force model employing automatic differentiation within the simsopt stellarator design software and use it in derivative-based coil optimization for a quasi-axisymmetric stellarator. We show that it is possible to significantly reduce point-wise forces throughout the coils, though this comes with trade-offs to fast particle losses and the minimum distance between coils and the plasma surface. The trade-off between magnetic forces and coil-surface distance is mediated by the minimum coil–coil distance for coils near the inboard side of the ‘bean’ cross-section of the plasma. The relationship between forces and fast particle losses is mediated by the normal field error. Coil forces can be lowered to a threshold with minimal deterioration to losses. Importantly, the magnet optimization approach here can be used also for tokamaks, other fusion concepts, and applications outside of fusion.

Hurwitz, Siena (ORCID:0000000166599659)↗

DOE OSTI2025

Performance of an X-γ ray detection system based on a thick silicon LGAD

We present the performance of a X-γ ray detection systems based on a 300 µm thick silicon low-gain avalanche diode (LGAD) and of an equivalent Diode structure without the gain layer, read-out by a custom-made low-noise charge amplifier. For the LGAD structure, the multiplication gains M s from 10.2 to 19.3 are measured, and the Equivalent Noise Charge (ENC) components have been studied in detail. As expected, a significant reduction of the white and 1/f voltage series and dielectric ENC components is observed with increasing gain, shortening the optimum peaking time to τ = 0.4÷8 μs depending on the gain value, with respect to the standard diode, while the parallel ENC component increases proportionally to the gain up to M s = 15. Setting the gain M s = 12.7, a minimum ENC = 34.5 electrons r.m.s. (298 eV FWHM) at τ = 1 μs is found. However, the electronic noise, evaluated on the pulser FWHM, is found to only marginally contribute to the width of the spectral lines of a radioactive source, which are dominated instead by the statistical noise of the charge multiplication within the LGAD structure. The spectral lines' excess width is found to be proportional to the signal multiplication gain M s and increases with the energy of the photons. A minimum line width of 850 eV FWHM at 13.9 keV is measured at room temperature, τ = 2 μs and a gain M s = 10.2.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Quantifying the Variability of “Fixed-Width” Buffers on Harvested Lands in Western Oregon and Washington

In contemporary forest management, buffers of unharvested trees are left along streams to protect riparian and aquatic ecosystems. Buffer regulations often focus on specific minimum width requirements, which aid in straightforward regulation and application, but minimum widths also suggest buffered edges are uniform and contain little variability. Conceptual papers suggest that alternative buffer configurations may offer greater flexibility in landscape-level protection, increase forest complexity, and enhance aquatic and riparian biodiversity. However, before considering alternatives to fixed-width buffers, it is necessary to quantify the inherent variability in current buffer practices present on the landscape. In this study, we used aerial imagery to quantify variability of buffer widths in two hundred randomly selected recently harvested units on managed land in Oregon and Washington with both fish-bearing and non-fish-bearing sections of stream. Wider buffers on larger streams had a greater magnitude of variability, but when normalized by stream size, variation was greater in smaller streams, and overall, variation ranged from 25% to 50% of the mean width. Despite local variability, buffer widths rarely fell below 9.14 m (30 ft). The variation quantified here provides an initial measure of variability to inform future management, given emerging interest in variable retention buffers.

58 GEOSCIENCES↗

DOE OSTI2025

The ability of Salmonella enterica subsp. enterica strains to form biofilms on abiotic surfaces and their susceptibility to selected essential oil components

Abstract The ability of Salmonella enterica subsp. enterica to persist and form biofilms on different surfaces can constitute a source of food contamination, being an issue of global concern. The objective of this study was to understand the biofilm formation profile of 14 S. enterica strains among different serovars and sources and to evaluate the ability of essential oil (EO) components (carveol, citronellol, and citronellal) to disinfect the biofilms formed on stainless steel and polypropylene surfaces. All the strains were able to form biofilms with counts between 5.34 to 6.78 log CFU cm−2. Then, the minimum inhibitory concentration (MIC) and minimum bactericidal concentration (MBC) of EO components were evaluated on two selected strains. All compounds inhibited the growth of Salmonella Typhimurium (strain 1; MIC = 800–1000 µg ml−1) and Salmonella Enteritidis (strain 5; MIC = 400–1000 µg ml−1) and only carveol showed bactericidal activity against strains 1 and 5 (MBC = 1200 µg ml−1). Biofilms were exposed to the EO components at 10 × MIC for 30 min and polypropylene surfaces were more difficult to disinfect showing reductions between 0.9 and <1.2 log CFU cm−2. In general, the S. enterica biofilms demonstrated a significant tolerance to disinfection, demonstrating their high degree of recalcitrance on food processing surfaces.

Pérez-Lavalle, Liliana (ORCID:0000000344825031)↗

DOE OSTI2025

Nonperturbative and perturbative dynamics of a light QCD axion: Dark matter and the strong 𝐶⁢𝑃 problem

Considerable theoretical efforts have gone into expanding the reach of the quantum chromodynamics (QCD) axion beyond its canonical mass–decay-constant relation. The 𝑍 𝒩 QCD axion model reduces the QCD axion mass naturally, by invoking a discrete 𝑍 𝒩 symmetry through which the axion field is coupled to 𝒩 copies of the Standard Model. Before the QCD phase transition at temperature 𝑇 QCD , the 𝑍 𝒩 potential has a minimum at misalignment angle 𝜃 = 𝜋. At 𝑇 QCD , 𝜃 = 𝜋 becomes a maximum; the axion potential becomes exponentially suppressed and develops 𝒩 minima—only one of which actually solves the strong 𝐶⁢𝑃 problem. Before 𝑇 QCD , 𝜃 relaxes toward 𝜋. After 𝑇 QCD , the axion field starts from around the hilltop and may have sufficient kinetic energy to overcome the newly suppressed potential barriers. Such a field evolution leads to nonperturbative effects via the self-interactions near the hilltop, which can cause the exponential growth of fluctuations and backreaction on the coherent motion. This behavior can influence the relic density of the field and the minimum in which it settles. We conduct the first lattice simulations of the 𝑍 𝒩 QCD axion using 𝒞osmoℒattice to accurately calculate dark matter abundances and find nonperturbative dynamics reduce the abundance by up to a factor of two. We furthermore find that the probability of solving the strong 𝐶⁢𝑃 problem tends to diverge considerably from the naïve expectation of 1/𝒩.

Axions↗

DOE OSTI2026

Observation of Suppressed Charged-Particle Production in Ultrarelativistic Oxygen-Oxygen Collisions

A hot and dense state of nuclear matter, known as the quark-gluon plasma, is created in collisions of ultrarelativistic heavy nuclei. Highly energetic quarks and gluons, collectively referred to as partons, lose energy as they travel through this matter, leading to suppressed production of particles with large transverse momenta (𝑝 T ). Conversely, high-𝑝 T particle suppression has not been seen in proton-lead collisions, raising questions regarding the minimum system size required to observe parton energy loss. Oxygen-oxygen (OO) collisions examine a region of effective system size that lies between these two extreme cases. The CMS detector at the CERN LHC has been used to quantify charged-particle production in inclusive OO collisions for the first time via measurements of the nuclear modification factor (𝑅 AA ). The 𝑅 AA is derived by comparing particle production to expectations based on proton-proton (𝑝⁢𝑝) data and has a value of unity in the absence of nuclear effects. The data for OO and 𝑝⁢𝑝 collisions at a nucleon-nucleon center-of-mass energy $\sqrt{s_{NN}}$ = 5.36 TeV correspond to integrated luminosities of 6.1 nb −1 and 1.02 pb −1 , respectively. The 𝑅 AA is below unity with a minimum of 0.69 ± 0.04 around 𝑝 T = 6 GeV. The data exhibit better agreement with theoretical models incorporating parton energy loss as compared to baseline models without energy loss.

FOS: Physical sciences↗

DOE OSTI2025

Precision Measurement of Net-Proton-Number Fluctuations in Au + Au Collisions at RHIC

We report precision measurements on cumulants (𝐶 𝑛 ) and factorial cumulants (𝜅 𝑛 ) of (net) proton number distributions up to fourth order in Au + Au collisions over center-of-mass energies $\sqrt{s_{NN}}$ = 7.7–27 GeV from phase II of the Beam Energy Scan program at RHIC. (Anti)protons are selected at midrapidity (|𝑦| < 0.5) within a transverse momentum range of 0.4 < 𝑝 𝑇 < 2.0 GeV/𝑐. Relative to various noncritical-point model calculations and peripheral collision 70%–80% data, the net proton 𝐶 4 /𝐶 2 measurement in 0%–5% collisions shows a minimum around 19.6 GeV for significance of deviation at ∼2–5⁢𝜎. A minimum in 𝐶 4 /𝐶 2 with respect to a noncritical baseline is expected to be a characteristic feature of the signature associated with a critical point in the QCD phase diagram. In addition, deviations from noncritical baselines around the same collision energy region are also seen in proton factorial cumulant ratios, especially in 𝜅 2 /𝜅 1 and 𝜅 3 /𝜅 1 . As a result, dynamical model calculations including a critical point are called for in order to understand these precision measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

New H erwig7 underlying event tune: From RHIC to LHC energies

We present parameter sets corresponding to new underlying event tunes for the H erwig7.3 Monte Carlo event generator. The existing H erwig tunes are in good agreement with LHC data, however, they are not typically designed for center-of-mass energies below $\sqrt{𝑠}$ = 300 GeV. The tunes presented in this study can describe midrapidity data collected at the nominal RHIC energy of $\sqrt{𝑠}$ = 200 GeV as well as higher center-of-mass energies utilized by experiments at the LHC and Tevatron. The base “New Haven” tune is developed by fitting minimum-bias simulations of proton-proton ($𝑝⁢𝑝$) collisions to midrapidity identified hadron and jet data from the STAR experiment. The “Nashville” tune includes a separate set of parameters developed by tuning to Tevatron proton-antiproton ($𝑝\bar{⁢𝑝}$) data at $\sqrt{𝑠}$ = 300, 900 and 1960 GeV from CDF, and LHC $𝑝⁢𝑝$ measurements from CMS at $\sqrt{𝑠}$ =7 TeV, in addition to the STAR measurements. Both new tunes demonstrate significant improvements over the recommended default tune currently included in the latest version of H erwig for minimum bias production. As such, we advocate using these tunes for future simulation studies at midrapidity by experimental collaborations at RHIC (STAR and sPHENIX) and the LHC (ALICE, ATLAS, and CMS).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Low-temperature magnetic structure and lattice response in SmCuAs 2

We investigated the structural and magnetic properties of single-crystalline SmCuAs 2 using high-resolution synchrotron x-ray diffraction and x-ray resonant magnetic scattering (XRMS) at the Sm 𝐿 2 and 𝐿 3 edges. Temperature-dependent diffraction measurements confirm that SmCuAs2 maintains its tetragonal symmetry from room temperature down to 8 K, with lattice parameters showing anomalous behavior below the resistivity minimum (𝑇≈30K). Notably, the c-axis lattice parameter exhibits a plateau and subsequent increase near the Néel temperature, indicating magnetoelastic coupling. XRMS measurements reveal a commensurate antiferromagnetic structure with a propagation vector 𝒒=(0,0,0.5). Our measurement shows that the Sm moments are aligned within the 𝒂𝒃 plane and arranged in a ++−− stacking along the 𝒄 axis. Comparison with related 𝑅⁢𝐸⁢CuAs 2 compounds (RE = Pr, Nd, and Gd) suggests that in-plane moment orientation and associated magnetic frustration play a key role in the emergence of the resistivity minimum. Differences in spin-orbit and magnetoelastic coupling across the series highlight their importance in governing low-temperature transport behavior.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Moduli axions, stabilizing moduli, and the large field swampland conjecture in heterotic M-theory

We compute the F- and D-term potential energy for the dilaton, complex structure, and Kähler moduli of realistic vacua of heterotic M-theory compactified on Calabi-Yau threefolds where, for simplicity, we choose ℎ 1,1 = ℎ 2,1 = 1. However, the formalism is immediately applicable to the “universal” moduli of Calabi-Yau threefolds with ℎ 1,1 = ℎ 1,2 > 1 as well. The F-term potential is computed using the nonperturbative complex structure, gaugino condensate and “world sheet instanton” superpotentials in theories in which the hidden sector contains an anomalous U⁡(1) structure group. The Green-Schwarz anomaly cancellation induces inhomogeneous “axion” transformations for the imaginary components of the dilaton and Kähler modulus—which then produce a D-term potential. V D is a function of the real components of the dilaton and Kähler modulus (s and t) that is minimized and precisely vanishes along a unique line in the s–t plane. Excitations transverse to this line have a mass m anom which is an explicit function of t. The F-term potential energy is then evaluated along the V D = 0 line. For values of t small enough that m anom ≳ M U —where M U is the compactification scale—we plot V F for a realistic choice of coefficients as a function of the Pfaffian parameter p. We find values of p for which V F has a global minimum at negative or zero vacuum density or a metastable minimum with positive vacuum density. In all three cases, the s, t, and associated “axion” moduli are completely stabilized. Finally, we show that, for any of these vacua, the large t behavior of the potential energy satisfies the “large scalar field” Swampland conjecture.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

First measurement of 𝐷 𝑠⁢1 ⁡(1 + )⁢(2536) + and 𝐷$^{*}_{𝑠⁢2}$⁡(2 + )⁢(2573) + production in proton-proton collisions at $\sqrt{s}$ = 13 TeV at the LHC

The production yields of the orbitally excited charm-strange mesons D s⁢1 ⁡(1 + )⁢(2536)+ and 𝐷$^{*}_{𝑠⁢2}$⁡(2 + )⁢(2573) + were measured for the first time in proton-proton (pp) collisions at a center-of-mass energy of $\sqrt{s}$ =13 TeV with the ALICE experiment at the LHC. The D$^{+}_{s⁢1}$ and D$^{*+}_{s⁢2}$ mesons were measured at midrapidity (|𝑦| < 0.5) in minimum-bias and high-multiplicity pp collisions in the transverse-momentum interval 2< 𝑝 T < 24 GeV/𝑐. Their production yields relative to the D$^{+}_{s}$ ground-state yield were found to be compatible between minimum-bias and high-multiplicity collisions, as well as with previous measurements in e ± ⁢p and e + ⁢e − collisions. The measured D$^{+}_{s⁢1}$/D$^{+}_{s}$ and D$^{*+}_{s⁢2}$/D$^{+}_{s}$ yield ratios are described by statistical hadronization models and can be used to tune the parameters governing the production of excited charm-strange hadrons in Monte Carlo generators, such as pythia 8.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Quantum Adiabatic Optimization with Rydberg Arrays: Localization Phenomena and Encoding Strategies

Quantum adiabatic optimization seeks to solve combinatorial problems using quantum dynamics, requiring the Hamiltonian of the system to align with the problem of interest. However, these Hamiltonians are often incompatible with the native constraints of quantum hardware, necessitating encoding strategies to map the original problem into a hardware-conformant form. While the classical overhead associated with such mappings is easily quantifiable and typically polynomial in problem size, it is much harder to quantify their overhead on the quantum algorithm, e.g., the transformation of the adiabatic timescale. In this work, we address this challenge on the concrete example of the encoding scheme proposed in [Nguyen , PRX Quantum , 010316 (2023)], which is designed to map optimization problems on arbitrarily connected graphs into Rydberg atom arrays. We consider the fundamental building blocks underlying this encoding scheme and determine the scaling of the minimum gap with system size along adiabatic protocols. Even when the original problem is trivially solvable, we find that the encoded problem can exhibit an exponentially closing minimum gap. We show that this originates from a quantum coherent effect, which gives rise to an unfavorable localization of the ground-state wave function. On the QuEra Aquila neutral atom machine, we observe such localization and its effect on the success probability of finding the correct solution to the encoded optimization problem. Finally, we propose quantum-aware modifications of the encoding scheme that avoid this quantum bottleneck and lead to an exponential improvement in the adiabatic performance. This highlights the crucial importance of accounting for quantum effects when designing strategies to encode classical problems onto quantum platforms. Published by the American Physical Society 2025

Bombieri, Lisa (ORCID:0009000950422897)↗

DOE OSTI2026

Shuttling Majorana zero modes in disordered and noisy topological superconductors

The braiding of Majorana zero modes (MZMs) forms the fundamental building block for topological quantum computation. Braiding protocols which involve the physical exchange of MZMs are typically envisioned on a network of topological superconducting wires. An important component of these protocols is the transport of MZMs, which can be performed by using electric gates to locally tune sections of the wire between topologically trivial and nontrivial phases. In this work, we numerically simulate this transport by tuning a single section of a superconducting wire which contains either disorder (uncorrelated and correlated) or noise. We focus on the impact of these additional effects on the diabatic error, which describes unwanted transitions between the ground state and excited states. We show that the behavior of the average diabatic error is predominantly controlled by the statistics of the minimum bulk energy gap which is suppressed in the presence of disorder. The increase in diabatic error can be several orders of magnitude and is most deleterious when the disorder correlation length is a finite fraction of the transport distance and negligible when these lengths are far apart. In the presence of noise, the diabatic error is significantly enhanced due to optical transitions which depend on the minimum bulk energy gap as well as the frequency modes present in the noise. The results presented here serve to further characterize the diabatic error in disordered and noisy settings, which are important considerations in practical implementations of physical braiding schemes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Minimizing the Electromechanical Stresses in Poloidal Field Coils by Optimizing their Numbers and Locations using FREDA Framework

Poloidal field (PF) and central solenoid (CS) coils play a crucial role in sustaining the equilibrium and preserving the shape of highly confined tokamak plasmas. Ensuring that PF coil current and mechanical stress stay within superconducting and structural limitations is an important check in the design assessment. Minimizing the PF coil currents and mechanical stresses influences reliability, cost, and performance. A free-boundary MHD equilibrium code—FreeGS is employed within the fusion reactor design and assessment (FREDA) whole facility modeling (WFM) framework to construct the plasma equilibrium based on the configuration and currents in the PF coils. Here, we present the capability of the FreeGS code to minimize the currents, forces, and electromagnetic stresses on the PF coils by optimizing their number, sizes, structures, and locations while maintaining an MHD stable plasma configuration with a large confinement factor. The workflow is initialized with a configuration of plasma parameters and coils’ locations from the 0-D tokamak build systems code in the FREDA framework. Then, FreeGS is called to calculate the initial equilibrium at the minimum total current in PF coils. Thereafter, FreeGS’s internal optimizer minimizes the currents and hoop and central forces on the PF coils while maintaining the reference equilibrium. Finally, the input configuration is updated with the optimized parameters for equilibria over the ramp-up phase of a burning-plasma operation. FREDA’s whole facility optimization capability, which includes all magnetic field coil systems, blanket, vacuum vessel (VV), first wall, divertor, etc., is under development and out of the scope for this study.

Hassan, Ehab [ORNL] (ORCID:0000000181060301)↗