SEARCH · Search NASA

Results for “theoretical developments”

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

SourceResearch connectionFrom yearThrough year

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

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Over the last few years, extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high degree of precision. An appealing and challenging route toward engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, i.e., to provide a unique framework that allows us to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory. This article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, such as the new theoretical framework of quantum electrodynamics density-functional formalism for the description of novel light–matter hybrid states. Those advances, and others being released soon as part of the Octopus package, will allow the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (quantum electrodynamical-materials).

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

Theoretical understanding of stability of the oxygen electrode in a proton-conductor based solid oxide electrolysis cell

The oxygen electrode in a proton-conductor based solid oxide cells is often a triple-conducting material that enables the transport and exchange of electrons (e - ), oxygen ions (O 2- ), and protons (H + ), thus expanding active areas to enhance the oxygen electrode activity. In this work, a theoretical model was developed to understand stability of tri-conducting oxygen electrode by studying chemical potentials of neutral species (i.e., μ o 2 , μ H 2 , and μ H 2 O ) as functions of transport properties, operating parameters, and cell geometry. Our theoretical understanding shows that (1): In a conventional oxygen-ion based solid oxide cell, a high μ o 2 (thus high oxygen partial pressure) exists in the oxygen electrode during the electrolysis mode, which may lead to the formation of cracks at the electrode/electrolyte interface. Further, while in a proton-conductor based solid oxide cell, the μ o 2 is reduced significantly, suppressing the crack formation, and resulting in improved performance stability (2). In a typical proton-conductor based solid oxide electrolyzer, the dependence of μ o 2 on the Faradaic efficiency is negligible. Hence, approaches to block the electronic current can improve the electrolysis efficiency while achieving stability (3). The difference of the μ o 2 (thus p o 2 ) between the oxygen electrode and gas phase can be reduced by using higher ionic conducting components and improving electrode kinetics, which lead to further improvement of electrode stability.

08 HYDROGEN↗

DOE OSTI2022

Experimental and kinetic modeling study of the low- temperature and high-pressure combustion chemistry of straight chain pentanol isomers: 1-, 2- and 3-Pentanol

Pentanols have received significant attention as a potential alternative fuel or fuel additive owing to their high energy densities and low vapor pressure. The development of robust chemical kinetic models for alternative fuels which can provide accurate and efficient predictions of combustion performance across a wide range of engine relevant conditions is important in developing cleaner, more efficient combustors. Although the high temperature oxidation kinetics of pentanol isomers has been researched considerably, their low temperature combustion chemistry needs further investigation. While previously proposed low temperature mechanisms for 1-pentanol based on analogy and rate rules need further refinement, the low temperature oxidation kinetics of 2-pentanol and 3-pentanol has not been studied previously by any means, experimentally or theoretically. A newly developed kinetic mechanism is presented in this work for the three straight chain pentanol isomers: 1-, 2- and 3-pentanol. Low temperature kinetics is based on a recent study by Lockwood et al., 2022 [20] involving theoretical calculations at the CCSD(T)/cc-pV∞Z level of theory for the oxidation pathways involving alcohol peroxy radicals. Rate of production analyses performed in this study highlight the importance of the newly added pressure-dependent reactions of the α-alcohol peroxy radical forming an RO 2 adduct. While the α-alcohol fuel radical reacts with O 2 to directly decompose via a chemically activated pathway at low pressures, the formation of the RO 2 adduct is favored at high pressures. Furthermore, the detailed model is comprehensively validated against new ignition experiments at low temperature and high pressure, together with the wide range of data available in the literature. Both qualitative and quantitative predictions of the experimental data using the proposed kinetic model are satisfactory for all three pentanol isomers studied here.

1-pentanol↗

DOE OSTI2021

Computation of linear MHD instabilities with the multi-region relaxed MHD energy principle

Over the last decade, a variational principle based on a generalisation of Taylor's relaxation, referred to as multi-region relaxed magnetohydrodynamics (MRxMHDs) has been developed. The numerical solutions of the MRxMHD equilibria have been constructed using the Stepped Pressure Equilibrium Code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 112502). In principle, SPEC could also be established to describe the MRxMHD stability of an equilibrium. In this work, a theoretical framework is developed to relate the second variation of the energy functional to the so-called Hessian matrix, enabling the prediction of MHD linear instabilities of cylindrical plasmas, and is implemented in SPEC. The negative and positive eigenvalues of the Hessian matrix predict the stability of an equilibrium. Verification studies of SPEC stability results with the M3D-C 1 code and the tearing mode $\Delta^{^{\prime}}$ criterion have been conducted for ideal and resistive MHD instabilities, respectively, in a pressureless cylindrical tokamak, and have shown good agreement. Overall, our stability analysis is a critical step towards understanding the MHD stability of three-dimensional MHDs where nested flux surfaces, magnetic islands and stochastic regions co-exist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Agent-Based Coordination Scheme for PV Integration (ABC4PV)

Renewables and especially photovoltaics (PV) have benefitted significantly from a host of incentives and policies targeted toward enhanced integration and adoption of specific energy technologies. However, with the push to move forward into a subsidy-free market framework, behind-the-meter residential PV applications have generally struggled to retain their value (unlike utility scale and commercial projects) [1]. This project focused on developing control-theoretic solutions aimed at improving the integration and interaction of behind-the-meter residential PV with other distribution system assets (controllable and non-controllable) to enhance the integrated value of residential PV. To this end, a suite of decentralized control methodologies have been developed to enable effective coordination and control of behind-the-meter residential load customers’ PV, battery storage systems (BSS), controllable loads and other similar assets within a distribution feeder. This interaction aims at procuring energy savings and, thus, energy bill savings. The main source of savings is drawn from reducing the effect of demand charge pricing and is realized at the feeder level, assuming community level interaction and management among the aforementioned assets. Optimal control of the assets is implemented with a distributed optimization methodology, leveraging consensus-based algorithms. The results gathered from the optimal control simulations demonstrates that the savings can be duly achieved and the algorithm decision times (to dynamically control asset set points, for example) are fast. As for the overall efficiency of PV+BSS systems, to procure energy savings from curtailment of the demand charge pricing effects, the optimal control is set up so as to minimize the variance of the load for all customers, throughout a feeder and throughout time in a rolling horizon scheduling with model predictive control. The control takes into account inter-temporal electrochemical storage (battery) degradation costs: specifically, we have developed a long-term lifetime model for the BSS that weighs in the effect of the degradation factor in the dispatch formulations, thus, a considerable operating cost that affects energy decision making. The levelized cost of energy (LCOE – redefined for the purpose of quantifying asset integration effectiveness through the customers’ energy cost) is shown to be below the threshold set for the combined PV+BSS topology of $ 0.14/kWh for multiple cases of PV penetration all the way up to 50%, provided that a policy of shared ownership of and savings is in place. Further, the LCOE calculated for the case before the deployment PV+BSS systems is also achievable, i.e. the deployment of PV+BSS, if planned and scheduled optimally. will have no effect on customers’ energy costs. From the control methodology viewpoint, the developed consensus-based algorithms are shown to converge for a wide range of problem cases (spanning normal operating scenarios and contingencies), guaranteeing dispatch solutions under forecasting errors, communication break-downs and cyber-security attacks. The proposed control solutions are scalable and real-time implementable, with dispatch computations and device set-point updates converging in less than 2s in most practical instances of the above events.

14 SOLAR ENERGY↗

DOE OSTI2021

Thermal transport in defective and disordered materials

With significant recent advancements in thermal sciences—such as the development of new theoretical and experimental techniques, and the discovery of new transport mechanisms—it is helpful to revisit the fundamentals of vibrational heat conduction to formulate an updated and informed physical understanding. The increasing maturity of simulation and modeling methods sparks the desire to leverage these techniques to rapidly improve and develop technology through digital engineering and multi-scale, electro-thermal models. With that vision in mind, this review attempts to build a holistic understanding of thermal transport by focusing on the often unaddressed relationships between subfields, which can be critical for multi-scale modeling approaches. For example, we outline the relationship between mode-specific (computational) and spectral (analytical) models. We relate thermal boundary resistance models based on perturbation approaches and classic transmissivity based models. We discuss the relationship between lattice dynamics and molecular dynamics approaches along with two-channel transport frameworks that have emerged recently and that connect crystal-like and amorphous-like heat conduction. Throughout, we discuss best practices for modeling experimental data and outline how these models can guide material-level and system-level design.

42 ENGINEERING↗

DOE OSTI2025

Photon–photon chemical thermodynamics of frequency conversion processes in highly multimode systems

Abstract Frequency generation in highly multimode nonlinear optical systems is inherently a complex process, giving rise to an exceedingly convoluted landscape of evolution dynamics. While predicting and controlling the global conversion efficiencies in such nonlinear environments has long been considered impossible, here, we formally address this challenge even in scenarios involving a very large number of spatial modes. By utilizing fundamental notions from optical statistical mechanics, we develop a universal theoretical framework that effectively treats all frequency components as chemical reactants/products, capable of undergoing optical thermodynamic reactions facilitated by a variety of multi-wave mixing effects. These photon–photon reactions are governed by conservation laws that directly determine the optical temperatures and chemical potentials of the ensued chemical equilibria for each frequency species. In this context, we develop a comprehensive stoichiometric model and formally derive an expression that relates the chemical potentials to the optical stoichiometric coefficients, in a manner akin to atomic/molecular chemical reactions. This advancement unlocks new predictive capabilities that can facilitate the optimization of frequency generation in highly multimode photonic arrangements, surpassing the limitations of conventional schemes that rely exclusively on nonlinear optical dynamics. Notably, we identify a universal regime of Rayleigh–Jeans thermalization where an optical reaction at near-zero optical temperatures can promote the complete and entropically irreversible conversion of light to the fundamental mode at a target frequency. Our theoretical results are corroborated by numerical simulations in settings where second-harmonic generation, sum-frequency generation and four-wave mixing processes can manifest.

Optics↗

DOE OSTI2025

Operation of helium sub-atmospheric multistage cryogenic centrifugal compressor trains: Part 1 – Steady state modeling and speed ratio selection

Helium cryogenic systems which can provide cooling below the normal boiling point of helium (approximately 4.2 K) are often required by superconducting radio-frequency niobium resonators utilized in modern high-energy particle accelerators. Achieving temperatures below 4.2 K generally involves operating a cryogenic vessel with liquid helium under sub-atmospheric conditions, thereby lowering the saturation pressure and corresponding saturation temperature. Over the last several decades, multi-stage cryogenic centrifugal compressor trains (CC’s) have been operated efficiently and reliably within large-scale cryogenic systems to continuously evacuate helium vapor generated by a device within the vessel, maintaining sub-atmospheric conditions in the vessel while pressurizing the return vapor to above atmospheric conditions. Traditionally, these CC systems have been operated using empirically derived control philosophies and insight gathered from previous operational experience. Recent efforts at the Facility for Rare Isotope Beams (FRIB) have been aimed at the development of a theoretical basis to characterize the operation of multi-stage cryogenic centrifugal compressor train and utilizing predictive model results to generate control parameters. The objective of this research was identifying operational points which adequately balance cryogenic system efficiency, stability, and overall ease of operation. Furthermore, this manuscript provides an overview of the predictive model development, characterization of the FRIB cryogenic centrifugal compressors and implementation of the predicted performance results during steady-state system operation.

Compressor train control↗

DOE OSTI2025

Unlocking New Regimes in Fractional Quantum Hall Effect with Quaternions

We demonstrate that formulating the composite-fermion theory of the fractional quantum Hall (FQH) effect in terms of quaternions greatly expands its reach and opens the door into many interesting issues that were previously not amenable to quantitative theoretical investigation. As an illustration, we explore the possibility of a nematic or a charge-density wave instability of the composite-fermion Fermi sea at half-filled Landau level and of the nearby FQH states by looking for a gap closing instability of the neutral magneto-roton excitation. As a result, our quaternion formulation of the FQH effect has been inspired by mathematical developments in the theoretical analyses of gravitational wave modes and cosmic microwave background radiation, where an important role is played by spin-weighted spherical harmonics that are nothing but monopole harmonics appearing in the spherical geometry for the FQH effect.

Composite fermions↗

DOE OSTI2023

Scaling hyporheic nitrogen cycling in large river alluvial aquifers

Research Element 1: Denitrification is a difficult flux to measure when fluxes are low and when the ecosystem is difficult to sample. These attributes describe Nyack floodplain, a large gravel-bedded aquifer adjacent to the Middle Fork Flathead River. Given that denitrification is simply respiration with nitrate as a terminal electron acceptor, we reasoned that N 2 production should negatively covary with oxygen respiration so long that oxygen is available. We sampled 16 wells 8 times throughout a year and analyzed these water samples on a membrane-inlet mass spectrometer for N 2 :Ar and O 2 . We fit multilevel models to N 2 :Ar as a function of dissolved oxygen. Within wells, there was no relationship between O 2 and N 2 :Ar. Among wells N 2 :Ar increased as O 2 declined. Based on the slope of this relationship, ~ -0.07, and estimated oxygen respiration in Nyack derived from prior studies, a possible range for denitrification is 0-8 mmol m -2 d -1 corresponding to a turnover rate for nitrate of 0-0.2/d. This estimated turnover rate is in the lower 20% of rates measured in the LINX 2 study, yet shows that even in a low nitrate (0.005 mmol/L) aquifer denitrification can be a substantial sink for dissolved N. Research Element 3 and 4: We developed a new theoretical approach to describe the relationship between water inflow, water outflow, and water age in porous media control volumes under dynamic flow conditions. Our approach had three distinct steps. First we described the application of well-established equations borrowed from Chemical Engineering to describe water age distributions in the hyporheic zone under steady state conditions. We published the associated results in the scientific journal PLOS ONE. Second, we developed a new theory ("flow spans") to describe how predicted steady-state water age distributions will respond to time varying imbalances in hyporheic discharge and recharge. Finally, we developed a proof-of-concept simulation model based on the new theory and used the model to examine the effects of flow regulation on a simulated river with characteristics representative of western U.S. montane rivers with expansive hyporheic zones.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Characterizing the impact of water uptake on light absorption by aerosol particles

Absorbing aerosols (AA) have an important impact on the global radiation budget. The composition and properties of AA can vary substantially throughout the atmosphere, depending on the particle source and the influence of chemical aging. Uncertainties associated with the radiative effects of AA remain substantial. A key contributor to this uncertainty is understanding the extent to which water uptake alters absorption by AA particles and how this depends on particle composition. We have used new and existing experimental tools to systematically characterize the relationship(s) between AA chemical and physical characteristics, relative humidity, and mixing of AA with other aerosol components through a combination of laboratory and field measurements. We have developed fundamental understanding through laboratory experiments that will enable interpretation of field observations at DOE ARM sites and facilitate process-based improvements in the simulation of absorption by carbonaceous aerosols in climate models. In the laboratory experiments we have considered a variety of different AA types, including fullerene soot (a surrogate for refractory black carbon, or BC) and nigrosine (a surrogate for partially soluble brown carbon, or BrC). By examining a range of AA types we will develop insight into how the chemical and physical properties of AA particles (e.g. whether they are refractory or soluble) impact the influence of water uptake on absorption. These AA particles will also be mixed—both internally and externally—with compounds having a wide range of hygroscopic properties to establish how the impact of water uptake on absorption will vary with photochemical processing. These laboratory observations have been compared with commonly used theoretical models, enabling development of empirical models that will facilitate improved representation of absorption by AA particles in regional and global climate models. Overall, our study has contributed to the mission of the Atmospheric System Research program by quantifying how interactions between aerosols and radiation depend on relative humidity. Through this, our work will improve understanding and model representation of aerosol processes as they affect the Earth’s radiation budget.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Dynamic density functional theory of polymers with salt in electric fields

Here we present a dynamic density functional theory for modeling the effects of applied electric fields on the local structure of polymers with added salt (polymer electrolytes). Time-dependent equations for the local electrostatic potential and volume fractions of polymer, cation, and anion of added salt are developed using the principles of linear irreversible thermodynamics. For such a development, a field theoretic description of the free energy of polymer melts doped with salts is used, which captures the effects of local variations in the dielectric function. Connections of the dynamic density functional theory with experiments are established by relating the three phenomenological Onsager’s transport coefficients of the theory to the mutual diffusion of electrolyte, ionic conductivity, and transference number of one of the ions. The theory is connected with a statistical mechanical model developed by Bearman and Kirkwood [J. Chem. Phys. 28, 136 (1958)] after relating the three transport coefficients to friction coefficients. The steady-state limit of the dynamic density functional theory is used to understand the effects of dielectric inhomogeneity on the phase separation in polymer electrolytes. The theory developed here provides not only a way to connect with experiments but also to develop multi-scale models for studying connections between local structure and ion transport in polymer electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

First-Principles Simulation of Beam-Induced Processes Underlying Atomic Manipulation in Electron Microscopes

The development of experimental methods and apparatuses capable of promoting atomically precise material manipulations holds great promise for realizing the ultimate limit of feature miniaturization in materials and devices. The ability to modify materials atom by atom is anticipated to usher in new technologies in areas as diverse as separation science, medicine, and quantum information science. Historically, scanning probe-based techniques have been the most prominent approaches in this space. However, these methods are best suited for the manipulation of surface-exposed regions of materials, as the strong perturbations required for bond scission are delivered most effectively to atoms in the near-proximity to the scanning probe. In contrast, convergent electron beams with energies tuned slightly below the threshold for inducing irreversible knock-on damage have recently been employed (within scanning transmission electron microscopy) to promote atomic-scale bond rearrangements in various beam-stable solids. Currently, however, the efficiency and selectivity of beam-induced atomic manipulation processes with focused electron beams are such that long irradiation times are required to induce a desired atomic rearrangement. With a better understanding of the underlying physics dictating the outcome of a given irradiation event, methods can be devised to improve the efficiency of these techniques so that their promise can be fully realized through widespread adoption.To this end, this Account details our recent efforts to develop and apply tractable first-principles simulation approaches for studying the response of materials to electric beam-like external electric potentials applied in real space. We briefly review the concepts and capabilities in the area of atomically precise materials manipulation and review the early demonstrations of accomplishments in this area, focusing on studies using scanned convergent electron beam probes in particular. We expound upon the depth of the challenge and identify critical shortcomings of theoretical methods that have previously been employed in the simulation of beam-induced processes. We then describe the computational methods that we have generalized from the concepts and tools most commonly applied to the study of molecular photochemistry and how our adaptations of these methods can be employed to capture the relevant dynamical phenomena for beam-induced processes ranging from the initial electron scattering to the ensuing multistate reactions. Here, we contextualize these methods within the current state of the art in this area, which has historically focused primarily on the simulation of inelastic image formation in the electron microscope for the purpose of interpreting the results of quantitative electron microscopy experiments. We demonstrate that the spatial distribution of state-specific excitation rates due to the presence of an external (probe) electric charge is inhomogeneous, such that irradiation at particular locations in materials can favor specific electronic transitions (and disallow others). In addition to the potential for excited-state reaction pathways to be accessed through the initial inelastic scattering of the tightly focused electron beam from the targeted atoms, we also identify favorable conditions for the electronically nonadiabatic evolution of the highly vibrationally excited system to open complex multistate reaction pathways. Implications of the early results for understanding the mechanisms and potential routes to improved efficiency and selectivity in beam-induced reactions are discussed. We conclude with a summary of the current state of theory and modeling capabilities in this area and provide our perspective on future directions for theoretical and experimental developments that we view as crucial to advancing the use of convergent electron beams in mode-specific, atomically precise platforms for direct-write materials modifications.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Summary: SEE4GEO

The seismoelectric effects technique (SEE) is a new and innovative approach for geothermal subsurface imaging and monitoring at reservoir scale. The objective of this project is to assess SEE in terms of data acquisition, cost and quality, and to determine its capability in comparison with classical imaging and monitoring techniques, particularly decoupled seismic and electromagnetic methods. This will be achieved by (1) development of a fast, true 3D numerical package, handling SEE imaging and subsurface properties characterization, including resistivity and permeability, (2) laboratory experiments performed in a controlled environment to define optimal deployment design, data quality, and inform field deployment, and (3) field surveys to ultimately test and draw lessons for practical use of SEE technology. There is a relatively extensive body of work in the literature on SEE, and members of this consortium have been involved in theoretical and numerical development of SEE modeling as well as laboratory experiments. Nevertheless, to our knowledge few, if any, documented efforts have been specifically targeting the use of SEE for geothermal subsurface imaging and monitoring. The strength and originality of our proposal rely on an integrated approach leveraging numerical, laboratory and field experiments, to properly document the practical use of SEE. Through this process, SEE in-hand technology for the geothermal industry will be able to progress from a TRL 1 to TRL 3.

58 GEOSCIENCES↗

DOE OSTI2022

Seismoelectric Effects for Geothermal Resources Assessment and Monitoring (SEE4GEO)

58 GEOSCIENCES↗

DOE OSTI2023

Laser-Induced Ultrafast Demagnetization and Spin Switching (Final Report)

The information technology revolution demands bigger and faster magnetic storage devices. The funded research focuses on a new frontier - Laser-induced ultrafast all-optical spin switching (AOS), where a single laser pulse can switch spins from one direction to another permanently, without a magnetic field, on a time scale of 1-10 picoseconds, 2-3 orders of magnitude faster than the existing technology. The funded research has successfully resolved a long-standing issue in ultrafast demagnetization whether the lattice contributes to the demagnetization. We developed a high level first-principles theory and employed a massively parallel computer to pin down the details of the physical mechanisms that are inaccessible to other methods. This unique feature allows us to compute the force on each atom in both bulk and monolayer systems. We find that it is indeed true there is a nonzero force on each atom, but the idea of vibration contribution is not crucial to our basic understanding of physics. This paves the way to future experimental, theoretical and technological developments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Game-Theoretic Approach for Grace-Period Policy in Supercomputers

Job scheduling at supercomputing facilities is important for achieving high utilization of these valuable resources while ensuring effective execution of jobs submitted by users. The jobs are scheduled according to their specified resource demands such as expected job completion times, and the available resources based on allocations. Jobs that overrun their allocated times are terminated, for example, after a grace-period. It is non-trivial and often very complex for users to accurately estimate the completion times of their jobs, and consequently they face a dilemma: underestimate the job time to have a higher priority and risk job termination due to overrun, or overestimate it to ensure its completion and risk its delayed execution. In this paper, we investigate whether providing grace-period can benefit facility performance by developing a game- theoretic model between a facility provider and multiple users for a simplified scheduling scenario based on job execution times. We present closed-form expressions for the provider’s and user’s best-response strategies to maximize their respective utility functions. We describe conditions under which offering a grace-period is advantageous to both facility provider and users by deriving the Nash equilibrium of the game.

He, Fei↗

DOE OSTI2024

Advanced Closures for Nonlinear Fluid Simulations of Plasmas (Final Report)

Advanced Closures for Nonlinear Fluid Simulations of Plasmas E. Held, Utah State University (Principal Investigator) J.-Y. Ji, Utah State University (Co-Investigator) Utah State University (USU) continued to develop and apply closures for hybrid fluid/kinetic simulations of magnetized fusion plasmas. The closures developed by USU's fusion theory and computation group use the drift ordering and are referred to as "continuum" and "general moment" approaches. They differ in their representation for the distribution functions which probabilistically describe the density of particles in 5 dimensional phase space. The two chosen velocity variables, needed in studies of kinetic physics for magnetized plasmas, are the speed and the component of a particle's velocity along the magnetic field. In magnetized plasmas, the charged particle trajectories are tied closely to the magnetic field lines. Both approaches account for binary interactions between particles using the Coulomb collision operator. Theoretical and computational development along both lines provided avenues for verification between the two approaches. This funding allowed USU to further develop physically accurate and numerically efficient models that advance our understanding of kinetic physics in experiments like the International Thermonuclear Experimental Reactor (ITER). The ITER device will likely prove the physical practicality of thermonuclear fusion by producing 500 megawatts of clean, abundant fusion power. This work helped to realize this goal by advancing our simulation capability using hybrid fluid/kinetic models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗