Fusion Nuclear Science Facility Research Articles

Analysis and design of fast flow liquid Li divertor for fusion nuclear science facility (FNSF) using coupled plasma boundary and LM MHD/heat transfer codes *

The SOLPS-ITER code is utilized to analyze the boundary plasma associated with a fast-flow lithium (Li) divertor configuration in the fusion nuclear science facility (FNSF) tokamak and identify operational regimes with acceptable divertor and core conditions. Plasma transport from the SOLPS-ITER code has been coupled with a liquid metal (LM) MHD/heat transfer code to model a Li open-surface divertor design and assess its impact on the scrape-off-layer (SOL) and core plasma performance. Simulations with only Neon (Ne) impurity seeding have been conducted to evaluate its impact on meeting FNSF design demands for the divertor and upstream plasma parameters. Simulation results indicate that Ne seeding significantly mitigates divertor heat flux but potentially reduces both upstream electron and main ion density due to fuel dilution. The combined application of Ne seeding and deuterium (D2) puffing is required to satisfy the FNSF design requirements on upstream density ( ne,sepOMP ∼1× 1020 m−3) and divertor energy flux ( q⊥,maxOdiv < 10 MW m−2). D2 puffing plays a role in counteracting upstream density drops and augmenting energy and momentum losses through atomic and molecular processes.The inlet Li flow velocity is systematically varied across a wide range to identify acceptable flows and corresponding LM surface temperatures. This comprehensive analysis identifies the acceptable Li flow parameters, LM surface temperature, and emitted Li fluxes necessary to meet the major design constraints. The emitted Li fluxes exhibit minimal impact on the main plasma at surface temperatures up to approximately ∼525 ∘C, corresponding emitted Li fluxes of up to φ Li ∼2 ×1023 atoms s−1. Uncertainties in the Li emission processes from the surface are also investigated, primarily influencing Li loss in the lower surface temperature range ( <525∘ C), with simulation results indicating a minor impact on the divertor and upstream plasma. Conversely, evaporation predominantly drives the Li loss processes at higher surface temperature ranges ( >525∘ C), contaminating both the divertor and upstream plasma.

Assessment of thermal and radiation induced creep in the dual cooled lead lithium blanket

The creep in recently designed blanket for the fusion nuclear science facility is described by decomposing it into thermal, irradiation and cavity swelling creep. Initial estimates of thermal creep using elastic high temperature rules showed that the recent heat transfer improvements at the first wall (FW) is feasible, and the FW can withstand heat flux of 0.35 MW/m2 at 3000 hrs without rupture. A viscoplastic model that is based on Norton's power law and relies on Multiphysics coupling of solid mechanics and heat transfer modules is used to capture the inelastic deformation in the blanket. Results showed that the design peak heat flux of 0.25 MW/m2 produced maximum thermal creep strain of 0.45 % and the relaxation of the thermal stress at the FW. Irradiation creep is prescribed to be proportional to the displacement damage dose and applied stress. The displacements from irradiation creep radially decrease from the FW to the back wall, but the maximum deformation is found at the back wall that is connected to the Helium manifold due to the high stresses at the region. Cavity swelling creep is dominant at steady state and combined with radiation creep to produce displacement of about 8 mm at the FW.

Magnetohydrodynamics in free surface liquid metal flow relevant to plasma-facing components

While flowing Liquid Metal (LM) Plasma-Facing Components (PFCs) represent a potentially transformative technology to enable long-pulse operation with high-power exhaust for fusion reactors, Magnetohydrodynamic (MHD) drag in the conducting LM will reduce the flow speed. Experiments have been completed in the linear open-channel LMX-U device [Hvasta et al 2018 Nucl. Fusion 58 01602] for validation of MHD drag calculations with either insulating or conducting walls, with codes similar to those used to design flowing LM PFCs for a Fusion Nuclear Science Facility [Kessel et al 2019 Fusion Sci. Technol. 75 886]. We observe that the average channel flow speed decreased with the use of conducting walls and the strength of the applied transverse magnetic field. The MHD drag from the retarding Lorentz force resulted in an increase of the LM depth in the channel that ‘piled up’ near the inlet, but not the outlet. As reproduced by OpenFOAM and ANSYS CFX calculations, the magnitude and characteristics of the pileup in the flow direction increased with the applied traverse magnetic field by up to 120%, as compared to the case without an applied magnetic field, corresponding to an average velocity reduction of ∼45%. Particle tracking measurements confirmed a predicted shear in the flow speed, with the surface velocity increasing by 300%, despite the 45% drop in the average bulk speed. The MHD effect makes the bulk flow laminarized but keeps surface waves aligned along the magnetic field lines due to the anisotropy of MHD drag. The 3D fringe field and high surface velocity generate ripples around the outlet region. It was also confirmed that the MHD drag strongly depends on the conductivity of the channel walls, magnetic field, and volumetric flow rate, in agreement with the simulations and a developed analytical model. These validated models are now available to begin to determine the conditions under which the ideal LM channel design of a constant flow speed and fluid depth could be attained.

Activation Analysis for the Inboard Region of FNSF Using SERPENT

This study describes an application of the SERPENT 2 code with the TENDL-2017 nuclear data library and the latest available model features of the Fusion Energy System Studies–Fusion Nuclear Science Facility (FNSF), to evaluate the activation of components after shutdown at 1, 10, and 100 years, assuming a plant lifetime of 8.5 full-power years. The primary parameters evaluated include the specific activity, decay heat, and waste disposal rating (WDR). The specific activity and decay heat are calculated with SERPENT 2 using a 360-deg model of the FNSF, while the WDR is calculated and classified based on the waste disposal limits established by the U.S. Nuclear Regulatory Commission under 10 CFR 61.55 as well as by using the Fetter approach. A python-based script developed for a previous high-level waste classification and analysis study was implemented and adapted to this research to calculate the WDR by comparing nuclide concentrations to the values established in 10 CFR 61.55 to generate a waste classification for each component surveyed. As only three short-lived isotopes have limitations for classifications beyond Class A, of which only 63Ni is present in appreciable quantities, there is a limit to the amount that short-lived isotopes contribute to the most significant waste analyzed here. In most cases, a handful of long-lived isotopes can be problematic, such as 59Ni and 94Nb, for example, which are solely responsible for multiple Class C classifications. The results herein reported heavily depend on the specific materials and mass/volume fractions in the specific model used in this study, which has changed and evolved since the inception of the FNSF concept and past studies. Therefore, the more significant contributions of this study may be the development of a modeling and simulation toolkit and a strategy to perform these calculations, so to help evaluate and optimize future fusion facilities.

Radiation Damage Analysis of FNSF Components Using McCad and MCNP

The Fusion Energy System Studies Fusion Nuclear Science Facility (FESS-FNSF) concept represents a transitional step between ITER and a commercial fusion power plant. The FNSF is a conceptualized D-T fueled tokamak with 518 MW of fusion power that has been extensively used to explore and optimize design features. The energetic 14.1-MeV neutrons can produce significant localized heating and activations, and can cause damage to plasma-facing components, which can determine maintenance/outage scheduling needs and also impact the lifetime of the device as a whole. This study illustrates a neutronics analysis that was conducted on a 22.5-degree symmetric sector of the FNSF with the goal of understanding the neutron heating and radiation damage that can be characterized by quantifying the displacements per atom (dpa). Concurrently, this study also focused on the development of analysis capabilities by converting a three-dimensional computer-aided design model of the FNSF into MCNP6.2 input using the McCad code. Accordingly, some confirmatory results on tritium production and the tritium breeding ratio (TBR) are provided to support model validation. The results produced by MCNP6.2 simulations showed that the highest heating and damage occurred in the outboard region, which concentrated approximately 290 MW of the total nuclear heating, in contrast to 97 MW within the inboard region. These results are consistent with previous studies that employed earlier versions of the FNSF concept and different modeling approaches. This study also provides additional details on neutron wall loading, as well as total heating from neutrons and gammas, results which show the total heating of the device (16 sectors) is approximately 477.83 ± 0.80% MW, indicating a neutron energy multiplication factor of 1.15. Additionally, the capability to calculate hydrogen and helium production, as well as dpa, is illustrated. Finally, the neutronics effects of using alternative materials to tungsten carbide were evaluated for the vacuum vessel, low-temperature shield, and structural ring components, which showed that compounds like YH2, Mg(BH4)2, and ZrH2 could reduce the total heating on the magnet and also reduce the TBR.

Multiphysics Modeling of Impurity Transport for FNSF Startup Scenario with ERO2.0

This study focuses on performing a multiphysics study using the ERO2.0 and UEDGE codes for two standard double null configurations for the Fusion Nuclear Science Facility: (a) 100% recycling and (b) 99% recycling. Results show that the main contributor to tungsten erosion along the divertor plates is impurities from the midplane waveguides. In addition, the standard high-recycling case (100% recycling) shows a significantly higher buildup of impurities along the divertor tiles during the startup phase, which can lead to a higher increase of energy loss in the plasma during steady-state operation. Last, for high recycling, anomalous diffusion can dominate over parallel field diffusion. The work performed in this study can be iteratively applied to a full operation scenario with additional physics such as those from neutrals, wall shaping, and additional external fields.

Numerical Analysis of Liquid Metal MHD Flow and Heat Transfer for Open-Surface Li Divertor in FNSF

Within the ongoing U.S.-based program on the development of liquid metal plasma-facing components, numerical simulations and analyses are performed to address the feasibility of the open-surface Li divertor. In the previous scoping studies (Smolentsev, 2021), heat-removal capabilities of the divertor were assessed using a simplified flow model for a slug-type velocity profile and constant flow thickness. Here, new analyses take into account forces acting on the flowing Li layer. Three reduced-order mathematical models are applied under the conditions of the U.S. Fusion Nuclear Science Facility (FNSF) to access magnetohydrodynamic (MHD) flow development effects, velocity distribution, and surface waves: 1) fully developed MHD flow; 2) quasi-2-D developing MHD flow; and 3) multiphase MHD flow. The obtained results for MHD flows and the surface heat flux computed with the plasma code scrape-off layer plasma simulation for ITER (SOLPS-ITER) are then used as input data to compute the temperature distribution in the divertor by solving the convection–diffusion energy equation.

Simulation of Liquid Lithium Divertor Geometry Using SOLPS-ITER

Plasma-facing component (PFC) geometries are evaluated for a Fusion Nuclear Science Facility (FNSF) design to find solutions that are compatible with flowing liquid lithium components while satisfying requirements from the perspective of plasma and neutral particle transport. Flowing liquid metal (LM) divertor systems offer important advantages due to continual regeneration of the surface material, but may require modifications from standard optimized tokamak divertor designs. Scrape Off Layer Plasma Simulator for the International Thermonuclear Experimental Reactor (SOLPS-ITER) simulations are used to compare standard vertical target, open, and balanced baffled geometries. It is found that the open geometry offers little control of neutral particles and has narrow operating windows where the divertor fluxes can be lowered to acceptable levels, while maintaining the required upstream density. The addition of neutral baffles allows greater control over the divertor conditions with less impact on the upstream density. The divertor neutral pressure in the baffled geometry is slightly higher than the open divertor, but not as large with a vertical target. A first-pass coupling of the SOLPS solution with sheath and implantation/erosion models to compute the lithium emission shows that the lithium is largely confined to the near-surface region due to the small ionization mean free path and large main ion flow toward the surface. The lithium source was varied by three orders of magnitude. At the highest rate, when the sourced lithium approaches the level of the main ion recycling flux, a high core lithium concentration (~10%) is observed with only moderate power dissipation to affect the divertor plasma.

Uncertainty Quantification of the Decay Gamma Source in Mesh-Based Shutdown Dose Rate Calculations

In fusion energy systems, part of the design effort is dedicated to the assessment of the shutdown dose rate (SDR) due to the decay photons that will be emitted from activated components. Monte Carlo transport codes are often used to obtain the neutron flux distribution in the problem domain. The neutron flux distribution is used in the rigorous 2-step (R2S) workflow to obtain the photon emission density distribution of decaying radionuclides. The photon emission density is then used as an input for a dedicated photon transport step to calculate the SDR. In this paper, the uncertainty of the decay gamma source due to the uncertainty of the neutron flux distribution in the R2S workflow is investigated. A scheme is developed to estimate the uncertainty of the decay gamma source, building on the concept of groupwise transmutation and using standard error propagation techniques. The applicability of the newly developed scheme is then demonstrated on one of the conceptual designs of the fusion nuclear science facility.

Divertor geometry modeling with the SOLPS-ITER code for reactor concepts with liquid metal divertors

• The SOLPS-ITER code has been applied to the analysis of the scrape-off layer (SOL) plasma associated with the fast flow Lithium (Li) divertor design for Fusion Nuclear Science Facility (FNSF). • Two flat divertor configurations are investigated and we found that simplified open geometry presents challenges for neutral control and impurity retention, but balanced geometry with baffling appears promising. • Lithium (Li) is sourced into the balanced configuration and scanning over a large range to understand what level of Li emission is needed to influence the divertor and upstream plasma conditions. • At low sourcing level ( ϕ Li ≤ 1 × 10 23 / s ), the Li shows a minor effect on the plasma and is well confined in the Private flux region (PFR) and divertor surfaces. • At high Li sourcing level ( ϕ Li > 1 × 10 24 / s ), the Li strongly affects the upstream and divertor plasma conditions and leads to strong detachment by increasing momentum and power losses. The SOLPS-ITER code has been applied to the analysis of the scrape-off layer (SOL) plasma associated with the fast flow Lithium (Li) divertor design for Fusion Nuclear Science Facility (FNSF). Two divertor configurations are investigated to determine puff levels that provide robust operational windows that meet the FNSF design requirements on upstream density and divertor flux with a flat divertor for flowing liquid metal (LM). Neon (Ne) is introduced to control the heat flux to the divertors while Deuterium (D 2 ) puffing is performed to maintain sufficiently high upstream density. We found that a simplified open geometry presents challenges for neutral control and impurity retention, but a balanced geometry with baffling is found to have acceptable puff operational windows. Lithium (Li) is sourced into the balanced configuration uniformly along the divertor surfaces and scanning over a large range to understand what level of Li emission is needed to influence the divertor and upstream plasma conditions. At low sourcing level ( ϕ Li ≤ 1 × 10 23 / s ), the Li shows a minor effect on the plasma and is well confined in the private flux region (PFR) and divertor surfaces. At moderate sourcing level ( 1 × 10 23 < ϕ Li ≤ 1 × 10 24 / s ), the Li starts affecting the upstream and divertor plasmas, and Li radiation and electron cooling are almost linearly increasing with the Li sourcing level. At high Li sourcing level ( ϕ Li > 1 × 10 24 / s ), the Li strongly affects the upstream and divertor plasma conditions and leads to strong detachment by increasing momentum and power losses. At high sourcing level, Li reaches above the x-point that can result in a substantial Li radiation in the core and sufficient dissipation to strongly affect the divertor and upstream plasmas.

Steady state thermo-mechanics and material property definition framework for analyzing DCLL blanket in the fusion nuclear science facility

A thermo-mechanics model that relies on creating the material property definition framework (MPDF) and multiphysics coupling of the heat transfer and the solid mechanics modules is developed to determine the structural integrity of the recently designed dual cooled lead lithium (DCLL) inboard blanket (IB) for the Fusion Nuclear Science Facility under steady state loads. The MPDF is called to supply fusion relevant neutron irradiation and temperature induced changes in material properties during multiphysics finite element runs, and PbLi temperature profiles are used to approximate Magnetohydrodynamics effect and the nuclear volumetric heating on the PbLi. Neutron irradiation and temperature induced reduction of the yield and ultimate strengths of F82H steel at the first wall (FW) are quantified for one year. A blanket in an assembly with gaps between blanket sectors and another blanket in an assembly with no gaps between blanket sectors, both exposed to radiation damage that lasted for one year are analyzed. Analysis using the elastic ITER structural design criteria for in-vessel components (ITER SDC-IC) design rules and a linear isotropic-hardening-type elastoplastic material model are used where most appropriate. The IB blanket with gaps between blanket sectors will withstand the steady state combined thermal and coolant loads for one year operational period but will fail if no gaps are allowed between blanket sectors. It is recommended that a gap of about 7.62 mm should be provided between IB blanket sectors during assembly which would close up during service, stop neutron streaming, reduce stresses and reduce bending of the FW into the scrape-off layer.

Plasma facing components with capillary porous system and liquid metal coolant flow

Liquid metal can create a renewable protective surface on plasma facing components (PFC), with an additional advantage of deuterium pumping and the prospect of tritium extraction if liquid lithium (LL) is used and maintained below 450 °C, the temperature above which LL vapor pressure begins to contaminate the plasma. LM can also be utilized as an efficient coolant, driven by the Lorentz force created with the help of the magnetic field in fusion devices. Capillary porous systems can serve as a conduit of LM and simultaneously provide stabilization of the LM flow, protecting against spills into the plasma. Recently, a combination of a fast-flowing LM cooling system with a porous plasma facing wall (CPSF) was investigated [A. Khodak and R. Maingi, Nucl. Mater. Energy 26, 100935 (2021)]. The system takes an advantage of a magnetohydrodynamics velocity profile as well as attractive LM properties to promote efficient heat transfer from the plasma to the LL at low pumping energy cost, relative to the incident heat flux on the PFC. In the case of a disruption leading to excessive heat flux from the plasma to the LM PFCs, LL evaporation can stabilize the PFC surface temperature, due to high evaporation heat and apparent vapor shielding. The proposed CPSF was optimized analytically for the conditions of a fusion nuclear science facility [Kessel et al., Fusion Sci. Technol. 75, 886 (2019)]: 10 T toroidal field and 10 MW/m2 peak incident heat flux. Computational fluid dynamics analysis confirmed that a CPSF system with 2.5 mm square channels can pump enough LL so that no additional coolant is needed.

Demonstration of RELAP5-3D for transient analysis of a dual coolant lead lithium fusion blanket concept

Our research aims to build a multiphysics framework for transient analysis of the Dual-Coolant Lead-Lithium (DCLL) blanket design of the proposed Fusion Nuclear Science Facility (FNSF) for all Anticipated Operational Occurrences (AOOs) using Reduced Order Models (ROMs). This paper presents an assessment of using RELAP5-3D for transient thermal-hydraulic blanket analysis of the DCLL design. Preliminary validation studies for the RELAP5-3D properties were conducted using a representative vertical flow loop model and a comparison study using experimental data from thermal convection corrosion loops at ORNL. These studies inspected the thermal hydraulic response of Lead-Lithium Eutectic (PbLi) systems within the RELAP5-3D code; mainly pressure drop and heat transfer, as compared to analytical calculations and experimental data. The Magnetohydrodynamic (MHD) pressure drop effect for liquid metal under the effects of a uniform magnetic field was implemented into our models using an equivalent forms loss coefficient based on the phasic momentum equations. We verified our implementation of MHD pressure drop within RELAP5-3D using a systematic approach consisting of two studies. Models were developed to perform comparison studies using experimental data from the Argonne Liquid metal EXperiment (ALEX) facility and the Magnetohydrodynamic PbLi Experiment (MaPLE) facility. Using the developed validation basis, we built a simplified DCLL channel model within RELAP5-3D, inspected MHD pressure drop and performed transient temperature analysis. The DCLL channel model used radial heating data from the International Thermonuclear Experimental Reactor (ITER) which was normalized to the FNSF design. A representative startup transient was performed using a relative power curve based on current PWR power ramp limits. This study demonstrated that RELAP5-3D ROMs are capable of accurate transient analysis of the DCLL blanket. We have since utilized MCNP to develop heating profiles based on the FNSF design and are working toward implementing them into our DCLL model.

Low cost, simpler HTS cable conductors for fusion energy systems

The cost and complexity of large, high-field superconducting magnet modules and related subsystems comprise 30% to 60% of the fusion reactor core capital cost. The strategic plan for the U.S. burning plasma research, the Fusion Energy Sciences Committee Report (FESAC) “Power the Future: Fusion and Plasmas’', and 2021 NASEM report “Key Goals and Innovations needed for a U.S. Fusion Pilot Plant” recommends that the U.S. pursue innovative science and technology to enable construction of a Fusion Pilot Plant (FPP) that produces net electricity from fusion at reduced capital cost. To achieve this, a novel combination of lower-cost high temperature superconductors (HTS) in cable configurations with co-wound reinforcement for higher current density are being investigated using a simplified construction strategy to produce compact stable coils. They would be capable of generating 20 T at up to 10-20 K. Small-scale, inexpensive test coils and prototypes will help develop each feature and validate cabled conductor design models. The near term goal is to validate engineering approaches, scientific models and fabrication capabilities applicable to fusion reactor development such as U.S. fusion nuclear science facility (FNSF), sustained high-power density tokamak facility (SHPD) and FPP designs. The design options include lower-cost, high-strength, quench resistant REBCO or Bi-2212 cables in an all metal coil design that simplifies HTS coil construction and quench protection system, with co-wound reinforcements that integrate stress management in HTS cable design and provides thermal mass to help prevent quench damage.

A cryostat for a 6 T conduction-cooled, no-insulation multi-pancake HTS solenoid

There is a growing requirement for high-field (>20 T) magnets capable of continuous operation, driven by the needs of both fundamental research and technological advance, particularly in application to an eventual pilot plant for magnetic confinement fusion. Even with HTS windings, such magnets will still require cryogenic cooling, and liquid helium (LHe) immersion, the typical solution to this problem, adds significantly to the operating expenses of such facilities. This reality makes cryogen-free cooling systems a necessity in future high-field magnet systems. The Princeton Plasma Physics Laboratory (PPPL) is exploring conduction-cooling systems of HTS pancake solenoids for a scanning tunneling microscopy (STM) facility at Princeton University, and potentially also for the central solenoid of the Fusion Nuclear Science Facility (FNSF). To these ends, PPPL is designing a cryostat to evaluate the thermal stability of a 5-6 T, 30 double-pancake (DP) REBCO insert coil of 40 mm ID / 70 mm OD, and smaller prototypes, operated in self-field with conduction cooling provided by a 2-stage GM cryocooler. The current design is expected to achieve 1st and 2nd stage temperatures of 44 K and 4-10 K, respectively, with the resistivity of DP-DP solder joints being the principal source of uncertainty in 2nd stage temperature predictions.

Design window for open-surface lithium divertor with helium-cooled substrate

This study for an open-surface flowing lithium (Li) divertor with a thin helium-cooled substrate includes a scaling analysis, numerical computations and an analytical solution to access the design window for a low vaporization divertor option where the maximum temperature at the Li free surface is limited to 450 °C. Three heat removal regimes were identified in the scaling analysis and then accessed in numerical computations: convection-dominated, diffusion-dominated and convection-diffusion. Of them, convection-dominated regime, in which all applied heat flux is removed by a high velocity Li layer, offers the highest heat removal capability. For a generic divertor, a high surface heat flux of 10 MW/m2 can be removed with the Li flowing at 7.3 m/s. Under the conditions of the US Fusion Nuclear Science Facility (FNSF), the required Li velocities are 6.5 m/s for the inner leg, and 15 m/s for the outer leg. The obtained analytical solutions for the diffusion-dominated regime, where the incident heat flux is removed by helium, suggests the lowest of the three regimes heat removal capability in the range of 0.5–1.2 MW/m2.

The Compactness and Inboard Radial Build of Fusion Nuclear Devices

Restraining the size of fusion power plants is considered an important avenue to make them a competitive energy source among other forms of energy production. The most critical contributor to the size of a tokamak is the inboard radial build, composed of multiple components with various functions. This build is the ultimate limit to size reduction. The Fusion Nuclear Science Facility is reviewed and each element of the inboard build is described, showing that the build, including breeding blanket, structural ring, vacuum vessel, low-temperature shield, and toroidal field and central solenoid (CS) coils, contributes 2.9 m of build, with 0.6 m of bore hole inside the CS coil, or 3.5 m to reach the plasma scrape-off layer. This implies that it would be challenging to make a significantly smaller build and simultaneously meet all the engineering requirements.

Design, construction, and testing of no-insulation small subscale solenoids for compact tokamaks

Fusion energy systems studies (FESS) for next-step devices based on the most promising magnetic configurations indicate that high magnetic fields and high current density for magnet coil systems may reduce device size and lower the cost. High current density and radiation resistant fusion magnets are particularly beneficial for low cost, low aspect ratio compact reactor designs such as Fusion Nuclear Science Facility (FNSF), Fusion Pilot Plant (FPP) of low-aspect ratio spherical tokamak (ST) or compact stellarators. Unlike typical high field magnets for nuclear magnetic resonance (NMR), magnetic resonance imaging (MRI), or high energy physics research applications, neutron irradiation damage to organic insulations in the coil winding pack is a critical issue for next generation fusion reactors where orders of magnitude higher neutron fluence than those in present experimental reactors such as ITER are expected. Moreover, a high coil winding pack current density is needed for high field magnets in low cost, compact radial build next step fusion reactors. The slow current charging time, however, is an issue in fully non-insulated coils. We present the design, construction and testing of subscale Nb3Sn solenoids, up to half the diameter of the central solenoid (CS) for National Spherical Torus Experiment (NSTX), with and without inter-layer insulation for enhancing radiation resistance while improving the winding pack current density and coil performance. The major radius for the NSTX/NSTX-Upgrade (reference for compact ST) is 0.854/0.934 meter, and 4.8 meters for FNSF. The magnetic field for NSTX/NSTX-U at the plasma center is 0.6/1 Tesla and 9 Tesla for FNSF. The inner diameter of the central solenoid coil is 0.2/0.4 meter for NSTX/NSTX-U and 1.2 meters for FNSF. The current charging and discharging behavior of the prototype solenoids was investigated to quantify performance advantages of a simplified coil fabrication without error-prone vacuum pressure impregnation (VPI) for the removal of epoxy organic insulation. Scalability of the no insulation concept for fusion can be addressed via engineered coating for novel insulation on contact resistance or intra-layer no insulations. Radiation resistance and coil winding efficiency were significantly improved in the no insulation coils for next step compact fusion reactors.

Status and Progress of Liquid Metal Thermofluids Modeling for the U.S. Fusion Nuclear Science Facility

At present, the U.S. Fusion Engineering Systems Study (FESS) considers several cooling/breeding concepts that utilize flowing liquid metals (LMs), Li, or eutectic PbLi alloy as working fluids for implementation of these concepts in the Fusion Nuclear Science Facility (FNSF). In this paper, we review recent modeling activities aimed at the investigation of LM flows and heat transfer relevant to the FESS-FNSF program. In particular, considerations are given to (1) development and validation & verification of computational magnetohydrodynamic (MHD) codes, (2) characterization of critical coupled MHD/heat transfer phenomena, and (3) design and analysis for selected LM applications in the FNSF. Under these three research thrusts, the reviewed topics including the MHD code HyPerComp Incompressible MHD solver for Arbitrary Geometries (HIMAG), MHD mixed-convection flows, MHD pressure drop in the blanket inlet/outlet manifolds, PbLi flows in a thermal convection loop, MHD PbLi flows in a dual-coolant lead-lithium blanket prototype, and a design window for a flowing Li divertor with He-cooled substrate.

Neutronics Calculations to Support the Fusion Evaluated Nuclear Data Library (FENDL)

In the design of fusion reactors, determining radiation levels due to neutrons and photons (gammas) throughout the reactor and its surroundings is important. Radiation transport codes need to have accurate cross-section libraries in order to produce accurate results. The Fusion Evaluated Nuclear Data Library (FENDL) is an international effort coordinated by the International Atomic Energy Agency, Nuclear Data Section, that assembles a collection of the best nuclear data for fusion applications. In the current FENDL-3.1d data library, neutron cross sections for 65 of the 180 isotopes present in the library come from ENDF/B-VII.1. Monte Carlo–based neutronics calculations using cross-section libraries from FENDL (versions 2.1 and 3.1d), ENDF/B (versions VII.1 and VIII.0), and candidate new evaluations for key structural elements/isotopes such as iron and chromium were performed. The calculations were performed in reactor-relevant models including a one-dimensional (1-D) cylindrical model of ITER, a three-dimensional (3-D) computer-aided design (CAD)–based model of ITER, and a 3-D CAD-based model of the U.S. Fusion Energy System Studies Fusion Nuclear Science Facility (FNSF). The results show that neutron fluxes calculated with different cross-section libraries can be as much as 12% higher and as much as 8% lower than those calculated with the reference cross-section library (FENDL-2.1). Nuclear heating calculated with different cross-section libraries can be as much as 14% higher and as much as 8% lower than those calculated with the reference cross-section library. Iron displacements per atom calculated with different cross-section libraries can be as much as 9% higher and as much as 9% lower than those calculated with the reference cross-section library. Helium production calculated with different cross-section libraries can be as much as 19% higher and as much as 2% lower than those calculated with the reference cross-section library. Tritium production in the ITER 1-D model’s nonbreeding regions calculated with different cross-section libraries can be as much as 246% higher and as much as 5% lower than those calculated with the reference cross-section library. The tritium breeding ratio in the FNSF 3-D model calculated with different cross-section libraries averaged 1% higher at the inboard and 1.4% higher at the outboard than those calculated with the reference cross-section library.