Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fusion experiments”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,279 records · Page 71

Investigation into the transport properties of planetary interiors through inelastic X-ray scattering experiments and quantum molecular dynamics (Final Technical Report)

This was a two-year research project for the period 08/15/2018 - 08/14/2020 conducted at the University of Nevada, Reno using Pronghorn, a newly built high-performance cluster located at the University. The short-term project goal was to provide computational support for the experimental campaigns at LCLS through data analysis of previous and upcoming scheduled experiments, along with performing state-of-the-art atomistic simulations of dense plasmas. The longer-term goal was to create a new computational high energy density physics group located at the University of Nevada, Reno. Guided by results from past and future high-resolution scattering experiments, we aimed to research and develop non-equilibrium and non-adiabatic atomistic simulations of dense plasmas that went beyond the Born-Oppenheimer approximation

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reconstruction of plasma equilibrium and eddy current induced error fields in the Lithium Tokamak eXperiment-β (Final Report)

This is the final report for the DOE award DE-SC0019239, titled “Reconstruction of plasma equilibrium and eddy current induced error fields in the Lithium Tokamak eXperiment-b” for the period July 1, 2018 – July 2, 2023. This award supported research into computational models for plasma equilibrium and 3D conducting structures in the Lithium Tokamak eXperiment-b (LTX-b) using the PSI-Tri/TokaMaker and PSI-Tet codes. Additional detail on progress in support of these tasks is provided below.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modelling the complex evaporated gas flow and its impact on particle spattering during laser powder bed fusion

The additive manufacturing (AM) of metals is becoming an increasingly important production process with the potential to replace traditional techniques such as casting. Laser Powder Bed Fusion (LPBF) is used in many applications to print metal parts from powder. The metal powder is heated locally with sufficient laser radiation that the liquid melt easily reaches its boiling temperature, which leads to a metallic vapour jet that can entrain both powder bed particles and molten droplets. The small size of laser-matter interaction site makes a detailed experimental analysis of the process challenging. Synchrotron X-ray imaging experiments are one of the few methods which can capture the dynamic melting and solidification processes. Comparing such experiments with computer simulations of the process is an important approach in order to better understand the manufacturing process and to analyse the influence of process parameters on the evaporated gas jet and the subsequent impact on particle ejection, leading to potentially reduced AM component quality. The melting and solidification of the metal powder is simulated using an Eulerian multiphase approach based on a control volume discretization of powder bed and substrate and a volume of liquid separation from melt and gas phase. The gas phase modelled as an ideal gas reaches velocities up to 100 m/s. Lagrangian particle tracking in the simulation demonstrates that the velocity fields calculated by the Eulerian multi-phase approach in combination with a standard drag-force model lead to particle accelerations in good agreement with those measured experimentally. In order to avoid numerical laborious Lagrangian calculations, a direct method to compare an Eulerian multiphase simulation with synchrotron X-ray experiments was introduced and validated. This approach is used to analyse the influence of process parameters including laser power and laser speed on the maximal acceleration of particles from the melt pool area. While the particle acceleration increases linearly with line energy in the conduction mode, a linear decrease of the acceleration with increasing line energy can be found in the transition mode before the accel-eration increases again with line energy in the keyhole mode.

36 MATERIALS SCIENCE↗

KEBLM: Knowledge-Enhanced Biomedical Language Models

Pretrained language models (PLMs) have demonstrated strong performance on many natural language processing (NLP) tasks. Despite their great success, these PLMs are typically pretrained only on unstructured free texts without leveraging existing structured knowledge bases that are readily available for many domains, especially scientific domains. As a result, these PLMs may not achieve satisfactory performance on knowledge-intensive tasks such as biomedical NLP. Comprehending a complex biomedical document without domain-specific knowledge is challenging, even for humans. Inspired by this observation, we propose a general framework for incorporating various types of domain knowledge from multiple sources into biomedical PLMs. We encode domain knowledge using lightweight adapter modules, bottleneck feed-forward networks that are inserted into different locations of a backbone PLM. For each knowledge source of interest, we pretrain an adapter module to capture the knowledge in a self-supervised way. We design a wide range of self-supervised objectives to accommodate diverse types of knowledge, ranging from entity relations to description sentences. Once a set of pretrained adapters is available, we employ fusion layers to combine the knowledge encoded within these adapters for downstream tasks. Each fusion layer is a parameterized mixer of the available trained adapters that can identify and activate the most useful adapters for a given input. Our method diverges from prior work by including a knowledge consolidation phase, during which we teach the fusion layers to effectively combine knowledge from both the original PLM and newly-acquired external knowledge using a large collection of unannotated texts. After the consolidation phase, the complete knowledge-enhanced model can be fine-tuned for any downstream task of interest to achieve optimal performance. Extensive experiments on many biomedical NLP datasets show that our proposed framework consistently improves the performance of the underlying PLMs on various downstream tasks such as natural language inference, question answering, and entity linking. These results demonstrate the benefits of using multiple sources of external knowledge to enhance PLMs and the effectiveness of the framework for incorporating knowledge into PLMs. Finally, while primarily focused on the biomedical domain in this work, our framework is highly adaptable and can be easily applied to other domains, such as the bioenergy sector.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Neural network for 3D inertial confinement fusion shell reconstruction from single radiographs

In inertial confinement fusion (ICF), x-ray radiography is a critical diagnostic for measuring implosion dynamics, which contain rich three-dimensional (3D) information. Traditional methods for reconstructing 3D volumes from 2D radiographs, such as filtered backprojection, require radiographs from at least two different angles or lines of sight (LOS). In ICF experiments, the space for diagnostics is limited, and cameras that can operate on fast timescales are expensive to implement, limiting the number of projections that can be acquired. To improve the imaging quality as a result of this limitation, convolutional neural networks (CNNs) have recently been shown to be capable of producing 3D models from visible light images or medical x-ray images rendered by volumetric computed tomography. Here, we propose a CNN to reconstruct 3D ICF spherical shells from single radiographs. We also examine the sensitivity of the 3D reconstruction to different illumination models using preprocessing techniques such as pseudo-flatfielding. To resolve the issue of the lack of 3D supervision, we show that training the CNN utilizing synthetic radiographs produced by known simulation methods allows for reconstruction of experimental data as long as the experimental data are similar to the synthetic data. We also show that the CNN allows for 3D reconstruction of shells that possess low mode asymmetries. Further comparisons of the 3D reconstructions with direct multiple LOS measurements are justified.

3D Reconstruction↗

Dynamic response of 17-4 stainless steel as a function of manufacturing method and heat treatment

We present a series of plate-impact experiments on 17-4 stainless steel to study the effect of manufacturing method and heat treatment on the Hugoniot elastic limit (HEL), Hugoniot, phase transformation stress, and spallation strength. Two traditional manufacturing methods were considered, wrought processing and casting, as well as two additive manufacturing methods, laser powder-bed fusion (LPBF) and wire-fed electron beam (EBAM). For both LPBF and EBAM 17-4 stainless steel variants, two billets were printed, enabling the application of two unique heat treatments. The HEL stress depended heavily on the thermal history, with the HEL increasing after the formation of Cu-rich precipitates via heat treatment. The Hugoniot response both below and above the phase transition was unaffected by the manufacturing method or heat treatment. The phase transition stress depended heavily on the thermal history, with its variation being attributed to the presence of various microstructural features. This is supported by a marked increase in the phase transition stress after precipitation hardening. These results suggest that the notion of the phase transition stress being dictated by bulk composition is an oversimplification and the stress fields generated by the meso-scale structure are a dominant force. The spallation strength was lower in the cast material compared to all other 17-4 stainless steel variants due to the presence of brittle δ-ferrite inclusions. Additionally, a drop in the tensile strain-rate was observed in the spallation response above the phase transition stress, which was hypothesized to stem from the kinetics of the reversion to the low-pressure phase during spall.

Compressive stress↗

Roadmap for warm dense matter physics

This roadmap presents the state-of-the-art, current challenges and near future developments anticipated in the thriving field of warm dense matter physics. Originating from strongly coupled plasma physics, high pressure physics and high energy density science, the warm dense matter physics community has recently taken a giant leap forward. This is due to spectacular developments in laser technology, diagnostic capabilities, and computer simulation techniques. Only in the last decade has it become possible to perform accurate enough simulations \& experiments to truly verify theoretical results as well as to reliably design experiments based on predictions. Consequently, this roadmap discusses recent developments and contemporary challenges that are faced by theoretical methods, and experimental techniques needed to create and diagnose warm dense matter. A large part of this roadmap is dedicated to specific warm dense matter systems and applications in astrophysics, inertial confinement fusion and novel material synthesis.

FOS: Physical sciences↗

Comprehensive new insights on the potential use of SiC as plasma-facing materials in future fusion reactors

Abstract The performance of silicon carbide as an alternative plasma facing material (PFM) was studied at various irradiation conditions relevant to ion energies and fluxes of a fusion reactor. This analysis involves detailed modeling of subsurface plasma/material interactions, sputtered particle transport above the surface and redeposition, and related changes in material composition and microstructure induced by steady-state and Edge Localized Mode ion fluxes. Transition of a crystalline SiC surface to semi-crystalline and amorphous phases was analyzed based on advanced modeling of DIII-D tokamak experiments where SiC was irradiated in single- and multiple- L-mode and H-mode discharges. This analysis shows that displacement damage, particle deposition/redeposition, and D accumulation on the SiC divertor surface can lead to significant microstructural changes that result in enhanced sputtering erosion in comparison with the original crystalline material. However, the resulting total net erosion rate for a full-coverage, advanced tokamak, SiC coated divertor may well be acceptably low. Moreover, the C sputtering yield from the evolved SiC surface can be seven times lower than from a pure graphite surface; this would imply significantly reduced tritium co-deposition rates in a D-T tokamak reactor, compared with a pure carbon surface. It was also determined that chemical sputtering of both C and Si should not result in any noticeable effect on the net erosion, for attached plasma regimes. Our results thus show encouraging results overall for use of SiC as a PFM in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of surface roughness on phase transition timing in megaampere pulsed-power–driven exploding conductors

An understanding of material phase transitions in megaampere pulsed-power–driven exploding conductors is important for predicting the growth of hydrodynamic instabilities in magneto-inertial fusion concepts. Here, this study analyzes phase transitions in electrical conductor explosions using 1D Lagrangian and 2D arbitrary Lagrangian–Eulerian resistive magnetohydrodynamic simulations to show that micrometer-scale surface roughness can lead to the electrothermal instability (ETI), a feedback effect that concentrates resistive heating and leads to early melting and ablation. Simulations of the Mykonos electrothermal instability II (METI-II) experiment show melting begins 19% sooner for machined rods with micrometer-scale surface roughness than for rods without these features. The surface magnetic field is 41 T around the initial region of melt, representing a lower magnitude than both the 86 T from 1D simulations and the 85 T threshold reported elsewhere. In 2D simulations with micrometer-scale surface roughness, temperature measurements indicate the critical point temperature of aluminum is reached 17% faster in comparison with 1D simulations. Values from 2D simulations with surface roughness align with predictions from ETI theory, and the observed temperature redistribution further supports the ETI as an underlying mechanism. Simulation results are validated against experimental photonic Doppler velocimetry data. This study shows 1D simulations are adequate to model conductors with sub-micrometer-scale surface roughness in this high-energy-density regime; however, 2D or 3D simulations are required to capture the full range of physics for accurately describing phase transitions in conductors with micrometer-scale or larger surface roughness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A heavy-ion production channel of 149 Tb via 63 Cu bombardment of 89 Y

We report the radionuclide 149 Tb (t 1/2 = 4.1 h) is a potential theranostic isotope which can simultaneously be used for targeted-alpha-particle therapy and positron-emission tomography. Feasibility experiments were performed to test a near-symmetric heavy-ion reaction of 63 Cu bombardment on monoisotopic 89 Y. The indirect reaction was studied to avoid isomer production. Offline gamma spectroscopy was used to quantify thick-target physical yields and experimental results show modest agreement to the fusion-evaporation code PACE4. A near-symmetric fission yield was also observed.

149Tb↗

An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks

We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.

97 MATHEMATICS AND COMPUTING↗

A theoretical approach for transient shock strengthening in high-energy-density laser compression experiments

In high-energy-density shock compression experiments, the desired state of compression is typically achieved by passing shock waves through a sequence of different materials. Here, a theoretical approach for transiently strengthening such shocks passing through interfaces of experimentally relevant impedance ratios is examined. A semi-analytical method based on characteristics analysis is used to solve the problem of a shock passing from one material to another through an intermediate region of non-uniform impedance; by appropriately designing this region, a greater shock strength can be achieved in the second material for a finite duration than in the absence of this region. When a shock passes into a material of higher impedance, the shock strengthening increases to a maximum before decreasing to an asymptotic value as the strength of the incident shock is increased. For shocks passing into materials of decreasing impedance, the shock strengthening increases monotonically as the strength of the incident shock is increased when the impedance ratio is above a critical threshold. Incorporating multiple intermediate materials can further increase the strength of the transmitted shock, with an exponential discretization of intermediate material impedances being the most effective distribution for strengthening strong shocks. The results suggest that up to 25% and 9% increases in pressure behind the leading shock can be achieved for materials of increasing and decreasing impedance, respectively. The technique is applied to the design of laser-driven dynamic compression experiments, and the results of the analysis are verified via comparison to simulations performed with the HYADES hydrodynamics code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of NBI-induced energetic ions on internal kink stability in HL-2M

A systematic investigation is carried out, studying the effect of the neutral beam injection induced energetic particles (EPs) on the n = 1 (n is the toroidal mode number) internal kink (IK) instability in the HL-2M tokamak, utilizing the MARS-F/K code [Liu et al., Phys. Plasmas 7, 3681 (2000) and 15, 112503 (2008)]. A high-beta sawteething HL-2M scenario, simulated by the TRANSP code [Breslau et al. Computer Software (2018)], is chosen for this study. Compared to the fluid model, non-perturbative magnetohydrodynamic (MHD)-kinetic hybrid computations with MARS-K show a generally stabilization effect on the IK, due to drift kinetic resonances associated with EPs. Further, the bounce resonance of trapped EPs has minor influence on the mode stability. In the absence of the plasma equilibrium flow and with the assumed particle pitch distribution, the transit resonance of co-current (countercurrent) passing EPs destabilizes (stabilizes) the IK. With plasma flow, both co- and countercurrent passing EPs tend to stabilize the mode, but the effect is stronger with the countercurrent particles. These modeling results provide useful guidance for interpreting MHD instabilities in the future high-performance experiments in HL-2M.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neural network based fast prediction of β N limits in HL-2M

Artificial neural networks (NNs) are trained, based on the numerical database, to predict the no-wall and ideal-wall β N limits, due to onset of the n = 1 ( n is the toroidal mode number) ideal external kink instability, for the HL-2M tokamak. The database is constructed by toroidal computations utilizing both the equilibrium code CHEASE (Lütjens et al 1992 Comput. Phys. Commun. 69 287) and the stability code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681). The stability results show that (1) the plasma elongation generally enhances both β N limits, for either positive or negative triangularity plasmas; (2) the effect is more pronounced for positive triangularity plasmas; (3) the computed no-wall β N limit linearly scales with the plasma internal inductance, with the proportionality coefficient ranging between 1 and 5 for HL-2M; (4) the no-wall limit substantially decreases with increasing pressure peaking factor. Furthermore, both the NN model and the convolutional neural network (CNN) model are trained and tested, producing consistent results. The trained NNs predict both the no-wall and ideal-wall limits with as high as 95% accuracy, compared to those directly computed by the stability code. Additional test cases, produced by the Tokamak Simulation Code (Jardin et al 1993 Nucl. Fusion 33 371), also show reasonable performance of the trained NNs, with the relative error being within 10%. The constructed database provides effective references for the future HL-2M operations. The trained NNs can be used as a real-time monitor for disruption prevention in the HL-2M experiments, or serve as part of the integrated modeling tools for ideal kink stability analysis.

Physics↗

Nonlinear modeling of the scaling law for the m/n = 3/2 error field penetration threshold

The scaling law for the m/n=3/2 error field (EF) penetration threshold is predicted numerically based on nonlinear single-fluid and two-fluid modeling using the TM1 code. The simulated penetration threshold of radial magnetic field b r at the plasma edge is scaled to the electron density n e , temperature T e , viscous time τ μ , toroidal field B t and the natural frequency ω in the form of b r /B t ∝n e αn T e αT τ μ αμ B t αB ω αω by scanning these parameters separately. Here, α n , α T , α μ , α B and α ω are the scaling coefficients on n e , T e , τ μ , B t and ω, respectively. Single-fluid modeling shows that the 3/2 EF threshold scales as b r /B t ∝n e 0.56 T e 0.6 τ μ -0.59 B t -1.15 ω, which is similar with the analytical scaling law in both the Rutherford and visco-resistive regimes. Yet, two-fluid modeling shows that the scaling law differs significantly in particular regarding the dependence on plasma rotation. In detail, the scaling coefficient α n on density decreases from 0.67 to 0.56 and α T on temperature decreases from 0.67 to 0.32, while α μ on viscous time is around -0.45 and α B on toroidal field decreases slightly from -1.15 to -1, when the ratio |ω E /ω *e | between plasma rotation frequency ω E and diamagnetic drift frequency ω *e varies from 0 to 10. Scans of the plasma rotation reveals that the penetration threshold linearly depends on the perpendicular rotation frequency (or natural frequency) ω ⊥e =ω E +ω *e , and there is a minimum in the required field amplitude when ω ⊥e 0. In addition, the enduring mystery of non-zero penetration threshold at zero plasma natural frequency in EF experiments is resolved by two-fluid simulations. We report that the very small island and smooth bifurcation in EF penetration near zero frequency is hard to detect in the experiment, leading to a finite penetration threshold within the capability of the experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neutral leakage, power dissipation and pedestal fueling in open vs closed divertors

The effect of divertor closure on neutral leakage, power dissipation and pedestal fueling is explored using the SOLPS code for different densities including attached and detached divertor conditions for DIII-D open and closed divertors. Here, a metric for divertor closure in terms of the percentage of neutrals escaping the divertor is presented. At the detachment onset the population of neutrals escaping the divertor is 12% for the closed divertor and 35% for the open divertor. This results in a lower contribution of the convective term to the total heat flux and increased power losses in the closed divertor. The comparison between the two geometries reveals that for the same upstream density, the upstream temperature is 23% higher in the closed divertor. The physical mechanisms for detachment at lower upstream density in the closed divertor are also identified with the reduction of the convective term directly resulting in 40% decrease and the increased power losses contributing with a factor of 2. The baffling also lowers the neutral flux refueling the pedestal resulting in further reduction of the core density for achieving detachment. The results introduced here show that the divertor recycling is the primary fueling mechanism in DIII-D. It is shown that fueling occurring in the region between the X-point and the outer midplane is at least four times lower in the closed divertor due the baffling structure which enhances target neutral trapping and consequently reduced upstream leakage. Due to the openness of the inner leg to the common flux and HFS regions both open and closed divertor experience almost identical fueling at the HFS highlighting the need to also consider the inboard divertor for shape optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of resonant magnetic perturbations on zonal flows and microturbulence

Results from a joint experimental and computational effort studying the effect of resonant magnetic perturbations (RMPs) on microturbulence levels and their connection to zonal flows in the DIII-D tokamak L-mode are presented. Beam emission spectroscopy measurements show a direct increase in density fluctuations at microturbulent scales with increasing RMP amplitude, suggesting that magnetic activity introduced by the RMP affects the regulation of microturbulence on DIII-D. This is analogous to how MHD-scale magnetic fluctuations arising from tearing modes have been observed in simulations to increase microturbulence levels in the reversed-field pinch (RFP). In the RFP, this is attributed to magnetic fluctuations eroding turbulence-limiting zonal flows; this work examines if a similar mechanism is present for DIII-D microturbulence. Gyrokinetic simulations find that the application of an RMP corresponds directly to a decrease in zonal flow levels, producing a similar increase of turbulent fluctuation levels over a range of RMP amplitudes as observed in the experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High Yield Muon Catalyzed Fusion

Fusion holds great promise as a clean and abundant energy source. However, traditional thermonuclear fusion encounters significant challenges due to the extreme temperatures required to overcome the coulomb barrier for two nuclei to fuse. In contrast, muon-catalyzed fusion presents an alternative approach that can surmount this barrier at significantly lower temperatures. Muons, with properties resembling those of electrons but 200 times heavier, can effectively reduce the atomic orbital radius, enabling central nuclei to overcome the coulomb force through the strong force. By introducing muons into a mixture of deuterium and tritium (two hydrogen isotopes), fusion is facilitated, releasing a 3.5MeV alpha particle and a 14.1MeV neutron. In the majority of cases, the muon is liberated and can initiate further fusions. However, approximately 0.8% of the time, it adheres to the alpha particle and remains bound until it either decays or undergoes reactivation through collisional ionization. To maximize the number of fusions per muon, it is crucial to enhance the cycling rate and reactivation fraction. Theoretical predictions and experimental data both suggest that the sticking rate decreases with increasing density. However, there exists a discrepancy between experimental observations and theoretical estimations regarding the extent of this decrease. To address these disparities, this experiment aims to investigate the cycling rate and sticking fraction under higher temperatures and pressures than previously explored.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗