Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three phase modeling”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Three-dimensional checkerboard spin structure on a breathing pyrochlore lattice

The standard approach to realize a spin-liquid state is through magnetically frustrated states, relying on ingredients such as the lattice geometry, dimensionality, and magnetic interaction type of the spins. While Heisenberg spins on a pyrochlore lattice with only antiferromagnetic nearest-neighbor interactions are theoretically proven disordered, spins in real systems generally include longer-range interactions. The spatial correlations at longer distances typically stabilize a long-range order rather than enhancing a spin-liquid state. Both states can, however, be destroyed by short-range static correlations introduced by chemical disorder. Here, using disorder-free specimens with a clear long-range antiferromagnetic order, we refine the spin structure of the Heisenberg spinel ZnFe 2 O 4 through neutron magnetic diffraction. The unique wave vector (1, 0, $\frac{1}{2}$) leads to a spin structure that can be viewed as alternatively stacked ferromagnetic and antiferromagnetic tetrahedra in a three-dimensional checkerboard form. Stable coexistence of these opposing types of clusters is enabled by the bipartite breathing pyrochlore crystal structure, leading to a second-order phase transition at 10 K. The diffraction intensity of ZnFe 2 O 4 is an exact complement to the inelastic scattering intensity of several chromate spinel systems which are regarded as model classical spin liquids. Our results challenge this attribution, and suggest instead of the six-spin ring mode, spin excitations in chromate spinels are closely related to the (1, 0, $\frac{1}{2}$) type of spin order and the four-spin ferromagnetic cluster locally at one tetrahedron.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Novel Phases and Dynamics of Multi-orbital Mott Insulators

The goal of this project is to gain theoretical insight into magnetism in the presence of spin-orbit coupling (SOC) for multi-orbital Mott insulators. Inspired by the exact solution of the Kitaev model that harbors a quantum spin liquid with novel excitations and struck by its relevance for materials with anisotropic orbital interactions, the PI will explore various new classes of 4d and 5d transition metal oxides. Starting with all electron Hamiltonians, the PI will derive minimal magnetic models. The aim will be to understand the role played by orbital frustration, even in the absence of any geometric frustration, in creating orbitally ordered as well as spin-orbital liquid phases. Orbital frustration arises primarily from the directional or anisotropic nature of d-orbitals in contrast to the isotropic nature of the spin degree of freedom. These models will be investigated by a variety of theoretical and numerical methods, including exact diagonalization, mean field theories, density matrix renormalization group and quantum Monte Carlo methods. Testable predictions for three experiments: nuclear magnetic resonance (NMR) and resonant x-ray scattering (RXS) to probe orbital ordering and pump-probe experiments to probe quantum dynamics will test the validity of these models for materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Volume IV. The DUNE far detector single-phase technology

The preponderance of matter over antimatter in the early universe, the dynamics of the supernovae that produced the heavy elements necessary for life, and whether protons eventually decay -- these mysteries at the forefront of particle physics and astrophysics are key to understanding the early evolution of our universe, its current state, and its eventual fate. DUNE is an international world-class experiment dedicated to addressing these questions as it searches for leptonic charge-parity symmetry violation, stands ready to capture supernova neutrino bursts, and seeks to observe nucleon decay as a signature of a grand unified theory underlying the standard model. Central to achieving DUNE's physics program is a far detector that combines the many tens-of-kiloton fiducial mass necessary for rare event searches with sub-centimeter spatial resolution in its ability to image those events, allowing identification of the physics signatures among the numerous backgrounds. In the single-phase liquid argon time-projection chamber (LArTPC) technology, ionization charges drift horizontally in the liquid argon under the influence of an electric field towards a vertical anode, where they are read out with fine granularity. A photon detection system supplements the TPC, directly enhancing physics capabilities for all three DUNE physics drivers and opening up prospects for further physics explorations. The DUNE far detector technical design report (TDR) describes the DUNE physics program and the technical designs of the single- and dual-phase DUNE liquid argon TPC far detector modules. Volume IV presents an overview of the basic operating principles of a single-phase LArTPC, followed by a description of the DUNE implementation. Each of the subsystems is described in detail, connecting the high-level design requirements and decisions to the overriding physics goals of DUNE.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Modeling the role of local crystallographic correlations in microstructures of Ti-6Al-4V using a correlated structure visco-plastic self-consistent polycrystal plasticity formulation

This paper presents a multi-level crystal plasticity-based simulation framework for modeling mechanical response and microstructure evolution of Ti-6Al-4V with α-lath/lamellar microstructures. Here, the model is a correlated structure visco-plastic self-consistent (CS-VPSC) formulation linking three scales: a single crystals micro-scale, a lath/lamellar colony meso-scale, and a lath/lamellar aggregate macro-scale. A selected hardening law for the evolution of critical resolved shear stress per slip system used in CS-VPSC is phenomenological. However, it adjusts the resistances of basal and prismatic slip systems based on the geometry of slip transfer between adjacent lamellae. Consistent with experimental evidences, the resolved shear stress on the pyramidal slip planes is dependent not only on the stress in the direction of slip but also on the two orthogonal shear stress components and the three normal stress components (non-Schmid effects). Electron backscatter diffraction (EBSD) data in conjunction with a procedure relying on α→β phase transformation is used to construct paired variants of α-laths/lamellae satisfying their local crystallographic correlations. The procedure fits volume fractions of individual laths/lamellae with the experimental EBSD data and selects a distribution of habit planes between adjacent variants with respect to the loading direction. The simulation framework is applied to interpret the deformation behavior in tension and compression along two sample directions of Ti-6Al-4V fabricated via laser powder bed fusion. Moreover, the model is used to simulate texture evolution during rolling of the material to large strains. It is demonstrated that the model is capable of predicting plastic anisotropy/asymmetry and the concomitant texture evolution. While the model reveals a significant effect of habit plane inclination with respect to the loading direction on yield stress, the comparison of the data and model predictions shows that a random distribution of habit planes fits the flow response. It is further inferred that the tension-compression asymmetry arises from the non-Schmid effects.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hidden symmetries, spin and charge of artificial magnetic monopoles

Here, we discuss the non-Abelian artificial magnetic monopoles associated with n-level energy crossing in quantum systems. We found that hidden symmetries reveal themselves as observables such as spin, hypercharge, and other physical degrees of freedom. We illustrated our results on concrete examples of two and three energy-level crossings. Our findings can be helpful for modeling various phenomena in physical systems such as SU(N) topological insulators and spintronic devices based on topological insulators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Experimental Investigation of Vapor Formation in Liquid CO2 Flow Through a Converging-Diverging Nozzle

Carbon dioxide is an attractive working fluid for many cycles, including for pumped thermal energy storage (PTES). A challenge with some proposed sCO2 PTES cycles is the operation of sCO2 machinery outside the typical bounds of experience, with local phase change from the liquid state being particularly unknown. Presently, there is insufficient data in the literature regarding multiphase CO2 to adequately design a multiphase-tolerant turbine, so generation of foundational data is required. This experimental study investigates the flow characteristics of sub-sonic liquid CO2 undergoing expansion and phase change in a converging-diverging nozzle. The nozzle is instrumented to measure static pressure, unsteady pressure, temperature, and density. The static pressure transducers are located at 27 axial locations to accurately characterize the pressure profile in the nozzle. High-accuracy RTDs are located at the entrance and exit of the nozzle, and three dynamic pressure transducers are strategically located to capture any unsteady phenomena. During testing, values of mass flow and nozzle inlet pressure are swept to vary the pressure drop and fluid properties. The measured total pressure drop in the nozzle is compared to a homogenous model and the Lockhart-Martinelli correlation method, with the latter predicting loss quite closely. The resulting data set is valuable for validating multiphase numerical models in a simple geometry before implementation of these models in turbomachinery design.

25 ENERGY STORAGE↗
DOE OSTI2021

Filamentation and focusing of electron beams due to interactions with plasma waves

Results of numerical modeling of the interaction of an electron beam propagating across relativistic plasma waves indicate that electron beam filamentation and focusing may occur under certain conditions. The model is based on solving the relativistic equation of motion in three dimensions for the individual electrons in a tenuous Gaussian beam, as they pass through a relativistic plasma wave. Several electron beam and plasma wave parameters were varied, and the results are summarized. One of the results is that the spacing of the electron beam filaments correlates with the wavelength of the plasma wave. The electron beam filaments appear as vertical slabs after the beam exits the plasma waves. The electron beam also compresses to a focus after it exits the plasma, and the focal distance depends on several parameters including the electron beam energy and phase velocity of the relativistic plasma wave. It is suggested that these focusing and filamentation phenomena may be the basis for diagnostics schemes for laser plasma interactions. The parameters used in the model were electron beam energies in the 5–50 keV range and plasma wave properties typical for the beat-wave produced by CO2 lasers, which correspond to the facilities available in our laboratory. The limitations of these results to lower energy density beam and plasma regimes and to higher energy density regimes will be discussed.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Significance of the structural configuration of B2 disorder in Co and Ti based Heusler alloys

We investigate here structural (at local and global levels) and transport properties for 𝑋 2 ⁢MnAl (𝑋= Co and Ti). Additionally, the magnetic properties were also studied for Ti 2 ⁢MnAl. Our x-ray diffraction results show that both the compounds stabilize in B2 disordered phase with cubic structure of 𝑃⁢𝑚⁢$\overline{3}$⁢𝑚 space group. Further, the structural configuration of the above disordered phase for both the compounds was identified using combined studies of x-ray absorption spectroscopy and multiple scattering calculations at the transition metal 𝐾 edges. Upon such identification, in the case of Co 2 ⁢Mn 1−𝑦⁢ Cr 𝑦 ⁢Al (𝑦= 0, 0.05, 0.1, 0.2) with change in 𝑦, we are able to establish a better connection quantitatively between the inverse of Mn-Co bonds and peak in the temperature-dependent resistivity. This highlights the crucial importance of a detailed understanding of the nature of B2 disorder. In the case of 𝑦=0, in the temperature range of study, the resistivity is driven by the functional form associated with (a) three-dimensional enhanced electron-electron Coulomb interaction scattering mechanism and (b) an unconventional one-magnon process. For Ti 2 ⁢MnAl, the transport shows metallic glasslike behavior at high temperature, while at low temperature it follows both the Cote-Meisel's model and quantum correction model. In this compound, the magnetic studies suggest the formation of superparamagnetic clusters in the paramagnetic matrix at low temperatures. Our density functional theory results are in line with the transport and magnetic properties. In literature, the spin polarization percentage (𝑃) for 𝑋= Co in B2 disordered phase is 76%. However, the present results emphasize the fact that in B2 disordered phase, the value of 𝑃 can range from 90% to 71% depending on the structural configuration introduced by swapping of the atomic positions of Mn and Al. For the compound under study, the value of percentage spin polarization obtained ranges between 82% to 85%. In addition, we also identify the origin of the difference in the shape of the Mn 3⁢𝑑 density of states for both the alloys. In conclusion, our results for 𝑋= Co alloy highlights the importance of identifying the specific structural configuration associated with a particular disorder category especially in the estimation of 𝑇 𝑐 and 𝑃 and for 𝑋= Ti, the physical properties can be tuned by varying the position of 𝐸 𝐹 and thereby its utilization in device applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Impact of titanium content on the thermo-mechanical and oxidation response of TiAlTa

This study aimed to analyze how the addition of titanium (Ti) to TiAlTa alloys affects their mechanical properties (Young’s modulus and hardness) and oxidation behavior. Alloys with three different Ti additions and equi-atomic Al and Ta were fabricated with nominal compositions of 33Ti-33Al-33Ta at%, 50Ti-25Al-25Ta at%, and 70Ti-15Al-15Ta at%. Phase identification of as-processed and compositionally homogenized alloys was conducted by coupling various electron microscopy, diffraction, and chemical analysis techniques. Composition and structures were modeled using CALPHAD-based phase predictions and compared with experimental results. Each of the alloys exhibited a unique set of microstructures with distinctly different ordered precipitates, dependent on cooling rate and chemical composition. Subsequently, nanoindentation tests were performed at temperatures of 25°C, 250°C, 500°C, and 750°C. Oxidation studies were conducted under static air at 750°C for up to 200 h. 33Ti-33Al-33Ta at% and 70Ti-15Al-15Ta at% alloys showed higher hardness and modulus than 50Ti-25Al-25Ta at% at room temperature. The 50Ti-25Al-25Ta at% alloy exhibited higher hardness and modulus values but lower oxidation resistance at 750°C when compared with the other specimens. The X-ray photoelectron spectroscopy experiments were conducted to further analyze the oxides present in the samples. Despite each sample displaying combinations of mixed oxide layers, TiO 2 was the preferred oxide forming in the 50Ti-25Al-25Ta composition, while the other samples preferred to form Al 2 O 3 . In conclusion, the findings from this systematic exploration of phase evolution in the TiAlTa alloy space provide insights for optimized material chemistries, processing, and properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ensemble Learning Based Convex Approximation of Three-Phase Power Flow

Though the convex optimization has been widely used in power systems, it still cannot guarantee to yield a tight (accurate) solution to some problems. To mitigate this issue, this paper proposes an ensemble learning based convex approximation for alternating current (AC) power flow equations that differs from the existing convex relaxations. The proposed approach is based on three-phase quadratic power flow equations in rectangular coordinates. To develop this data-driven convex approximation of power flows, the polynomial regression (PR) is first deployed as a basic learner to fit convex relationships between the independent and dependent variables. Then, ensemble learning algorithms such as gradient boosting (GB) and bagging are introduced to combine learners to boost model performance. Based on the learned convex approximation of power flow, optimal power flow (OPF) is formulated as a convex quadratic programming problem. The simulation results on IEEE standard cases of both balanced and unbalanced systems show that, in the context of solving OPF, the proposed data-driven convex approximation outperforms the conventional semi-definite programming (SDP) relaxation in both accuracy and computational efficiency, especially in the cases that the conventional SDP relaxation fails

Convex approximation↗
DOE OSTI2022

Assessment of four strain energy decomposition methods for phase field fracture models using quasi-static and dynamic benchmark cases

Abstract Strain energy decomposition methods in phase field fracture models separate strain energy that contributes to fracture from that which does not. However, various decomposition methods have been proposed in the literature, and it can be difficult to determine an appropriate method for a given problem. The goal of this work is to facilitate the choice of strain decomposition method by assessing the performance of three existing methods (spectral decomposition of the stress or the strain and deviatoric decomposition of the strain) and one new method (deviatoric decomposition of the stress) with several benchmark problems. In each benchmark problem, we compare the performance of the four methods using both qualitative and quantitative metrics. In the first benchmark, we compare the predicted mechanical behavior of cracked material. We then use four quasi-static benchmark cases: a single edge notched tension test, a single edge notched shear test, a three-point bending test, and a L-shaped panel test. Finally, we use two dynamic benchmark cases: a dynamic tensile fracture test and a dynamic shear fracture test. All four methods perform well in tension, the two spectral methods perform better in compression and with mixed mode (though the stress spectral method performs the best), and all the methods show minor issues in at least one of the shear cases. In general, whether the strain or the stress is decomposed does not have a significant impact on the predicted behavior.

Zhang, Shuaifang↗
DOE OSTI2020

Transient Coupled Chemical-Thermal-Fluid Field Simulation for Sealed Aluminum-clad Spent Nuclear Fuel Storage Canister

As the first step toward developing three-dimensional (3D) multi-physics computational fluid dynamics (CFD) model for unsealed and vented canister storage system, a 3D CFD model coupled with bulk gas radiolysis reactions was developed first for sealed DOE standard canisters filled with inert gas and trace amount of air and water. The workflow for constructing canister-scale 3D CFD models and coupling with gas phase radiolysis reactions were established, which can be readily extended to unsealed, vented canister storage system. This interim milestone report documents the theory of the model, workflow to establish radiolysis reaction network, and initial simulations of the evolutions of thermal fields and hydrogen gas concentrations within sealed DOE standard canisters over long period of time. In addition, a mesh refinement test was done to show that increasing the models mesh refinement had negligible impact upon the temperature profiles.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Status and Plan of DUNE at Fermilab

The discovery and the solid confirmation of neutrino flavor oscillation have been made throughout late 1990 to 2000s. Such oscillation is the clear proof of non-zero neutrino mass since this is caused by the fact that the neutrino mass eigenstates differ from that of the flavor. The neutrinos in the Standard Model, however, are massless leptons, and thus the neutrino flavor oscillation necessitates either the Standard Model to be significantly modified or a new theoretical paradigm to describe the universe better. In order to provide experimental information, it is essential for the properties of the neutrino oscillation be well understood through precision measurements of oscillation parameters, including the mixing angles, determination of the mass hierarchy between the three mass eigenstates and determination of the CP phase (Delta_CP) between neutrino and anti-neutrino oscillations. Among the several next generation neutrino experiments, this presentation covers the Deep Underground Neutrino Experiment (DUNE) which, with the high intensity neutrino beam facility, will help accomplish these goals. In this report, I will discuss the current status and plan of the DUNE experiment, including the results of the detector performance studies from the large scale prototype detectors operating at CERN and the timeline for construction and data taking of the experiment, as well as its physics prospects.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Building MCP-native hierarchical AI scientist ecosystems: a perspective on scaling multi-agent scientific discovery

Large language models (LLMs) are evolving from chatbots with limited tool-using capabilities to agentic AI systems that can perform deep research, assist in proposing hypotheses, help design experiments, automate data analysis, and draft scientific reports. However, there are currently two bottlenecks limiting LLMs' real-world impact on the broader scientific research community beyond academic demonstrations: lack of interoperability (repetitive manual tool-integration is required across scenarios) and the need for scalable coordination (unstructured communication and memory become brittle as the number of agents grows). In this Perspective, we argue that the next phase of agentic scientific discovery requires the development of an ecosystem of protocol-native agents and tools organized through hierarchies inspired by human society, beyond the current paradigm of a single monolithic “AI scientist”. We use Model Context Protocol (MCP) as a concrete example of an emerging interoperability layer for scientific tool and context exchange, and we propose three complementary pathways to increase the scaling capabilities of an MCP-native scientific ecosystem by addressing the composability issues: (1) MCP servers for high-value scientific tools maintained by domain experts, (2) automated transformation of existing code repositories into MCP services, and (3) autonomous invention and evolution of new agents and workflows. Finally, we provide a practical roadmap for scaling AI-driven scientific discovery by expanding tool supply and coordination in MCP-native scientific ecosystems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

From Electronic Structure to Ion Transport: Photoelectron Spectroscopy and Molecular Dynamics Simulations Reveal the Role of Anions in Lithium Battery Electrolytes

Electrolyte anions are pivotal for lithium battery performance, yet their fundamental electronic structural properties are not well understood. In this work, we employ a combination of negative-ion photoelectron spectroscopy (NIPES), ab initio calculations, and molecular dynamics (MD) simulations to investigate the electronic structures of three representative electrolyte anions. This multiscale approach enables us to elucidate how their intrinsic electronic properties govern anion–solvent interactions in gas-phase clusters, as well as lithium-ion (Li + ) solvation structures and ion transport behavior in the condensed phase. NIPES reveals that difluoro(oxalato)borate (DFOB – ), bis(fluorosulfonyl)imide (FSI – ), and bis(oxalato)borate (BOB – ) all exhibit high electron binding energies, with vertical/adiabatic detachment energies increasing from DFOB – (6.09/5.70 eV) to FSI – (6.80/6.10 eV) to BOB – (6.82/6.40 eV), correlating with enhanced oxidation stability. Ab initio calculations reveal that DFOB – /FSI – –solvent complexes bind Li + ∼ 10 kcal/mol stronger than BOB – series, aligning with the strength of a Li + –anion model. DFOB – exhibits pronounced charge localization on both oxygen and fluorine atoms, enabling their involvement in Li + coordination. In contrast, fluorine atoms in FSI – are largely electron-depleted and remain excluded from direct Li + binding. MD simulations further demonstrate that LiDFOB and LiFSI systems exhibit Li + diffusion coefficients three and five times higher than those of LiBOB across four common solvents. Notably, LiFSI salt in acetonitrile (AN) exhibits the fastest Li + diffusion among 12 electrolyte systems, highlighting the synergistic effect of FSI – and AN in promoting ion mobility. In conclusion, these findings provide a molecular-level understanding of the critical roles of anion and its microsolvation in optimizing Li + diffusion dynamics, once again emphasizing the positioning of FSI – and DFOB – as prime candidates for next-generation electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2023

Modified Virtual Oscillator Based Operation of Grid Forming Converters with Single Voltage Sensor

Virtual oscillator based control architecture for grid forming converter with inner current control loop for additional objective accomplishments is proposed in this paper. In general, virtual oscillator is based on Lienard type of oscillator equation with a cubic nonlinearity. Virtual oscillator alone cannot accommodate unbalanced or harmonic mitigation while operating converters in grid forming mode. Therefore, to accomplish better tracking performance, incorporating the provision for nonlinear and/or unbalanced loads, an inner current control loop is utilized which is based upon Lyapunov energy function type control architecture due to its advantages. To reduce the overall cost of implementation and ensure proper phase difference in the generated voltages, the implemented virtual oscillator operates with only one voltage feedback from the point of common coupling as presented in this paper. To verify the effectiveness of the proposed approach, the overall system has been modeled in MATTLAB/Simulink and PLECS domain and case studies showing the successful production of low harmonic load voltages under nonideal loading condition along with other important results are presented in this paper.

dual second order generalized integrator↗
DOE OSTI2020

Complex Disorder in Type-I Clathrates: Synthesis and Structural Characterization of A8GaxSn46−x (A = Rb, Cs; 6.9 < x < 7.5)

Exploratory studies in the systems Rb–Ga–Sn and Cs–Ga–Sn yielded the cubic type-I clathrates with refined compositions Rb8GaxSn46−x and Cs8GaxSn46−x (6.9 < x < 7.5). Nearly single-phase materials with good crystallinity were obtained from stoichiometric reactions of the elements. The structures were characterized by means of single-crystal X-ray diffraction methods. Both Rb8GaxSn46−x and Cs8GaxSn46−x represents cases, where a Group 13 element randomly substitutes a Group 14 element in the structure. The extent of Ga/Sn mixing is apparently governed by the drive of the system to achieve an optimal valence electron count, and hence, Rb8GaxSn46−x and Cs8GaxSn46−x (x ≈ 8) can be regarded as Zintl phases. This notion is supported by structure refinements on a multitude of single-crystal X-ray diffraction data, which also confirm that both types of cages in the cubic type-I structure are fully occupied by Rb and Cs atoms. The open-framework, comprised of 46 nodes per formula unit, adapts to the incorporation of nearly eight Ga atoms within the matrix of Sn, whereby small, short-range distortions result. The exact nature of these effects is still unclear, as so far, the structural variations could only be modeled as both positional and occupational disorder at one of three framework sites. Since vacancies in the structures of the binary type-I clathrates A8Sn46−x☐x (A = Rb, Cs; ☐ = missing Sn atom) are also known to cause local distortions, the latter were also synthesized with the same protocols used for the synthesis of A8GaxSn46−x and structurally re-analyzed. The results from the latter studies confirm that homogeneity issues abound, and that the final structures/compositions are an intricate function of the experimental conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Center for Tokamak Transients Simulations

The major goals for the Center for Tokamak Transients Simulations (CTTS) were to develop detailed numerical simulations of disruptive transients in tokamak experiments and to use them in theoretical studies of the thermal-quench (TQ) and current-quench (CQ) phases of disruption. The University of Wisconsin-Madison component of CTTS contributed relevant development for the NIMROD code (https://nimrodteam.org) and applied it to disruptions that involve global vertical displacement. Code-development contributions include boundary conditions that represent sheath effects which occur when tokamak plasma is in contact with the surrounding vessel wall during vertical displacement events (VDEs). We also implemented a reduced model of energetic runaway electrons (REs) to simulate the effect of these energetic electrons on the CQ phase of tokamak disruptions. Our computations of VDEs include 2D and 3D verification through benchmarks with two other widely used macroscale plasma simulation codes, M3D-C1 and JOREK. Over nonlinear evolution to the final termination of plasma current, the 2D simulation results on magnetic axis location, toroidal current, and halo current track each other well. With the 3D benchmark computation, all three codes found the same qualitative behavior of the asymmetric instabilities having increasing growth rate after contact with the wall, followed by nonlinear excitation of other toroidal harmonics. Our NIMROD computations for an idealized configuration show agreement with an analytical assessment by V. D. Pustovitov [Nucl. Fusion 55, 113032 (2015)] that the plasma and non-ideal wall remain in force-balance such that net forces only result from magnetic stresses over the outer surface of the wall. Our study of the influence of boundary conditions on VDEs shows that extended-MHD simulations of VDEs are sensitive to electron energy transport to the wall and that boundary conditions on flow are only significant to the extent that they influence the electron energy transport. Resistive-MHD NIMROD simulations of the tokamak discharges in the Madison Symmetric Torus show that magnetic perturbations of poloidal wavenumber m=3 lead to chaotic magnetic topology over the edge region, which helps deconfine energetic electrons. A second area of RE study shows that the coupled resistive MHD/reduced RE model reproduces the resistive-hose particle-beam instability. With respect to educational opportunities, this effort involved and supported four graduate students at the University of Wisconsin-Madison and one postdoctoral associate. The results of this study have been disseminated through journal publications, conference presentations, technical reports, and PhD dissertations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗