Engineering topics

Chen, J. L.

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

12 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↗

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.↗