Engineering Papers⌕ Search

Engineering topics

Russell, D. A.

Publications and source records attributed to Russell, D. A..

At least 19 records

Correlation between the relative blob fraction and plasma parameters in NSTX

Here, we present a simple analysis of the fraction of blobs observed in the scrape-off layer of NSTX and calculates the correlation between this blob fraction and various plasma parameters. The measurements were made using the gas puff imaging (GPI) diagnostic on a database of 103 shots from 2010 with neutral beam power P NBI = 0–6 MW. The blob fraction ranged from f b = 0.1% to 4.8%, where f b is defined as the fraction of time the GPI signal was larger than three times its average value within 0–7.5 cm outside the separatrix. The blob fraction generally decreased with increasing neutral beam injection power and was lower on average in H-mode plasmas than Ohmic and L-mode plasmas. The blob fraction had its highest correlation with the measured poloidal turbulence velocity and its radial gradient just inside the separatrix. This is in part consistent with a model in which the blob generation rate increases with the poloidal velocity shear.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reduced-model scrape-off layer turbulence (nSOLT) simulations comparing three fueling scenarios

The 2D scrape-off-layer turbulence code (nSOLT) includes 1D Boltzmann neutral–plasma interactions, a model of divertor recycling (introduced here), and a fixed source of plasma concentrated at the core-side boundary. Three fueling methods are considered herein: (1) neutral injection in the far-SOL is accomplished by specifying the density of Franck–Condon distributed neutrals streaming in from the boundary. (2) Divertor recycling is modeled by injecting a fraction of the particle parallel flux in the scrape-off layer (SOL) back into the edge region as a source of plasma. (3) A constant source fuels the edge plasma from the core-side boundary to model pellet injection. For machine parameters (B, R m , and L // ) illustrative of the MAST-U device, and for a deuterium plasma, turbulent equilibria are obtained that share the same plasma fueling rate for each of the three fueling methods, with only one of the sources on in each case. In the presence of self-consistent turbulence, quasi-steady plasma and neutral (deuterium) profiles, fueling efficiencies, SOL transparencies, and heat flux widths are compared. Characteristics of the turbulent fluctuations, including skewness, cross-phases, and power spectra, are described. The calculated fueling efficiencies, SOL transparencies to neutral penetration, and many of the turbulent properties are remarkably similar for all three fueling methods despite significant differences in the plasma profiles. The nonlinear states of the three cases are dominated by separatrix-spanning vortex cells that control particle and heat losses into the SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reduction of blob-filament radial propagation by parallel variation of flows: Analysis of a gyrokinetic simulation

Data from the XGC1 gyrokinetic simulation is analyzed to understand the three-dimensional spatial structure and the radial propagation of blob-filaments generated by quasisteady turbulence in the tokamak edge pedestal and scrape-off layer plasma. Spontaneous toroidal flows vary in the poloidal direction and shear the filaments within a flux surface resulting in a structure that varies in the parallel direction. Here, this parallel structure allows the curvature and grad-B induced polarization charge density to be shorted out via parallel electron motion. As a result, it is found that the blob-filament radial velocity is significantly reduced from estimates which neglect parallel electron kinetics, broadly consistent with experimental observations. Conditions for when this charge shorting effect tends to dominate blob dynamics are derived and compared with the simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Blob wakes in NSTX

Transient small-scale structures were identified in the wake of blobs moving poloidally through the scrape-off layer of high-powered H-mode plasmas in NSTX, using the gas puff imaging (GPI) diagnostic. These blob wakes had a poloidal wavelength in the range of λ pol = 3.5 ± 0.7 cm, which is significantly smaller than the average blob scale of L pol ~12 cm, and the wakes had a poloidal velocity of V pol = 1.5 ± 1.0 km/s in the electron diamagnetic direction, which is opposite to the blob poloidal velocity in these shots. These wakes were radially localized 0–4 cm outside the separatrix and occurred within ~50 μs after the passage of a blob through the GPI field of view. The clearest wakes were seen when the GPI viewing angle was well aligned with the local B field line, as expected for such small-scale structures given the diagnostic geometry. In conclusion, a plausible theoretical interpretation of the wakes is discussed: the observed wakes share some features of drift waves and/or drift-Alfvén waves which could be excited by the blobs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A reduced model of neutral-plasma interactions in the edge and scrape-off-layer: Verification comparisons with kinetic Monte Carlo simulations

The 2D scrape-off-layer turbulence code (SOLT) is extended to include neutral-plasma interactions. A Boltzmann equation is derived for the evolution of the bi-normally-averaged neutral distribution function, G( x ,vx,t), in the radial dimension, and this evolution is included in the new code (nSOLT). Neutral-plasma interactions are mediated by charge-exchange (CX) and ionization rates based on poloidally-averaged plasma density and temperature. Good agreement is obtained between asymptotically stationary neutral density profiles from nSOLT simulations and those previously obtained from the Monte Carlo neutral transport code DEGAS 2, for time-averaged NSTX H-mode plasma profiles. The sensitivity of the nSOLT neutral profiles to atomic physics parameters, with and without CX physics, is included in the comparison. In addition, nSOLT simulations that evolve the plasma in 1D, using radial diffusion as a proxy for turbulent (blob) transport, illustrate the convergence to a self-consistent neutral-plasma equilibrium sustained by a neutral source at the far-SOL boundary and plasma heating in the core; equilibria consistent with typical NSTX Ohmic L-mode plasmas are described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mean flows and blob velocities in scrape-off layer (SOLT) simulations of an L-mode discharge on Alcator C-Mod

Two-dimensional scrape-off layer turbulence (SOLT) code simulations are compared with an L-mode discharge on the Alcator C-Mod tokamak [M. Greenwald, et al., Phys. Plasmas 21, 110501 (2014)]. Density and temperature profiles for the simulations were obtained by smoothly fitting Thomson scattering and mirror Langmuir probe (MLP) data from the shot. Simulations differing in turbulence intensity were obtained by varying a dissipation parameter. Mean flow profiles and density fluctuation amplitudes are consistent with those measured by MLP in the experiment and with a Fourier space diagnostic designed to measure poloidal phase velocity. Blob velocities in the simulations were determined from the correlation function for density fluctuations, as in the analysis of gas-puff-imaging (GPI) blobs in the experiment. In the simulations, it was found that larger blobs moved poloidally with the ExB flow velocity, vE , in the near-SOL, while smaller fluctuations moved with the group velocity of the dominant linear (interchange) mode, vE + 1/2 vdi, where vdi is the ion diamagnetic drift velocity. Comparisons are made with the measured GPI correlation velocity for the discharge. The saturation mechanisms operative in the simulation of the discharge are also discussed. It is found that neither sheared flow nor pressure gradient modification can be excluded as saturation mechanisms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Blob-hole correlation model for edge turbulence and comparisons with NSTX gas puff imaging data [Blob-hole correlation model for edge turbulence and comparisons with NSTX GPI data]

We report that gas puff imaging (GPI) observations made in NSTX have revealed two-point spatial correlations of edge and scrape-off layer turbulence in the plane perpendicular to the magnetic field. A common feature is the occurrence of dipole-like patterns with significant regions of negative correlation. In this paper, we explore the possibility that these dipole patterns may be due to blob-hole pairs. Statistical methods are applied to determine the two-point spatial correlation that results from a model of blob-hole pair formation. It is shown that the model produces dipole correlation patterns that are qualitatively similar to the GPI data in several respects. Effects of the reference location (confined surfaces or scrape-off layer), a superimposed random background, hole velocity and lifetime, and background sheared flows are explored and discussed with respect to experimental observations. Additional analysis of the experimental GPI dataset is performed to further test this blob-hole correlation model. A time delay two-point spatial correlation study did not reveal inward propagation of the negative correlation structures that were postulated to correspond to holes in the data nor did it suggest that the negative correlation structures are due to neutral shadowing. However, tracking of the highest and lowest values (extrema) of the normalized GPI fluctuations shows strong evidence for mean inward propagation of minima and outward propagation of maxima, in qualitative agreement with theoretical expectations. Finally, other properties of the experimentally observed extrema are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theory based scaling of edge turbulence and implications for the scrape-off layer width

Turbulence and plasma parameter data from the National Spherical Torus Experiment (NSTX) is examined and interpreted based on various theoretical estimates. In particular, quantities of interest for assessing the role of turbulent transport on the midplane scrape-off layer heat flux width are assessed. Because most turbulence quantities exhibit large scatter and little scaling within a given operation mode, this paper focuses on length and time scales and dimensionless parameters between operational modes including Ohmic, low (L), and high (H) modes using a large NSTX edge turbulence database. These are compared with theoretical estimates for drift and interchange rates, profile modification saturation levels, a resistive ballooning condition, and dimensionless parameters characterizing L and H mode conditions. It is argued that the underlying instability physics governing edge turbulence in different operational modes is, in fact, similar, and is consistent with curvature-driven drift ballooning. Saturation physics, however, is dependent on the operational mode. Five dimensionless parameters for drift-interchange turbulence are obtained and employed to assess the importance of turbulence in setting the scrape-off layer heat flux width λ q and its scaling. An explicit proportionality of the width λ q to the safety factor and major radius (qR) is obtained under these conditions. Lastly, quantitative estimates and reduced model numerical simulations suggest that the turbulence mechanism is not negligible in determining λ q in NSTX, at least for high plasma current discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mean flows and blob velocities in scrape-off layer (SOLT) simulations of an L-mode discharge on Alcator C-Mod

Two-dimensional scrape-off layer turbulence (SOLT) code simulations are compared with an L-mode discharge on the Alcator C-Mod tokamak [M. Greenwald, et al., Phys. Plasmas 21, 110501 (2014)]. Density and temperature profiles for the simulations were obtained by smoothly fitting Thomson scattering and mirror Langmuir probe (MLP) data from the shot. Simulations differing in turbulence intensity were obtained by varying a dissipation parameter. Mean flow profiles and density fluctuation amplitudes are consistent with those measured by MLP in the experiment and with a Fourier space diagnostic designed to measure poloidal phase velocity. Blob velocities in the simulations were determined from the correlation function for density fluctuations, as in the analysis of gas-puff-imaging (GPI) blobs in the experiment. In the simulations, it was found that larger blobs moved poloidally with the ExB flow velocity, v E , in the near-SOL, while smaller fluctuations moved with the group velocity of the dominant linear (interchange) mode, v E + 1/2 v di , where v di is the ion diamagnetic drift velocity. Comparisons are made with the measured GPI correlation velocity for the discharge. The saturation mechanisms operative in the simulation of the discharge are also discussed. In conclusion, it is found that neither sheared flow nor pressure gradient modification can be excluded as saturation mechanisms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Blob structure and motion in the edge and SOL of NSTX

Here, the structure and motion of discrete plasma blobs (a.k.a. filaments) in the edge and scrape-off layer of NSTX is studied for representative Ohmic and H-mode discharges. Individual blobs were tracked in the 2D radial versus poloidal plane using data from the gas puff imaging diagnostic taken at 400 000 frames s -1 . A database of blob amplitude, size, ellipticity, tilt, and velocity was obtained for ~45 000 individual blobs. Empirical relationships between various properties are described, e.g. blob speed versus amplitude and blob tilt versus ellipticity. The blob velocities are also compared with analytic models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Turbulent transport regimes and the scrape-off layer heat flux width

Understanding the responsible mechanisms and resulting scaling of the scrape-off layer (SOL) heat flux width is important for predicting viable operating regimes in future tokamaks, and for seeking possible mitigation schemes. Here in this paper, we present a qualitative and conceptual framework for understanding various regimes of edge/SOL turbulence and the role of turbulent transport as the mechanism for establishing the SOL heat flux width. Relevant considerations include the type and spectral characteristics of underlying instabilities, the location of the gradient drive relative to the SOL, the nonlinear saturation mechanism, and the parallel heat transport regime. We find a heat flux width scaling with major radius R that is generally positive, consistent with previous findings [J. W. Connor et al., Nucl. Fusion 39, 169 (1999)]. The possible relationship of turbulence mechanisms to the neoclassical orbit width or heuristic drift mechanism in low (L) -mode and high (H) -mode plasmas is considered, together with implications for future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiple Beam Two-Plasmon Decay: Linear Threshold to Nonlinear Saturation in Three Dimensions

The linear stability of multiple coherent laser beams with respect to two-plasmon-decay instability in an inhomogeneous plasma in three dimensions has been determined. Cooperation between beams leads to absolute instability of long-wavelength decays, while shorter-wavelength shared waves are shown to saturate convectively. The multibeam, in its absolutely unstable form, has the lowest threshold for most cases considered. Finally, nonlinear calculations using a three-dimensional extended Zakharov model show that Langmuir turbulence created by the absolute instability modifies the convective saturation of the shorter-wavelength modes, which are seen to dominate at late times.

42 ENGINEERING↗

Development of a Practical Methodology for Elastic-Plastic and Fully Plastic Fatigue Crack Growth

A practical engineering methodology has been developed to analyze and predict fatigue crack growth rates under elastic-plastic and fully plastic conditions. The methodology employs the closure-corrected effective range of the J-integral, delta J(sub eff) as the governing parameter. The methodology contains original and literature J and delta J solutions for specific geometries, along with general methods for estimating J for other geometries and other loading conditions, including combined mechanical loading and combined primary and secondary loading. The methodology also contains specific practical algorithms that translate a J solution into a prediction of fatigue crack growth rate or life, including methods for determining crack opening levels, crack instability conditions, and material properties. A critical core subset of the J solutions and the practical algorithms has been implemented into independent elastic-plastic NASGRO modules. All components of the entire methodology, including the NASGRO modules, have been verified through analysis and experiment, and limits of applicability have been identified.

McClung, R. C.↗

Guidelines for Proof Test Analysis

These guidelines integrate state-of-the-art elastic-plastic fracture mechanics (EPFM) and proof test implementation issues into a comprehensive proof test analysis procedure in the form of a road map which identifies the types of data, fracture mechanics based parameters, and calculations needed to perform flaw screening and minimum proof load analyses of fracture critical components. Worked examples are presented to illustrate the application of the road map to proof test analysis. The state-of-the art fracture technology employed in these guidelines is based on the EPFM parameter, J, and a pictorial representation of a J fracture analysis, called the failure assessment diagram (FAD) approach. The recommended fracture technology is validated using finite element J results, and laboratory and hardware fracture test results on the nickel-based superalloy Inconel 718, the aluminum alloy 2024-T3511, and ferritic pressure vessel steels. In all cases the laboratory specimens and hardware failed by ductile mechanisms. Advanced proof test analyses involving probability analysis and multiple-cycle proof testing (MCPT) are addressed. Finally, recommendations are provided on how to account for the effects of the proof test overload on subsequent service fatigue and fracture behaviors.

Chell, G. G.↗

A Comparison of Single-Cycle Versus Multiple-Cycle Proof Testing Strategies

Single-cycle and multiple-cycle proof testing (SCPT and MCPT) strategies for reusable aerospace propulsion system components are critically evaluated and compared from a rigorous elastic-plastic fracture mechanics perspective. Earlier MCPT studies are briefly reviewed. New J-integral estimation methods for semielliptical surface cracks and cracks at notches are derived and validated. Engineering methods are developed to characterize crack growth rates during elastic-plastic fatigue crack growth (FCG) and the tear-fatigue interaction near instability. Surface crack growth experiments are conducted with Inconel 718 to characterize tearing resistance, FCG under small-scale yielding and elastic-plastic conditions, and crack growth during simulated MCPT. Fractography and acoustic emission studies provide additional insight. The relative merits of SCPT and MCPT are directly compared using a probabilistic analysis linked with an elastic-plastic crack growth computer code. The conditional probability of failure in service is computed for a population of components that have survived a previous proof test, based on an assumed distribution of initial crack depths. Parameter studies investigate the influence of proof factor, tearing resistance, crack shape, initial crack depth distribution, and notches on the MCPT versus SCPT comparison. The parameter studies provide a rational basis to formulate conclusions about the relative advantages and disadvantages of SCPT and MCPT. Practical engineering guidelines are proposed to help select the optimum proof test protocol in a given application.

McClung, R. C.↗

Development of a Practical Methodology for Elastic-Plastic and Fully Plastic Fatigue Crack Growth

A practical engineering methodology has been developed to analyze and predict fatigue crack growth rates under elastic-plastic and fully plastic conditions. The methodology employs the closure-corrected effective range of the J-integral, (Delta)J(sub eff), as the governing parameter. The methodology contains original and literature J and (Delta)J solutions for specific geometries, along with general methods for estimating J for other geometries and other loading conditions, including combined mechanical loading and combined primary and secondary loading. The methodology also contains specific practical algorithms that translate a J solution into a prediction of fatigue crack growth rate or life, including methods for determining crack opening levels, crack instability conditions, and material properties. A critical core subset of the J solutions and the practical algorithms has been implemented into independent elastic-plastic NASGRO modules. All components of the entire methodology, including the NASGRO modules, have been verified through analysis and experiment, and limits of applicability have been identified.

McClung, R. C.↗

NASA Contractor Report: Guidelines for Proof Test Analysis

These Guidelines integrate state-of-the-art Elastic-Plastic Fracture Mechanics (EPFM) and proof test implementation issues into a comprehensive proof test analysis procedure in the form of a Road Map which identifies the types of data, fracture mechanics based parameters, and calculations needed to perform flaw screening and minimum proof load analyses of fracture critical components. Worked examples are presented to illustrate the application of the Road Map to proof test analysis. The state-of-the-art fracture technology employed in these Guidelines is based on the EPFM parameter, J, and a pictorial representation of a J fracture analysis, called the Failure Assessment Diagram (FAD) approach. The recommended fracture technology is validated using finite element J results, and laboratory and hardware fracture test results on the nickel-based superalloy IN-718, the aluminum alloy 2024-T351 1, and ferritic pressure vessel steels. In all cases the laboratory specimens and hardware failed by ductile mechanisms. Advanced proof test analyses involving probability analysis and Multiple Cycle Proof Testing (MCPT) are addressed. Finally, recommendations are provided on to how to account for the effects of the proof test overload on subsequent service fatigue and fracture behaviors.

Chell, G. G.↗

A Comparison of Single-Cycle Versus Multiple-Cycle Proof Testing Strategies

Single-cycle and multiple-cycle proof testing (SCPT and MCPT) strategies for reusable aerospace propulsion system components are critically evaluated and compared from a rigorous elastic-plastic fracture mechanics perspective. Earlier MCPT studies are briefly reviewed. New J-integral estimation methods for semi-elliptical surface cracks and cracks at notches are derived and validated. Engineering methods are developed to characterize crack growth rates during elastic-plastic fatigue crack growth (FCG) and the tear-fatigue interaction near instability. Surface crack growth experiments are conducted with Inconel 718 to characterize tearing resistance, FCG under small-scale yielding and elastic-plastic conditions, and crack growth during simulated MCPT. Fractography and acoustic emission studies provide additional insight. The relative merits of SCPT and MCPT are directly compared using a probabilistic analysis linked with an elastic-plastic crack growth computer code. The conditional probability of failure in service is computed for a population of components that have survived a previous proof test, based on an assumed distribution of initial crack depths. Parameter studies investigate the influence of proof factor, tearing resistance, crack shape, initial crack depth distribution, and notches on the MCPT vs. SCPT comparison. The parameter studies provide a rational basis to formulate conclusions about the relative advantages and disadvantages of SCPT and MCPT. Practical engineering guidelines are proposed to help select the optimum proof test protocol in a given application.

McClung, R. C.↗