Search NASA⌕ Search

SEARCH · Search NASA

Results for “Jump conditions”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Evolution of coupled weakly driven waves in a dissipative plasma

The nonlinear collisional dynamics of coupled driven plasma waves in the presence of background dissipation is studied analytically within kinetic theory. Sufficiently near marginal stability, phase space correlations are poorly preserved and time delays become unimportant. The system is then shown to be governed by two first-order coupled autonomous differential equations of cubic order for the wave amplitudes and two complementary first-order equations for the evolution of their phases. That system of equations can be decoupled and further simplified to a single second-order differential equation of Liénard's type for each amplitude. Numerical solutions for this equation are obtained in the general case, while analytic solutions are obtained for special cases in terms of parameters related to the spacing of the resonances of the two waves in frequency space, e.g., wave lengths and oscillation frequencies. These parameters are further analyzed to find classes of quasi-steady saturation and pulsating scenarios. To classify equilibrium points, local stability analysis is applied, and bifurcation conditions are determined. When the two waves saturate at similar amplitude levels, their combined signal is shown to invariably exhibit amplitude beating and phase jumps of nearly π. In conclusion, the obtained analytical results can be used to benchmark simulations and to interpret eigenmode amplitude measurements in fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Modeling low cycle fatigue (LCF) of additively manufactured Hastelloy X using An accelerated crystal plasticity fatigue damage model

This paper presents a microstructure-based model for low cycle fatigue (LCF) behavior and life of Nickel-based alloy Hastelloy X manufactured using laser-powder bed fusion (L-PBF) additive manufacturing (AM). AM Hastelloy X, a solution-strengthened alloy, is tested at elevated temperature under fully reversed LCF conditions at different strain levels. A generalized plane strain finite element model is generated from electron backscatter diffraction (EBSD) characterization. The constitutive behavior of the material under fatigue is modeled using crystal plasticity and calibrated with both monotonic tensile and cyclic stress–strain data. The fatigue micro-crack initiation and propagation in the microstructure is modeled using a modified Chaboche fatigue damage model. An embedded boundary condition with a homogenous medium is used to apply the cyclic deformation and prevent numerically introduced over-constraints during fatigue simulation. A ‘cycle-jump’ method is used to accelerate the fatigue simulation and reduce the computational cost. The simulation results are compared to LCF experiments, showing satisfactory matches in cyclic stress behavior and number of cycles to macro-crack initiation for all applied strain ranges. In addition, the model illustrates the potential for quantifying microscale fatigue life impacting factors such as microstructure and surface roughness, which is needed to accurately quantify the reliability of AM components in service.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermal characterization of the build chamber in electron beam melting

Electron beam powder bed fusion, commonly termed electron beam melting (EBM), offers great versatility in multiple-part processes and can produce high quality as-built components due to low residual stresses. The EBM process is complex and requires careful thermal management, including uniform and consistent preheating of the powder bed, in order to ensure quality and consistency of the product. However, most of the simulations in the literature focus on the selective melting stage of the process. As of today, optimal conditions for initial pre-heating temperatures are only available for specific materials. Thus, in order to extend the EBM technology to other desired build materials, a much better understanding of the pre-heating stages is required. In this work, numerical and experimental approaches are combined in order to investigate the effects and sensitivities of heat removal from the build plate during EBM pre-heating stages. For this purpose, a carefully reconstructed numerical model of the build chamber of an ARCAM Q20 + machine is developed. It includes all main parts of the chamber and all relevant heat transfer mechanisms, whereas special attention is paid to radiation heat exchange between various bodies. In order to validate the model, dedicated experiments are performed, in which a system of thermocouples is installed in the build chamber, allowing direct measurement of the local temperatures of the start-plate and heat shields. A good agreement between the simulation and experimental findings is achieved, leading to a better understanding of the thermal processes characteristic to the pre-heating stages. This basic analysis is followed by a representative pre-heating case, where a powder bed is also considered. The energy required to obtain the desired pre-heating temperatures is evaluated, and the role of the powder bed in heat transfer within the chamber is assessed. The pre-heating stage, simulated in the present work, is supposed to create proper conditions for sintering and consequent melting of the powder, which are highly dependent on the local temperatures and heat transfer features. Thus, the findings of the reported work present a step towards a better understanding of the thermal processes that characterize EBM. The reported model can be further used to provide realistic boundary condition inputs for other meso- or macro-scale models as a function of time and geometry. The model can serve also for verification of machine settings, i.e., jump-safe and melt-safe ones which actually provide desired preheat, and for development of settings for new powders.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The essential character of the neutrino mechanism of core-collapse supernova explosions

Abstract Calibrating with detailed 2D core-collapse supernova (CCSN) simulations, we derive a simple CCSN explosion condition based solely upon the terminal density profiles of state-of-the-art stellar evolution calculations of the progenitor massive stars. This condition captures the vast majority of the behaviour of the one hundred 2D state-of-the-art models we performed to gauge its usefulness. The goal is to predict, without resort to detailed simulation, the explodability of a given massive star. We find that the simple maximum fractional ram pressure jump discriminant we define works well ∼90 per cent of the time and we speculate on the origin of the few false positives and false negatives we witness. The maximum ram pressure jump generally occurs at the time of accretion of the silicon/oxygen interface, but not always. Our results depend upon the fidelity with which the current implementation of our code F ornax adheres to Nature and issues concerning the neutrino–matter interaction, the nuclear equation of state, the possible effects of neutrino oscillations, grid resolution, the possible role of rotation and magnetic fields, and the accuracy of the numerical algorithms employed remain to be resolved. Nevertheless, the explodability condition we obtain is simple to implement, shows promise that it might be further generalized while still employing data from only the unstable Chandrasekhar progenitors, and is a more credible and robust simple explosion predictor than can currently be found in the literature.

Astronomy & Astrophysics↗
DOE OSTI2025

Efficient Quantum Gibbs Samplers with Kubo–Martin–Schwinger Detailed Balance Condition

Lindblad dynamics and other open-system dynamics provide a promising path towards efficient Gibbs sampling on quantum computers. In these proposals, the Lindbladian is obtained via an algorithmic construction akin to designing an artificial thermostat in classical Monte Carlo or molecular dynamics methods, rather than being treated as an approximation to weakly coupled system-bath unitary dynamics. Recently, Chen, Kastoryano, and Gilyén (arXiv:2311.09207) introduced the first efficiently implementable Lindbladian satisfying the Kubo–Martin–Schwinger (KMS) detailed balance condition, which ensures that the Gibbs state is a fixed point of the dynamics and is applicable to non-commuting Hamiltonians. This Gibbs sampler uses a continuously parameterized set of jump operators, and the energy resolution required for implementing each jump operator depends only logarithmically on the precision and the mixing time. In this work, we build upon the structural characterization of KMS detailed balanced Lindbladians by Fagnola and Umanità, and develop a family of efficient quantum Gibbs samplers using a finite set of jump operators (the number can be as few as one), akin to the classical Markov chain-based sampling algorithm. Compared to the existing works, our quantum Gibbs samplers have a comparable quantum simulation cost but with greater design flexibility and a much simpler implementation and error analysis. Moreover, it encompasses the construction of Chen, Kastoryano, and Gilyén as a special instance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Micropinch formation dynamics in X pinches

High temporal resolution x-ray streak camera studies of micropinch formation in Cu hybrid x pinches reveal key plasma conditions. Analysis of Ne-like Cu lines indicate an average electron temperature of about 200 eV and 4.5×10 28 m -3 electron density. Here, the spectra suggest that the electron temperature jumps to about 1 keV, inferred from the continuum and the postcontinuum line emission that includes Li-like Cu lines. There is no sign of a rapid temperature change or a substantial surge in radiation emission during the 200 ps precontinuum x-ray burst, suggesting that the radiative collapse process does not play a major role in micropinch formation. Two-dimensional extended Magnetohydrodynamic (MHD) simulations, coupled to a collisional-radiative spectral analysis code, suggest the significance of the rapid radial implosion of high-temperature, low-density plasma, the axial outflow, and the dynamic plasma pressure in micropinch formation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

The Dependence of Electron-only Magnetic Reconnection on Plasma Beta and Magnetic Shear: Theory and Observations

We investigate the onset conditions for electron-only magnetic reconnection through an analysis of 86 reconnection events observed by the Magnetospheric Multiscale mission in magnetosheath and magnetospheric current sheets. Our study presents the first observational evidence that the occurrence of electron-only reconnection depends on both the magnetic-field shear angle (θ) across the current sheet and the cross-sheet electron beta jump, quantified either by the inflow electron beta difference (Δβ e ) or by a modified electron beta jump (Δβ$^*_e$) defined by the cross-sheet electron pressure difference normalized by the average of the inflow magnetic pressure in the two inflow regions. We extend the diamagnetic-drift suppression mechanism from standard reconnection to electron-only reconnection by deriving relations applicable to electron-scale current sheets. Using scaling estimates, the electron-only reconnection suppression relation is derived as Δβ$^*_e$ ≳ 2(L/d e ) tan (θ/2), while a natural extension of the ion-coupled reconnection suppression criterion is Δβ e ≳ 2(L/d e ) tan (θ/2). At low Δβ e and Δβ$^*_e$ values, reconnection is detected across a large range of magnetic shear (from low to high), whereas high Δβ e and Δβ$^*_e$ conditions restrict reconnection to large shear configurations: shear angles as small as θ ≈ 5° occur when Δβ e (or Δβ$^*_e$) ≲0.1. However, only events exceeding θ ≳ 70° are found when Δβ e ≳ 1.0; this apparent restriction weakens when employing Δβ$^*_e$, as the corresponding data points shift toward smaller values of the modified electron beta jump. These findings quantitatively support theoretical predictions of reconnection suppression in high-β e plasmas at low shears, attributed to X-line drift induced by pressure gradients across electron-scale current sheets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A cycle-jump acceleration method for the crystal plasticity simulation of high cycle fatigue of the metallic microstructure

For metallic materials, high-cycle fatigue life is sensitive to underlying microstructure features including secondary phases, textures, grains morphology, etc. The traditional, data-based safe-life approaches for modeling fatigue don’t explicitly consider the microstructure and can’t guide study in microstructure modification for improved fatigue property. Crystal plasticity-based simulation provides increased model fidelity at the expense of immense computation time, making it inapplicable for high cycle fatigue. In this work, an acceleration method based on cycle-jump approach (Lesne and Savalle, 1989) was developed for microstructure-based high-cycle FE simulation using crystal plasticity constitutive-model. We report this method demonstrated high efficiency in benchmark tests of various conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Cause and impact of low-frequency chirping modes in DIII-D hybrid discharges

Significant variations in MHD activity and fast-ion transport are observed in the DIIID high-beta, steady-state hybrid discharges with a mixture of Electron Cyclotron (EC) waves and Neutral Beam Injection (NBI). When Electron Cyclotron Heating (ECH) or Current Drive (ECCD) is applied, Alfven Eigenmodes (AEs) are usually suppressed and replaced by low-frequency bursting modes. The analysis of a recently compiled database of hybrid discharges suggests that the change of the fast-ion pressure especially the perpendicular pressure is the main factor responsible for the instability transition. The lower ratio of fast-ion injection speed v inj to Alfven speed v alfven and slight drop of safety factor q min during ECCD may also facilitate the transition. The database shows that AEs mainly occur when the fast-ion fraction P f /P total is less than 0.53 and v inj /v alfven is greater than 0.50, while low-frequency bursting modes appear in the opposite regime. In this research, P f and P total are the central fast-ion pressure from classical prediction and total plasma pressure, respectively. The reason why the instability transition correlates with P f /P total and v inj /v alfven is that they can significantly modify the drive of low-frequency bursting modes and AEs. The explanation is supported by the observation that low-frequency bursting modes are rarely seen in the hybrids with NBI only, with EC waves and counter-NBI, or with high plasma density. A careful check of the low-frequency bursting modes suggests that they are mainly chirping (neoclassical) tearing modes [referred to as chirping (N)TM], i.e. the mode frequency firstly jumps up from the steady (N)TM frequency, then chirps down, and finally returns to the steady (N)TM frequency. Sometimes, the (N)TMs are fully stabilized and replaced with pure fishbones. Lastly, the resonance condition calculation and ’Kick’ model simulations suggest that (N)TMs and fishbones can interact through modification of the fast ion distribution in phase space, which influences the drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Open Source Fault-tolerant Grid Frequency Measurement for Solar Inverters

The Discrete Fourier transform (DFT) based measurement algorithms are one of the most common measurement algorithms for grid parameter estimation such as rms, phase angle, frequency. Over the past few years, many DFT based algorithms have been developed to enhance its measurement accuracy under steady-state and/or dynamic grid conditions. For example, an adaptive band-pass filter utilizing exponential modulation filter has been proposed to reduce measurement errors at the presence of large frequency deviations. Measurement accuracy of different algorithms including FIR filter, extended Kalman filtering (EKF), and enhanced DFT method have been compared in detail under different grid conditions. Two artificial signals that have 90-degree phase difference were constructed by the Clarke transformation to address the frequency spectrum leakage of DFT. A multi-module approach was developed to enhance both steady-state and dynamic measurement accuracies, in which each module was developed to eliminate some specific errors. Besides DFT-based measurement algorithms, some signal model-based algorithms have been developed to further improve the accuracy under dynamic conditions. However, a key drawback of the state-of-the-art algorithms is that they cannot perform measurements accurately during system transient faults. In the Blue Cut Fire event, there was a phase angle jump of about 26 degrees in the voltage waveform during the transient fault. The phase angle jump fault will cause waveform discontinuity, and these algorithms will fail to provide reliable measurements during this period because they typically assume the waveform to be measured is continuous, no matter what method (DFT, PLL, EKF, FIR, or Taylor WLS) is used for estimation. In fact, the measurement errors during the system transient faults like phase-jump is not required in the IEEE Standard. As a result, although a measurement instrument can pass the strict IEEE Standard, it could still be the source of the problem in the future if we have similar system transient faults, which could happen again. Therefore, developing the fault-tolerant measurement technology is the key to solve the problem.

14 SOLAR ENERGY↗
DOE OSTI2021

Pore-Scale Dynamics of Liquid CO2–Water Displacement in 2D Axisymmetric Porous Micromodels Under Strong Drainage and Weak Imbibition Conditions: High-Speed μPIV Measurements

Resolving pore-scale transient flow dynamics is crucial to understanding the physics underlying multiphase flow in porous media and informing large-scale predictive models. Surface properties of the porous matrix play an important role in controlling such physics, yet interfacial mechanisms remain poorly understood, in part due to a lack of direct observations. This study reports on an experimental investigation of the pore-scale flow dynamics of liquid CO 2 and water in two-dimensional (2D) circular porous micromodels with different surface characteristics employing high-speed microscopic particle image velocimetry (μPIV). The design of the micromodel minimized side boundary effects due to the limited size of the domain. The high-speed μPIV technique resolved the spatial and temporal dynamics of multiphase flow of CO 2 and water under reservoir-relevant conditions, for both drainage and imbibition scenarios. When CO 2 displaced water in a hydrophilic micromodel (i.e., drainage), unstable capillary fingering occurred and the pore flow was dominated by successive pore-scale burst events (i.e., Haines jumps). When the same experiment was repeated in a nearly neutral wetting micromodel (i.e., weak imbibition), flow instability and fluctuations were virtually eliminated, leading to a more compact displacement pattern. Energy balance analysis indicates that the conversion efficiency between surface energy and external work is less than 30%, and that kinetic energy is a disproportionately smaller contributor to the energy budget. This is true even during a Haines jump event, which induces velocities typically two orders of magnitude higher than the bulk velocity. These novel measurements further enabled direct observations of the meniscus displacement, revealing a significant alteration of the pore filling mechanisms during drainage and imbibition. While the former typically featured burst events, which often occur only at one of the several throats connecting a pore, the latter is typically dominated by a cooperative filling mechanism involving simultaneous invasion of a pore from multiple throats. This cooperative filling mechanism leads to merging of two interfaces and releases surface energy, causing instantaneous high-speed events that are similar, yet fundamentally different from, burst events. Finally, pore-scale velocity fields were statistically analyzed to provide a quantitative measure of the role of capillary effects in these pore flows.

Li, Yaofa↗
DOE OSTI2023

Nano-infrared imaging of metal insulator transition in few-layer 1T-TaS 2

Abstract Among the family of transition metal dichalcogenides, 1T-TaS 2 stands out for several peculiar physical properties including a rich charge density wave phase diagram, quantum spin liquid candidacy and low temperature Mott insulator phase. As 1T-TaS 2 is thinned down to the few-layer limit, interesting physics emerges in this quasi 2D material. Here, using scanning near-field optical microscopy, we perform a spatial- and temperature-dependent study on the phase transitions of a few-layer thick microcrystal of 1T-TaS 2 . We investigate encapsulated air-sensitive 1T-TaS 2 prepared under inert conditions down to cryogenic temperatures. We find an abrupt metal-to-insulator transition in this few-layer limit. Our results provide new insight in contrast to previous transport studies on thin 1T-TaS 2 where the resistivity jump became undetectable, and to spatially resolved studies on non-encapsulated samples which found a gradual, spatially inhomogeneous transition. A statistical analysis suggests bimodal high and low temperature phases, and that the characteristic phase transition hysteresis is preserved down to a few-layer limit.

42 ENGINEERING↗
DOE OSTI2020

Effect of Solute Elements on Boron Segregation in Boron-Containing Steels

The addition of boron to steel alloys results in an increase in both hardenability and casting defects. The casting difficulties are predicted to stem from a metatectic reaction, δ + γ → L + γ, where a fully solidified material begins to locally remelt as the temperature decreases. Another possible source of casting defects is a boride-rich phase that is predicted to remain liquid at low temperatures. To experimentally determine which reaction is the likely source of the casting defects, the predicted reactions and the effect of solute elements on those reactions are investigated. Levitation zone melting is used to control segregation in a ternary Fe–C–B alloy and a commercial 22MnB5 alloy. Carbon segregation and a peritectic reaction result in a peritectic jump during directional solidification where the first directionally solidified (DS) zone undergoes δ-bcc solidification followed by a peritectic jump to steady state planar solidification of γ-fcc in the second DS zone. The presence of other solute elements in the zone melted 22MnB5 alloy lead to a breakdown in the planar solidification front before steady state solidification could be achieved in the second DS zone. With a cellular solid/liquid interface, boron-rich intercellular liquid formed low melting iron boro-carbide particles. The controlled solidification conditions in a levitation zone melter were unable to prevent ~0.003 wt% boron from segregating to high enough levels to form boride particles. Therefore, it is likely that during commercial casting, the formation of the low melting boride phase from interdendritic segregation is a key source of the casting issues.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Measurement bias in self-heating x-ray free electron laser experiments from diffraction studies of phase transformation in titanium

X-ray self-heating is a common by-product of X-ray Free Electron Laser (XFEL) techniques that can affect targets, optics, and other irradiated materials. Diagnosis of heating and induced changes in samples may be performed using the x-ray beam itself as a probe. However, the relationship between conditions created by and inferred from x-ray irradiation is unclear and may be highly dependent on the material system under consideration. Here, we report on a simple case study of a titanium foil irradiated, heated, and probed by a MHz XFEL pulse train at 18.1 keV delivered by the European XFEL using measured x-ray diffraction to determine temperature and finite element analysis to interpret the experimental data. We find a complex relationship between apparent temperatures and sample temperature distributions that must be accounted for to adequately interpret the data, including beam averaging effects, multivalued temperatures due to sample phase transitions, and jumps and gaps in the observable temperature near phase transformations. The results have implications for studies employing x-ray probing of systems with large temperature gradients, particularly where these gradients are produced by the beam itself. Finally, this study shows the potential complexity of studying nonlinear sample behavior, such as phase transformations, where biasing effects of temperature gradients can become paramount, precluding clear observation of true transformation conditions.

Crystallography↗
DOE OSTI2022

MPC solution for optimal load shifting for buildings with ON/OFF staged packaged units: Experimental demonstration, and lessons learned

Small and medium-sized commercial buildings (SMCB) are significant demand response resources, and it is important to develop grid-responsive control algorithms that exploit those resources and create financial benefits for building owners and HVAC service providers. Furthermore, unlike large-sized commercial buildings, there is an opportunity to have universally applicable control solutions for many SMCBs since those buildings have a consistent HVAC system configuration: SMCBs are commonly served by multiple-staged air conditioning units controlled by their own thermostats. Despite the demand response potential and scalability, however, very few control solutions are available for SMCBs. Typical model predictive control (MPC) and heuristic control approaches for cooling load shifting that lower thermostat setpoints before an electric price jump are suitable mainly for large-sized commercial buildings where a continuous capacity modulation is possible, e.g., via dampers in variable air volume terminal units. However, those approaches can cause undesired, high peaks for SMCBs due to the nature of ON/OFF unit staging and narrow thermostat deadbands. This could discourage the use of advanced grid-responsive controls for SMCBs due to the concern of high demand charges, and has to be resolved. This paper presents a MPC solution that overcomes this challenge. It has a hierarchical MPC structure where an upper level MPC is responsible for electrical load shifting in response to an electric price signal while a lower level MPC is responsible for coordinating compressor stages to eliminate unnecessary peaks and follows the setpoints determined by the upper level MPC. In this work, two one-month, comprehensive laboratory tests have been carried out to demonstrate load shifting and cost savings for the algorithm. Interesting trade-offs between energy efficiency and load flexibility were observed and are discussed, and lessons learned for applying MPCs for SMCBs are also presented.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

A Large-Scale Hardware Experiment Demonstration of Operating High Inverter-Based Resource Power Systems With Grid-Forming Inverters

This paper presents experimental hardware results from a microgrid system as the penetration level of grid-forming (GFM) inverters increases. The experiment aims to showcase the advantages of GFM inverters in enhancing system stability and to investigate the operational challenges in systems relying entirely on inverter-based resources (IBRs). Five test scenarios are devised to progressively increase GFM inverter penetration levels: 0% (S1), 19% (S2), 37% (S3), 68% (S4), and 100% (S5). To ensure consistency, identical loading conditions and dynamic events are applied across all scenarios. The conducted tests include load step changes, output variation of grid-following (GFL) inverters during islanded mode, transition operations (such as synchronization to the grid and islanding), and rateof- change-of-frequency (ROCOF) and voltage jump tests during grid-connected mode. Key findings from the tests are summarized as follows: (1) Scenario S1 proved most challenging for load steps, and GFL variations because of insufficient GFM capacity, leaving the diesel generator unable to handle transient ridethrough; (2) Scenario S5 was the most difficult for transition operations, because the lack of a diesel generator to maintain stiff bus voltage resulted in unexpected reactive power flows during synchronization, causing the point of common coupling (PCC) breaker to trip; (3) Scenarios S4 and S5 were particularly challenging for ROCOF tests because of minimal system inertia, leading to large transients that triggered breaker trips; and (4) Scenario S4 exhibited the strongest voltage recovery capability because of the presence of the diesel generator and a higher number of GFM inverters, both of which possess the highest capacity for reactive power injection.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Mechanisms of Retention and Transport of Fission Products in Virgin and Irradiated Nuclear Graphite

To fulfill the design and licensing requirements for advanced reactors, the U.S. Department of Energy (DOE) is specifically interested in targeted research to quantify the mechanisms of fission product transport and retention in graphitic grades that will likely be used in next-generation very/high temperature reactors, fluoride salt-cooled high temperature reactors and molten salt reactors. We have examined an integrated experimental-computational approach to determine the diffusivities of selected fission products – ruthenium (Ru) and silver (Ag) in five nuclear graphite grades – POCO AXF-5Q, POCO ZXF-5Q, PCEA, IG110, and NBG-18. Experimental investigation for diffusion of other fission products (iodine, cesium, and krypton) in graphite could not be carried out due to extreme evaporation due to high vapor pressure, despite several attempts incorporating various experimental schemes. Experimental investigation with irradiated graphite could not be carried out because necessary approval to transport irradiated graphite could not be obtained due to procedural reasons at ORNL and national policy change in UK. An effort was made to “produce” by ion sputtering, but it was discontinued because the outcome did not adequately represent irradiated graphite. Thin film and “cup-and-cap” methods were employed to deposit fission products, diffusion annealed were carried out in an argon atmosphere, and concentration profiles were determined by dynamic secondary ion mass spectroscopy depth profiling technique. Regardless of graphite type, anomalous diffusion behavior of Ru (i.e., lower than expected based on Arrhenius temperature dependence) at an intermediate temperature, around 600 to 700 °C, was observed. Diffusion behavior of silver on the other hand agreed well to the Arrhenius temperature dependence although the excessive evaporation of silver hampered the accurate determination of diffusivities at temperature higher than 900 °C. In general, silver was observed to diffuse faster than Ru in the temperature range compared from 500 to 900 °C, had consistently had higher pre-exponential factor corresponding to higher jump frequency. To understand the abnormal diffusional behavior of Ru in graphite in intermediate temperature range, a detailed microstructural analysis was carried out. Clustering of Ru into particles on the surface of graphite was observed around 700 °C corresponding to de-wetting and spheroidization. This change would influence the boundary condition of diffusion, although features of clustering was much smaller (<1 mm) than the SIMS raster size (120 x 120 mm). Raman spectroscopy also demonstrated that Raman peak intensity for Ru decreased significantly for samples annealed higher than 700 °C and presence of Ru in graphite altered the crystallinity of graphite. To further elucidate the abnormal diffusion behavior observed in intermediate temperature investigated, the vibrational properties of graphene, bilayer graphene and Kr atoms in bilayer graphene were investigated with molecular dynamics. Krypton, which is a fission product, was introduced between a bilayer graphene system to evaluate (i) whether the phonon dispersion curves change with high temperature in bilayer graphene (ii) whether the addition of Kr atoms between the bilayer introduces any significant change in the vibrational properties, and (iii) whether any change in the vibrational properties can be correlated to the diffusive behavior of Kr atoms. The results show that the Kr atoms indeed exhibit an anomalous diffusive behavior at 550°C. Phonon analysis shows that a few phonon frequencies are enhanced as much as 20 to 30 % with the addition of just 0.2 at.% Kr. So, this may explain the unexpected self-diffusivity of fission products.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Analysis of nonlinear features associated with radio-frequency sheaths

An analytic study is conducted of nonlinear features that arise when radio-frequency (RF) sheath boundary conditions are solved self-consistently with the electromagnetic fields in a magnetized plasma. An instability associated with RF sheath–plasma interactions is identified and shown to be closely linked with sheath–plasma resonance and the occurrence of multiple solutions. It is proposed that this instability serves as a physical mechanism that enables jumping from one solution branch to another, providing a physical basis for the prediction of an RF sheath hysteresis in the time domain. It is further shown that this instability can lead to non-intuitive features such as spontaneous symmetry breaking in an otherwise completely symmetric problem. In plane-stratified geometry, two analytically solvable examples are considered for the cases of an unbounded half-space and a bounded domain

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗