Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “average atom”

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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Simulations of plasma current induced by toroidal field ramping down on tokamaks

In tokamak experiments, the BT ramping down can drive off-axis and parallel inductive current. This approach leads to a decrease in q95 and an increase in normalized beta, βN. The off-axis inductive current also broadens the current profile. Typically, the 1.5D transport code is used to simulate the time evolution of the plasma current profile, which is based on the flux-surface-averaged Faraday’s law. The ONETWO code is one of those transport codes. However, this code cannot simulate the situation of the evolving toroidal field BT. This article proposes a new time-dependent model to take into account the BT ramping down situation. A modified formula of flux-surface-averaged Faraday’s law was derived to consider the effect of the BT ramping rate on the current evolution, and it was implemented in the ONETWO code. Then, the modified ONETWO code was used to simulate the BT ramping down experiment on DIII-D. The simulation result of the plasma current evolution with BT ramping down shows a broader current profile with smaller ohmic current induced by the poloidal field, compared to that without BT ramping down.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Recrystallization, cracking, and erosion of dispersoid-strengthened tungsten materials during exposure to divertor plasmas

In this study, we investigated the effects of combined intense particle and heat flux exposure on advanced tungsten plasma-facing materials within the DIII-D fusion facility. Our test matrix included two types of dispersoid-strengthened tungsten (containing either 100 nm diameter TiO 2 or Ni particles), along with high-purity polycrystalline tungsten as a reference. This experiment relied on a sample geometry angled at 15° relative to the divertor surface, thereby allowing the surfaces to intercept steady-state perpendicular heat fluxes (q ⟂ ) ranging from 10.1 to 19.6 MW/m 2 . During each shot, the samples were exposed to 42 Hz edge-localized modes (ELMs), allowing us to test the material response to transient heating. We correlated the exposure conditions with extensive post-test surface composition analysis and microscopy to determine how the plasma modified each surface. The angled specimens closest to the strike point received the highest combined heat and particle flux and melted midway through the experiment. EBSD analysis revealed they were completely recrystallized throughout, with an average grain size >100 µm. On the other hand, the specimens that received a lower steady state heat flux survived with more superficial surface damage. Whereas the high-purity polycrystalline tungsten exhibited a higher surface roughness, the dispersoid-strengthened material exhibited more extensive shallow inter-granular cracking. In addition, the surface was depleted of dispersoids following plasma exposure, possibly because of evaporation and/or sputtering. The results described here provide insights into the performance of these materials in a fusion environment which can guide further optimization for use in long-pulse devices.

Kolasinski, Robert D. [Sandia National Laboratorie↗
NASA NTRS1993

Characterization of Magnetoacoustic Emission Related to Structural Properties of Ferromagnets

An extensive study of magnetoacoustic emission (MAE) properties has been performed over the past several years. As a result, the dependence of the spectral characteristics of MAE on certain microstructural variations and uniaxially applied stress in a particular type of low carbon steel are now well known. The embrittlement-causing concentration of certain atomic species, e. g., tin, sulphur, phosphorous etc., at the grain boundaries of this steel creates strong potential barriers resisting the motion of non-180 domain walls which is the source of MAE. (Since the only type of non-180 domain walls in this material are 90 domain walls, the term 90 domain wall will be used throughout this paper in place of non-180 domain wall.) An MAE burst produced during one-half cycle of a hysteresis loop at a low AC magnetic field frequency (e. g., 0.7 Hz) shows two sub-peaks; the leading peak is usually sharp and short-lived, while the trailing peak is usually smooth and long-lasting. It has been shown that the enhanced domain wall-defect interaction, due to the strengthened potential barriers, causes an increase in the asymmetry of the MAE signal by suppressing the leading sub-peak and amplifying the trailing sub-peak. This phenomena is due to the delayed motion of the 90 domain walls. The effect of a tensile stress applied parallel to the external AC magnetic field is to diminish the MAE. On the other hand, the amplitude of the MAE burst has been shown to be a non-monotonic function of the stress amplitude. Recently, our study has concentrated on obtaining quantitative values for parameters computed from the MAE spectra averaged over a sufficient number of cycles to achieve statistical stability. Nevertheless, certain fundamental elements of the MAE characteristics remain unexplained.

Namkung, M.↗
DOE OSTI2022

Evaluation of the GaAs displacement damage metric using updated nuclear data

The emerging use of the physics-based athermal recombination-corrected displacement per atom (arc-dpa) model for the displacement damage efficiency has motivated a re-evaluation of the historical empirically-derived GaAs damage response function with the purpose of highlighting needs for future analytical and experimental work. The 1-MeV neutron damage equivalence methodology used in the ASTM E-722 standard for GaAs has been re-evaluated using updated nuclear data. This yielded a higher fidelity representation of the GaAs displacement kerma and, through the use of the refined PKA recoil energy-dependent damage efficiency model, an updated 1-MeV(GaAs) displacement damage function. This re-evaluation included use of the Norgett-Robinson-Torrens (NRT) model for an updated threshold treatment, rather than the sharp-threshold Kinchin-Pease model used in the current ASTM standard. The underlying nuclear data evaluations have been updated to use the ENDF/VIII.0 {sup 75}As and TENDL-2019 {sup 71}Ga/{sup 69}Ga evaluations. The displacement kerma and 1-MeV-equivalent damage responses were calculated using a modified NJOY-2016 code which allowed for refinements in some of the damage models. This paper shows that an updated displacement damage function, based upon the latest nuclear data, is consistent with the experimental data used to develop the current ASTM E-722 GaAs standard. Using a double ratio approach to compare the available experimental data with the calculated response, the average legacy double ratio was found to be 0.97 ± 0.05 and the average updated double ratio was found to be 0.94 ± 0.05. (authors)

36 MATERIALS SCIENCE↗
DOE OSTI2026

Achievement of a high-density, high-confinement, and high-beta tokamak plasma regime in DIII-D, and implications for a lower-current path for ITER and FPP

Experiments on DIII-D have demonstrated a density-confinement synergy that enables sustainment of high performance in a previously unattained parameter regime of simultaneous very high energy confinement quality (H 98y2 ≥ 1.5), very high line-average density Greenwald fraction (ƒ Gr = πa 2 < n >/I P ≥ 1.4), and high toroidal beta (β T ≥ 3%). Tokamak operation in this regime is essential for a compact steady-state FPP, as well as for Q=10 with 500 MW of fusion power in ITER at I P << 15 MA. These experiments leveraged the knowledge that, in the high-poloidal-beta (β P ) regime, impurity and density gradients can enhance turbulence stabilization caused by high α MHD (α MHD ~(dβ P )⁄dr). This was described by theoretical predictions and gyrokinetic transport simulations [M.T. Kotschenreuther et al, 2024 Nucl. Fusion, 64 076033], and later confirmed by experiments on DIII-D [S. Ding et al, 2024 Nature 629 555]. To increase both β P and β T , the new experiments increased the ideal-wall stability β N -limit by using a smaller plasma-outer wall distance and higher triangularity in the plasma cross section (top/bottom average δ~0.9), enabled by the recent “shape & volume rise” (SVR) modification to the DIII-D divertor. The higher triangularity also contributed to achieving higher ƒ Gr by enabling higher pedestal density. At high density, the pedestal is ballooning limited and exhibits small and frequent ELMs, while the divertor is near detachment even without any impurity seeding. High plasma performance was attained and sustained reproducibly, with the eventual terminations brought about by an MHD mode destabilized as the current profile slowly continued to evolve. A path to stationary fully noninductive operation might include ECH injection to reduce both core impurity accumulation and the electron collisionality, thus increasing the bootstrap current. These experiments provide the first experimental demonstration of the ƒ Gr , H 98y2 , and β T values required simultaneously for ITER Q = 10 at I P < 10 MA, pointing to practical ways to improve the energy confinement in a fusion reactor.

Garofalo, Andrea M. [General Atomics, San Diego, C↗
DOE OSTI2020

The MFiX Particle-in-Cell Method (MFiX-PIC) Theory Guide

MFiX (Multiphase Flow with Interphase eXchanges) is an open-source multiphase flow solver developed at the National Energy Technology Laboratory. Within the code, users have access to a single phase or interpenetrating continua-based multiphase two-fluid model (TFM), a discrete element model (DEM), and a particle-in-cell model (PIC). TFM, DEM, and PIC can all be used to create multiphase simulations that include hydrodynamics, chemical reactions, and heat transfer.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Assessing high fidelity multi-component models to facilitate safeguards at Gas Centrifuge Enrichment Plants

We report that the International Atomic Energy Agency (IAEA) inspectors routinely carry out environmental sampling (ES) as a verification method. Collection of environmental swipe samples at various locations in Gas Centrifuge Enrichment Facilities (GCEPs) is an important process in detecting misuse of a declared facility and possibly the existence of undeclared nuclear material. These samples are measured for isotopic composition in uranium containing particles by Thermal Ionization Mass Spectrometry (TIMS) or Inductively Coupled Plasma Mass Spectrometry (ICP-MS). Even though ES is highly effective in detecting the absolute value of enrichments and their deviations from the declared values, it cannot explain the cause of those changes. Several potential explanations can serve as possibilities for particles detected above or other than the declared enrichment. These include normal and non-malicious events such as the design of the enrichment cascades, unintentional failures of the machines, or overshoot during startup of a cascade. It can also be the result of deliberate misuse by the facility operators. The primary objective of this work is to understand how these factors affect the enrichments produced by a cascade and quantify anticipated multi-isotopic concentrations for each case. The following methodology is employed to determine signatures at a particular facility. 1) Utilize a new two-dimensional multi-component diffusion code to obtain centrifuge performance data and use that information to design and perform cascade analysis. Compare and contrast the results with previous 1-D radially averaged solutions from the Pancake code. 2) Design a GCEP cascade with the production goal of 19.75% 235 U for each set of machine data above. Investigate two cascade scenarios that include enrichment of natural uranium (NU) feed to 19.75% 235U in a single cascade compared to a two-step process of NU to 5% and then 5% to 19.75%. 3) Simulate the intentional vs. unintentional off-normal scenarios in the cascades to assess the differences in isotopic concentrations. A non-ideal squared-off cascade model developed at the University of Virginia is used to calculate flow rates and isotopic concentrations of the process gas. The analysis is performed using the Rome machine model operated at 600 m/s rotor speed. The upper and lower bounds of normal and abnormal enrichments in a typical facility are used in conjunction with ES results to understand the root causes of such observations.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2022

Elucidation of the physical factors that control activated transport of penetrants in chemically complex glass-forming liquids

Understanding the activated transport of penetrant or tracer atoms and molecules in condensed phases is a challenging problem in chemistry, materials science, physics, and biophysics. Many angstrom- and nanometer-scale features enter due to the highly variable shape, size, interaction, and conformational flexibility of the penetrant and matrix species, leading to a dramatic diversity of penetrant dynamics. Based on a minimalist model of a spherical penetrant in equilibrated dense matrices of hard spheres, a recent microscopic theory that relates hopping transport to local structure has predicted a novel correlation between penetrant diffusivity and the matrix thermodynamic dimensionless compressibility, S 0 ( T ) (which also quantifies the amplitude of long wavelength density fluctuations), as a consequence of a fundamental statistical mechanical relationship between structure and thermodynamics. Moreover, the penetrant activation barrier is predicted to have a factorized/multiplicative form, scaling as the product of an inverse power law of S 0 ( T ) and a linear/logarithmic function of the penetrant-to-matrix size ratio. This implies an enormous reduction in chemical complexity that is verified based solely on experimental data for diverse classes of chemically complex penetrants dissolved in molecular and polymeric liquids over a wide range of temperatures down to the kinetic glass transition. The predicted corollary that the penetrant diffusion constant decreases exponentially with inverse temperature raised to an exponent determined solely by how S 0 ( T ) decreases with cooling is also verified experimentally. Our findings are relevant to fundamental questions in glassy dynamics, self-averaging of angstrom-scale chemical features, and applications such as membrane separations, barrier coatings, drug delivery, and self-healing.

Science & Technology - Other Topics↗
DOE OSTI2026

Divertor detachment characterization in negative triangularity discharges in DIII-D via 2D edge-plasma transport modeling

Edge fluid modeling of the first divertor-plasma detachment experiments in negative triangularity (NT) discharges in the DIII-D tokamak is presented using the 2D multifluid edge transport code UEDGE, including cross-field particle drifts. Experiments were performed where the lower single-null magnetic equilibrium had a strong NT (δ≈−0.5), that is, where the magnetic X-point is at a larger major radius than the core magnetic axis. Divertor-plasma detachment was induced by increasing the core plasma density in DIII-D via intrinsic gas puffing. Here density scans are performed with UEDGE to reach a detached plasma and to quantitatively recover the experimental rollover of the ion saturation current on the outer divertor target plate. The simulations cover experiments with both signs of the toroidal magnetic field, B T , where the ion magnetic Grad-B drifts are directed into (forward B T ) and out of (reverse B T ) the divertor region. Consistent with experiments with neutral beam power injection, the NT simulations reproduce: 40% higher density is needed to reach detachment onset with forward B T compared with reverse B T , and the absence of deep detachment is found with reverse B T . Similarly, comparison between Ohmic discharges in NT and positive triangularity (PT) shaping confirms that a substantially higher density is needed to achieve detachment in NT than in PT, with NT requiring an line-average density of at least the Greenwald density or higher. Simulation results suggest that higher densities are needed to reach detachment in negative compared to PT because these discharges have a shorter midplane-to-target distance along the total magnetic field B, a shorter outer divertor poloidal leg length (0.06 m vs 0.2 m), and reduced radial transport [near-scrape-off layer (SOL) D ⊥ /χ ⊥ =0.3/0.5 vs D ⊥ /χ ⊥ =1.0/1.0, all in (m 2 /s)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1999

Interstellar PAH emission in the 11-14 micron region: new insights from laboratory data and a tracer of ionized PAHs

The Ames infrared spectral database of isolated, neutral and ionized polycyclic aromatic hydrocarbons (PAHS) shows that aromatic CH out-of-plane bending frequencies are significantly shifted upon ionization. For solo- and duet-CH groups, the shift is pronounced and consistently toward higher frequencies. The solo-CH modes are blueshifted by an average of 27 cm-1 and the duet-CH modes by an average of 17 cm-1. For trio- and quartet-CH groups, the ionization shifts of the out-of-plane modes are more erratic and typically more modest. As a result of these ionization shifts, the solo-CH out-of-plane modes move out of the region classically associated with these vibrations in neutral PAHS, falling instead at frequencies well above those normally attributed to out-of-plane bending, vibrations of any type. In addition, for the compact PAHs studied, the duet-CH out-of-plane modes are shifted into the frequency range traditionally associated with the solo-CH modes. These results refine our understanding of the origin of the dominant interstellar infrared emission feature near 11.2 microns, whose envelope has traditionally been attributed only to the out-of-plane bending of solo-CH groups on PAHS, and provide a natural explanation for the puzzling emission feature near 11.0 microns within the framework of the PAH model. Specifically, the prevalent but variable long-wavelength wing or shoulder that is often observed near 11.4 microns likely reflects the contributions of duet-CH units in PAH cations. Also, these results indicate that the emission between 926 and 904 cm-1 (10.8 and 11.1 microns) observed in many sources can be unambiguously attributed to the out-of-plane wagging, of solo-CH units in moderately sized (fewer than 50 carbon atom) PAH cations, making this emission an unequivocal tracer of ionized interstellar PAHS.

Astronomy↗
DOE OSTI2025

Re-evaluating the prompt fission neutron spectrum of spontaneously fissioning 252 Cf

The prompt fission neutron spectrum (PFNS) of spontaneously fissioning 252 Cf is a Neutron Data Standards observable. Nearly all fission spectra of actinides were measured relative to it, using efficiencies derived from it, or analyzed with simulations validated by it. The current Standards evaluation was published by W. Mannhart in 1987. It could not be updated because the evaluation input, experimental mean values and covariances, were lost. First, we attempt to reproduce it. However, Mannhart’s evaluation can only be reproduced within its one-σ uncertainties as some of its aspects (e.g., experimental covariances, rejected data points) remain unknown. Therefore, a new evaluation is presented: We revisit all existing experimental 252 Cf(sf) PFNS data, including those published after the release of the current Standards evaluation, and re-estimate associated covariances. The newly evaluated 252 Cf(sf) PFNS differs distinctly from Mannhart’s below 300 keV and extends it to lower and higher outgoing neutron energies (500 eV–25 MeV). The new evaluated uncertainties are larger from 3–9 MeV and smaller otherwise. Spectrum averaged cross sections of importance to the International Reactor Dosimetry and Fusion File community calculated with the new spectrum are close to those calculated with Mannhart’s evaluation and agree with experimental values well within their uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

LANL's Contribution to the Design and Preparation of the ANL Bubble Experiment #2

In collaboration with Argonne National Laboratory (ANL), Los Alamos National Laboratory (LANL) is assisting in the design and development of portions of the second ANL Bubble Experiment to be performed in late 2020 at ANL. The ANL Bubble Experiment, as called in this report, is a series of direct electron irradiations of a uranyl sulfate solution to produce radiolysis-induced gas bubbles of hydrogen and oxygen. The gas bubbles formed in the solution enhance mixing and heat transfer. The study of the dynamics, shape, and size of radiolysis-induced gas bubbles is of great importance to understand and characterize the thermal and fluid behavior of the solution, especially for the solution based, neutron-induced fission production technique for Mo-99. During the first experiment performed in 2014, ANL’s 35 MeV electron linear accelerator provided average powers of 6, 12, and 15 kW using a rastered beam to homogeneously heat the 15 x 15 x 80-cm uranyl-sulfate solution. Gas bubble size, shape, velocity, solution temperature, and hydrogen and oxygen concentrations were recorded during the irradiations. The details of the experimental setup and results for the first experiment are described in the ANL reports, Design and Construction of Experiment for Direct Electron Irradiation of a Uranyl Sulfate Solution: Bubble Formation and Thermal Hydraulics Studies and Experimental Results for Direct Electron Irradiation of a Uranyl Sulfate Solution: Bubble Formation and Thermal Hydraulics Studies. The purpose of this report is to describe the issues and lessons learned associated with the first experiment and elaborate on designs to improve and obtain more accurate experimental results for the upcoming experiment.

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

OrbNet Denali: A machine learning potential for biological and organic chemistry with semi-empirical cost and DFT accuracy

We present OrbNet Denali, a machine learning model for an electronic structure that is designed as a drop-in replacement for groundstate density functional theory (DFT) energy calculations. The model is a message-passing graph neural network that uses symmetryadapted atomic orbital features from a low-cost quantum calculation to predict the energy of a molecule. OrbNet Denali is trained on a vast dataset of 2.3 × 10 6 DFT calculations on molecules and geometries. This dataset covers the most common elements in biochemistry and organic chemistry (H, Li, B, C, N, O, F, Na, Mg, Si, P, S, Cl, K, Ca, Br, and I) and charged molecules. OrbNet Denali is demonstrated on several well-established benchmark datasets, and we find that it provides accuracy that is on par with modern DFT methods while offering a speedup of up to three orders of magnitude. For the GMTKN55 benchmark set, OrbNet Denali achieves WTMAD-1 and WTMAD-2 scores of 7.19 and 9.84, on par with modern DFT functionals. For several GMTKN55 subsets, which contain chemical problems that are not present in the training set, OrbNet Denali produces a mean absolute error comparable to those of DFT methods. For the Hutchison conformer benchmark set, OrbNet Denali has a median correlation coefficient of R 2 = 0.90 compared to the reference DLPNO-CCSD(T) calculation and R 2 = 0.97 compared to the method used to generate the training data (ωB97X-D3/def2- TZVP), exceeding the performance of any other method with a similar cost. Similarly, the model reaches chemical accuracy for noncovalent interactions in the S66x10 dataset. For torsional profiles, OrbNet Denali reproduces the torsion profiles of ωB97X-D3/def2-TZVP with an average mean absolute error of 0.12 kcal/mol for the potential energy surfaces of the diverse fragments in the TorsionNet500 dataset.

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

Mid-infrared trace detection with parts-per-quadrillion quantitation accuracy: Expanding frontiers of radiocarbon sensing

Detection sensitivity is a critical characteristic to consider during selection of spectroscopic techniques. However, high sensitivity alone is insufficient for spectroscopic measurements in spectrally congested regions. Two-color cavity ringdown spectroscopy (2C-CRDS), based on intra-cavity pump–probe detection, simultaneously achieves high detection sensitivity and selectivity. This combination enables mid-infrared detection of radiocarbon dioxide ( 14 CO 2 ) molecules in room-temperature CO 2 samples, with 1.4 parts-per-quadrillion (ppq, 10 − 15 ) sensitivity (average measurement precision) and 4.6-ppq quantitation accuracy (average calibrated measurement error for 21 samples from four separate trials) demonstrated on samples with 14 C/C up to ∼ 1.5 × natural abundance ( ∼ 1,800 ppq). These highly reproducible measurements, which are the most sensitive and quantitatively accurate in the mid-infrared, are accomplished despite the presence of orders-of-magnitude stronger, one-photon signals from other CO 2 isotopologues. This is a major achievement in laser spectroscopy. A room-temperature-operated, compact, and low-cost 2C-CRDS sensor for 14 CO 2 benefits a wide range of scientific fields that utilize 14 C for dating and isotope tracing, most notably atmospheric 14 CO 2 monitoring to track CO 2 emissions from fossil fuels. The 2C-CRDS technique significantly enhances the general utility of high-resolution mid-infrared detection for analytical measurements and fundamental chemical dynamics studies.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Challenging excited states from adaptive quantum eigensolvers: subspace expansions vs. state-averaged strategies

The prediction of electronic structure for strongly correlated molecules represents a promising application for near-term quantum computers. Significant attention has been paid to ground state wavefunctions, but excited states of molecules are relatively unexplored. In this work, we consider the adaptive, problem-tailored (ADAPT)-variational quantum eigensolver (VQE) algorithm, a single-reference approach for obtaining ground states, and its state-averaged generalization for computing multiple states at once. We demonstrate for both rectangular and linear H4, as well as for BeH2, that this approach, which we call multistate-objective, Ritz-eigenspectral (MORE)-ADAPT-VQE, can make better use of small excitation manifolds than an analogous method based on a single-reference ADAPT-VQE calculation, q-sc-EOM. In particular, MORE-ADAPT-VQE is able to accurately describe both avoided crossings and crossings between states of different symmetries. In addition to more accurate excited state energies, MORE-ADAPT-VQE can recover accurate transition dipole moments in situations where traditional ADAPT-VQE and q-sc-EOM struggle. These improvements suggest a promising direction toward the use of quantum computers for difficult excited state problems.

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

Preliminary Neutronics Design and Analysis of the Fast Modular Reactor

General Atomics is developing a new 100-MW(thermal) fast modular reactor (FMR) that provides safe, carbon-free electricity and is capable of incremental capacity additions. The modular design allows it to be factory built and assembled onsite to keep the capital cost low, while the use of dry cooling facilitates siting to complement renewables in nearly any location. The FMR uses high-assay low-enriched uranium-dioxide fuel encapsulated by recognized irradiationresistant silicon carbide composite (SiGA®) cladding that is derisked in the current accident-tolerant fuel program. The FMR fuel assembly is a hexagonal fuel bundle of 120 fuel rods. The total length of the fuel assembly is less than 4 m, with an active fuel length of 1.8 m. The fuel assemblies are configured in an annular core that is located and supported by the reactor internals. The coolant material is helium at a normal operating pressure of 7 MPa. The core is surrounded by zirconium silicide (Zr 3 Si 2 ) and graphite reflector blocks. The fuel, coolant, internals, and reflectors are contained within a reactor pressure vessel. Here, the preliminary nuclear design and analysis established the arrangement of the active core and reflector blocks. The nuclear design analyses of the FMR defined the design parameters, such as fuel enrichments, excess reactivity, fueling scheme, fuel cycle, power distribution, and control rod worth. The preliminary conceptual design determined the three-batch fueling scheme with the allowable total power peaking factor of 1.5. The average discharge burnup is 100 GW days per ton of uranium.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Creep-enhanced vortex pinning revealed through nonmonotonic relaxation of the Campbell length

Here, we study the effects of flux creep on the linear AC response of the vortex lattice in single crystals Ca 3 ⁢Ir 4 ⁢Sn 13 by measuring the Campbell penetration depth, 𝜆 C ⁡(𝑇,𝐻,𝑡) . Thermal fluctuations release vortices from shallow pinning sites, only for them to become re-trapped by deeper potential wells, causing an initial increase of the effective Labusch parameter, which is proportional to the pinning well curvature. This effect cannot be detected in conventional magnetic relaxation measurements but is revealed by our observation of a nonmonotonic time evolution of 𝜆 C ⁡(𝑇,𝐻,𝑡) , which directly probes the average curvature of the occupied pinning centers. The time evolution of 𝜆 C ⁡(𝑇,𝐻,𝑡) was measured at different temperatures in samples with different densities of pinning centers produced by electron irradiation. The curves can be collapsed together when plotted on a logarithmic time scale 𝑡 → 𝑇 ⁢ln ⁡(𝑡/𝑡 0 ) confirming that the time evolution is driven by flux creep. The 𝜆 C ⁡(𝑇,𝐻,𝑡) is hysteretic with a noticeable nonmonotonic relaxation in the presence of a vortex density gradient (after zero-field cooling), but is monotonic after field cooling, where the vortex density is uniform. This result quantitatively corroborates the novel picture of vortex creep based on the strong pinning theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Solution of the linear wave-particle kinetic equation for global modes of arbitrary frequency in a tokamak

The linear response of a plasma to perturbations of arbitrary frequency and wavelength is derived for any axisymmetric magnetized toroidal plasma. An explicit transformation to action-angle coordinates is achieved using orthogonal magnetic coordinates and the Littlejohn Lagrangian, establishing the validity of this result to arbitrary order in normalized Larmor radius. The global resonance condition for compressional modes is clarified in more detail than in previous works, confirming that the poloidal orbit-average of the cyclotron frequency gives the desired result at lowest order in Larmor radius. The global plasma response to the perturbation at each resonance is captured by a poloidal and gyroaverage of the perturbing potential. A “global gyroaveraging” of the potential is a natural by-product of this analysis which takes into account the changing of the magnetic field over an orbit. The resonance condition depends on two arbitrary integers which completely separately capture the effects poloidal non-uniformity and finite Larmor radius in generating sidebands. We learn that poloidal sidebands generated for compressional modes are dominated by the change in gyrofrequency over the orbit, which is very different to shear modes where the gyrofrequency only contributes via a finite Larmor radius effect. This increases the number of bounce harmonics required to compute the linear drive, giving a more complicated resonance map. An example calculation is given comparing resonance of shear and compressional modes in a published DIII-D case.

Compressional↗