SEARCH · Engineering Papers

Results for “Fusion reactors”

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 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2025

MX precipitate behavior in an irradiated advanced Fe-9Cr steel: Self-ion irradiation effects on phase stability

In an effort to optimize Fe-9Cr reduced activation ferritic/martensitic (RAFM) steels and to inform the design and operation of fusion reactors, this work represents the first in a series of cohesive studies dedicated to the evolution of MX-TiC precipitates under accelerated single and dual ion irradiations. This study investigates CNA9, a simplified Fe-9Cr RAFM steel featuring initial MX-TiC precipitate densities of (2.3±0.3)×10²¹ m⁻³. This material was subjected to single self-ion irradiation at damage levels ranging from 1 to 100 displacements per atom (dpa) over a temperature range of 300 to 600°C, with a nominal dose rate of 7×10⁻⁴ dpa/s. Irradiation-induced coarsening was observed, as evidenced by statistically significant increases in mean diameter sizes, at 15 dpa at both 500°C and 600°C, whereas no coarsening was noted at 300°C or 400°C. Further, complete dissolution of precipitates occurred at damage levels of 50 and 100 dpa across the two temperatures tested (300°C and 500°C) while no significant changes were observed at any doses below 15 dpa at 500°C. Experimentally parameterized recoil resolution modeling suggests that the observed radiation stability of MX-TiC precipitates is intricately linked to diffusional changes of solutes resulting from the co-evolution of microstructural features within the experiments. The findings align with current theoretical perspectives on radiation-induced precipitate stability in complex alloys.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Innovative laser-based methods for monitoring fuel retention in ITER

This paper addresses the challenge of tritium inventory management in ITER and future fusion reactors, highlighting the importance of accurate tritium measurement and its spatial distribution within the vacuum vessel. Given ITER’s operational constraints, especially the limit on tritium retention, precise measurement is essential for both safety and regulatory compliance. To tackle these questions, the paper presents the T-monitor diagnostic system developed by Forschungszentrum Jülich, which uses Laser-Induced Desorption (LID) in combination with Diagnostic Residual Gas Analysis (DRGA) to measure hydrogen isotope concentrations on the surface of divertor tiles. The system integrates a high-power laser, advanced optical components, and a Fast Scanning Mirror Unit (FSMU) for accurate laser spot positioning with rapid response. Designed to measure in situ tritium retention, the diagnostic provides high-resolution spatial mapping, vital for evaluating detritiation strategies. The laser heating process increases the divertor surface temperature to 1600 K within the laser spot, promoting hydrogen isotope desorption. Accurate measurements require the precise control of laser parameters, including pulse duration and spot size, with a target relative accuracy of 20%. The optical design includes both in-vessel and ex-vessel components, such as durable high-reflectivity mirrors made of gold and copper, selected not only for their infrared performance but also for their transmission of visible wavelengths for observation purposes. To protect optical components from contamination, a pneumatic shutter is used.

Hydrogen isotopes↗

DOE OSTI2024

Density drop at the divertor target in the prototype material plasma exposure eXperiment (Proto-MPEX)

The steady-state linear device “Material Plasma Exposure eXperiment” (MPEX) is currently under construction at Oak Ridge National Laboratory with the goal of enabling Plasma-Material Interaction studies at future fusion reactor relevant plasma conditions. In this work, a newly in-house developed hybrid Particle-In-Cell code-PICOS++ is applied to understand the experimental results obtained from the prototype of MPEX referred to as the “Proto-MPEX” during its helicon-only and helicon with ion cyclotron resonance heating (ICRH) experiments. This study explains the physics of the experimentally observed plasma density-drop at the divertor target in Proto-MPEX device during ICRH. In contrast to previous work on ICRH in MPEX, this study demonstrates that the mirror force plays a central role in the Proto-MPEX plasma transport during ICRH, which has new features not previously explored. Force balance analyses reveal that the temperature anisotropy produced by ICRH leads to a significant increase in the mirror force downstream of the resonance where the magnetic field is diverging. This force accelerates ions toward the target and leads to a drop in plasma density to ensure conservation of particle flux. Simulations with ICRH where the magnetic field divergence downstream of the resonance has been removed, do not produce plasma acceleration nor density drop at the target despite efficient ion heating at the resonance. Moreover, simulation results demonstrate that for a given ICRH power, lowering the source rate produces ions with increased perpendicular energy which interact with the mirror force to produce higher plasma acceleration which increases the strength of the density-drop at the target. The strength of the density drop appears to reach an asymptotic limit at a certain threshold ICRH power. Simulations show that this threshold power increases with increasing particle source rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Turbulence link to L-mode, I-mode, and H-mode confinement in the DIII-D tokamak

Understanding the physics of low-confinement (L-), improved-confinement (I-), and high-confinement (H-) modes is critical for fusion reactors. The finding herein reports observations of two types of turbulence coexisting near the L-mode edge, one magnetohydrodynamic (MHD)-like and another micro-tearing mode (MTM)-like, linked to the H-mode and I-mode confinement in the DIII-D tokamak. Ion-scale magnetic and density turbulence is measured using a Faraday-effect radial-interferometer-polarimeter and beam-emission-spectroscopy (BES). Broadband turbulence spectra of up to ~600 kHz are observed in two discharges where transitions between L-mode, I-mode, and H-mode occurs. Turbulence is found to be inversely correlated with confinement, meaning lower turbulence power at higher confinement. Distinctively, the high-frequency (HF, >~100 kHz) magnetic turbulence power changes by the most (55%) during transitions primarily involving energy confinement change, whereas the low-frequency (LF, <~100 kHz) magnetic and density turbulence power changes by the most (80%) during transitions primarily involving particle confinement change. The LF turbulence amplitude oscillates with and leads to deuterium-alpha emission oscillations before an H-mode. These results imply that HF turbulence mainly affects energy confinement whereas LF turbulence can affect particle confinement. The magnetic and density turbulence exhibits coherence up to 0.6 and cross-phase magnitude close to π/2 in most cases, suggesting they have a common origin in both the LF and HF ranges. BES suggests that LF turbulence resides at the edge (ρ = 0.95) and HF turbulence can be at the outer core (ρ = 0.8) or edge (ρ = 0.95). Comparisons of measurements, theory, and gyrokinetic simulations suggest that HF turbulence is MTM-like in all cases, whereas LF turbulence is more consistent with MHD-like modes and the exact instability might change during transitions—except that a drift-wave origin is possible in a low collisionality H-mode. These results suggest that the H-mode involves suppressed MHD-like turbulence, whereas the I-mode mitigates MTM-like turbulence along with largely unchanged MHD-like turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Taming harmful bursts and heat flux in high-confinement tokamak plasmas

A major challenge in tokamak fusion research is first-wall erosion caused by steady heat loads and sudden energy bursts known as edge-localized modes. Divertor detachment reduces steady-state heat flux, while resonant magnetic perturbations can suppress these instabilities. However, integrating the two has been difficult because they require conflicting operating conditions. Here we demonstrate simultaneous achievement of resonant magnetic perturbations mitigated small edge-localized modes and impurity seeded partial divertor detachment in plasmas with an ITER-similar shape on the DIII-D tokamak. Experiments and simulations show that resonant magnetic perturbations facilitate detachment by redistributing particles, lowering the core density and increasing the scrape-off layer density, thereby reducing the amount of injected gas required. Cooling-gas injection eliminates the secondary heat-flux peak created by three-dimensional magnetic lobes, while edge cooling weakens the plasma response to the applied magnetic fields. These advances illustrate a viable pathway for integrating edge stability control with power exhaust in future fusion reactors.

magnetically confined plasmas↗

DOE OSTI2025

Understanding the formation of a low-pressure pedestal in the presence of a strong internal transport barrier in DIII-D high β p plasmas

As a promising scenario for fusion reactors, the high poloidal-beta (β P ) scenario is characterized by a strong large radius internal transport barrier (ITB), which significantly enhances the overall confinement quality and the bootstrap current fraction for fully non-inductive operation. It is frequently observed that in the presence of a strong ITB, the pedestal height is lower and is accompanied by small edge localized modes (ELMs), which further improves the compatibility of a high performance core with an edge solution. A mechanism for the formation of the low pedestal is proposed in this paper. It is found that the strong ITB creates an off-axis bootstrap current to clamp the local safety factor q, and thus the magnetic shear in the outer core/pedestal region is increased. Gyrokinetic simulations with the CGYRO code show that the higher magnetic shear brings the experimental profiles into the range where the growth rate of drift-wave instabilities and thus transport is higher, and therefore a lower pedestal gradient is expected. Here, the combination of low pedestal and high magnetic shear further enhances the turbulent transport across the whole pedestal, consistent with power balance analysis. Such a positive feedback mechanism ultimately results in a lower pressure pedestal as observed in experiments. Under such a low pedestal, linear simulations with BOUT++ predict the growth rates of peeling–ballooning modes to be lower across the whole toroidal mode number spectra, and the nonlinear BOUT++ simulation exhibits lower saturated fluctuation intensity as well, consistent with the experimentally observed lower ELM size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Quantum Information for Fusion Energy Sciences (Final Technical Report)

The simulation of plasma dynamics is a critical area of Fusion Energy Sciences (FES) due to it’s usefulness in predicting, controlling, and confining plasmas in the context of potential fusion reactors. The simulation of plasmas is a computationally difficult problem in both classical and quantum physics, motivating investigation into the potential of quantum computers to simulate these systems. This project took several concrete steps towards this goal by developing tools for improving the control, characterization, and calibration of quantum gates on a superconducting quantum computer, developing error suppression and mitigation tools to reduce errors on the quantum computer, and utilizing these advancements to simulate reduced models of plasma dynamics on the quantum computer. In order to efficiently simulate plasma physics, an optimal control method which synthesizes, directly at the pulse level, any quantum gate on qubit and qutrit systems was developed. Using four superconducting transmon quantum processors at Rigetti and LLNL, it was demonstrated that any arbitrary quantum gate on qubits and qutrits could be implemented with high fidelity, leading to a significantly reduced length of a gate sequence. A problem of interest in FES is the nonlinear optical process of laser pulse compression within a plasma. Since quantum physics is linear, simulating nonlinear operations is not naturally feasible on a quantum computer, however it is possible to simulated a quantized version of the nonlinear process. A quantization approach to convert nonlinear wave-wave interaction problems to Hamiltonian simulation problems was developed and demonstrated using two qubits on a Rigetti device. In this experiment, a number of error suppression and mitigation techniques were investigated to determine how best to utilize the finite quantum resources. This study provides an example of how plasma problems may be solved on near-term, noisy quantum computing platforms and identified a promising set of techniques. Building on the insights of these experiments, the investigation turned to linear electron-plasma wave physics. A connection was identified between a local one-dimensional lattice spin model and linear wave phenomena, allowing a plasma physics problem to be efficiently mapped to the quantum computer. In this framework, reflection and transmission of plasma waves at a sharp boundary was studied, as well as the propagation of waves through an inhomogeneous plasma medium. In addition to the suite of error suppression and mitigation techniques developed, this experiment introduced the use of a digital-analog gate scheme designed to efficiently simulate the plasma Hamiltonian. With hardware available at the conclusion of the project, simulation at the scale of 9 qubits and 15 timesteps (60 entangling layers) was achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Assessing Stainless Steel Compatibility in Flowing FLiBe Salt

To address the considerable interest in LiF-BeF 2 (FLiBe) compatibility for fission and fusion reactor applications, static and flowing compatibility experiments were conducted to assess the compatibility with type 316H stainless steel. In static testing at 550° and 650 °C, small mass changes were measured and posttest characterization of the FLiBe showed increased levels of Fe, Cr, Ni and Mn in the salt. Adding Be in the static salt test reduced the dissolution of Fe and Ni. Further, an initial assessment of mass transfer in flowing FLiBe without a Be addition was conducted using a monometallic 316H thermal convection loop (TCL) operated for 1000 h with a peak temperature of 650 °C. Similar to prior results in flowing FLiNaK salt, the 316H specimens exhibited small mass losses in the hot leg. Posttest characterization of the 316H specimens suggested Cr surface depletion in the hot and cold legs and possibly Fe deposition in the cold leg. To further understand this behavior, Cr and Fe dissolution was measured in static FLiBe at 550–650 °C.

316H stainless steel↗

DOE OSTI2025

Deuterium trapping mechanisms in reduced activation ferritic martensitic steels and their correlation with mechanical strengthening

Development of high-strength materials often involves introduction of additional strengthening microstructures that also serve as tritium trapping sites. Such additions in fusion material development could degrade the fuel efficiency in fusion reactors and raise radiological concerns. The contribution of individual microstructure features in hydrogen trapping must be evaluated to ensure fuel efficiency and radiological safety. This study explores the mechanistic origins of deuterium trapping in reduced-activation ferritic–martensitic steels and its correlation to mechanical strengthening. A series of model alloys and engineering steels were fabricated and subjected to different heat treatments to control deuterium trapping site density. Deuterium retention was evaluated using D 2 gas charging and thermal desorption spectroscopy, focusing on the role of grain boundary, dislocation, M 23 C 6 precipitates, and TiC precipitates. Multiscale microstructure characterization and synchrotron X-ray diffraction were performed to characterize microstructure, which was correlated to the deuterium retention property. Results show that TiC precipitates exhibit the highest deuterium trapping capacity, followed by M 23 C 6 precipitates. Dislocation and grain boundary demonstrate the lowest and similar efficiencies. Furthermore, the relationship of trapping quantity and mechanical strengthening of these microstructure features was quantified, demonstrating that TiC precipitates offer highest deuterium trapping per unit of mechanical strengthening.

Retention↗

DOE OSTI2025

Vacuum vessel design with lofted toroidal surfaces for a QHS configuration

A vacuum vessel in between the hot plasma and the magnetic coils is crucial for a fusion reactor. In this study, a tool is presented to perform the initial design of such a vessel automatically. Starting from the last closed flux surface, an algorithm is developed which lofts this manifold out as much as possible. At the same time, the algorithm makes sure that no intersection with the coils takes place. In the presented work, this algorithm is applied to the HSX stellarator and it has been verified that the resulting vessel wall can be built from an engineering perspective. Using field line tracing, the strikepoints on the vessel wall indicate the initial locations for a divertor.

HSX↗

DOE OSTI2025

Implementation of high-speed data acquisition at DIII-D

Research at the DIII-D National Fusion Facility in San Diego focuses on short pulse plasma discharges that specialize on various shaping profiles. High-speed data collection is a critical component for the operation of many of DIII-D’s diagnostics and is fundamental for capturing high-resolution data used in experimental data analysis. Differing techniques enable the plasma control system (PCS) to perform complex real-time feedback control on microsecond time scales. This work presents a comprehensive overview of data acquisition, focusing on the hardware and software used in reliable data acquisition at DIII-D. The robust nature of the data acquisition system allows for various techniques to coexist seamlessly. However, as modern systems capable of nanosecond resolution become more common, existing architectures will need to be modified. Here, by addressing the key challenges of high-speed data acquisition, DIII-D is able to provide real-time data used in plasma operation and has the ability to acquire high fidelity data needed for future experimental fusion reactors, such as ITER.

Control↗

DOE OSTI2026

Corroborating VNA and thermal measurements of transmission loss on the DIII-D ECH waveguide system

Electron cyclotron heating (ECH) and current drive (ECCD) will play a large role in tokamak-based fusion reactors. At the DIII-D tokamak, 110 GHz microwaves injected into the plasma can provide core heating and current drive as well as impurity control, neoclassical tearing mode mitigation, and breakdown assistance. Understanding the physics of these processes relies on accurate estimates of injected ECH power. DIII-D’s ECH system consists of six MW-class Microwave Power Products (MPP) gyrotron microwave sources. Operating the gyrotrons far from the tokamak removes them from magnetic field interference, so 31.75 mm inner-diameter corrugated waveguides transmit the microwave power the 80 m from the gyrotrons to steerable launchers in the tokamak chamber. Estimates of injected power rely on knowing the generated power at the source and then subtracting transmission loss. Conventional transmission loss measurements based on calorimetric dummy loads are onerous and only possible during extended maintenance periods. This work examines two tools that provide more flexibility for the transmission loss measurements. Furthermore, a resistive temperature detector (RTD) array installed along a waveguide measures heat lost to the transmission line, and low power time domain reflectometry (TDR) measurements with a vector network analyzer (VNA) allows loss measurements without burdensome hardware modifications.

ECH↗

DOE OSTI2024

Changes in dislocation punching behavior due to hydrogen-seeded helium bubble growth in tungsten

The accumulation of gas atoms in tungsten is a topic of long-standing interest to the plasma-facing materials community due the metal's use as a divertor material in some tokamak fusion reactors. The nucleation and growth of He/H gas bubbles (along with their isotopes) can result from impinging fluxes of these gases which give rise to damage at the W divertor surface. The inclusion of He or H in W has been studied extensively by the community, finding that He bubbles modify the surface through periodic dislocation punching and bursting mechanisms while H bubbles impact the metal through plastic-strain induced material failure. However, the mechanisms which are present during the combined flux of both He and H is not well-studied atomistically. Motivated by this, an atomistic modeling study is conducted using molecular dynamics to assess the behavior of mixed concentration He:H bubbles in W. Here we find that the introduction of H into a growing He bubble results in a dramatic change in the nature and presence of dislocation loops which are typically generated via dislocation punching in over-pressurized He bubbles. Most notably, at high H concentrations, there is a switchover in energetic favorability from glissile 1/2<111> dislocations to sessile <100> dislocations. This thermodynamic crossover could imply a significant reduction in W surface morphology changes than with pure He bubbles and, additionally, we show this to have implications on the trapping of H in the bubbles and their associated dislocations.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Interpreting experimental measurements of helium bubbles using stochastic cluster dynamics models of heterogeneous nucleation and growth in irradiated ferritic alloys

Among a number of other advantageous features, ferritic/martensitic steels are being considered as fusion reactor structural materials due to their low intrinsic swelling under irradiation. However, under high-energy neutron irradiation, He produced through (n, α) reactions stabilizes vacancy clusters, which then act as seeds for bubble formation and growth, negating the intrinsic swelling resistance of these alloys. Standard models of irradiation damage accumulation and microstructural evolution consider homogeneous nucleation as the basis for bubble formation and growth. However, this generally does not explain the large bubble densities and sizes observed experimentally under a number of different conditions. Here, we propose a more realistic physical model of bubble nucleation, formation, and growth designed to capture recent experimental measurements of He-bubble formation and evolution during co-implantation of He and Fe ions in model ferritic alloys. We find that experimental results are explained only when the following three features are considered simultaneously: (i) heterogeneous nucleation of He-vacancy bubbles at defect sinks (e.g., dislocations, grain boundaries, and second-phase precipitates), (ii) nucleation and growth of bubbles via the ‘trap mutation’ mechanism (i.e., spontaneous production of Frenkel pairs due to absorption of He atoms), and (iii) transition from a growth-limited, He-stabilized bubble regime to a ‘runaway’ void-type growth scenario. The model is implemented into a stochastic cluster dynamics framework that takes advantage of cluster size grouping methods to accelerate the simulations, allowing us to reach 10 dpa of simulated irradiated dose, and to capture cluster sizes in excess of 20 nm. Further, a careful extrapolation exercise conducted assuming classical nucleation theory leads to excellent agreement with the experimental measurements at 50 dpa of irradiation.

36 MATERIALS SCIENCE↗

DOE OSTI2025

The interaction of edge dislocations with hydrogen-helium bubbles in tungsten

In fusion reactors, plasmas-facing components undergo degradation due to displacement damage and injection of gas impurities. Significant amounts of hydrogen (H) and helium (He) impurities can be introduced into materials through plasma exposure and nuclear transmutation, and their synergistic effects lead to the formation of mixed He/H bubbles under irradiation, changing their hardening mechanisms. In this study, molecular dynamics (MD) models were developed to study the interaction between ½<111> edge dislocations with mixed He/H-bubbles in W. The critical resolved shear stress (CRSS) that is required for dislocation breakaway increases with He/H-to-vacancy ratio at 600 K, indicating strong pinning effects of mixed He/H bubbles. However, as the temperature increases up to 1400 K, the mixed He/H bubbles become unstable as the H atoms are increasingly emitted from the mixed He/H bubbles and migrate into W matrix, which lowers bubble pressure and CRSS. Overall, these results highlight the synergistic effects of He and H on dislocation behavior in irradiated W during deformation.

36 - MATERIALS SCIENCE↗

DOE OSTI2026

Molecular dynamics simulations of reflection and sputtering behavior of boron under deuterium ion irradiation

Boronization is a commonly used method of wall conditioning in fusion reactors. The application of boron films to the plasma-facing materials results in enhanced plasma performance due to the reduction of intrinsic impurities. This is primarily driven by a reduction in oxygen content that is chemically trapped in the boron film. The reactive nature of these boron films also raises questions concerning interactions with hydrogen isotopes. In this work, boron-deuterium interactions were studied using molecular dynamics (MD). Reactive force field potentials were used to model the chemical interactions between B and D. An amorphous boron substrate was irradiated by D atoms at varying incident energies, 10 eV < E i < 150 eV and angles, 0° < α < 85°. The reflection probability was calculated and compared to results from the commonly-used binary collision approximation (BCA) method. This comparison found that the BCA underestimated the reflection probability at E i < 35 eV and α > 45 ∘ . The source of this discrepancy was found to be the surface binding energy model. The BCA calculation with an isotropic surface binding energy model was more closely aligned to the MD result. This, in combination with a correction function based on the MD results allows for corrections to the reflection probability of deuterium impinging on boron surfaces. The sputtering of the substrate material was also studied. While this study did not contain sufficient events to quantitatively describe the sputtering behavior, some qualitative results emerged: namely, chemical sputtering of B and D-containing molecules (BD, BD 2 , BD 3 ) at low ( < 20 eV) incident deuterium energies. This result suggests that chemical sputtering could be a significant factor in limiting boron coating lifetime when exposed to lower ion energies, such as those in detached plasmas. The results show that chemical interactions should be taken into account when modeling the interactions between D ions and B surfaces.

Boron↗

DOE OSTI2026

In situ studies on heavy ion irradiation and partial oxidation induced stacking faults and nanograins in tungsten nanolaminates

Understanding the microstructural evolution of tungsten (W) under extreme irradiation environments is critical for its application as a plasma-facing material in future fusion reactors. Heavy ion irradiation study offers expedited damage accumulation and thus allows irradiation tolerance property prediction in a short period of time. In this study we explore in situ Kr ion irradiation of W nanolaminates at 800 °C and uncover a complex interplay of irradiation-induced transformations, including the emergence of stacking faults, nanograin formation, and generation of thickened grain boundaries defined as GB regions with the thickness of up to 15 nm. High-resolution transmission electron microscopy studies reveal the formation of extended planar faults and a metastable hexagonal close-packed (hcp) phase within the body-centered cubic (bcc) matrix. These structural transitions are facilitated by irradiation-induced shear. They are further stabilized by the presence of oxygen, which promotes stacking fault formation and vacancy trapping. Elevated temperatures enhance defect mobility, enabling dynamic recrystallization into ultrafine grains and the thickening of GBs due to defect absorption and impurity segregation. These findings unravel non-equilibrium phase transformation pathways in heavy ion irradiated W and highlight the critical role of impurity-mediated defect dynamics in governing its radiation tolerance properties.

Wazeer, Adil [Purdue University]↗

DOE OSTI2025

Effect of alloying on intrinsic ductility in WTaCrV high entropy alloys

Tungsten (W) exhibits desirable properties for extreme applications, such as the divertor in magnetic fusion reactors, but its practicality remains limited due to poor formability and insufficient irradiation resistance. In this work, we study the intrinsic ductility of body-centered cubic WTaCrV based high entropy alloys (HEAs), which are known to exhibit excellent irradiation resistance. The ductility evaluations are carried out using a criterion based on the competition between the critical stress intensity factors for emission (K Ie ) and cleavage (K Ic ) in the {110} slip planes and {110} crack planes, which are evaluated within the linear elastic fracture mechanics framework and computed using density functional theory. The results suggest that increasing the alloying concentrations of V and reducing the concentrations of W can significantly improve the ductility in these HEAs. The elastic anisotropy for these HEAs is analyzed using the Zener anisotropy ratio and its correlation with the concentration of W in the alloys is studied. Results indicate that these alloys tend to be fairly isotropic independently from the concentration of W in them. The computed data for the elastic constants of these HEAs is also compared against available experimental data. The results are in good agreement, validating the robustness and accuracy of the computational methods. Multiple phenomenological ductility metrics were also computed and analyzed against the analytical model. Some metrics, mainly the surrogate D parameter, show a good correlation with the Rice model. The potential of these empirical metrics to serve as surrogate screening models for optimizing the compositional space is also discussed.

Anisotropy↗