Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface profiling”

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 307 records · Page 17

Seismic Monitoring at the Farnsworth CO2-EOR Field Using Time-Lapse Elastic-Waveform Inversion of 3D-3C VSP Data

During the Development Phase of the U.S. Southwest Regional Partnership on Carbon Sequestration, supercritical CO2 was continuously injected into the deep oil-bearing Morrow B formation of the Farnsworth Unit in Texas for Enhanced Oil Recovery (EOR). The project injected approximately 94 kilotons of CO2 to study geologic carbon storage during CO2-EOR. A three-dimensional (3D) surface seismic dataset was acquired in 2013 to characterize the subsurface structures of the Farnsworth site. Following this data acquisition, the baseline and three time-lapse three-dimensional three-component (3D-3C) vertical seismic profiling (VSP) data were acquired at a narrower surface area surrounding the CO2 injection and oil/gas production wells between 2014 and 2017 for monitoring CO2 injection and migration. With these VSP datasets, we inverted for subsurface velocity models to quantitatively monitor the CO2 plume within the Morrow B formation. We first built 1D initial P-wave (Vp) and S-wave (Vs) velocity models by upscaling the sonic logs. We improved the deep region of the Vp and Vs models by incorporating the deep part of a migration velocity model derived from the 3D surface seismic data. We improved the shallow region of 3D Vp and Vs models using 3D traveltime tomography of first arrivals of VSP downgoing waves. We further improved the 3D baseline velocity models using elastic-waveform inversion (EWI) of the 3D baseline VSP upgoing data. Our advanced EWI method employs alternative tomographic and conventional gradients and total-variation-based regularization to ensure the high-fidelity updates of the 3D baseline Vp and Vs models. We then sequentially applied our 3D EWI method to the three time-lapse datasets to invert for spatiotemporal changes of Vp and Vs in the reservoir. Our inversion results reveal the volumetric changes of the time-lapse Vp and Vs models and show the evolution of the CO2 plume from the CO2 injection well to the oil/gas production wells.

42 ENGINEERING↗

The Effects of Contaminant-Aging on Decontamination Efficacy for Rapid Remediation of Concrete Surfaces - 20412

Understanding the relationship between decontamination efficacy and contaminant aging is imperative for developing effective remediation strategies following a large-scale contamination event. Non-destructive decontamination of concrete following a contamination incident becomes more difficult with increased time between contamination and decontamination because contaminants may chemically bind to concrete and/or penetrate into the concrete subsurface. In this work, we evaluated the decontamination efficacy of two decontamination methods on concrete samples contaminated with soluble Cesium-137 (Cs-137) and silica particles with 0.5 μm and 2 μm diameters. Concrete samples were aged between 1 and 59 days, with half of the coupons receiving 1 mL of artificial rainwater about once every three days. After aging, coupons were either decontaminated or analyzed to determine contaminant penetration depths. Coupons were decontaminated with 0.1 M Potassium Chloride (KCI) solution applied either by flowing solution across the contaminated coupon face or non-destructive power washing. Contaminant depth profiles were created by removing the top surface of the coupon at least 15 times, measuring the activity in each removed layer, normalizing the total activity removed to the measured decontamination efficacy following the last layer removal, and determining the penetration thickness using the coupon dimensions, bulk density, and mass of material removed for each layer. Pressurized washing effectively removed particles, but not soluble cesium, despite the majority of Cs-137 being located within the first millimeter of the subsurface. Soluble cesium removal by flowing 0.1 M KCI solution across the contaminated surface dropped to less than 12% within the first ten days. This drop in Cs-137 removal efficacy occurred within the first five days for coupons experiencing artificial rainfall events. Flowing 0.1 M KCI solution across the contaminated face was ineffective at removing particles, likely because particles had localized in surface depressions. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Antarctic permafrost processes and antiphase dynamics of cold-based glaciers in the McMurdo Dry Valleys inferred from 10 Be and 26 Al cosmogenic nuclides

Abstract. Soil and sediment mixing and associated permafrost processes are not widely studied or understood in the McMurdo Dry Valleys of Antarctica. In this study, we investigate the stability and depositional history of near-surface permafrost sediments to ∼ 3 m depth in the Pearse and lower Wright valleys using measured cosmogenic 10Be and 26Al depth profiles. In Pearse Valley, we estimate a minimum depositional age of ∼ 74 ka for the active layer and paleoactive-layer sediments (< 0.65 m). Combined depth profile modelling of 10Be and 26Al gives a depositional age for near-surface (< 1.65 m) permafrost in Pearse Valley of 180 +20/-40 ka, implying that the deposition of permafrost sediments predates MIS 5 advances of Taylor Glacier. Deeper permafrost sediments (> 2.09 m) in Pearse Valley are thus inferred to have a depositional age of > 180 ka. At a coastal, lower-elevation site in neighbouring lower Wright Valley, 10Be and 26Al depth profiles from a second permafrost core exhibit near-constant concentrations with depth and indicate the sediments are either vertically mixed after deposition or sufficiently young so that post-depositional nuclide production is negligible relative to inheritance. 26Al/10Be concentration ratios for both depth profiles range between 4.0 and 5.2 and are all lower than the nominal surface production rate ratio of 6.75, indicating that prior to deposition, these sediments experienced complex, yet similar, exposure–burial histories. Assuming a single-cycle exposure–burial scenario, the observed 26Al/10Be ratios are equivalent to a total minimum exposure–burial history of ∼ 1.2 Myr. In proximity to the depth profile core site, we measured cosmogenic 10Be and 26Al in three granite cobbles from thin, patchy drift (Taylor 2 Drift) in Pearse Valley to constrain the timing of retreat of Taylor Glacier. Assuming simple continuous exposure, our minimum, zero-erosion exposure ages suggest Taylor Glacier partially retreated from Pearse Valley no later than 65–74 ka. The timing of retreat after 65 ka and until the Last Glacial Maximum (LGM) when Taylor Glacier was at a minimum position remains unresolved. The surface cobble ages and permafrost processes reveal Taylor Glacier advances during MIS 5 were non-erosive or mildly erosive, preserving the underlying permafrost sediments and peppering boulders and cobbles upon an older, relict surface. Our results are consistent with U/Th ages from central Taylor Valley and suggest changes in moisture delivery over Taylor Dome during MIS 5e, 5c, and 5a appear to be associated with the extent of the Ross Ice Shelf and sea ice in the Ross Sea. These data provide further evidence of antiphase behaviour through retreat of a peripheral lobe of Taylor Glacier in Pearse Valley, a region that was glaciated during MIS 5. We suggest a causal relationship of cold-based glacier advance and retreat that is controlled by an increase in moisture availability during retreat of sea ice and perhaps the Ross Ice Shelf, as well as, conversely, a decrease during times of sea ice and Ross Ice Shelf expansion in the Ross Sea.

58 GEOSCIENCES↗

In Situ Studies of Copper-Zirconia-Zinc Oxide Model Catalysts for CO 2 Hydrogenation

A combination of in situ X-ray photoelectron spectroscopy (XPS) and infrared reflection absorption spectroscopy (IRAS) was used to investigate the formation of surface intermediates from CO 2 hydrogenation on copper-zirconia-zinc oxide model catalysts under reaction conditions. Copper clusters with different numbers of atoms (n = 1, 4, 13) were deposited onto bare and ZrO 2 -modified ZnO powder supports to systematically examine the effects of Cu cluster size and the synergy between metal and metal-oxide components. Under low pressure CO 2 hydrogenation conditions (CO 2 :H 2 = 1:9, 0.4 mbar, 300–600 K), XPS and IRAS identify the most prominent intermediates as carbonate (CO 3 *), formate (HCOO*) and methoxy (CH 3 O*), which is the final surface-bound intermediate leading to methanol. The temperature profiles are consistent with a mechanism in which CO 2 is adsorbed as carbonate species (HCO 3 *, CO 3 *) followed by hydrogenation reactions to formate (HCOO*) and methoxy (CH 3 O*), but the relative yields strongly depend on surface composition. Specifically, the presence of ZrO 2 promotes CO 2 adsorption and activation and improves the thermal stability of the Cu clusters against loss of surface area. Moreover, the ternary Cu 4 /ZrO 2 /ZnO surface is significantly more active than Cu 4 /ZnO and ZrO 2 /ZnO surfaces for the formation of methoxy (CH 3 O*), indicating that Cu–ZrO 2 interfaces promote the formation of key intermediates leading to methanol. Finally, the yields of intermediates are similar for all Cu cluster sizes (Cu 1 , Cu 4 and Cu 13 ), indicating that the primary role of Cu is to provide H atoms via H 2 dissociation and spillover. In conclusion, these molecular-level insights provide a fundamental understanding of the enhanced efficiency of ternary Cu–ZrO 2 –ZnO catalysts and establish design principles for developing improved catalysts for CO 2 conversion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis of position-dependent cavity parameters in irradiated metals to obtain insight on fundamental defect migration phenomena

Motion of point defects is a fundamental process that governs microstructure and properties of materials. Here, we examine the near-surface and grain boundary cavity swelling depth profiles in neutron- and ion-irradiated simple metals (Cu, Ni, and Fe-Cr) and investigate diffusional broadening of implanted Ni ions in Fe-Cr alloys. Vacancy migration energies and radiation-enhanced diffusion were experimentally estimated. Cavity denuded zone widths near planar sinks are shown to be dependent on temperature, damage rate, balance of point defects and sinks, and vacancy migration energies. An enhanced cavity swelling zone adjacent to the void-denuded zone was observed in specimens with low to moderate sink strength and interpreted as evidence of 1D gliding interstitial clusters. Radiation-enhanced diffusional broadening of implanted Ni ions (at 400–550 °C) in Fe and Fe-Cr was up to 250 nm toward the surface, which leads to a much broader near-peak-damage region where cavity swelling is suppressed. Diffusional broadening of implanted ions is calculated to be similarly pronounced for self-ion irradiations, particularly near the peak and higher swelling temperature regimes. Adequately high ion energies (~8–15 MeV) are recommended for ion-irradiation studies to provide a sufficiently broad midrange safe analysis region with minimized surface and injected ion effects.

36 MATERIALS SCIENCE↗

Retrieved Number Concentration of Cloud Condensation Nuclei (RNCCN) Profile Value-Added Product Report

The cloud condensation nuclei (CCN) concentration at cloud base is the most relevant measure of the aerosol that influences droplet formation in clouds. Since the CCN concentration depends on supersaturation, a more general measure of the CCN concentration is the CCN spectrum (values at multiple supersaturations). The CCN spectrum is now measured at the surface at several U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility observatories and by the ARM Mobile Facility (AMF) but is not measured at the cloud base. Rather than rely on expensive aircraft measurements for all studies of aerosol effects on clouds, a way to project CCN observations made at the surface to cloud base is needed. Remote sensing of aerosol extinction provides information about the vertical profile of the aerosol but cannot be directly related to the CCN concentration because the aerosol extinction is strongly influenced by humidification, particularly near cloud base. Ghan and Collins (2004) and Ghan et al. (2006) propose a method to remove the influence of humidification from the extinction profiles and tie the “dry extinction” retrieval to the surface CCN concentration, thus estimating the CCN profile. This methodology has been implemented as ARM’s Retrieved Number Concentration of CCN (RNCCN) Profile Value-Added Product (VAP).

54 ENVIRONMENTAL SCIENCES↗

Chemical Mixing Induced by Internal Gravity Waves in Intermediate-mass Stars

Internal gravity waves can cause mixing in the radiative interiors of stars. We study this mixing by introducing tracer particles into 2D hydrodynamic simulations. Following the work of Rogers & McElwaine, we extend our study to different masses (3, 7, and 20 M ⊙ ) and ages (ZAMS, midMS, and TAMS). The diffusion profiles of these models are influenced by various parameters such as the Brunt–Väisälä frequency, density, thermal damping, the geometric effect, and the frequencies of waves contributing to these mixing profiles. We find that the mixing profile changes dramatically across age. In younger stars, we noted that the diffusion coefficient increases toward the surface, whereas in older stars the initial increase in the diffusion profile is followed by a decreasing trend. We also find that mixing is stronger in more massive stars. Hence, future stellar evolution models should include this variation. In order to aid the inclusion of this mixing in 1D stellar evolution models, we determine the dominant waves contributing to these mixing profiles and present a prescription that can be included in 1D models.

79 ASTRONOMY AND ASTROPHYSICS↗

Rational surfaces, flows and radial structure in the TJ-II stellarator

Abstract In this work, we report on the results obtained by measuring several turbulent quantities well inside the plasma edge by means of a Langmuir probe during dynamical rotational transform scans in the TJ-II stellarator, while applying a radial electric field to the edge plasma using a biasing probe. By calculating the intermittence parameter from floating potential measurements, we are able to identify a major low order rational surface and hence relate the probe measurements to the local value of the rotational transform. Based on the former, we are able to show that the poloidal plasma velocity (and hence radial electric field) has a significant radial structure that is clearly related to the rotational transform profile and in particular the lowest order rational surfaces in the range studied. The poloidal velocity is also affected by the edge biasing. The particle flux Γ was also found to exhibit a radial pattern, as did the flow shear suppression term ω E × B , but the relation of the former to the low-order rational surfaces was less clear. We surmise that this lack of direct correspondence is due to an unknown term in the turbulence evolution equation: the instability growth rate, γ . We make use of a reduced Magnetohydrodynamic turbulence model to interpret the results. Overall, a picture is obtained in which the plasma self-organizes towards a state with a clear radial pattern of the radial electric field, in line with expectations from some numerical studies describing the spontaneous formation of an ‘ E × B staircase’, consisting of alternating layers with fast and slow radial transport. In this state, the radial profiles of various quantities (density, temperature, pressure) will not be smooth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shape fluctuations and optical transition of He 2 * excimer tracers in superfluid He 4

Metastable He$_{2}^{*}$ excimer molecules have been utilized as tracer particles of the normal component in superfluid $^4$He (He II) which can be imaged via laser-induced fluorescence. These excimer molecules form tiny bubbles in He II and can bind to quantized vortices at sufficiently low temperatures, thereby allowing for direct visualization of vortex dynamics in an inviscid superfluid. However, the $a^{3}\Sigma^+_{u}$${\rightarrow}$$c^{3}\Sigma^+_{g}$ optical absorption line, which is responsible for the fluorescence imaging of the He$_{2}^{*}$ molecules, is controlled by fluctuations on the bubble shape, and its exact line profile is not known at low temperatures. In this paper, we present a bubble model for evaluating the surface fluctuation eigenmodes of the excimers in He II. The line profile of the $a^{3}\Sigma^+_{u}{\rightarrow}c^{3}\Sigma^+_{g}$ transition is calculated at different temperatures by considering both the zero-point and thermal fluctuations on the bubble shape. We show that, as the temperature drops from 2~K to 20 mK, the peak absorption strength is enhanced by a factor of about five, accompanying a blueshift of the peak location by about 2 nm. A double-peak line profile due to the rotational levels of the molecular core can be resolved. This bubble model also allows us to evaluate the stiffness of the He$_{2}^{*}$ bubbles and hence their diffusion constant in He II due to scattering off thermal phonons. Our results will aid the design of future experiments on imaging quantized vortices in He II using He$_{2}^{*}$ tracers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Engineering Redox Flow Battery Electrodes with Spatially Varying Porosity Using Non‐Solvent‐Induced Phase Separation

Redox flow batteries (RFBs) are a promising electrochemical platform for efficiently and reliably delivering electricity to the grid. Within the RFB, porous carbonaceous electrodes facilitate electrochemical reactions and distribute the flowing electrolyte. Tailoring electrode microstructure and surface area can improve RFB performance, lowering costs. Electrodes with spatially varying porosity may increase electrode utilization and provide surface area in reaction‐limited zones; however, the efficacy of such designs remains an open area of research. Herein, a non‐solvent‐induced phase‐separation (NIPS) technique that enables the reproducible synthesis of macrovoid‐free electrodes with well‐defined across‐thickness porosity gradients is described. The monotonically varying porosity profile is quantified and the physical properties and surface chemistries of porosity‐gradient electrodes are compared with macrovoid‐containing electrode, also synthesized by NIPS. Then, the electrochemical and fluid dynamic performance of the porosity‐gradient electrodes is evaluated, exploring the effect of changing the direction of the porosity gradient and benchmarking against the macrovoid‐containing electrode. Lastly, the performance is examined in a vanadium RFB, finding that the porosity‐gradient electrode outperforms the macrovoid electrode, is independent of gradient direction, and performs favorably compared to advanced electrodes in the contemporary literature. It is anticipated that the approach motivates further exploration of microstructurally tailored electrodes in electrochemical systems.

25 ENERGY STORAGE↗

Residual stresses and microstructure within Allvac 718Plus laser powder bed fusion bars

The residual stresses within Allvac 718Plus bars built using LBPF with different laser powers, speeds and table displacements were measured using diffraction and mechanical methods and modelled with the baseplate attached and then removed. The residual stress profiles within the bars from the top surface down through the bulk were all quite similar, becoming less tensile due to a strain hardening mechanism. Table displacement has the greatest impact on residual stress, decreasing with increasing displacement/powder layer thickness. There was good agreement amongst the modeling and measurements. The microstructures were examined and varied slightly with energy density with higher densities having larger grains and enhanced post solidification diffusion. We found energy density had minimal impact on the residuals stresses within the parameters to produce dense material.

36 MATERIALS SCIENCE↗

Chromium-coated cladding analysis under simulated LOCA burst conditions

Of the near-term accident-tolerant fuel concepts identified, chromium coatings have been shown to slow cladding oxidation without altering fuel system geometry or greatly affecting neutronic performance. To identify how coated cladding tubes perform under high-temperature accident conditions, pressurized-tube burst tests were conducted in the Severe Accident Test Station at Oak Ridge National Laboratory. To analyze how the coating affects cladding behavior, coated and uncoated cladding burst tests were simulated by using the BISON fuel performance code. Thermocouple data from these tests were fit into axial and azimuthal profiles and combined to generate 3D cladding surface temperatures, and pressure transducer data were compared until cladding failure. Additionally, the cladding temperatures at failure and the internal gas pressure evolution show relatively good agreement between the simulation and experiment results. Simulations were then extrapolated to demonstrate the coatings effectiveness to increasing the cladding burst temperature by using a parametric evaluation of the initial tube gas pressure. This work demonstrates the possibility of an increased cladding failure margin under transient conditions due to the addition of chromium coating, and more pertinently, an increased cladding failure margin from more accurate experimental characterization.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Initial hydration process of calcium silicates in Portland cement: A comprehensive comparison from molecular dynamics simulations

As the main components of Portland cement, calcium silicates show substantial differences in their hydration reactivity which have not been fully explained. A comprehensive comparison of the initial hydration process of calcium silicates, namely dicalcium silicate (C{sub 2}S) and tricalcium silicate (C{sub 3}S), was conducted using molecular dynamics simulations. The initial hydration process was divided into three stages using cut-off times of 0.001 ns and 3 ns. The hydration of M{sub 3}-C{sub 3}S (010) was more evident than that of β-C{sub 2}S (100), supported by the hydroxylation degree, radial distribution function, atomic density profile, etc. The coordination number of the surface Ca atoms might be the underlying reason for such a difference. Interactions between cement surfaces and water molecules were mainly characterised by solid OH bonding and Ca-water O bonding. Dissolution of Ca atoms was observed, although quite scarce, while no dissolution of Si atoms was observed.

36 MATERIALS SCIENCE↗

Nonlinear control of safety factor gradient in tokamaks using spatially variable electron cyclotron current drives

Active control of plasma properties may be necessary to achieve stable operation of next-generation tokamaks over large time scales. Such control algorithms can regulate the plasma properties to avoid the onset of magnetohydrodynamic (MHD) instabilities. For instance, the global and local properties of the safety factor profile are linked to the onset of neoclassical tearing modes (NTMs). This work proposes a model-based control approach for deterring/delaying NTMs through active modulation of one of the safety factor properties - the gradient of the safety factor profile at a particular rational safety factor surface. In particular, a novel control-oriented model for the local safety factor gradient is developed. The nonlinear control model is governed by a nonautonomous ordinary differential equation that accounts for a given rational safety factor surface’s spatial variation over time. Further, to improve the controllability of the spatially evolving parameter, the control model treats ECH&CD positions, along with noninductive powers, as controllable variables. A nonlinear control algorithm based on feedback linearization with optimization is synthesized to achieve the objective of regulating the safety factor gradient around a given target. The proposed algorithm allocates optimal ECH&CD positions, in addition to auxiliary powers, at each time instant as the rational safety factor surface drifts to locations with a low control authority. Stability guarantees of the proposed control law are also discussed in this work. The proposed algorithm is tested for a DIII-D tokamak scenario in nonlinear simulations carried out using the Control Oriented Transport SIMulator (COTSIM). Both fixed and moving ECH&CD cases are studied, and their outcomes are compared. Simulation results demonstrate that enthusiastic regulation of the safety factor gradient can be achieved during the ramp-up and flat-top phases of tokamak operation in both fixed and moving ECH&CD cases. However, real-time updates of ECH&CD positions can prevent the saturation of auxiliary powers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ion Size Effects on the Thermodynamic, Kinetic, and Mechanical Properties during Ion Exchange in Solid-State Electrolytes

Ion exchange offers a pathway to impose residual compressive stresses to mitigate the electro-chemo-mechanical cracking of solid-state electrolytes such as lithium lanthanum zirconium oxide. This study uses a coupled multiscale framework (integrating density functional theory (DFT), molecular dynamics (MD), and continuum modeling) to examine how exchange ion size influences stress, diffusion, fracture toughness, and electronic properties. Larger isovalent ions (Na + , Ag + , K + ) were exchanged with Li + , with DFT confirming their preference for octahedral 96h sites and a linear relationship between ion size and chemical free expansion coefficient. MD simulations reveal stress and concentration effects on exchange ion diffusivity at elevated temperatures, with Na + and Ag + maintaining favorable mobility while K + showing concentration-dependent clustering. Continuum modeling predicts the range of fracture strength improvements and the required ion exchange concentration profile. It was shown that a 5% surface exchange concentration can induce ∼0.6 GPa of surface compressive stress using Na + and ∼1.0 GPa of surface compressive stress using Ag + . On the other hand, larger ion exchange species may penalize Li + diffusivity by increasing the activation volume and activation energy. Interestingly, Na + has a negligible penalty on Li-ion diffusivity. The room temperature Li + ion diffusivity is reduced by ∼40% with Ag + ion exchange. Electronic band structure analysis shows no size-dependent change in the bandgap, though Ag + introduces localized defect states near the valence band maximum. This study highlights ion size as a key factor in optimizing LLZO properties, offering a framework to improve the solid-state battery performance.

Jagad, Harsh D. [Brown Univ., Providence, RI (Unit↗

The scaling of kinetic and transport behaviors in the solution-phase chemistry of a plasma–liquid interface

The reactions at a plasma–liquid interface often involve species such as the solvated electron or the hydroxyl radical, which initiate the reduction or oxidation of solution-phase reactants (so-called scavengers) or are consumed by their own second-order recombination. Here, the mathematical scaling of the reaction–diffusion equations at the interface is used to obtain a characteristic time that can be used to determine the transition from highly efficient scavenger reduction or oxidation to lower efficiencies due to transport limitations. The characteristic time (tc) is validated using numerical solutions of the reaction–diffusion equations. When the scavenger kinetics are faster than second-order recombination, this characteristic transition time scales proportionally with the scavenger diffusivity (Ds) and the square of the scavenger bulk concentration (SB) and inversely proportional to the electron flux (J) squared; that is, t c = D s S B 2 F 2 /J 2 , where F is Faraday’s constant. However, when the scavenger kinetics are comparable or slower than second-order recombination, this scaling does not hold. Extending this analysis to three-dimensional systems shows that the profile of the electron flux on the surface affects the spatial location where reactions are most effective. Finally, the assessment of the implications of these behaviors for the reactor design highlights how effectively controlling the electron flux and solution transport may be necessary to improve the efficiency of scavenger reactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Radiation induced non-linear oscillations in ITER baseline scenario plasmas in DIII-D

Abstract This work shows how the radiation brought about by metals or metal-equivalent radiators such as Kr and Xe produces non-linear dynamics on otherwise stationary β N flattops of DIII-D ITER Baseline Scenario demonstration discharges. The Kr and Xe gases are used to reproduce the radiative loss rates of W in present machines that operate at core temperatures much lower than the expected ITER temperature. Experiments on DIII-D with injection of Kr and Xe, as well as with sources of intrinsic metals reach the range of radiated fraction values expected in the ITER core and experience slow oscillations in temperature and radiated power. In many cases of high radiated fraction, the core temperature decreases enough for the safety factor profile to rise above the 1/1 rational surface, naturally eliminating sawteeth and occasionally producing a persistent helical core. The oscillations can be reproduced by a modified Lotka–Volterra system for temperature and radiated fraction if diffusion and noise are included, which indicates that the interplay between temperature and radiation can be the main cause of the cyclic nature of the system. A new physics based model which includes equations for temperature, density and input power can also reproduce the oscillations observed in the experiments. The present results suggest that the non-linearity of the system can be increased by the inclusion of the inherently non-linear alpha heating term, which is proportional to ∼ n e 2 T i 2 , and obtains oscillations in the model when added to an otherwise more stationary system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗