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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Diffusivity Measurements Made Instant and Easy

A compact common path interferometer (CPI) system has been developed to measure the diffusivity of liquid pairs. The CPI is an optical technique that can be used to measure changes in the gradient of the refraction index of transparent materials. It uses a shearing interferometer that shares the same optical path from a laser light source to the final imaging plane. The molecular diffusion coefficient of liquids can be determined from the physical relations between changes in the optical path length and liquid phase properties. When the data obtained by using the CPI have been compared with similar results from other techniques, the instrument has been demonstrated to be far superior to other instruments for measuring the diffusivity of miscible liquids while staying very compact and robust (ref. 1). Because of its compactness and ease of use, the CPI has been adopted for use in studies of interface dynamics as well as other diffusion-controlled process applications (ref. 2). This progress will permit experiments in microgravity that can quantitatively answer basic science questions about mass and thermal diffusion and their effect in transport processes. This instrument is a spinoff of a diagnostic development for microgravity fluid physics experiments at the NASA Glenn Research Center that has used optics and electronics existing in the fluid physics laboratory for feasibility studies. Optical diagnostic techniques have become an integral part of many areas of measurement applications in industrial and research laboratories. Many types of interferometers and their phase-shifted versions have been used as instruments for measuring optical wave fronts for lens testing and combustion and fluid flow diagnostics. One of these, the point diffraction interferometer, is considered to be robust (see, for example, ref. 3) because it has a common-path design. The point diffraction interferometer is difficult to align and has a limited measurement range for liquid-phase applications. Interferometry and schlieren techniques have been widely used for many years for gas-flow-visualization. On the other hand, the CPI, which is basically a Wollaston prism Polarizer in combination with an analyser (crystalline quartz retardation plate) is used instead of a point diffraction interferometer unit. The advantage of using the CPI over other optical techniques is that it can make quantitative measurements in liquids with the large index of refraction variations that often occur in interface dynamics studies. This can be simply accomplished by using different prism polarizers for each particular experimental condition.

Rashidnia, Nasser↗
NASA NTRS2002

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗
NASA NTRS2004

Non-Additive Voltametric Currents From a Mixture of Two, Three and Four Redox-Active Compounds and Electroanalytical Implications

We have published recently the effect of dissimilar diffusion coefficients on the size of the voltammetric waves from a mixture of two redox-active compounds. Similarly, at the potential range where three redox-active species, decamethylferrocene (dMeFc), ferrocene (Fc) and N-methylphenothiazine (MePTZ), are oxidized simultaneously with rates controlled by linear diffusion, electrogenerated radicals diffusing outwards from the electrode react with the original species diffusing towards the electrode from the bulk; thus, Fc(+) reacts with dMeFc producing Fc and dMeFc(+), while MePTZ(+) reacts both with dMeFc producing MePTZ and dMeFc(+), and with Fc producing MePTZ and Fc(+). These reactions replace dMeFc with Fc at the second plateau, and both dMeFc and Fc with MePTZ at the third plateau. Since the diffusion coefficients of the three species are not equal, the mass-transfer limited currents of the second and the third oxidation wave plateaus change by approx. 10%. Numerical simulations of the experimental voltamograms support this mechanism. Similar results were also obtained for a mixture of four redoxactive compounds. The implications of this non-additive nature of currents on: (a) the use of internal voltammetric standards for quantitative analysis of a mixture of redox-active compounds; and, (b) the half wave potentials (E1/2) of the 2nd, 3rd and 4th waves for qualitative analysis, will be discussed.

Dass, Amala↗
NASA NTRS2005

The SMART Theory and Modeling Team: An Integrated Element of Mission Development and Science Analysis

When targeting physical understanding of space plasmas, our focus is gradually shifting away from discovery-type investigations to missions and studies that address our basic understanding of processes we know to be important. For these studies, theory and models provide physical predictions that need to be verified or falsified by empirical evidence. Within this paradigm, a tight integration between theory, modeling, and space flight mission design and execution is essential. NASA's Magnetospheric MultiScale (MMS) mission is a pathfinder in this new era of space research. The prime objective of MMS is to understand magnetic reconnection, arguably the most fundamental of plasma processes. In particular, MMS targets the microphysical processes, which permit magnetic reconnection to operate in the collisionless plasmas that permeate space and astrophysical systems. More specifically, MMS will provide closure to such elemental questions as how particles become demagnetized in the reconnection diffusion region, which effects determine the reconnection rate, and how reconnection is coupled to environmental conditions such as magnetic shear angles. Solutions to these problems have remained elusive in past and present spacecraft missions primarily due to instrumental limitations - yet they are fundamental to the large-scale dynamics of collisionless plasmas. Owing to the lack of measurements, most of our present knowledge of these processes is based on results from modern theory and modeling studies of the reconnection process. Proper design and execution of a mission targeting magnetic reconnection should include this knowledge and have to ensure that all relevant scales and effects can be resolved by mission measurements. The SMART mission has responded to this need through a tight integration between instrument and theory and modeling teams. Input from theory and modeling is fed into all aspects of science mission design, and theory and modeling activities are tailored to SMART needs during mission development and science analysis. In this presentation, we will present an overview of SMART theory and modeling team activities. In particular, we will provide examples of science objectives derived from state-of-the art models, and of recent research results that continue to be utilized in SMART mission development.

Hesse, Michael↗
NASA NTRS1987

Method of measuring field funneling and range straggling in semiconductor charge-collecting junctions

Electric-field funneling length is measured while irradiating a semiconductor charge-collecting junction with electron-hole-pair generating charged particles at a first junction bias voltage. The bias voltage is then reduced to a second level in order to reduce the depth of the depletion region such that the total charge can no longer be collected by drift and measured in the energy band previously displayed in the multichannel analyzer. This is representative of the maximum electric field funnelling length which may be calculated by measuring the difference at the second bias voltage level of the depletion width and the ion penetration range. The bias voltage is further lowered to a third level at which the particles are collected over a spread of energy levels while at least some of the particles are still collected at the selected energy level. From this the different depths of penetration of the particles are determined while additional effects due to diffusion are minimized.

Zoutendyk, John A.↗
NASA NTRS2009

Unsteady Extinction of Opposed Jet Ethylene/Methane HIFiRE Surrogate Fuel Mixtures vs Air

A unique idealized study of the subject fuel vs. air systems was conducted using an Oscillatory-input Opposed Jet Burner (OOJB) system and a newly refined analysis. Extensive dynamic-extinction measurements were obtained on unanchored (free-floating) laminar Counter Flow Diffusion Flames (CFDFs) at 1-atm, stabilized by steady input velocities (e.g., U(sub air)) and perturbed by superimposed in-phase sinusoidal velocity inputs at fuel and air nozzle exits. Ethylene (C2H4) and methane (CH4), and intermediate 64/36 and 15/85 molar percent mixtures were studied. The latter gaseous surrogates were chosen earlier to mimic ignition and respective steady Flame Strengths (FS = U(sub air)) of vaporized and cracked, and un-cracked, JP-7 "like" kerosene for a Hypersonic International Flight Research Experimentation (HIFiRE) scramjet. For steady idealized flameholding, the 100% C2H4 flame is respectively approx. 1.3 and approx.2.7 times stronger than a 64/36 mix and CH4; but is still 12.0 times weaker than a 100% H2-air flame. Limited Hot-Wire (HW) measurements of velocity oscillations at convergent-nozzle exits, and more extensive Probe Microphone (PM) measurements of acoustic pressures, were used to normalize Dynamic FSs, which decayed linearly with pk/pk U(sub air) (velocity magnitude, HW), and also pk/pk P (pressure magnitude, PM). Thus Dynamic Flame Weakening (DFW) is defined as % decrease in FS per Pascal of pk/pk P oscillation, namely, DFW = -100 d(U(sub air)/U(sub air),0Hz)/d(pkpk P). Key findings are: (1) Ethylene flames are uniquely strong and resilient to extinction by oscillating inflows below 150 Hz; (2) Methane flames are uniquely weak; (3) Ethylene / methane surrogate flames are disproportionately strong with respect to ethylene content; and (4) Flame weakening is consistent with limited published results on forced unsteady CFDFs. Thus from 0 to approx. 10 Hz and slightly higher, lagging diffusive responses of key species led to progressive phase lags (relative to inputs) in the oscillating flames, and caused maximum weakening. At 20 to 150 Hz, diffusion-rate-limited effects diminished, causing flames to "regain strengnth," and eventually become completely insensitive beyond 300 Hz. Detailed mechanistic understanding is needed. Overall, ethylene flames are remarkably resilient to dynamic extinction by oscillating inflows. They are the strongest, with the notable exception of H2. For HIFiRE tests, the 64%/36% surrogate disproportionally retains the high dynamic FS of ethylene, so the potential for loss of scramjet flameholding (flameout) due to low frequency oscillations is significantly mitigated.

Vaden, Sarah N.↗
NASA NTRS2009

Development of High Temperature Dissimilar Joint Technology for Fission Surface Power Systems

NASA is developing fission surface power (FSP) system technology as a potential option for use on the surface of the moon or Mars. The goal is to design a robust system that takes full advantage of existing materials data bases. One of the key components of the power conversion system is the hot-side Heat Exchanger (HX). One possible design for this heat exchanger requires a joint of the dissimilar metals 316L stainless steel and Inconel 718, which must sustain extended operation at high temperatures. This study compares two joining techniques, brazing and diffusion bonding, in the context of forming the requisite stainless steel to superalloy joint. The microstructures produced by brazing and diffusion bonding, the effect of brazing cycle on the mechanical tensile properties of the alloys, and the strength of several brazed joints will be discussed.

Locci, Ivan E.↗
NASA NTRS2011

Advancing Solar Irradiance Measurement for Climate-Related Studies: Accurate Constraint on Direct Aerosol Radiative Effect (DARE)

Earth's climate is driven primarily by solar radiation. As summarized in various IPCC reports, the global average of radiative forcing for different agents and mechanisms, such as aerosols or CO2 doubling, is in the range of a few W/sq m. However, when solar irradiance is measured by broadband radiometers, such as the fleet of Eppley Precision Solar Pyranometers (PSP) and equivalent instrumentation employed worldwide, the measurement uncertainty is larger than 2% (e.g., WMO specification of pyranometer, 2008). Thus, out of the approx. 184 W/sq m (approx.263 W/sq m if cloud-free) surface solar insolation (Trenberth et al. 2009), the measurement uncertainty is greater than +/-3.6 W/sq m, overwhelming the climate change signals. To discern these signals, less than a 1 % measurement uncertainty is required and is currently achievable only by means of a newly developed methodology employing a modified PSP-like pyranometer and an updated calibration equation to account for its thermal effects (li and Tsay, 2010). In this talk, we will show that some auxiliary measurements, such as those from a collocated pyrgeometer or air temperature sensors, can help correct historical datasets. Additionally, we will also demonstrate that a pyrheliometer is not free of the thermal effect; therefore, comparing to a high cost yet still not thermal-effect-free "direct + diffuse" approach in measuring surface solar irradiance, our new method is more economical, and more likely to be suitable for correcting a wide variety of historical datasets. Modeling simulations will be presented that a corrected solar irradiance measurement has a significant impact on aerosol forcing, and thus plays an important role in climate studies.

Tsay, Si-Chee↗
NASA NTRS2016

Materials Outgassing Rate Decay in Vacuum at Isothermal Conditions

As a laboratory for scientific research, the International Space Station has been in Low Earth Orbit for nearly 20 years and is expected to be on-orbit for another 10 years. The ISS has been maintaining a relatively pristine contamination environment for science payloads. Materials outgassing induced contamination is currently the dominant source for sensitive surfaces on ISS and modeling the outgassing rate decay over a 20 to 30 year period is challenging. Materials outgassing is described herein as a diffusion-reaction process using ASTM E 1559 rate data. The observation of -1/2 (diffusion) or non-integers (reaction limited) as rate decay exponents for common ISS materials indicate classical reaction kinetics is unsatisfactory in modeling materials outgassing. Non-randomness of reactant concentrations at the interface is the source of this deviation from classical reaction kinetics. A diffusion limited decay was adopted as the result of the correlation of the contaminant layer thicknesses on returned ISS hardware, the existence of high outgassing silicone exhibiting near diffusion limited decay, and the confirmation of non-depleted material after ten years in the Low Earth Orbit.Keywords: Materials Outgassing, ASTM E 1559, Reaction Kinetics, Diffusion, Space Environments Effects, Contamination

Huang, Alvin Y.↗
NASA NTRS2016

Extension of the Viscous Collision Limiting Direct Simulation Monte Carlo Technique to Multiple Species

There are many flows fields that span a wide range of length scales where regions of both rarefied and continuum flow exist and neither direct simulation Monte Carlo (DSMC) nor computational fluid dynamics (CFD) provide the appropriate solution everywhere. Recently, a new viscous collision limited (VCL) DSMC technique was proposed to incorporate effects of physical diffusion into collision limiter calculations to make the low Knudsen number regime normally limited to CFD more tractable for an all-particle technique. This original work had been derived for a single species gas. The current work extends the VCL-DSMC technique to gases with multiple species. Similar derivations were performed to equate numerical and physical transport coefficients. However, a more rigorous treatment of determining the mixture viscosity is applied. In the original work, consideration was given to internal energy non-equilibrium, and this is also extended in the current work to chemical non-equilibrium.

Liechty, Derek S.↗
NASA NTRS2016

Electron Energization and Mixing Observed by MMS in the Vicinity of an Electron Diffusion Region During Magnetopause Reconnection

Measurements from the Magnetospheric Multiscale (MMS) mission are reported to show distinct features of electron energization and mixing in the diffusion region of the terrestrial magnetopause reconnection. At the ion jet and magnetic field reversals, distribution functions exhibiting signatures of accelerated meandering electrons are observed at an electron out-of-plane flow peak. The meandering signatures manifested as triangular and crescent structures are established features of the electron diffusion region (EDR). Effects of meandering electrons on the electric field normal to the reconnection layer are detected. Parallel acceleration and mixing of the inflowing electrons with exhaust electrons shape the exhaust flow pattern. In the EDR vicinity, the measured distribution functions indicate that locally, the electron energization and mixing physics is captured by two-dimensional reconnection, yet to account for the simultaneous four-point measurements, translational invariant in the third dimension must be violated on the ion-skin-depth scale.

Chen, Li-Jen↗
NASA NTRSDate not supplied

Numerical Investigation of Film Coefficient Engineering Methodology for Dissociated, Chemically Reacting Boundary Layers

Aerothermal analysis of spacecraft planetary entry is heavily dependent on heritage engineering models. The film coefficient heat transfer model examined in this paper estimates the convective heating to the vehicle for a laminar, dissociated, chemically reacting boundary layer for an Earth atmosphere. This model requires information about the vehicle and flow field for a given trajectory point and estimates a proportional relationship between enthalpy potential and convective heat flux. In practice it is the aerothermal engineer who must decide which assumptions are appropriate for their application. This work looks at numerous CFD simulations for an arbitrary, axisymmetric flight vehicle to analyze the relative importance of both the mass and energy constraints imposed at the wall boundary, as well as the effect of various diffusion models. The sensitivity of the engineering model is assessed based on the impact of the CFD assumptions as an effort to make it available to design engineers across industry.

numerical boundary layer↗
NASA NTRS2023

Numerical Investigation of Film Coefficient Approximation for Chemically Reacting Boundary-Layer Flows

Aerothermal analysis of spacecraft planetary entry is heavily dependent on heritage engineering models. The film coefficient heat transfer model examined in this paper estimates the convective heating to the vehicle for a laminar, dissociated, chemically reacting boundary layer for an Earth atmosphere. This model requires information about the vehicle and flowfield for a given trajectory point and estimates a proportional relationship between enthalpy potential and convective heat flux. In practice it is the aerothermal engineer who must decide which assumptions are appropriate for his/her application. This work looks at numerous CFD simulations for an arbitrary, axisymmetric flight vehicle to analyze the relative importance of both the mass and energy constraints imposed at the wall boundary, as well as the effect of various diffusion models. Within the subset of tested energy boundary conditions, it is found that the most desirable energy boundary condition is the radiative equilibrium boundary condition, which permits conservative estimates of convective heat flux, but also generates flowfield-dependent spatial thermal distributions along the surface. Other key findings are presented in an effort to make the film coefficient engineering model readily available to design engineers across industry.

Film coefficient↗
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 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↗