Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TRANSPORT THEORY”

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

Refining the two-band model for highly compensated semimetals using thermoelectric coefficients

In studying compensated semimetals, the two-band model has proven extremely useful in capturing electrical conductivity under magnetic field, as a function of density and mobility of electronlike and holelike carriers. However, it rarely offers practical insight into magnetothermoelectric properties. Here, in this work, we report the field dependence of thermoelectric (TE) coefficients in a highly compensated semimetal NbSb 2 , where we find the Seebeck (𝑆 𝑥⁢𝑥 ) and Nernst (𝑆 𝑥⁢𝑦 ) coefficients increase quadratically and linearly with applied magnetic field, respectively. Such field dependence was predicted in previous work that studied a system of two parabolic bands, within semiclassical Boltzmann transport theory when the following two conditions are simultaneously met: 𝜔 𝑐 ⁢𝜏 ≫ 1 and tan ⁡𝜃 𝐻 ≪ 1, where 𝜔 𝑐 and 𝜏 refer to the cyclotron frequency and relaxation time, respectively, and 𝜃 𝐻 is the Hall angle. Under these conditions, we find the field dependence of the TE coefficients directly provides a relation between the electronlike (𝑛 𝑒 ) and holelike (𝑛 ℎ ) carrier densities, which in turn can be used to refine two-band model fitting. With this, we find the compensation factor ($\frac{|Δ⁢𝑛|}{𝑛_{𝑒}}$, where Δ⁢𝑛 = 𝑛 𝑒 −𝑛 ℎ ) of NbSb 2 is two orders of magnitude smaller than what was found in unrestricted fitting, resulting in a larger saturation-field scale for magnetoresistance. Within the same framework of the semiclassical theory, we also deduce that the thermoelectric Hall angle tan⁡ 𝜃 𝛾 = $\frac{𝑆_{𝑥⁢𝑦}}{𝑆_{𝑥⁢𝑥}}$ can be expressed as ($\frac{Δ⁢𝑛}{𝑛_{𝑒}}$ ×𝜔 𝑐 ⁢𝜏) −1 , which serves as a parameter to predict the degree of compensation. Our findings offer crucial insights into identifying empirical conditions for field-induced enhancement of TE performance and into engineering-efficient thermoelectric devices based on semimetallic materials.

electrical conductivity↗
DOE OSTI2026

The Physics Imposed on a Streaming Operator by Spherical Transport Problems

The streaming operator, which generates a displacement of a particle on a straight line at a constant speed in transport theory, is derived algebraically from a spherical coordinate formulation of Newton’s second law. This derivation leads to an operator that has more partial derivatives than a Cartesian coordinate formulation of the operator. The additional partial derivatives, which are with respect to the normalized velocity variables of a particle, take into account the intrinsic curvature of a ball. Moreover, these partial derivatives mitigate ray effects, which arise when a finite number of normalized velocities (also called directions or discrete ordinates) are used to simulate a continuous S 2 sphere of directions, by rotating the polar axis of the S 2 sphere into the radial direction of the coordinate system. As a consequence of this rotation, the number of actual discrete ordinates is greatly amplified to an enormous number of effective discrete ordinates by a multiplier that is equal to the number of patches that partitions a spherical surface. In addition to the derivation of the streaming operator, we provide in closed form a solution to the system of characteristic equations that is equivalent to the streaming operator. Furthermore, the solution to the system of characteristic equations enables the construction of an integral operator that is the inverse to the streaming operator. Examples in which ray effects are immensely mitigated by spherical coordinates are presented.

integral operator↗
NASA NTRS2024

Improved Modelling of SEP Event Onset Within the WSA–Enlil–SEPMOD Framework

Multi-spacecraft observations of solar energetic particle (SEP) events not only enable a deeper understanding and development of particle acceleration and transport theories but also provide important constraints for model validation efforts. However, because of computational limitations, a given physics-based SEP model is usually best suited to capture a particular phase of an SEP event, rather than its whole development from onset through decay. For example, magnetohydrodynamic (MHD) models of the heliosphere often incorporate solar transients only at the outer boundary of their so-called coronal domain – usually set at a heliocentric distance of 20–30 R⊙. This means that particle acceleration at coronal mass ejection (CME)-driven shocks is also computed from this boundary onwards, leading to simulated SEP event onsets that can be many hours later than observed, since shock waves can form much lower in the solar corona. In this work, we aim to improve the modelled onset of SEP events by inserting a “fixed source” of particle injection at the outer boundary of the coronal domain of the coupled WSA–Enlil 3D MHD model of the heliosphere. The SEP model that we employ for this effort is Solar Energetic Particle MODel (SEPMOD), a physics-based test-particle code based on a field line tracer and adiabatic invariant conservation. We apply our initial tests and results of SEPMOD’s fixed-source option to the 2021 October 9 SEP event, which was detected at five well-separated locations in the inner heliosphere – Parker Solar Probe, STEREO-A, Solar Orbiter, BepiColombo, and near-Earth spacecraft.

Solar eruptions↗
DOE OSTI2023

Synthesis and Observation of Emergent Phenomena in Epitaxial Heusler Compound Heterostructures

The proposal was on the synthesis and observation of emergent phenomena in epitaxial Heusler compound heterostructures. The large range of properties and number of Heusler compounds opens up a wide number of potential compounds that will exhibit emergent phenomena. The similarity, large range of relatively inexpensive, large area, high crystal quality, III-V bulk substrates, lattice parameters and the ability to tune the lattice parameters through ternary or quaternary III-V compound semiconductor epitaxial growth, makes III-V semiconductors an ideal choice for substrates for epitaxial growth of Heusler compounds. A number of half Heusler compounds have been predicted to exhibit band inversion, making them topological and are therefore expected to exhibit spin-momentum locked topological surface states with linear dispersion. Others are predicted to be semimetals with Weyl points and others semiconducting and magnetic. During the course of this grant, emphasis has been on investigating Heusler compounds with emergent phenomena and demonstrating the ability to tune their properties through alloying and strain. We have grown toplogical semimetal (PtLuSb, PtMnBi), Weyl (Co 2 MnAl, Co 2 TiGe), half metal (PtMnSb, Co 2 MnSi, Co 2 MnAl x Si 1-x , Co2FeAl), and semiconducting (CoTiSb, NiTiSn) and tuned their properties through alloying and epitaxial strain. We also investigated the closely related materials of rare-earth monopnictide, some of which have also been predicted to be topological. During the attempts to grow the PtMnBi, it was discovered that Bi, another predicted topological material when ultrathin, could be grown epitaxially on InSb, results for which are also reported here. The main focus for this effort has been on using variable photon energy and spin-dependent angle resolved photoemission (ARPES) to determine bulk band structure and surface states of pristine epitaxial films grown on III-V semiconductor and MgO substrates and correlate results with theory and transport measurements. Theory has been critical to interpretation of experimental results and has been essential in guiding experiments. The research benefited from several strong collaborations between the PIs and the beamline scientists at the Advanced Light Source at Lawrence Berkeley Laboratory, the Stanford Linear Accelerator Center (SLAC) at Stanford and at the Max Lab at Lund University in Sweden. The strong experiment - theory collaboration between the PI’s groups, the Palmstrøm group at UCSB and the Janotti group at the University of Delaware, has been critical for interpreting the experimental ARPES and magnetotransport measurements results and making predictions to guide experiments. Weekly interactive Zoom meetings made this work well. A collaboration between the Palmstrøm group and Dr. Alexei Fedorov at the Advanced Light Source (ALS) resulted in significant modifications to his end chamber to accommodate the vacuum suitcase that was designed and constructed in the Palmstrøm group at UCSB. In collaboration with beamline scientists, Drs. Makoto Hashimoto and Donghui Lu at SLAC, Palmstrøm made modifications to the vacuum suitcase and developed special sample holders that allowed samples to be grown in the Palmstrøm MBE systems at UCSB and transported in the UHV vacuum suitcase to SLAC for ARPES measurements. The development of the vacuum suitcase was essential for this grant as it has allowed variable photon energies to be used to identify surface versus bulk states on samples that could not be capped and decapped using As- or Sb-capping layers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Effects of aluminum content on thermoelectric performance of Al x CoCrFeNi high-entropy alloys

Introducing a non-regular distribution in the mass and bonding by including distinctly different elements can reduce the phonon transport even within structurally well-ordered materials. These distributions are a quality of all high-entropy alloys (HEAs), however the inclusion of aluminum in Al x CoCrFeNi is particularly impactful due to the large mismatch in atomic mass with other components. The resultant low phonon conductivity is a requirement for high thermoelectric performance, motivating the investigation of the effects of Al content on phonon transport as well as other thermoelectric properties. This work examines the phonon and electron transport and thermoelectric conversion properties with various Al contents (0 ≤ x Al ≤ 2) in this Cantor alloy using first-principles calculations, molecular dynamics, and semi-classical Boltzmann transport theory. The calculated phonon density of states and thermoelectric properties present reasonable agreements with experiments, including neutron scattering. A large reduction of phonon conductivity (k L ) is observed even with low x Al s, which we attribute to effective phonon scatterings by the large mass mismatch. However, its temperature dependence is not significant, demonstrating a minor contribution of interphonon scattering. In contrast, electrical conductivity (σ) and Seebeck coefficient (S) increase with temperature at higher x Al s with body-centered cubic structures. Therefore, the thermoelectric figure of merit (ZT) of Al x CoCrFeNi HEAs is enhanced by increasing the Al content mainly due to the increase of the thermoelectric power factor (σS 2 ) at high temperatures, while at low temperatures the phonon-scattering enhancement by mass mismatch is also important.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Connection Between Mass Flux Transport and Eddy Diffusivity in Convective Atmospheric Boundary Layers

Turbulence parameterizations for convective boundary layer in coarse-scale atmospheric models usually consider a combination of the eddy-diffusive transport and a non-local transport, typically in the form of a mass flux term, such as the widely adopted eddy-diffusivity mass-flux (EDMF) approach. These two types of turbulent transport are generally considered to be independent of each other. Using results from large-eddy simulations, here, we show that a Taylor series expansion of the updraft and downdraft mass-flux transport can be used to approximate the eddy-diffusivity transport in the atmospheric surface layer and the lower part of the mixed layer, connecting both eddy-diffusivity and mass-flux transport theories in convective conditions, which also quantifies departure from the Monin-Obukhov similarity (MOS) in the surface layer. This study provides a theoretical support for a unified EDMF parameterization applied to both the surface layer and mixed layer and highlights important correction required for surface models relying on MOS

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Collisionless cooling of perpendicular electron temperature in the thermal quench of a magnetized plasma

Thermal quench of a nearly collisionless plasma against an isolated cooling boundary or region is an undesirable off-normal event in magnetic fusion experiments, but an ubiquitous process of cosmological importance in astrophysical plasmas. Parallel transport theory of ambipolar-constrained tail electron loss is known to predict rapid cooling of the parallel electron temperature $T_{e\Vert}$ although $T_{e\Vert}$ is difficult to diagnose in actual experiments. Instead direct experimental measurements can readily track the perpendicular electron temperature $T_{e\bot}$ via electron cyclotron emission. The physics underlying the observed fast drop in $T_{e\bot}$ requires a resolution. Here two collisionless mechanisms, dilutional cooling by infalling cold electrons and wave-particle interaction by two families of whistler instabilities, are shown to enable fast $T_{e\bot}$ cooling that closely tracks the mostly collisionless crash of $T_{e\Vert}$. These findings motivate both experimental validation and reexamination of a broad class of plasma cooling problems in laboratory, space, and astrophysical settings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Selective permeation up a chemical potential gradient to enable an unusual solvent purification modality

The need for energy-efficient recovery of organic solutes from aqueous streams is becoming more urgent as chemical manufacturing transitions toward nonconventional and bio-based feedstocks and processes. In addition to this, many aqueous waste streams contain recalcitrant organic contaminants, such as pharmaceuticals, industrial solvents, and personal care products, that must be removed prior to reuse. We observe that rigid carbon membrane materials can remove and concentrate organic contaminants via an unusual liquid-phase membrane permeation modality. Surprisingly, detailed thermodynamic calculations on the chemical potential of the organic contaminant reveal that the organic species has a higher chemical potential on the permeate side of the membrane than on the feed side of the membrane. This unusual observation challenges conventional membrane transport theory that posits that all permeating species move from high chemical potential states to lower chemical potential states. Based on experimental measurements, we hypothesize that the organic is concentrated in the membrane relative to water via favorable binding interactions between the organic and the carbon membrane. The concentrated organic is then swept through the membrane via the bulk flow of water in a modality known as “sorp-vection.” We highlight via simplified nonequilibrium thermodynamic models that this “uphill” chemical potential permeation of the organic does not result in second-law violations and can be deduced via measurements of the organic and water sorption and diffusion rates into the carbon membrane. Moreover, this work identifies the need to consider such nonidealities when incorporating unique, rigid materials for the separations of aqueous waste streams.

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

Space-Dependent Calculation of the Multiplicity Moments for Shells With the Inclusion of Scattering

In recent work, we extended the methodology of multiplicity counting in nuclear safeguards by elaborating the one-speed stochastic transport theory of the calculation of the so-called multiplicity moments, i.e., the factorial moments of the number of neutrons emitted from a fissile item, following a source event from an internal neutron source [spontaneous fission and (α; n) reactions]. Calculations were made for solid spheres and cylinders, with the source being homogeneously distributed within the item. Recent measurements of the Rocky Flats Shells during the Measurement of Uranium Subcritical and Critical (MUSIC) campaign conducted by Los Alamos National Laboratory and assisted by the University of Michigan inspired us to extend the model to spherical shell geometry with a point source in the middle of the central cavity. Comparison of the calculated results with the experimental ones indicated that accounting for fission as the only neutron reaction (the standard procedure in the point model, adapted also in our work so far) was not sufficient for reaching good agreement with measurements. The model was therefore extended to include elastic scattering into the one speed formalism, whereas the effect of inelastic scattering was accounted for in an empirical way. After these extensions, good agreement was found between the calculated and the measured values. The paper describes the extension of the theory and provides concrete quantitative results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Prehydrodynamic evolution and its signatures in final-state heavy-ion observables

We investigate the effects of prehydrodynamic evolution on final-state observables in heavy-ion collisions using state-of-the art event simulations coupled to different prehydrodynamic scenarios, which include the recently developed effective kinetic transport theory evolution model KøMPøST. Differential flow observables are found to be mostly insensitive to the details of prehydrodynamic evolution. The main effect we observe is in the p T spectra, particularly the mean transverse momentum. However, at least part of this effect is a consequence of the underlying conformal invariance assumption currently present in such approaches, which is known to be violated in the temperature regime probed in heavy-ion collisions. This assumption of early time conformal invariance leads to an artificially large out-of-equilibrium bulk pressure when switching from (conformal) prehydrodynamic evolution to hydrodynamics (using the nonconformal QCD equation of state), which in turn increases the transverse momentum. Our study indicates that a consistent treatment of prehydrodynamic evolution in heavy-ion collisions requires the use of nonconformal models of early time dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Phase-space methods for neutrino oscillations: Extension to multibeams

The phase-space approach (PSA), which was originally introduced in Lacroix [] to describe neutrino flavor oscillations for interacting neutrinos emitted from stellar objects is extended to describe arbitrary numbers of neutrino beams. The PSA is based on mapping the quantum fluctuations into a statistical treatment by sampling initial conditions followed by independent mean-field evolution. A new method is proposed to perform this sampling that allows treating an arbitrary number of neutrinos in each neutrino beams. We validate the technique successfully and confirm its predictive power on several examples where a reference exact calculation is possible. We show that it can describe many-body effects, such as entanglement and dissipation induced by the interaction between neutrinos. Due to the complexity of the problem, exact solutions can only be calculated for rather limited cases, with a limited number of beams and/or neutrinos in each beam. The PSA approach considerably reduces the numerical cost and provides an efficient technique to accurately simulate arbitrary numbers of beams. Examples of PSA results are given here, including up to 200 beams with time-independent or time-dependent Hamiltonians. We anticipate that this approach will be useful to bridge exact microscopic techniques with more traditional transport theories used in neutrino oscillations. It will also provide important reference calculations for future quantum computer applications where other techniques are not applicable to classical computers. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Shifting and Splitting of Resonance Lines due to Dynamical Friction in Plasmas

We report a quasilinear plasma transport theory that incorporates Fokker-Planck dynamical friction (drag) and pitch angle scattering is self-consistently derived from first principles for an isolated, marginally unstable mode resonating with an energetic minority species. It is found that drag fundamentally changes the structure of the wave-particle resonance, breaking its symmetry and leading to the shifting and splitting of resonance lines. In contrast, scattering broadens the resonance in a symmetric fashion. Comparison with fully nonlinear simulations shows that the proposed quasilinear system preserves the exact instability saturation amplitude and the corresponding particle redistribution of the fully nonlinear theory. Even in situations in which drag leads to a relatively small resonance shift, it still underpins major changes in the redistribution of resonant particles. This novel influence of drag is equally important in plasmas and gravitational systems. In fusion plasmas, the effects are especially pronounced for fast-ion-driven instabilities in tokamaks with low aspect ratio or negative triangularity, as evidenced by past observations. The same theory directly maps to the resonant dynamics of the rotating galactic bar and massive bodies in its orbit, providing new techniques for analyzing galactic dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Experimental Measurements of Ion Diffusion Coefficients and Heating in a Multi-Ion-Species Plasma Shock

Collisional plasma shocks generated from supersonic flows are an important feature in many astrophysical and laboratory high-energy-density plasmas. Compared to single-ion-species plasma shocks, plasma shock fronts with multiple ion species contain additional structure, including interspecies ion separation driven by gradients in species concentration, temperature, pressure, and electric potential. We present time-resolved density and temperature measurements of two ion species in collisional plasma shocks produced by head-on merging of supersonic plasma jets, allowing determination of the ion diffusion coefficients. Here our results provide the first experimental validation of the fundamental inter-ion-species transport theory. The temperature separation, a higher-order effect reported here, is valuable for advancements in modeling HED and ICF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Invariant Forms of Dissolution Fingers

Dissolution of fractured and porous media introduces a positive feedback between fluid flow and reactant transport, leading to the emergence of pronounced, fingerlike channels. We investigate the formation of these structures using a microfluidic Hele-Shaw cell with a soluble bottom. Our experiments show that the shape of dissolution fingers is invariant and reveals itself over time as the fingers extend into the system. Here, by combining reactive-transport theory and conformal mapping techniques we derive these invariant forms. We relate these results to natural dissolution fingers in karst landscapes, and illustrate how to determine the groundwater flow rate responsible for their formation based on the finger shape.

Flows in porous media↗
DOE OSTI2020

LA17-V-InverseRadTransportProblem-PD2Kb. Final Project Report

In the inverse radiation transport problem, a radiation signal emanating from an unknown object is used to reconstruct geometric and material parameters of the object , e.g., the mass of special nuclear material (SNM), the isotopic composition of the SNM (e.g., the enrichment of 235U), the thickness of shielding, etc. The inverse radiation transport problem is one of the most fundamental and important in the emergency response mission, and has been addressed in previous NA-22 and NA-42/84 investments by using derivative-free optimization methods. Our project was aimed at (i) consolidating previous work by applying polynomial chaos expansion (PCE) for reducing significantly the number of forward model simulations required by derivative-free optimization methods; (ii) developing a general-purpose framework and stand-alone software for solving the inverse radiation transport problem by applying a new “predictive modeling” methodology (predictive modeling of coupled multi-physics systems, PM_CMPS) which enables the fusing of experimental and computational information; and (iii) developing several benchmarks (for unscattered gamma rays) for demonstrating the quantification of non-Gaussian features (e.g., skewness and kurtosis, which characterize asymmetries and long tails) of response distributions by deterministically and efficiently propagating parameter uncertainties (covariances and higher-order moments, where available) using high-order sensitivities computed by the recently developed comprehensive adjoint sensitivity analysis methodology (C-ASAM). With regard to (i), application of PCE, the most significant accomplishment of this project was an implementation of PCE into INVERSE that significantly reduced computation time for the DiffeRential Evolution Adaptive Metropolis (DREAM) and Differential Evolution (DE) methods. For inverse transport problems with measurements of neutron leakage, neutron multiplication, and neutron-induced gamma-rays, the solution times for the DREAM and DE methods were reduced by factors of 30-300, reducing run times from hours to minutes while maintaining the same accuracy as the standard (non-PCE) DREAM and DE. Speedups were further improved by exploring sparse grid quadrature sets using the Toolkit for Adaptive Stochastic Modeling and Non-Intrusive ApproximatioN (Tasmanian) toolset. An initial exploration of sparse grids indicated that further speedups factors of 5-10 are obtainable. The implementation of PCE in INVERSE is not in a form that users can easily access. With regard to (ii), solving the inverse problem using a predictive modeling paradigm, the most significant accomplishment of this project were the derivation and implementation of adjoint-based second derivatives in neutron and photon transport theory. There are two codes that compute these derivatives. SENSPG computes them for uncollided gamma rays and SENSMG computes them for neutron leakage. SENSMG is open-source code available on github. Other significant accomplishments included a stand-alone code, MULTI-PRED, that implements the predictive modeling paradigm, along with a detailed manual. The code was implemented in the existing INVERSE code for inverse problems, but a head-to-head comparison of inverse methods was not achieved. We also derived and implemented the capability to compute sensitivities for multiplicity measurements into the SENSMG code, and this capability was used in a data assimilation application, but it was not used for solving inverse problems. With regard to (iii), developing benchmarks to demonstrate the importance of high-order sensitivity and uncertainty analysis, the most significant accomplishment of this project were two test problems that showed how important second-order sensitivities are.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Characteristics of Ultralow Aspect Ratio Toroidal Plasmas and Non-Solenoidal Startup and Edge Stability Studies at Near-Unity Aspect Ratio (Final Report)

This grant supported a long-term research and development program to explore the physics and technology of spherical tokamak (ST) plasma behavior at near-unity aspect ratio ( A ). The PEGASUS Toroidal Experiment is a university-scale ( R ) ~ 0.40 m, α ~ .35 m, I p ≤ 0.3 MA, B T ≤ 0.15 T) experiment that accesses A → 1, to provide large values of the normalized plasma current, I N = I p /α B T , or equivalently the toroidal field utilization factor I p / I TF . The research activities focused on: 1) developing the technique of Local Helicity Injection (LHI) as a viable means of initiating and growing an ST plasma without the use of a central solenoid; 2) exploring the equilibrium and stability properties of spherical tokamak plasmas at extremely low aspect ratio; and 3) obtaining new insights in the properties of the enhanced confinement H-mode regime at A ~ 1. The initial (Phase I) facility employed a unique high-stress central solenoid Ohmic for current drive and heating with simple magnet current waveforms. A subsequent major upgrade (Phase II) offered more flexible control of the plasma evolution along with a powerful LHI system for plasma initiation, heating and sustainment. The LHI technique consists of injection of high-density current streams into vacuum toroidal and vertical magnetic fields. These streams become unstable and support magnetic reconnection to allow the system to relax into a tokamak-like magnetic configuration. This program focused on: 1) development of electron current sources that provide sufficiently high injection currents in the challenging edge plasma region; 2) validation of theoretical limits to achievable I p arising from magnetic helicity and energy conservation plus absolute limits arising from relaxation to a minimal energy states; 3) tests of a power balance model to describe the evolution of the plasma current I p (t); and 4) study of magnetic fluctuations and their role in the I p evolution. This culminated in the production of toroidal currents of I p ≥ 0.2 MA with only ~8 kA of injected current. Operation with Ohmic induction alone in Phase I created plasmas with soft limits on the achievable I p and field utilization ( I p / I TF ~ 1) due to the onset of large low-order (2/1 or 3/2) internal resonant tearing modes. With the expanded capabilities in Phase II operation, methods employed to mitigate those modes and achieve I p / I TF > 1 were: 1) I TF ramp-downs; 2) non-solenoidal startup via LHI; and 3) extremely low I TF Ohmic operations assisted by LHI current injector pre-ionization. Extensive use of LHI-driven startup and drive for non-solenoidal operation produced favorable current profiles to stably operate tokamak plasmas with unprecedentedly high I N = I p / aB T ~ 15 and I p / I TF > 2. With strong reconnection-driven anomalous ion heating, this in turn provided access to tokamak plasmas with B t up to 100%, a world record 2 times higher than that previously achieved in advanced or spherical tokamaks. This high B t regime had a large fraction (~50%) of the plasma volume in a region of an absolute minimum- B configuration. The MHD stability limit in this regime was the ideal external kink mode without a stabilizing wall. These results provided the first demonstration of access to tokamak plasmas with very high B t , high i N , high elongation κ, and low internal inductance ℓ i , and as such served as a proof-ofprinciple for access to this long-sought ST regime. The naturally low B T needed for stability at A ≤ 1.2 reduces the L-H mode power threshold, P LH , so that H-mode plasmas were achieved with Ohmic heating alone. The power threshold for the transition to the H-mode regime was found to be 10× the accepted values from the ITER international scaling, confirming initial trends indicated by experiments on MAST and NSTX, and suggesting that the physics behind this transition is not yet fully understood. The density dependence of the transition and the observed insensitivity to whether the plasma is bounded by a limiter or magnetic divertor are consistent with the recent FM 3 model. Multiple unstable toroidal magnetic modes were measured during ELM (Edge Localized Mode) crashes, with generally lower toroidal mode numbers than seen at higher A , which was attributed to the larger peeling mode instability drive at low A . First-ever edge current profile measurements through an ELM crash showed current-carrying filaments were generated, as seen in ELM simulations and in electromagnetic blob transport theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

EvaluateFlux, VARPOW, VARPEAK, and DIF3DtoVTK Utility Programs for DIF3D

The DIF3D-VARIANT solver is part of the DIF3D code which was developed at ANL for engineering level predictive analysis of nuclear reactor cores. Over the course of the last 20 years, the DIF3D-VARIANT option of DIF3D has become increasingly more used because of improvements in computer memory capabilities and performance. This is beneficial to the users of DIF3D as DIF3D-VARIANT is able to apply both p- and h- spatial mesh refinement strategies in diffusion and transport theory. Given its greater use, utility programs were created to assist in post-processing the DIF3D output. This document discusses the basic execution process of DIF3D and the output files it creates but is primarily focused on being the manual for the utility programs EvaluateFlux, VARPOW, VARPEAK, and DIF3DtoVTK.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Rational sub-nanometer manipulation of polymer morphology for efficient chemical separations

The ability to control polymer morphology on the sub-nanometer length scale has broad implications for chemical separations. To achieve such control on easily processable systems, this proposal focuses on the synthesis and characterization of polymers containing appended labile moieties that are easily detached by thermolysis or UV irradiation deep within the glassy state. Once liberated, these moieties can diffuse from the polymer matrix as gaseous products, leaving behind templated pathways for selective diffusion and sorption of small molecules. With a specific target of creating polymeric membrane materials with unprecedented diffusion and (ad)sorption characteristics for chemical separations, synthesis of new materials will be complemented with advanced metrologies, simulations, and evaluation of thermodynamic and transport theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗