Engineering topics

Chen, J. L.

Publications and source records attributed to Chen, J. L..

14 records

DOE OSTI2023

Long-pulse high-performance H-mode plasmas achieved on EAST

A record duration of a 310 s H-mode plasma (H 98y2 ~ 1.3, n e /n GW ~ 0.7, f BS > 50%) has been recently achieved on experimental advanced superconducting tokamak (EAST) with metal walls, exploiting the device's improved long-pulse capabilities. The experiment demonstrates good control of tungsten concentration, core/edge MHD stability, and particle and heat exhaust with an ITER-like tungsten divertor and zero injected torque, establishing a milestone on the path to steady-state long-pulse high-performance scenarios in support of ITER and CFETR. Important synergistic effects are leveraged toward this result, which relies purely on radio frequency (RF) powers for heating and current drive (H&CD). On-axis electron cyclotron heating enhances the H&CD efficiency from lower hybrid wave injection, increasing confinement quality and enabling fully non-inductive operation at high density (n e /n GW ~ 70%) and high poloidal beta (β P ~ 2.5). A small-amplitude grassy edge localized mode regime facilitates the RF power coupling to the H-mode edge and reduces divertor sputtering/erosion. The high energy confinement quality (H 98y2 ~ 1.3) is achieved with the experimental and simulated results pointing to the strong effect of Shafranov shift on turbulence. Transport analysis suggests that trapped electron modes dominate in the core region during the record discharge. The detailed physics processes (RF synergy, core-edge integration, confinement properties, etc.) of the steady-state operation will be illustrated in the content. In the future, EAST will aim at accessing more relevant dimensionless parameters to develop long-pulse high-performance plasma toward ITER and CFETR steady-state advanced operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Characterization of beam ion loss in high poloidal beta regime on EAST

A critical issue for achieving the integrated operation of steady-state long-pulse high-confinement (H-mode) plasmas on experimental advanced superconducting tokamak (EAST) is to improve beam ion population confinement during neutral beam injection (NBI). To study the characterization of beam ion loss and improve beam ion confinement, the steady-state long pulse scenario discharges were conducted on EAST (β p $\geqslant$ 2.0, β N $\geqslant$ 1.7, q 95 $\geqslant$ 6.7 and H 98y2 $\geqslant$ 1.1) with NBI heating. Based on neutron yield, the beam voltage and line-averaged electron density were adjusted from 50 kV to 60 kV and 4.4 × 10 19 m –3 to 5.0 × 10 19 m –3 , respectively. The results show that the dominant mechanisms of beam ion loss are shine-through loss, prompt loss, and stochastic ripple loss. The shine-through loss fraction is determined by initial velocity, flight time and entire beam path. The change in prompt loss fraction is caused by the change in the deposition of beam ions. The change in stochastic ripple loss fraction is caused by the change in the initial fraction of trapped-confined ions. Detailed physics shows that the prompt loss fraction during counter-Ip injections (~45%) is far larger than during co-Ip injections (~5%) due to the finite orbit width. The lost ions are mainly deposited on the lower divertor or below the midplane since the direction of magnetic drift is vertical down. The orbit types of prompt loss during counter-Ip injections are mainly trapped-lost and ctr-passing lost. To minimize the prompt loss fraction during counter-Ip injections, a reversed Ip configuration (rev-Ip) discharge #94758 was conducted. The result suggests that the beam ion wall load fraction during counter-Ip tangential injection (~3%) is far lower than that in normal Ip configuration (nor-Ip) discharge #94820. It is also found that the confinement of beam ion population in the counter-Ip injection #94758 was greatly improved when compared to #94820. Finally, this study can provide unique support for the improvement of beam ion population confinement and for the performance evaluation of the NBI system on EAST and future tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Edge-localized-mode simulation in CFETR steady-state scenario

We report the EPED1 model and self-consistent core-pedestal coupling in integrated modeling are used to design the pedestal structure of the China Fusion Engineering Testing Reactor (CFETR) steady-state scenario. The key parameters, such as β p and q 95 , are based on the grassy edge-localized-mode (ELM) experimental database. In this work, we use the BOUT++ six-field two-fluid code to simulate the onset of the ELM in the CFETR steady-state scenario. The ELM size is around 0.2% in nonlinear simulations, which is in the experimental range of the grassy ELM discharges, 0.1%–1% observed in multiple tokamak devices. Linear and nonlinear simulations show that the dominant high-n ballooning modes peak around n = 40. Compared to type-I ELM crashing dynamics, grassy ELM crashing has a smaller initial crash and is then followed by three phases of turbulence spreading, which are dominated by multi-modes, a high-n mode of n = 45 and low-n mode of n = 5, respectively. In contras to type-I ELM, the perturbation of the high-n mode has a narrow width around ψ = 0.95, and magnetic island formation and reconnection occur only beyond ψ = 0.95, leading to a small initial crash. Mode–mode interaction in the multi-mode coexistence stage stops the growth of individual modes and reduces the transport of particles and heat, and these are the two reasons why the ELM size is small. In–out asymmetry of transient heat flux with a ratio of E out /E in = 3.5 is found during grassy ELM crash. The rise and delay times of the heat flux match the calculation from the free-streaming model. To evaluate the erosion of the divertor target, the energy fluence at the outer divertor target is calculated, which is 0.029 MJ m -2 , 5.5 times smaller than the tungsten melting limit 0.16 MJ m -2 . The calculated energy fluency still follows the experimental scaling law from type-I ELM experiments. The fluctuation eddies in the toroidal direction show a filament structure at the outer mid-plane. Parallel heat flux patterns with a toroidal mode number n = 10 are found at the outer divertor with an amplitude of 680 MW m -2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Theoretical analysis of key factors achieving reversed magnetic shear q-profiles sustained with lower hybrid waves on EAST

The phase space analysis technique is applied to provide new insights into the fully non-inductive lower hybrid current driven (LHCD) discharges on EAST (Garofalo et al 2017 Nucl. Fusion 57 076037). The analysis shows that there are bounded and unbounded topologies of lower hybrid (LH) waves in phase space. For typical parameters on EAST, the propagation domain for 4.6 GHz LH waves is bounded, while unbounded for 2.45 GHz LH waves and one of the conditions is recognized to be good for achieving an off-axis current profile driven by 4.6 GHz LH waves on EAST. The parametric analysis on the potential power deposition (PPD) region for those experiments demonstrates that the reversed magnetic shear dominates the wave behavior, and confines the LH power absorption to the far off-axis region as long as it occurs in the 4.6 GHz dominated LHCD discharge. GENRAY/CQL3D simulations also confirm this effect of reversed shear on the power absorption of 4.6 GHz waves. A bounded propagation domain and a positive feedback loop between magnetic shear reversal and off-axis LHCD profile could explain the LH-only-sustained strong reversed magnetic shear observed in experiments on multiple machines. In contrast, there is no such effect for 2.45 GHz waves on EAST since the PPD region is less sensitive to the change of reversed shear due to the unbounded propagation domain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Advances in physics understanding of high poloidal beta regime toward steady-state operation of CFETR

Experimental and modeling investigations of the high βp scenarios on the DIII-D and EAST tokamaks show advantages in high energy confinement, avoidance of n = 1 MHD, and core-edge integration with reduced heat flux, making this scenario an attractive option for CFETR steady-state operation. Experiments show that plasmas with high confinement and high density can be achieved with neutral beam injection on DIII-D (β p ~ 2.2, β N ~ 3.5, f BS ~ 50%, f Gw ~ 1.0, H 98y2 ~ 1.5) and pure RF power on EAST (β P ~ 2.0, β N ~ 1.6, f BS ~ 50%, f Gw ~ 0.8, H 98y2 > 1.3). By tailoring the current density profile, a q-profile with local (off-axis) negative shear is achieved that yields improved confinement and MHD stability. Transport analysis and simulation suggest that the combination of high density gradient and high Shafranov shift allows turbulence stabilization and higher confinement. Using on-axis ECH injection, tungsten accumulation is avoided on EAST, and this is reproduced in modeling. Reduced heat flux (by >40%) and maintenance of high core confinement is achieved with active feedback control of the radiated divertor, an important result for long pulse operation in tokamaks. In conclusion, the improved physics understanding and validated modeling tools are used to design a 1GW steady-state scenario for CFETR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Grassy ELM regime at low pedestal collisionality in high-power tokamak plasma

Spontaneous mitigation of giant ELMs and appearance of grassy ELMs have been observed repeatedly at low pedestal collisionality (${\nu }_{\mathrm{e},\mathrm{p}\mathrm{e}\mathrm{d}}^{{\ast}}\sim 0.15$) in the high-power (P inj > 13 MW) hybrid scenario in DIII-D. Higher β p and higher q 95 appear to be beneficial to achieving the grassy ELM regime. The grassy ELM H-mode plasma shows high energy confinement performance (H 98y2 up to over 1.5) under the conditions of high neutral beam torque and high core rotation. The pedestal width appears to exceed the EPED1.0 model prediction by more than 50%. Pedestal stability analysis performed with the ELITE code indicates that the stability against low-n kink/peeling modes is improved with increased auxiliary heating power and the operational point in the grassy ELM regime is located near the ballooning boundary. The pedestal stability characteristics during the grassy ELM crashes have been investigated in comparison with the giant ELM crashes based on plasma profiles experimentally measured with high resolution and accuracy. It has been found that the underlying mechanism for the observed small-amplitude ELM crashes is mainly the expansion of the ballooning stability boundary induced by an initial radially localized collapse in the pedestal, which helps to stop the growth of instabilities and further collapse of the pedestal. The effect of electron density pedestal on mitigating edge localized instabilities has been analyzed by numerical simulation, suggesting that the electron density pedestal characterized by high n e,sep /n e,ped and low density gradient helps to stabilize peeling-ballooning modes because of a low pressure pedestal gradient and to lower the ballooning boundary mainly because of a low ion diamagnetic frequency in the pedestal region, thus triggering ballooning instabilities and producing the intrinsic grassy ELMs. Numerical simulation of the Chinese fusion engineering test reactor (CFETR) with the SOLPS code indicates that the separatrix density might be insensitive to the electron diffusivity in the pedestal region and increase with the power flowing from the core region to the edge region. Furthermore, pedestal stability analysis suggests that the flat density pedestal with high separatrix density obtained in the high-power plasma in CFETR would make the operational point close to the ballooning boundary, which is considered to help destabilize ballooning instabilities and facilitate the access to the grassy ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Progress in extending high poloidal beta scenarios on DIII-D towards a steady-state fusion reactor and impact of energetic particles

To prepare for steady-state operation of future fusion reactors (e.g. the International Thermonuclear Experimental Reactor and China Fusion Engineering Test Reactor (CFETR)), experiments on DIII-D have extended the high poloidal beta ( β P ) scenario to reactor-relevant edge safety factor q 95 ~ 6.0, while maintaining a large-radius internal transport barrier (ITB) using negative magnetic shear. Excellent energy confinement quality ( H 98y2 > 1.5) is sustained at high normalized beta ( β N ~ 3.5). This high-performance ITB state with Greenwald density fraction near 100% and q min ≥ 3 is achieved with toroidal plasma rotation V tor ~ 0 at ρ ≥ 0.6. This is a key result for reactors expected to have low V tor . At high β P (≥1.9), large Shafranov shift can stabilize turbulence leading to a high confinement state with a low pedestal and an ITB. At lower β P (<1.9), negative magnetic shear in the plasma core contributes to turbulence suppression and can compensate for reduced Shafranov shift to continue to access a large-radius ITB and excellent confinement with low V tor , consistent with the results of gyrofluid transport simulations. These high- β P cases are characterized by weak/no Alfvén eigenmodes (a.e.) and classical fast-ion transport. At high density, the fast-ion deceleration time decreases and Δ β fast is lower; these reduce a.e. drive. The reverse-shear Alfvén eigenmodes are weaker or stable because the negative magnetic shear region is located at higher radius, away from the peaked fast-ion profile. Resistive wall modes can be a limitation at simultaneous high β N , low internal inductance, and low rotation. Analysis suggests that additional off-axis external current drive could provide a more stable path at reduced q 95 . Based on a DIII-D high- β P plasma with large-radius ITB, two scenarios are proposed for CFETR Q = 5 steady-state operation with ~1 GW fusion power: a lower-$l_i$($l_i$ ~ 0.66) and a higher-$l_i$($l_i$ ~ 0.75) case. Using a Landau closure model, multiple energetic particle (EP) effects on the a.e. stability are analyzed modifying the growth rate of the a.e.s triggered by the neutral-beam-injection EPs and alpha particles, although the stabilizing/destabilizing effect is weak for the cases analyzed. The stabilizing effects of the combined EP species β , energy, and density profile in CFETR need further investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2006

Revised Atmospheric Angular Momentum Series Related to Earth's Variable Rotation under Consideration of Surface Topography

The atmospheric angular momentum is closely related to variations in the Earth rotation. The atmospheric excitation function (AEF), or namely atmospheric effective angular momentum function, is introduced in studying the atmospheric excitation of the Earth's variable rotation. It may be separated into two portions, i.e, the "wind" terms due to the atmospheric motion relative to the mantle and the "pressure" terms due to the variations of atmospheric mass distribution evident through surface pressure changes. The AEF wind terms during the period of 1948-2004 are re-processed from the NCEP/NCAR (National Centers for Environmental Prediction-National Center for Atmospheric Research) reanalysis 6-hourly wind and pressure fields. Some previous calculations were approximate, in that the wind terms were integrated from an isobaric lower boundary of 1000 hPa. To consider the surface topography effect, however, the AEF is computed by integration using the winds from the Earth's surface to 10 hPa, the top atmospheric model level, instead of from 1000 hPa. For these two cases, only a minor difference, equivalent to approx. 0.004 milliseconds in length-of-day variation, exists with respect to the axial wind term. However, considerable differences, equivalent to 5-6 milliarcseconds in polar motion, are found regarding equatorial wind terms. We further compare the total equatorial AEF (with and without the topographic effect) with the polar motion excitation function (PMEF) during the period of 1980-2003. The equatorial AEF gets generally closer to the PMEF, and improved coherences are found between them when the topography effect is included. Keywords: Atmospheric angular momentum, Atmospheric excitation function, Earth rotation, Topography, Wind, Pressure.

Zhou, Y. H.↗

NASA NTRS2006

Terrestrial Water Mass Load Changes from Gravity Recovery and Climate Experiment (GRACE)

Recent studies show that data from the Gravity Recovery and Climate Experiment (GRACE) is promising for basin- to global-scale water cycle research. This study provides varied assessments of errors associated with GRACE water storage estimates. Thirteen monthly GRACE gravity solutions from August 2002 to December 2004 are examined, along with synthesized GRACE gravity fields for the same period that incorporate simulated errors. The synthetic GRACE fields are calculated using numerical climate models and GRACE internal error estimates. We consider the influence of measurement noise, spatial leakage error, and atmospheric and ocean dealiasing (AOD) model error as the major contributors to the error budget. Leakage error arises from the limited range of GRACE spherical harmonics not corrupted by noise. AOD model error is due to imperfect correction for atmosphere and ocean mass redistribution applied during GRACE processing. Four methods of forming water storage estimates from GRACE spherical harmonics (four different basin filters) are applied to both GRACE and synthetic data. Two basin filters use Gaussian smoothing, and the other two are dynamic basin filters which use knowledge of geographical locations where water storage variations are expected. Global maps of measurement noise, leakage error, and AOD model errors are estimated for each basin filter. Dynamic basin filters yield the smallest errors and highest signal-to-noise ratio. Within 12 selected basins, GRACE and synthetic data show similar amplitudes of water storage change. Using 53 river basins, covering most of Earth's land surface excluding Antarctica and Greenland, we document how error changes with basin size, latitude, and shape. Leakage error is most affected by basin size and latitude, and AOD model error is most dependent on basin latitude.

Seo, K.-W.↗

NASA NTRS1998

New Estimates of Hydrological and Oceanic Excitations of Variations of Earth's Rotation, Geocenter and Gravitational Field

Hydrological mass transport in the geophysical fluids of the atmosphere-hydrosphere-solid Earth surface system can excite Earth's rotational variations in both length-of-day and polar motion. These effects can be computed in terms of the hydrological angular momentum by proper integration of global meteorological data. We do so using the 40-year NCEP data and the 18-year NASA GEOS-1 data, where the precipitation and evapotranspiration budgets are computed via the water mass balance of the atmosphere based on Oki et al.'s (1995) algorithm. This hydrological mass redistribution will also cause geocenter motion and changes in Earth's gravitational field, which are similarly computed using the same data sets. Corresponding geodynamic effects due to the oceanic mass transports (i.e. oceanic angular momentum and ocean-induced geocenter/gravity changes) have also been computed in a similar manner. We here compare two independent sets of the result from: (1) non-steric ocean surface topography observations based on Topex/Poseidon, and (2) the model output of the mass field by the Parallel Ocean Climate Model. Finally, the hydrological and the oceanic time series are combined in an effort to better explain the observed non-atmospheric effects. The latter are obtained by subtracting the atmospheric angular momentum from Earth rotation observations, and the atmosphere- induced geocenter/gravity effects from corresponding geodetic observations, both using the above-mentioned atmospheric data sets.

Chao, Benjamin F.↗

NASA NTRS1994

Micromechanical characterization of nonlinear behavior of advanced polymer matrix composites

Due to the presence of curing stresses and oriented crystalline structures in the matrix of polymer matrix fiber composites, the in situ nonlinear properties of the matrix are expected to be rather different from those of the bulk resin. A plane stress micromechanical model was developed to retrieve the in situ elastic-plastic properties of Narmco 5260 and Amoco 8320 matrices from measured elastic-plastic properties of IM7/5260 and IM7/8320 advance composites. In the micromechanical model, the fiber was assumed to be orthotropically elastic and the matrix to be orthotropic in elastic and plastic properties. The results indicate that both in situ elastic and plastic properties of the matrices are orthotropic.

Gates, Thomas S.↗

NASA NTRS1988

Redundant manipulators for momentum compensation in a micro-gravity environment

This paper is concerned with the implementation and assessment of joint motion management strategies for kinematically redundant robotic manipulators operating in the micro-gravity environment of Space Station. These robots must be capable of conducting experiments and manufacturing processes without disturbing the micro-gravity environment through base reactions/motions. The redundant degrees of freedom of the manipulator permit the inverse kinematic problem to be solved simultaneously with the minimization of a cost function. The cost function in this case is the weighted sum of the squares of the base forces and moments and is minimized over discrete time segments. The Generalized Inverse Method and Rayleigh Ritz technique are used to solve the combined optimization/inverse kinematics problem. Numerical examples include various robotic configurations and degrees of manipulator redundancy.

Quinn, R. D.↗

NASA NTRS1988

Modeling creep behavior of fiber composites

A micromechanical model for the creep behavior of fiber composites is developed based on a typical cell consisting of a fiber and the surrounding matrix. The fiber is assumed to be linearly elastic and the matrix nonlinearly viscous. The creep strain rate in the matrix is assumed to be a function of stress. The nominal stress-strain relations are derived in the form of differential equations which are solved numerically for off-axis specimens under uniaxial loading. A potential function and the associated effective stress and effective creep strain rates are introduced to simplify the orthotropic relations.

Chen, J. L.↗

NASA NTRS1986

Optimum structural sizing of conventional cantilever and joined wing configurations using equivalent beam models

The present paper describes an approach for the optimum sizing of single and joined wing structures that is based on representing the built-up finite element model of the structure by an equivalent beam model. The low order beam model is computationally more efficient in an environment that requires repetitive analysis of several trial designs. The design procedure is implemented in a computer program that requires geometry and loading data typically available from an aerodynamic synthesis program, to create the finite element model of the lifting surface and an equivalent beam model. A fully stressed design procedure is used to obtain rapid estimates of the optimum structural weight for the beam model for a given geometry, and a qualitative description of the material distribution over the wing structure. The synthesis procedure is demonstrated for representative single wing and joined wing structures.

Hajela, P.↗