Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective diffusivity”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Equilibrium solute segregation to matrix- θ' precipitate interfaces in Al-Cu alloys from first principles

Particular combinations of solute atoms segregated to the interface of the Al and θ ' - Al 2 Cu in Al-Cu alloys can help stabilize θ ' precipitates at high temperatures. Stabilization of such precipitates is determined by a combination of thermodynamics (including driving forces for coarsening and transformation and solute segregation tendencies) and kinetic effects (including solute diffusion and interfacial mobility in the presence of interfacial solute segregation). For some alloys such as recent Al-Cu-based alloys, multiple solutes segregate in significant quantities to interfaces, and solute-solute interactions at the interface are important, with multiple types of solutes competing for similar interfacial sites. To treat this situation, we develop and apply a statistical mechanics approach to calculate the temperature-dependent equilibrium solute atoms distribution near the coherent and semicoherent interfaces between the Al matrix and the θ ' - Al 2 Cu precipitates. The developed approach is applied to the investigation of Si, Mn, and Zr segregation at the interface, as particular combinations of these elements affect the thermal stability of the θ ' precipitates. We demonstrate that because Si and Mn atoms segregate on the same semicoherent interface, the presence of Si reduces the concentration of Mn solutes at the interface. Si atoms preferably occupy the first layer of the interface and compete with Mn atoms for one type of particular sites in the layer. Mn atoms preferably occupy the second layer of the semicoherent interface, and the Mn-Mn interaction plays an important role in their distribution. Zr atoms mostly segregate on one of the two nonequivalent sites of the second layer of the coherent interface. Due to symmetry properties of the coherent interface, the calculations show that the segregation Zr of atoms to this interface will likely lead to the formation of L 1 2 ordered Al 3 Zr layer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Plants mediate soil organic matter decomposition in response to sea level rise

Abstract Tidal marshes have a large capacity for producing and storing organic matter, making their role in the global carbon budget disproportionate to land area. Most of the organic matter stored in these systems is in soils where it contributes 2–5 times more to surface accretion than an equal mass of minerals. Soil organic matter ( SOM ) sequestration is the primary process by which tidal marshes become perched high in the tidal frame, decreasing their vulnerability to accelerated relative sea level rise ( RSLR ). Plant growth responses to RSLR are well understood and represented in century‐scale forecast models of soil surface elevation change. We understand far less about the response of SOM decomposition to accelerated RSLR . Here we quantified the effects of flooding depth and duration on SOM decomposition by exposing planted and unplanted field‐based mesocosms to experimentally manipulated relative sea level over two consecutive growing seasons. SOM decomposition was quantified as CO 2 efflux, with plant‐ and SOM ‐derived CO 2 separated via δ 13 CO 2 . Despite the dominant paradigm that decomposition rates are inversely related to flooding, SOM decomposition in the absence of plants was not sensitive to flooding depth and duration. The presence of plants had a dramatic effect on SOM decomposition, increasing SOM ‐derived CO 2 flux by up to 267% and 125% (in 2012 and 2013, respectively) compared to unplanted controls in the two growing seasons. Furthermore, plant stimulation of SOM decomposition was strongly and positively related to plant biomass and in particular aboveground biomass. We conclude that SOM decomposition rates are not directly driven by relative sea level and its effect on oxygen diffusion through soil, but indirectly by plant responses to relative sea level. If this result applies more generally to tidal wetlands, it has important implications for models of SOM accumulation and surface elevation change in response to accelerated RSLR .

Mueller, Peter↗
DOE OSTI2019

Electronic Structures and Magnetic Properties of Transition Metal Doped CsPbI{sub 3} Perovskite Compounds by First-Principles Calculation

Transition metal doped cesium lead halide (CsPbI{sub 3}) perovskite compounds were studied for application in photovoltaic solar cells. Electronic structures, chemical shifts of {sup 207}Pb and {sup 127}I-NMR, vibration modes in infrared and Raman spectra of transition metals (Mn{sup 2+}, Fe{sup 2+} or Cu{sup 2+})-doped CsPbI{sub 3} perovskite compounds were studied by the first-principles calculation using density functional theory. The CsPb(Fe)I{sub 3} perovskite crystals had a slight perturbation of crystal field in the coordination structure. The electron density distribution was delocalized on the 5p orbital of I atom, the 3d orbital of Fe atom and the 6p orbital of Pb atom. The first excited process was based on ligand metal charge transfer from the 5p orbital on I atom to the 3d orbital of Fe atom. The chemical shifts of {sup 127}I-NMR were associated with the electron correlation of electron-nuclear spin interaction and nuclear quadrupole interactions based on electron field graduate. The asymmetric vibrations of Pb–I bonds stretching mode related to electron conductivity with scattering of the carrier diffusion as phonon effectiveness. The slight perturbation of the coordination structure in the CsPb(Fe)I{sub 3} perovskite crystal will improve the photovoltaic and optical properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Modeling the Mesoscale Transport of Lithium-Magnetite Electrodes Using Insight from Discharge and Voltage Recovery Experiments

A multi-scale mathematical model, which accounts for mass transport on the crystal and agglomerate length-scales, is used to investigate the electrochemical performance of lithium-magnetite electrochemical cells. Experimental discharge and voltage recovery data are compared to three sets of simulations, which incorporate crystal-only, agglomerate-only, or multi-scale transport effects. Mass transport diffusion coefficients are determined by fitting the simulated voltage recovery times to experimental data. In addition, a further extension of the multi-scale model is proposed which accounts for the impact of agglomerate size distributions on electrochemical performance. The results of the study indicate that, depending on the crystal size, the low utilization of the active material is caused by transport limitations on the agglomerate and/or crystal length-scales. For electrodes composed of small crystals (6 and 8 nm diameters), it is concluded that the transport limitations in the agglomerate are primarily responsible for the long voltage recovery times and low utilization of the active mass. In the electrodes composed of large crystals (32 nm diameter), the slow voltage recovery is attributed to transport limitations on both the agglomerate and crystal length-scales.

25 ENERGY STORAGE↗
DOE OSTI2018

PLUMECALC, Version 2.3.x

PLUMECALC is a post processor to random walk particle tracking simulations of subsurface solute transport. PLUMECALC takes results of a single particle tracking run without retention processes and incorporates the effects of matrix diffusion or other retention processes for a user-specified solute source.

Dash, Zora↗
DOE OSTI2017

Simulations and Data analysis for the 35 ton Liquid Argon detector as a prototype for the DUNE experiment

The Deep Underground Neutrino Experiment (DUNE) is a next-generation neutrino experiment which will be built at the Sanford Underground Research Facility (SURF), and will receive a wide-band neutrino beam from Fermilab, 1300~km away. At this baseline DUNE will be able to study many of the properties of neutrino mixing, including the neutrino mass hierarchy and the value of the CP-violating complex phase ($\delta_{CP}$). DUNE will utilise Liquid Argon (LAr) Time Projection Chamber (TPC) (LArTPC) technology, and the Far Detector (FD) will consist of four modules, each containing 17.1~kt of LAr with a fiducial mass of around 10~kt. Each of these FD modules represents around an order of magnitude increase in size, when compared to existing LArTPC experiments. \\ The 35 ton detector is the first DUNE prototype for the single (LAr) phase design of the FD. There were two running periods, one from November 2013 to February 2014, and a second from November 2015 to March 2016. During t he second running period, a system of TPCs was installed, and cosmic-ray data were collected. A method of particle identification was developed using simulations, though this was not applied to the data due to the higher than expected noise level. A new method of determining the interaction time of a track, using the effects of longitudinal diffusion, was developed using the cosmic-ray data. A camera system was also installed in the detector for monitoring purposes, and to look for high voltage breakdowns. \\ Simulations concerning the muon-induced background rate to nucleon decay are performed, following the incorporation of the MUon Simulations UNderground (MUSUN) generator into the DUNE software framework. A series of cuts which are based on Monte Carlo truth information is developed, designed to reject simulated background events, whilst preserving simulated signal events in the $n \rightarrow K^{+} + e^{-}$ decay channel. No background events are seen to survive the app lication of these cuts in a sample of 2~$\times$~10$^9$ muon! s, representing 401.6~years of detector live time. This corresponds to an annual background rate of <~0.44~events$\cdot$Mt$^{-1}\cdot$year$^{-1}$ at 90\% confidence, using a fiducial mass of 13.8~kt.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Predictive Large-Eddy Simulation of Supercritical-Pressure Reactive Flows in the Cold Ignition Regime

This report describes a one-year study performed under DOE sponsorship, continuing the investigation of high-pressure turbulent reactive flows. The interest is in the effect of the chemical species distribution in high-pressure turbulent flows in the presence of strong temperature gradients as would occur during reactions in realistic flows. The prime example of such flows are boundary layers in which the wall is at a lower temperature than that of the fluid, as would be the case in Diesel engines. Previous DOE BES work in the program further highlighted the importance of the boundary layer as a configuration for fundamental studies: soot formation in boundary layers is still a problem poorly understood and depends on the availability of particular chemical species at that location, and when boundary layer Large Eddy Simulations results were compared with measurements, the agreement was unfavorable, showing that this important ‘unit’ problem is not well understood. To understand this unit problem, modeling and Direct Numerical Simulations of this unit problem were conducted for the simplest possible multispecies system, that is a binary-species system. The results discovered a new phenomenon, that is, Soret effect induced uphill diffusion. The far-reaching implication is that through an imposed wall/fluid temperature difference it is possible to control the distribution of the species in the boundary layer. Because the results have been documented in one paper published in the refereed literature, and also in conference papers, this final report is in the form of an Executive Summary succinctly describing the results and putting them in perspective with respect to existing information. The refereed and conference papers published are individually listed as Appendices and attached to this report. One manuscript is still in progress and is thus not listed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Crystal Growth in Gels from the Mechanisms of Crystal Growth to Control of Polymorphism: New Trends on Theoretical and Experimental Aspects

A gel can be considered to be a two-phase (liquid and solid) system, which lacks flow once it reaches a stationary state. The solid phase is usually a tridimensional polymeric mesh, while the liquid phase is usually found in three forms: contained in great cavities, retained in the capillary pores between micelles, or adsorbed on the surface of a micelle. The influence of the use of gels in crystal growth is diverse and depends on the type of gel being used. A decrease in solubility of any solute in the liquid may occur if the solvent interacts extensively with the polymeric section, hence, the nucleation in gels in these cases apparently occurs at relatively low supersaturations. However, if the pore size is small enough, there is a possibility that a higher supersaturation is needed, due to the compartmentalization of solvents. Finally, this may also represent an effect in the diffusion of substances. This review is divided into three main parts; the first evaluates the theory and practice used for the obtainment of polymorphs. The second part describes the use of gels into crystallogenesis of different substances. The last part is related to the particularities of protein crystal polymorphism, as well as modern trends in gel growth for high-resolution X-ray crystallography.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A PDE Model of Breast Tumor Progression in MMTV-PyMT Mice

The evolution of breast tumors greatly depends on the interaction network among different cell types, including immune cells and cancer cells in the tumor. This study takes advantage of newly collected rich spatio-temporal mouse data to develop a data-driven mathematical model of breast tumors that considers cells’ location and key interactions in the tumor. The results show that cancer cells have a minor presence in the area with the most overall immune cells, and the number of activated immune cells in the tumor is depleted over time when there is no influx of immune cells. Interestingly, in the case of the influx of immune cells, the highest concentrations of both T cells and cancer cells are in the boundary of the tumor, as we use the Robin boundary condition to model the influx of immune cells. In other words, the influx of immune cells causes a dominant outward advection for cancer cells. We also investigate the effect of cells’ diffusion and immune cells’ influx rates in the dynamics of cells in the tumor micro-environment. Sensitivity analyses indicate that cancer cells and adipocytes’ diffusion rates are the most sensitive parameters, followed by influx and diffusion rates of cytotoxic T cells, implying that targeting them is a possible treatment strategy for breast cancer.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

EXTENSION OF THE MURAM RADIATIVE MHD CODE FOR CORONAL SIMULATIONS

We present a new version of the MURaM radiative magnetohydrodynamics (MHD) code that allows for simulations spanning from the upper convection zone into the solar corona. We implement the relevant coronal physics in terms of optically thin radiative loss, field aligned heat conduction, and an equilibrium ionization equation of state. We artificially limit the coronal Alfvén and heat conduction speeds to computationally manageable values using an approximation to semi-relativistic MHD with an artificially reduced speed of light (Boris correction). We present example solutions ranging from quiet to active Sun in order to verify the validity of our approach. We quantify the role of numerical diffusivity for the effective coronal heating. We find that the (numerical) magnetic Prandtl number determines the ratio of resistive to viscous heating and that owing to the very large magnetic Prandtl number of the solar corona, heating is expected to happen predominantly through viscous dissipation. We find that reasonable solutions can be obtained with values of the reduced speed of light just marginally larger than the maximum sound speed. Overall this leads to a fully explicit code that can compute the time evolution of the solar corona in response to photospheric driving using numerical time steps not much smaller than 0.1 s. Numerical simulations of the coronal response to flux emergence covering a time span of a few days are well within reach using this approach.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Global Simulations of the Inner Regions of Protoplanetary Disks with Comprehensive Disk Microphysics

The gas dynamics of weakly ionized protoplanetary disks (PPDs) are largely governed by the coupling between gas and magnetic fields, described by three non-ideal magnetohydrodynamical (MHD) effects (Ohmic, Hall, ambipolar). Previous local simulations incorporating these processes have revealed that the inner regions of PPDs are largely laminar and accompanied by wind-driven accretion. We conduct 2D axisymmetric, fully global MHD simulations of these regions (∼1–20 au), taking into account all non-ideal MHD effects, with tabulated diffusion coefficients and approximate treatment of external ionization and heating. With the net vertical field aligned with disk rotation, the Hall-shear instability strongly amplifies horizontal magnetic field, making the overall dynamics dependent on initial field configuration. Following disk formation, the disk likely relaxes into an inner zone characterized by asymmetric field configuration across the midplane, which smoothly transitions to a more symmetric outer zone. Angular momentum transport is driven by both MHD winds and laminar Maxwell stress, with both accretion and decretion flows present at different heights, and modestly asymmetric winds from the two disk sides. With anti-aligned field polarity, weakly magnetized disks settle into an asymmetric field configuration with supersonic accretion flow concentrated at one side of the disk surface, and highly asymmetric winds between the two disk sides. In all cases, the wind is magneto-thermal in nature, characterized by a mass loss rate exceeding the accretion rate. More strongly magnetized disks give more symmetric field configuration and flow structures. Deeper far-UV penetration leads to stronger and less stable outflows. Implications for observations and planet formation are also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

BETO 2021 Peer Review - Cell Free and Immobilization Technologies (CFIT)

Today, several key factors negatively impact the production of fuels and chemicals from renewable sources. Common hindrances in the biological production of biochemicals are: (1) end-product or intermediate toxicity to the microbial biocatalyst, (2) the diversion of carbon to biomass formation, and (3) co-production of undesired byproducts. A particularly attractive alternative is to eliminate the biocatalyst entirely and instead operate the desired metabolic pathways in isolation, thus circumventing the roadblocks of biological toxicity, lower yields, and lack of specificity. However, cell-free enzyme systems still suffer from low productivities owing in part to the effects of free diffusion of intermediates, lack of long term enzyme stability, cofactor cost or inefficient recycling rates, and finally, the cost of enzyme production/purification. This project represents a new effort to propose innovative and cost competitive routes to producing biochemicals from a variety of feedstocks using cell free approaches. These routes will help reduce the current risk and cost associated with classical cell free production. Cell free technologies show promise for application to the production of toxic/inhibitory products or products difficult to separate from microbial growth media and can help reduce the production barriers in multiple areas of biological conversion of feedstocks to biochemicals.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2017

Exploring the Use of Muon Momentum for Detection of Nuclear Material Within Shielded Spent Nuclear Fuel Dry Casks

The inspection of the storage containers that house used nuclear fuel is an important issue facing the nuclear industry; currently, there are limited options available to provide for even minimal inspections. During the long storage life of the fuel, intermediate handling and transportation activities may take place that can create concerns regarding fuel integrity (safety) and continuity of knowledge (safeguards and security). Accidents and natural disasters also contribute to this concern. The reliable and safe management of used nuclear fuel is among the critical tasks to be addressed for the advancement of fission-based nuclear energy in the United States and elsewhere. Since the early 1950's, when the first nuclear power plant began to produce electricity, approximately 65,000 metric tons of used nuclear fuel have been generated in the United States. The postponement of the opening of a permanent repository until 2048 will lead to increased used fuel accumulation at reactor sites. A rapid increase in storage installations with more than 9,000 casks storing {approx}420,000 assemblies is anticipated by 2050 even if no new nuclear power plants are built. After used nuclear fuel has been placed inside a cask, the cask is sealed, not allowing for visual inspection. Conventional methods, e.g., x-rays, are limited in that they cannot penetrate dense, well-shielded objects. More sophisticated techniques, such as penetrating neutrons or proton radiography, require the use of an expensive accelerator. Efforts by Ziock et al. using neutron and gamma-ray imagers concluded that '...the signatures are deemed insufficient to uniquely identify individual casks of the same type...'. It was found that the gamma-ray and neutron images were dominated by diffusion and thermalization effects resulting from the large amount of shielding. One promising mode of inspection would use cosmic-ray muons. These muons are naturally generated, are able to penetrate dense materials, and are incident from all directions above the horizon. Tomography performed using these muons could potentially provide efficient, inexpensive, remote assessment of used nuclear fuel. However, current muon applications are limited in that they have not included momentum measurement. If momentum measurement could be incorporated, that would represent a big step forward because it would allow more accurate estimation of the muon observables, e.g., scattering variance. This approach has the potential to reconstruct the contents of a cask with improved resolution and identify partial defects, fuel assemblies, and individual fuel rods.

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

Prediction of divertor heat flux width for ITER using BOUT++ transport and turbulence module

Investigation of turbulent transport dynamics in scrape-off-layer (SOL) and divertor heat flux width prediction is performed for ITER. Both BOUT++ transport and BOUT++ turbulence codes are applied to capture the physics on different temporal scales. Simulations start with an ITER 15MA baseline scenario profile generated by CORSICA (Kim et al 2015 Paper ITER_D_R9T8J9 v1.1). In BOUT++ transport code, the plasma parameters ( n i , T i , T e ) and radial electric () profiles are evolved to steady state. Here, the initial plasma profiles inside the separatrix are taken from CORSICA scenario studies. Transport coefficients are calculated by inverting the plasma profiles inside the separatrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Reduced diffusion and enhanced retention of multiple radionuclides from pore structure characterization of barrier materials for enhanced repository performance

Fluid flow and chemical transport in porous media are the macroscopic consequences of pore structure, which integrates geometry (e.g., pore size and surface area, pore-size distribution) and topology (e.g., pore connectivity). Low-permeability geological media whose pores are poorly interconnected will exhibit the characteristics of anomalous diffusion and sample size-dependent effective porosity, which will strongly impact long-term net diffusion and retention of radionuclides in geological repository settings involving different host rocks and barrier materials. A suite of innovative and complementary experimental approaches is utilized to study the microscopic pore structure and macroscopic fluid flow & chemical transport for a range of host rocks and barrier materials, in addition to standard clay minerals and reference rocks. With a particular focus on quantifying the presence and magnitude of “isolated” pores for a reduced effective porosity in low-permeability geomedia, the integrated methodologies for basic properties and pore structure characterization of these geomedia include X-ray diffraction, thin section petrography, grain size distribution, water immersion porosimetry after vacuum-pulling for full saturation, mercury intrusion porosimetry, nitrogen physisorption, scanning electron microscopy, X-ray computed tomography, and (ultra-)small angle neutron (X-ray) scattering. In addition, custom-designed gas diffusion, tracer recipe involving a range of anionic and cationic chemicals with subsequent analyses by laser ablation and inductively coupled plasma-mass spectrometry, along with batch sorption, column transport, and imbibition tests were conducted for coupled effects of pore structure and chemical retention/transport. From the perspectives of pore structure in conjunction with multiple and complementary approaches to examining a range of sample sizes under different observational scales, we find that the poor pore connectivity is prevalent in low-permeability media (mudstone and crystalline rock) that is related to geological processes (e.g., compaction, diagenesis and thermal maturation). For example, the deep and organic matter-rich mudstones have a much smaller effective porosity than the total porosity (as a result of poor pore connectivity) and associated diffusion coefficient, and the effective porosity & diffusion coefficients are also dependent upon the sample sizes used in the measurement. Similarly, most of the pore space in the shallow mudstone is also controlled by pore-throat diameters in the 5-50 nm range of intergranular pore types from its fine-grained nature, but with an overall good pore connectivity. However, the nm-sized pore space (physically pore-network architecture) and strong sorption capacities (chemical retention from clay minerals) of both shallow and deep mudstones lead to the synergistic retention of cationic radionuclides and their utilities as effective host rocks and barrier materials. Our unique approaches of studying how the micro-scale pore structure affect macro-scale fluid flow, diffusion & retention, and chemical transport produce improved mechanistic understanding, and realistic quantification, of diffusion and retention of typical radionuclides in a range of generic host rocks and barrier materials (clay/shale, salt, crystalline rock, and tuff), with the overall results leading to scientifically-based understanding of enhanced isolation (from both diffusion and retention) of radionuclides and improved confidence on the long-term performance of geological repository to store high-level radioactive wastes. In addition to the training of 25 undergraduates, graduates, and postdocs of UTA, the scientists (organizations) involved in performing this work (e.g., discussion, sample sharing, and operation of SANS and SAXS instruments) include Ed Matteo, Yifeng Wang, and Kristopher Kuhlman (Sandia National Laboratories), Jens Birkholzer, Liange Zheng, Tim Kneafsey, and Sharon Borglin (Lawrence Berkeley National Laboratory), Mavrik Zavarin (Lawrence Livermore National Laboratory), Yukio Tachi and Yuta Fukatsu (Japan Atomic Energy Agency), Mieke de Craen (Euridice, Belgium), Markus Bleuel (NIST), Wei-Ren Chen, Gergely Nagy, Changwoo Do, William Heller, Larry Anovitz, and Kenneth Littrell (ORNL), as well as Jan Illvsky, Ivan Kuzmenko, Ju-Sang Park and Jon Almers (ANL). Key deliverables include a total of 13 peer-reviewed journal articles (nine published and three under review), 23 presentations at scientific conferences (AAPG, AAPG Southwest Section, AGU, Asian Clay Conference, GSA, GSA South-Central Section, IHLRWM, InterPore, International Conference on Chemistry and Migration Behavior of Actinides and Fission Products in the Geosphere, International Conference on Coupled Processes in Fractured Geological Media: Observation, Modeling and Application), and academic institutions (UTA, New Mexico State University; University of Poitiers, France; University of Helsinki, Finland; Uppsala University, Sweden; Istanbul Technical University, Turkey) and other organizations (Andra, France; Posiva Oy, Finland).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Disentangling the effects of non-adiabatic interactions upon ion self-diffusion within warm dense hydrogen

Warm dense matter is a material state in the region of parameter space connecting condensed matter to classical plasma physics. In this intermediate regime, we investigate the significance of non-adiabatic electron-ion interactions upon ion dynamics. To disentangle non-adiabatic from adiabatic electron-ion interactions, we compare the ion self-diffusion coefficient from the non-adiabatic electron force field computational model with an adiabatic, classical molecular dynamics simulation. A classical pair potential developed through a force-matching algorithm ensures the only difference between the models is due to the electronic inertia. Here, we implement this new method to characterize non-adiabatic effects on the self-diffusion of warm dense hydrogen over a wide range of temperatures and densities. Ultimately we show that the impact of non-adiabatic effects is negligible for equilibrium ion dynamics in warm dense hydrogen.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Dust-void evolution driven by turbulent dust flux can induce runaway migration of Earth-mass planets

Torques from asymmetric dust structures (so-called dust-void and filamentary structures) formed around low-mass planets embedded in a nonturbulent dust-gas disk can exceed the torques produced by the gas disk component and then go on to dominate the planet’s orbital dynamics. Here, we investigate how these structures (hence the dust torque) change when the effect of turbulent dust diffusion and dust feedback are included, along with the direct implications on the migration of Earth-like planets. Using the FARGO3D code, we performed 2D and 3D multifluid hydrodynamic simulations, focusing on a non-migrating planet with a mass of M p = 1.5 M ⊕ in 2D and on migrating planets with M p ∈ [1.5, 12] M ⊕ in 3D. We varied the δ-dimensionless diffusivity parameter in the range [0, 3 × 10 −3 ] and considered three different Stokes numbers, St = {0.04, 0.26, 0.55}, which are representative of the gas-dominated, the transitional, and the gravity-dominated regimes, respectively. In our 2D models, we find that turbulent diffusion of dust prevents the formation of the dust-void and filamentary structures when δ > 3 × 10 −4 . Otherwise, dust structures survive turbulent diffusion flow. However, dust and total torques become positive only in transitional and gravity-dominated regimes. In our 3D models, we find that the dust-void is drastically modified and the high-density ring-shaped barrier delineating the dust-void disappears if δ ≳ 10 −4 , due to the effect of dust turbulent diffusion along with the back-reaction of the dust. For all values of δ, the filament in front of the planet is replaced by a low-density trench. Remarkably, as we have allowed the planets to migrate, the evolving dust-void can drive either runaway migration or outward (inward) oscillatory-torque migration. Our study thus suggests that low-mass Earth-like planets can undergo runaway migration in dusty disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Simulations of self- and Xe diffusivity in uranium mononitride including chemistry and irradiation effects

A combination of density functional theory and empirical potential atomic scale simulations have been used to determine a model for defect stability and mobility in uranium mononitride (UN), as a function of temperature (T) and N 2 partial pressure (p N 2 ). Using the model, predictions of hypo-stoichiometry under U-rich conditions compare favorably to CALPHAD calculations using the TAF-ID database. Furthermore, our predictions of U and N self-diffusivity are in good agreement with experiments carried out as a function of T at specific partial pressures under thermal equilibrium. The validated atomic scale data have then been implemented within a cluster dynamics method to simulate irradiation-enhanced defect concentrations. All defects and clusters studied have significantly enhanced concentrations, with respect to thermal equilibrium, as T is lowered. The irradiation-enhanced Xe diffusivity is compared to post-irradiation annealing and in-pile experiments. In conclusion, the contributions of various defects and clusters to non-stoichiometry, self-diffusivity, and Xe diffusivity are discussed.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗