Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Configuration interaction”

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

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

At least 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Thermodynamic theory of crystal plasticity: Formulation and application to polycrystal fcc copper

We present a thermodynamic description of crystal plasticity. Our formulation is based on the Langer-Bouchbinder-Lookman thermodynamic dislocation theory (TDT), which asserts the fundamental importance of an effective temperature that describes the state of configurational disorder and therefore the dislocation density of the crystalline material. We extend the TDT description from isotropic plasticity to crystal plasticity with many slip systems. Finite-element simulations show favourable comparison with experiments on polycrystal fcc copper under uniaxial compression, tension, and simple shear. Here, the thermodynamic theory of crystal plasticity thus provides a thermodynamically consistent and physically rigorous description of dislocation motion in crystals. We also discuss new insights about the interaction of dislocations belonging to different slip systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Asymmetric equilibrium core structures of pyramidal-II < c + a > dislocations in ten hexagonal-close-packed metals

The structures of pyramidal-II < c + a > dislocations, one of the most important defects in structural hexagonal-close-packed (HCP) metals, have not been fully characterized for many of the HCP metals in use today. Here, we employ ab initio informed phase-field dislocation dynamics to determine the minimum energy structure of pyramidal {1¯1¯22} < 11¯23 > dislocations in ten HCP metals, including Be, Co, Mg, Re, Ti, Zn, Cd, Hf, Y, and Zr. As input for the simulations, we calculate, using first-principles density functional theory, the {1¯1¯22} generalized stacking fault energy (GSFE) curves for all ten metals. From these calculations, it is found that magnetism in Co is necessary for achieving a local minimum in the GSFE curve. We observe in simulations that edge and screw character dislocations split into two partials separated by a low-energy intrinsic stacking fault. The splitting distance is shown to scale inversely with the local minimum energy normalized by the product of its shear modulus and Burgers vector. Interestingly, some HCP metals exhibit an asymmetric structure, with either unequal partial Burgers vectors or widths, in contrast to the symmetric configuration expected from linear elastic dislocation theory. We explain these structures by properties of the local maxima in their GSFE curves. Metals with larger degrees of elastic anisotropy result in dislocations with larger splitting distances than would be expected under the commonly used assumption of elastic isotropy. Furthermore, these findings on the sizes and asymmetry in the structures of pyramidal-II < c + a > dislocations are fundamental to understanding how these dislocations glide and interact or react with other defects when these metals are mechanically strained.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Sensitivity of Wind Turbine Array Downstream Effects to the Parameterization Used in WRF

The Weather Research and Forecasting (WRF) Model has been extensively used for wind energy applications, and current releases include a scheme that can be applied to examine the effects of wind turbine arrays on the atmospheric flow and electricity generation from wind turbines. Herein we present a high-resolution simulation using two different wind farm parameterizations: 1) the “Fitch” parameterization that is included in WRF releases and 2) the recently developed Explicit Wake Parameterization (EWP) scheme. We compare the schemes using a single yearlong simulation for a domain centered on the highest density of current turbine deployments in the contiguous United States (Iowa). Pairwise analyses are applied to diagnose the downstream wake effects and impact of wind turbine arrays on near-surface climate conditions. On average, use of the EWP scheme results in small-magnitude wake effects within wind farm arrays and faster recovery of full WT array wakes. This in turn leads to smaller impacts on near-surface climate variables and reduced array–array interactions, which at a systemwide scale lead to summertime capacity factors (i.e., the electrical power produced relative to nameplate installed capacity) that are 2%–3% higher than those from the more commonly applied Fitch parameterization. It is currently not possible to make recommendations with regard to which wind farm parameterization exhibits higher fidelity or to draw inferences with regard to whether the relative performance may vary with prevailing climate conditions and/or wind turbine deployment configuration. However, the sensitivities documented herein to the wind farm parameterization are of sufficient magnitude to potentially influence wind turbine array siting decisions. Thus, our research findings imply high value in undertaking combined long-term high-fidelity observational studies in support of model validation and verification.

17 WIND ENERGY↗
DOE OSTI2021

Mechanistic understanding of support effect on the activity and selectivity of indium oxide catalysts for CO 2 hydrogenation

Herein we present a mechanistic study on the support effect (ZrO 2 and CeO 2 ) of In 2 O 3 catalysts in CO 2 hydrogenation by a combined experimental and computational approach. Kinetic experiments and surface characterization suggested that the activity of In 2 O 3 catalysts cannot be simply correlated with the abundance of surface oxygen vacancies (O v ) formed by either H 2 -reduction or thermal treatment, which has been frequently invoked in previous studies. The support effect should originate from the electronic interactions between In 2 O 3 and the support oxide, rather than geometric factors or the difference in the particle size of In 2 O 3 . Theoretical modelling revealed that surface O v facilitate the formation and stabilization of the formate (HCOO*) intermediate. While a carbonate-like structure is favored for CO 2 adsorption on CeO 2 -supported or unsupported In 2 O 3 catalysts, CO 2 tends to bind strongly in a bent configuration on the O v site at the In 2 O 3 -ZrO 2 interface. The distinct CO 2 adsorption structures on different supported In 2 O 3 catalysts may account for the different reaction energy profiles in the subsequent hydrogenation reactions, especially the rate-limiting step, i.e., hydrogenation of HCOO* to CH 2 O* and methoxy (CH 3 O*). The relatively higher methanol selectivity of In 2 O 3 catalyst supported on ZrO 2 with respect to that on CeO 2 are suggested to stem from the greater energy difference (Δ$E_a$) between the parallel hydrogenation and C-O bond cleavage of HCOO*, which leads to the formation of methanol and CO, respectively. We report this study underlines the important role of metal-oxide-interface in determining the catalytic behavior of oxide-supported In 2 O 3 catalysts in CO 2 conversion.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Machine learning based algorithms for uncertainty quantification in numerical weather prediction models

Complex numerical weather prediction models incorporate a variety of physical processes, each described by multiple alternative physical schemes with specific parameters. The selection of the physical schemes and the choice of the corresponding physical parameters during model configuration can significantly impact the accuracy of model forecasts. There is no combination of physical schemes that works best for all times, at all locations, and under all conditions. It is therefore of considerable interest to understand the interplay between the choice of physics and the accuracy of the resulting forecasts under different conditions. This paper demonstrates the use of machine learning techniques to study the uncertainty in numerical weather prediction models due to the interaction of multiple physical processes. The first problem addressed herein is the estimation of systematic model errors in output quantities of interest at future times, and the use of this information to improve the model forecasts. The second problem considered is the identification of those specific physical processes that contribute most to the forecast uncertainty in the quantity of interest under specified meteorological conditions. In order to address these questions we employ two machine learning approaches, random forests and artificial neural networks. The discrepancies between model results and observations at past times are used to learn the relationships between the choice of physical processes and the resulting forecast errors. Numerical experiments are carried out with the Weather Research and Forecasting (WRF) model. The output quantity of interest is the model precipitation, a variable that is both extremely important and very challenging to forecast. The physical processes under consideration include various micro-physics schemes, cumulus parameterizations, short wave, and long wave radiation schemes. The experiments demonstrate the strong potential of machine learning approaches to aid the study of model errors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Machine learning force field model for kinetic Monte Carlo simulations of itinerant Ising magnets

Here, we present a scalable machine learning (ML) framework for large-scale kinetic Monte Carlo (kMC) simulations of itinerant electron Ising systems. As the effective interactions between Ising spins in such itinerant magnets are mediated by conducting electrons, the calculation of energy change due to a local spin update requires solving an electronic structure problem. Such repeated electronic structure calculations could be overwhelmingly prohibitive for large systems. Assuming the locality principle, a convolutional neural network (CNN) model is developed to directly predict the effective local field and the corresponding energy change associated with a given spin update based on Ising configuration in a finite neighborhood. As the kernel size of the CNN is fixed at a constant, the model can be directly scalable to kMC simulations of large lattices. Our approach is reminiscent of the ML force field models widely used in first-principles molecular dynamics simulations. Applying our ML framework to a square-lattice double-exchange Ising model, we uncover unusual coarsening of ferromagnetic domains at low temperatures. Our work highlights the potential of ML methods for large-scale modeling of similar itinerant systems with discrete dynamical variables.

machine learning↗
DOE OSTI2025

Symmetry breaking forms split-off flat bands in quantum oxides controlling metal versus insulator phases

The crystal structure used as input to electronic structure calculation of conventionally bonded solids consists primarily of the standard crystallographic degrees of freedom. Quantum solids sometimes have additional microscopic degrees of freedom (m-DOF) nested within the crystallographic structure. These might include local motifs including positional (Peierls dimers, deformed octahedra, Jahn-Teller distortions), magnetic (local moment configurations) and dipolar (local ferroelectric configurations). Such motifs can be observed experimentally via local probes that avoid averaging, and theoretically as distinct total energy lowering features relative to more simplified “average crystallographic structures”. Here we examine the ability of electronic structure methods independent of strong correlation physics to explain the broad phenomenology of metal-insulator selectivity in quantum oxides by energy-lowering symmetry breaking. We do this by avoiding the restriction of considering only (i) electron-electron interactions in a fixed unresponsive lattice—as done in Mott strong correlation explanations—allowing, however, (ii) local positional and magnetic motifs that coexist within a crystallographic landscape. We find in a broad range of quantum oxides with different symmetry breaking modes the formation of split-off flat bands that can also control metallic versus insulating characteristics. The split-off band effect is common to magnetic and non-magnetic materials, including both binary and ternary oxides. One finds that when such local symmetry breaking motifs are considered in mean-field-like (e.g. density functional theory) electronic structure calculations of quantum oxides, they explain many of the observed trends in metal vs insulator phases discussed otherwise in prior literature via strong correlation in symmetry unbroken view.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Report on Initial Tribological Studies of Graphite in Dry Argon and Molten Salt Environment

This report documents completion of the Advanced Reactor Technologies (ART) Level 2 Milestone M2TG-24OR0501081: “Complete report on initial tribological studies within molten salt environment”, due January 31, 2024. This milestone is a carryover milestone from FY23 M2TG-23OR0501091. The report summarizes the initial studies on tribological characterization of the graphite–graphite sliding interface in a dry argon and molten FLiNaK salt environment. The experiments were conducted on a high-temperature multifunctional tribometer placed in a glovebox. This configuration enables probing the wear and friction properties at high temperatures in an inert environment. The worn surfaces were analyzed using high-resolution microscopy and white-light profilometry to determine the wear modes and wear rates. This report discusses the previous wear and friction studies on 316H stainless steel (SS) sliding against graphite in molten salt for different temperatures, sliding speeds, and salt amount. Moreover, key tribological parameters of pebble-on-pebble interaction in gas-cooled and molten salt reactors (MSRs) are analyzed and used to design an experimental test matrix for future studies. The outcomes of this project could provide key information about the tribological behavior of graphite pebbles in gas-cooled and MSRs, thereby contributing to safer operation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

The proton radius (puzzle?) and its relatives

We review determinations of the electric proton charge radius from a diverse set of low-energy observables. We explore under which conditions it can be related to Wilson coefficients of appropriate effective field theories. This discussion is generalized to other low-energy constants. This provides us with a unified framework to deal with a set of low-energy constants of the proton associated with its electromagnetic interactions. Unambiguous definitions of these objects are given, as well as their relation with expectation values of QCD operators. We show that the proton radius obtained from spectroscopy and lepton–proton scattering (when both the lepton and proton move with nonrelativistic velocities) is related to the same object of the underlying field theory with O ( α ) precision. The model dependence of these analyses is discussed. The prospects of constructing effective field theories valid for the kinematic configuration of present, or near-future, lepton–proton scattering experiments are explored.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

In situ characterization of laser-generated melt pools using synchronized ultrasound and high-speed X-ray imaging

Metal additive manufacturing is a fabrication method that forms a part by fusing layers of powder to one another. An energy source, such as a laser, is commonly used to heat the metal powder sufficiently to cause a molten pool to form, which is known as the melt pool. The melt pool can exist in the conduction or the keyhole mode where the material begins to rapidly evaporate. The interaction between the laser and the material is physically complex and difficult to predict or measure. Here, high-speed X-ray imaging was combined with immersion ultrasound to obtain synchronized measurements of stationary laser-generated melt pools. Furthermore, two-dimensional and three-dimensional finite-element simulations were conducted to help explain the ultrasonic response in the experiments. In particular, the time-of-flight and amplitude in pulse-echo configuration were observed to have a linear relationship to the depth of the melt pool. These results are promising for the use of ultrasound to characterize the melt pool behavior and for finite-element simulations to aid in interpretation.

42 ENGINEERING↗
DOE OSTI2025

TRUST Sensors in Environments: Thermocouples (SE-TC), Release FY25

The Delivery Environments Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) project is a broad project intended to analyze simplified problems experimentally and with modeling and simulation. The purpose of analyzing these simplified problems is to extend solution methods to more complex problems, as well as understand deficiencies and gaps in knowledge of methods currently used in more complex analyses. The TRUST project encompasses several smaller testbeds intended to isolate individual phenomena. The testbed under consideration in this report is the Sensors in Environments: Thermocouples testbed. In previous years, the purpose of this testbed was to quantify uncertainty of thermocouple sensors. To accomplish this, an aluminum plate was placed in a thermal chamber and subject to various types of thermal loading. Thermocouples were placed in various locations on the aluminum plate in various configurations (e.g., embedded in the plate, placed under Kapton tape), and an effort was made to quantify uncertainty in these measurements. Finite element simulations were performed to investigate how sensitive these measurements were to parameters such as the boundary conditions on the plate and material properties. However, a fundamental source of uncertainty in this analysis was the convective heat transfer from the plate. Convective heat transfer is a complex physical phenomenon comprised of a number of interacting sub-processes and is difficult to predict accurately a priori. As such, the main purpose of this testbed in FY25 was to better understand, both experimentally and numerically, the convective heat transfer from the plate. This is a highly applicable problem to several more complex problems, as convective heat transfer occurs in nearly all problems where a body is moving through air. Numerically, this required a two-step approach. First, the air flow in the thermal chamber was in vestigated using computational fluid dynamics. The commercial solver Fluent was used to perform these simulations. From these simulations, a heat transfer coefficient over the surface of the plate was calculated. This heat transfer was then used as boundary conditions for finite element heat transfer simulations within the plate, which were performed using Abaqus. Significant effort was devoted to automating the handoff between these two solvers. Experimentally, previous thermocouple results in the plate were used to validate the time-dependent thermal profiles produced from Abaqus. Further experimental efforts were performed both to help validate the Fluent simulations and to inform its boundary conditions. For example, hot-wire anemometers were used to measure the velocity in the chamber, which would be particularly useful in understanding the chamber inlet velocity. Thermocouple measurements were also taken in the chamber, instead of only on the plate, to serve as validation evidence for the Fluent simulations. Numerical results showed that the Fluent to Abaqus workflow matched previous plate thermocouple measurements well. This type of handoff is useful for more complex experiments, or those that are not able to be examined in as great of detail as this testbed, as it was performed without any experimental input. Experimental results, however, were more mixed. The anemometers proved unreliable, with inconsistent measurements across all anemometers, even at locations that were nearly identical. On the other hand, the thermocouples provided a relatively rich view of the temperature field in the chamber.

42 ENGINEERING↗
DOE OSTI2024

Collaborative Research Proposal: Time Resolved Optical Emissions Spectroscopy and Laser Induced Fluorescence Spectroscopy of Nanosecond Pulsed Discharges in a Gas-Liquid Water Film Reactor

Electrical discharge plasma formed in contact with liquid water is of interest for a wide range of applications in chemical, biomedical, agricultural, electrical, and materials science and engineering. Such plasma reactors are of very timely importance since they also have significant disinfection capability by inactivating bacteria, viruses and other pathogens. Many types of plasma sources including those driven by AC, DC, RF, microwave, and pulsed electrical power supplies coupled to a wide range of different electrode configurations and reactor designs have been developed and explored which contact the plasma with liquid water. Recent roadmaps have recommended that further work is needed to develop our understanding of the fundamental chemical and physical process which occur at the interface of non-thermal plasma with liquid water solutions in order to advance this large diversity of applications which ultimately depend upon efficient production of key reactive chemical species. There is a wide range of interacting factors that affect the chemical reactions that occur in the plasma, in the liquid phase, and at the interface. These factors ultimately govern the key reactive chemical species formed and used in the various applications and they include a) the reactor design and input parameters, b) the discharge and transport processes, c) the plasma properties, and d) the resulting chemical reactions. For pulsed discharges, the power supply design and output parameters control the applied voltage, frequency, rise time, and width (duration) of the applied pulses. The reactor design involves specification of the gas-liquid contacting methods, electrode gap distance, reactor volume and shape, and gas and liquid flow rates achievable. The gas and liquid compositions as well as the liquid properties such as pH and conductivity are also of key importance in determination of the resulting chemical reactions. In addition, the important plasma properties include plasma gas temperature, electron density, electron energy (distribution), and size of the plasma channels which are all affected by the discharge and transport properties. Many studies have focused on specific aspects of these various processes. Of particular important and relevance to the proposed work is the utilization of nanosecond pulses to generate plasma in gas-liquid systems, liquid bubbles, underwater, and in gases. Recent advances in nanosecond pulses provide significant advantages in utilization with liquid water. For example, fast rise time and short pulses are less sensitive to water conductivity, such short pulses may provide advantages in fast temporal quenching of the plasma, and fundamental analysis of pulse properties, including pulse shape and width, may be facilitated by investigation of single filamentary fast pulses. Many studies have also dealt with the role of the gas composition on formation of reactive oxygen species (ROS) (i.e., hydroxyl radicals – ·OH, hydrogen peroxide – H 2 O 2 , various atomic oxygen species, ozone-O 3 , hydroperoxyl radicals HO 2 ·) and reactive nitrogen species (RNS) (i.e., nitrogen oxides – NO, NO 2 , N 2 O - collectively termed NO X , nitrite-NO 2 - , nitrate-NO 3 - , peroxynitrite-ONOO - ). The hydroxyl radical is the critical species in many chemical oxidation reactions for chemical degradation of toxic compounds in water and gases and for synthesis of some compounds. The mixture of various nitrogen oxide species is important for many biochemical and biological processes involved in biomedical and agricultural applications including disinfection and fertilizer production. The present proposal focuses on determination of the effects of time resolved electron density and hydroxyl radicals on plasma chemical reactions through collaboration with the Princeton Collaborative Low Temperature Plasma Research Facility (PCRF) at the Princeton Plasma Physics Laboratory (PPPL). In order to further investigate the role of the plasma generated electrons and hydroxyl radicals on the overall formation of the key species including hydrogen peroxide, hydroxyl radicals, and nitrogen oxides, the proposed work, thus seeks to determine high resolution time resolved electron density by optical emissions spectroscopy and time resolved hydroxyl radicals using laser induced fluorescence in the nanosecond discharge reactor. The combination of data on electron density and hydroxyl radical concentration will be utilized in the present work to more fully characterize the chemical reaction processes in this system and to advance the design, development, and operation of such chemical reactors for a wide range of applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

General Analysis of Data Collected from DOE Lidar Buoy Deployments Off Virginia and New Jersey

Pacific Northwest National Laboratory (PNNL) operates two AXYS WindSentinel lidar buoys for the U.S. Department of Energy’s Wind Energy Technologies. The purpose of these buoys is to collect hub-height winds and supporting meteorological and oceanographic information to facilitate the development of offshore wind energy in the U.S. In general, each buoy is deployed for a year or more at a given location in order to capture at least a full annual cycle of weather conditions. The initial deployment for one buoy was off the coast of Virginia beginning in 2015, and the other buoy was first deployed off the coast of New Jersey beginning in 2016. Over the last two years, PNNL has had an opportunity to analyze the data collected during these first two deployments. This report describes a substantial analysis of data collected by the two lidar buoys operated off the coasts of Virginia and New Jersey. The centerpiece instrument for each buoy as deployed off Virginia and New Jersey is a lidar system, which is designed to measure the horizontal wind vector from approximately 40 m to 200 m above the sea surface with vertical resolution of 40 m. Since the initial deployments, the original lidars have been replaced with more powerful Leosphere 866 v2 systems. The analyses in this report will apply to the original Vindicator systems. In addition to the wind profiles from the lidars, the buoys collect near-surface measurements of wind speed and direction, air temperature, relative humidity, barometric pressure, and solar irradiance. Oceanographic variables measured include the two-dimensional wave spectrum, water temperature and conductivity, and ocean current vectors to a depth of 90 m. An assessment of overall data recovery and a basic analysis of the data collected was provided in a previous report This report substantially extends that analysis. The various sections describe the development and application of an inertial measurement unit (IMU) data recovery scheme for the New Jersey deployment; a climatological analysis of winds at hub height and at the surface together with thermodynamic variables measured at the surface for the full deployment periods; an analysis of oceanographic observations describe sea state; the development of a refinement for NOAA’s WaveWatch III model to allow its application to near-shore areas; a basic climatology of ocean currents observed from the buoys; an analysis of observed winds in the framework of Monin-Obukhov Similarity Theory; and the development and evaluation of techniques to extract turbulence intensity and turbulence kinetic energy from the lidars. The analyses contained in this report provide a great deal of new information about offshore conditions on the U.S. East Coast. In addition, the experience gained will inform both configurations and analysis of data from future deployments of these lidar buoy systems.

17 WIND ENERGY↗
DOE OSTI2020

High broadband photoconductivity of few-layered MoS 2 field-effect transistors measured using multi-terminal methods: effects of contact resistance

Among the layered two dimensional semiconductors, molybdenum disulfide (MoS 2 ) is considered to be an excellent candidate for applications in optoelectronics and integrated circuits due to its layer-dependent tunable bandgap in the visible region, high ON/OFF current ratio in field-effect transistors (FET) and strong light–matter interaction properties. In this study, using multi-terminal measurements, we report high broadband photocurrent response ( R ) and external quantum efficiency (EQE) of few-atomic layered MoS 2 phototransistors fabricated on a SiO 2 dielectric substrate and encapsulated with a thin transparent polymer film of Cytop. The photocurrent response was measured using a white light source as well as a monochromatic light of wavelength λ = 400 nm–900 nm. We measured responsivity using a 2-terminal configuration as high as R = 1 × 10 3 A W -1 under white light illumination with an optical power P opt = 0.02 nW. The R value increased to 3.5 × 10 3 A W -1 when measured using a 4-terminal configuration. Using monochromatic light on the same device, the measured values of R were 10 3 and 6 × 10 3 A W -1 under illumination of λ = 400 nm when measured using 2- and 4-terminal methods, respectively. The highest EQE values obtained using λ = 400 nm were 10 5 % and 10 6 % measured using 2- and 4-terminal configurations, respectively. The wavelength dependent responsivity decreased from 400 nm to the near-IR region at 900 nm. The observed photoresponse, photocurrent–dark current ratio (PDCR), detectivity as a function of applied gate voltage, optical power, contact resistances and wavelength were measured and are discussed in detail. The observed responsivity is also thoroughly studied as a function of contact resistance of the device.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ecological Assembly Processes Are Coordinated between Bacterial and Viral Communities in Fractured Shale Ecosystems

The ecological drivers that concurrently act upon both a virus and its host and that drive community assembly are poorly understood despite known interactions between viral populations and their microbial hosts. Hydraulically fractured shale environments provide access to a closed ecosystem in the deep subsurface where constrained microbial and viral community assembly processes can be examined. Here, we used metagenomic analyses of time-resolved-produced fluid samples from two wells in the Appalachian Basin to track viral and host dynamics and to investigate community assembly processes. Hypersaline conditions within these ecosystems should drive microbial community structure to a similar configuration through time in response to common osmotic stress. However, viral predation appears to counterbalance this potentially strong homogeneous selection and pushes the microbial community toward undominated assembly. In comparison, while the viral community was also influenced by substantial undominated processes, it assembled, in part, due to homogeneous selection. When the overall assembly processes acting upon both these communities were directly compared with each other, a significant relationship was revealed, suggesting an association between microbial and viral community development despite differing selective pressures. These results reveal a potentially important balance of ecological dynamics that must be in maintained within this deep subsurface ecosystem in order for the microbial community to persist over extended time periods. More broadly, this relationship begins to provide knowledge underlying metacommunity development across trophic levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Global particle buildup simulations with gas puff scan: application to WEST discharge

This paper deals with the distribution of sources, transport, and exhaust of particles in a tokamak. Knowledge and understanding of all the physical phenomena involved in the global particle buildup are necessary to study and predict density regimes and subsequently to develop optimized scenarios for tokamak operation in order to control heat and particle exhaust. Neutral particles and their interactions with plasma are central in this perspective. This paper discusses the impact of varying the intensity of particle fueling in 2D transport simulations of a WEST discharge. Simulations are performed with an updated version of SOLEDGE-HDG that allows a more realistic transport of neutrals using a self-consistent diffusive model based on charge exchange and ionization processes. New code capabilities allow the entire WEST poloidal cross section to be simulated in a realistic configuration for both geometry and the range of control parameters. A gas puff scan illustrates the main features of the sheath-limited, high-recycling, and detached regimes, such as the buildup of the temperature gradient and the pressure drop in the scrape-off layer (SOL), the target temperature falling to 1 eV, and the ionization source moving away from the targets, as well as the particle flux rollover. A crude estimate of wall erosion is also provided, showing the respective role of each plasma wall component in each of these regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Thermodynamics and its prediction and CALPHAD modeling: Review, state of the art, and perspectives

Thermodynamics is a science concerning the state of a system, whether it is stable, metastable, or unstable, when interacting with its surroundings. The combined law of thermodynamics derived by Gibbs about 150 years ago laid the foundation of thermodynamics. In Gibbs combined law, the entropy production due to internal processes was not included, and the 2nd law was thus practically removed from the Gibbs combined law, so it is only applicable to systems under equilibrium, thus commonly termed as equilibrium or Gibbs thermodynamics. Gibbs further derived the classical statistical thermodynamics in terms of the probability of configurations in a system in the later 1800's and early 1900's. With the quantum mechanics (QM) developed in 1920's, the QM-based statistical thermodynamics was established and connected to classical statistical thermodynamics at the classical limit as shown by Landau in the 1940's. In 1960's the development of density functional theory (DFT) by Kohn and co-workers enabled the QM prediction of properties of the ground state of a system. On the other hand, the entropy production due to internal processes in non-equilibrium systems was studied separately by Onsager in 1930's and Prigogine and co-workers in the 1950's. In 1960's to 1970's the digitization of thermodynamics was developed by Kaufman in the framework of the CALculation of PHAse Diagrams (CALPHAD) modeling of individual phases with internal degrees of freedom. CALPHAD modeling of thermodynamics and atomic transport properties has enabled computational design of complex materials in the last 50 years. Our recently termed zentropy theory integrates DFT and statistical mechanics through the replacement of the internal energy of each individual configuration by its DFT-predicted free energy. The zentropy theory is capable of accurately predicting the free energy of individual phases, transition temperatures and properties of magnetic and ferroelectric materials with free energies of individual configurations solely from DFT-based calculations and without fitting parameters, and is being tested for other phenomena including superconductivity, quantum criticality, and black holes. Those predictions include the singularity at critical points with divergence of physical properties, negative thermal expansion, and the strongly correlated physics. Furthermore, those individual configurations may thus be considered as the genomic building blocks of individual phases in the spirit of the materials genome®. This has the potential to shift the paradigm of CALPHAD modeling from being heavily dependent on experimental inputs to becoming fully predictive with inputs solely from DFT-based calculations and machine learning models built on those calculations and existing experimental data through newly developed and future open-source tools. Furthermore, through the combined law of thermodynamics including the internal entropy production, it is shown that the kinetic coefficient matrix of independent internal processes is diagonal with respect to the conjugate potentials in the combined law, and the cross phenomena that the phenomenological Onsager flux and reciprocal relationships are due to the dependence of the conjugate potential of a molar quantity on nonconjugate molar quantities and other potentials, which can be predicted by the zentropy theory and CALPHAD modeling.

42 ENGINEERING↗
DOE OSTI2024

Simulations of vertical displacement oscillatory modes and global Alfvén Eigenmodes in JET geometry

Vertical Displacement Oscillatory Modes (VDOM), with frequency in the Alfvén range, are natural modes of oscillation of magnetically confined laboratory plasmas with elongated cross-section. These axisymmetric modes arise from the interaction between the plasma current, which is in equilibrium with currents flowing in external coils, and perturbed currents induced on a nearby conducting wall. The restoring force exerted by these perturbed currents on the vertical motion of the plasma column leads to its oscillatory behavior. An analytic model for VDOM was proposed based on an idealized 'straight tokamak' equilibrium with uniform equilibrium current density. This article introduces the first numerical simulations of VDOM in a realistic JET tokamak configuration, using the extended-MHD code NIMROD and drawing comparisons with Global Alfvén Eigenmodes (GAE). The results show qualitative agreement with analytic predictions regarding mode frequency and radial structure, supporting the identification of VDOM as a fundamental oscillation mode in tokamak plasmas. VDOM and GAE are modeled in a representative JET discharge, where axisymmetric perturbations with toroidal mode number n = 0 driven unstable by fast ions were observed. The two modes are examined separately using a forced oscillator within the NIMROD code, which enables a comparison of their characteristics and helps identify the experimentally observed mode possibly as a GAE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗