Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability regime”

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 19 records

Mode spectrum characteristics and onset of the low-shear MHD stability regime

Equilibria with extended regions of weak magnetic shear, including some tokamak scenarios and stellarators, can be susceptible to pressure-driven internal MHD instabilities even though there is no mode rational surface in the plasma. Nonresonant modes, in particular, can have properties that are unattractive for confinement, including displacing substantial volumes of the plasma and leading to more efficient pressure gradient flattening in the nonlinear regime. The onset and linear properties of the low shear stability regime are examined using the initial-value, extended-MHD code M3D-C. For monotonic q-profiles, we demonstrate a clear correlation between the convergents associated with the continued fraction representation of q 0 and the spectrum of unstable modes. Nonresonant modes are observed to be destabilized preferentially to any other resonant instability with the same toroidal mode number when n > 1. Using the observed connection between the spectrum associated with q 0 and the overall equilibrium stability properties, we suggest a technique for reducing the uncertainty on both q 0 and magnetic shear in the core region, obtained either via measurement or through the analysis and reconstruction of experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Behavior and mechanisms of Doppler wind lidar error in varying stability regimes

Abstract. Wind lidars are widespread and important tools in atmospheric observations. An intrinsic part of lidar measurement error is due to atmospheric variability in the remote-sensing scan volume. This study describes and quantifies the distribution of measurement error due to turbulence in varying atmospheric stability. While the lidar error model is general, we demonstrate the approach using large ensembles of virtual WindCube V2 lidar performing a profiling Doppler-beam-swinging scan in quasi-stationary large-eddy simulations (LESs) of convective and stable boundary layers. Error trends vary with the stability regime, time averaging of results, and observation height. A systematic analysis of the observation error explains dominant mechanisms and supports the findings of the empirical results. Treating the error under a random variable framework allows for informed predictions about the effect of different configurations or conditions on lidar performance. Convective conditions are most prone to large errors (up to 1.5 m s−1 in 1 Hz wind speed in strong convection), driven by the large vertical velocity variances in convective conditions and the high elevation angle of the scanning beams (62∘). Range-gate weighting induces a negative bias into the horizontal wind speeds near the surface shear layer (−0.2 m s−1 in the stable test case). Errors in the horizontal wind speed and direction computed from the wind components are sensitive to the background wind speed but have negligible dependence on the relative orientation of the instrument. Especially during low winds and in the presence of large errors in the horizontal velocity estimates, the reported wind speed is subject to a systematic positive bias (up to 0.4 m s−1 in 1 Hz measurements in strong convection). Vector time-averaged measurements can improve the behavior of the error distributions (reducing the 10 min wind speed error standard deviation to <0.3 m s−1 and the bias to <0.1 m s−1 in strong convection) with a predictable effectiveness related to the number of decorrelated samples in the time window. Hybrid schemes weighting the 10 min scalar- and vector-averaged lidar measurements are shown to be effective at reducing the wind speed biases compared to cup measurements in most of the simulated conditions, with time averages longer than 10 min recommended for best use in some unstable conditions. The approach in decomposing the error mechanisms with the help of the LES flow field could be extended to more complex measurement scenarios and scans.

17 WIND ENERGY↗

Observation of a new pedestal stability regime in MAST Upgrade H-mode plasmas

Abstract The first pedestal stability and structure analysis on the new MAST Upgrade (MAST-U) spherical tokamak H-mode plasmas is presented. Our results indicate that MAST-U pedestals are close to the low toroidal mode number ( n ) peeling branch of the peeling-ballooning instability, in contrast with MAST H-mode pedestals which were deeply in the high- n ballooning branch. This offers the possibility of reaching the ELM-free quiescent H-mode (Burrell et al 2005 Plasma Phys. Control. Fusion 47 B37–B52) or high-performance super H-mode (Snyder et al 2015 Nucl. Fusion 55 083026; Snyder et al 2019 Nucl. Fusion 59 086017) regimes. In addition, the coupling between the peeling and ballooning branches is weak in MAST-U, suggesting that a path to very high pedestal pressure gradient at high density may exist with sufficient heating power. A possible explanation for the differences between MAST and MAST-U pedestal stability is given in terms of plasma shaping parameters, in particular squareness and elongation, as well as the pedestal top temperature and collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Linear Stability Regime Analysis of the Compressible Reacting Mixing Layer

Previous investigations have shown that a compressible reacting mixing layer can develop two peaks in the mean density weighted vorticity profile. Linear stability analyses show that at these peaks two distinct 'outer' instability modes appear in addition to the more common central mode, which exists unaccompanied in incompressible nonreacting flows. The present study parametrically analyzes the effects of compressibility, heat release, stoichiometry, and density ratio on the amplification rate and obliquity of each stability mode. The mean profiles used in the spatial stability calculation are generated by self-similar solutions of the compressible boundary layer equations combined with the assumption of infinitely fast chemistry. It is shown that the influence of stoichiometry and density ratio on the peaks of the density weighted vorticity profile determines which modes will dominate. Of particular interest are the conditions where two modes are equally amplified, causing the mixing layer to develop into a 'colayer' structure.

Day, M. J.↗

Machine-learning identification of the variability of mean velocity and turbulence intensity for wakes generated by onshore wind turbines: Cluster analysis of wind LiDAR measurements

Light detection and ranging (LiDAR) measurements of isolated wakes generated by wind turbines installed at an onshore wind farm are leveraged to characterize the variability of the wake mean velocity and turbulence intensity during typical operations, which encompass a breadth of atmospheric stability regimes and rotor thrust coefficients. The LiDAR measurements are clustered through the k-means algorithm, which enables identifying the most representative realizations of wind turbine wakes while avoiding the imposition of thresholds for the various wind and turbine parameters. Considering the large number of LiDAR samples collected to probe the wake velocity field, the dimensionality of the experimental dataset is reduced by projecting the LiDAR data on an intelligently truncated basis obtained with the proper orthogonal decomposition (POD). The coefficients of only five physics-informed POD modes are then injected in the k-means algorithm for clustering the LiDAR dataset. The analysis of the clustered LiDAR data and the associated supervisory control and data acquisition and meteorological data enables the study of the variability of the wake velocity deficit, wake extent, and wake-added turbulence intensity for different thrust coefficients of the turbine rotor and regimes of atmospheric stability. Furthermore, the cluster analysis of the LiDAR data allows for the identification of systematic off-design operations with a certain yaw misalignment of the turbine rotor with the mean wind direction.

17 WIND ENERGY↗

Flow-shear destabilization of multiscale electron turbulence

The impact of sheared ${\mathbf{E} \hskip -1pt \times \hskip -1pt \mathbf{B}}$ flow on multiscale turbulence is studied with nonlinear gyrokinetic simulations. Simulations are based on DIII-D-like, high-confinement mode (H-mode) pedestal parameters in the regime of low ion temperature gradient drive, where there is a broad spectrum of electron temperature gradient (ETG)-driven turbulence. It is found that ${\mathbf{E} \hskip -1pt \times \hskip -1pt \mathbf{B}}$ shear can have a significant effect on ETG-driven electron transport, with an unexpected transition from a turbulence stabilization regime at moderate to large shearing rates $γ_E$ to a novel turbulence destabilization regime at low levels of $γ_E$. In the turbulence stabilization regime, the electron energy flux decreases monotonically with $γ_E$, even when $γ_E$ is small compared to electron mode growth rates. The stabilizing effect comes dominantly from the electron, not ion, gyrokinetic equation. In the novel destabilization regime, reduction of zonal energy results from the interaction of $γ_E$-modulated nonlinear drive in the zonal ion gyrokinetic equation, increasing the electron transport over a broad range of shearing rates. Finally, neither of these effects has been observed in previous electron-scale simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Implementation and Validation of a Generalized Actuator Disk Parameterization for Wind Turbine Simulations Within the FastEddy Model

ABSTRACT Fast and accurate large‐eddy simulation (LES) of the atmospheric boundary layer plays a crucial role in advancing wind energy research. Long‐duration wind farm studies at turbine‐resolving scales have become increasingly important to understand the intricate interactions between large wind farms and the atmospheric boundary layer. However, the prohibitive computational cost of these turbulence‐ and turbine‐resolving simulations has precluded such modeling to be exercised on a regular basis. To that end, we implement and validate the generalized actuator disk (GAD) model in the computationally efficient, graphics processing unit (GPU)–resident, LES model FastEddy. We perform single‐turbine simulations under three atmospheric stabilities (neutral, unstable, and stable) and compare them against observations from the Scaled Wind Farm Technology (SWiFT) facility and other LES codes from the recent Wakebench turbine wake model benchmark. Our idealized LES results agree well with observed wake velocity deficit and downstream recovery across stability regimes. Turbine response in terms of rotational speed, generated power, torque, and thrust coefficient are well predicted across stability regimes and are consistent with the LES results from the benchmark. The FastEddy simulations are found to be at least two orders of magnitude more efficient than the traditional CPU‐based LES models, opening the door for realistic LES simulations of full wind plants as a viable standard practice.

17 WIND ENERGY↗

Thermodynamic and kinematic drivers of atmospheric boundary layer stability in the central Arctic during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC)

Observations collected during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) provide a detailed description of the impact of thermodynamic and kinematic forcings on atmospheric boundary layer (ABL) stability in the central Arctic. This study reveals that the Arctic ABL is stable and near-neutral with similar frequencies, and strong stability is the most persistent of all stability regimes. MOSAiC radiosonde observations, in conjunction with observations from additional measurement platforms, including a 10 m meteorological tower, ceilometer, microwave radiometer, and radiation station, provide insight into the relationships between atmospheric stability and various atmospheric thermodynamic and kinematic forcings of ABL turbulence and how these relationships differ by season. We found that stronger stability largely occurs in low-wind (i.e., wind speeds are slow), low-radiation (i.e., surface radiative fluxes are minimal) environments; a very shallow mixed ABL forms in low-wind, high-radiation environments; weak stability occurs in high-wind, moderate-radiation environments; and a near-neutral ABL forms in high-wind, high-radiation environments. Surface pressure (a proxy for synoptic staging) partially explains the observed wind speeds for different stability regimes. Cloud frequency and atmospheric moisture contribute to the observed surface radiation budget. Unique to summer, stronger stability may also form when moist air is advected from over the warmer open ocean to over the colder sea ice surface, which decouples the colder near-surface atmosphere from the advected layer, and is identifiable through observations of fog and atmospheric moisture.

54 ENVIRONMENTAL SCIENCES↗

A Real Time Investigation of Morphological Evolution During Solidification of Different Alloy Systems

Solidification phenomenon which occur at the solid/liquid (s/I) interface play a major role in the determination of structure and hence the technologically important properties of a casting. However, metals being opaque, conclusions related to several important phenomenon such as boundary layer thickness, morphological evolution, and eutectic and cell spacing are deduced from quenching experiments and subsequent post solidification metallographic analysis. Consequently, limited information is obtained about the dynamics of the process. This paper will discuss the recent efforts at the Space Science Laboratory, NASA Marshall Space Flight Center, to view and quantify in-situ and in real time the dynamics of the solidification process and to measure interfacial undercooling. First, a high resolution x-ray transmission microscope (XTM) has been developed to monitor fundamental interfacial phenomena during directional solidification of metals and alloys. The XTM operates in the range of 10-100 KeV and through projection is capable of achieving magnification of up to 16OX. Secondly, an innovative collapsible furnace has been designed to quantify interfacial undercooling by measuring the temperature of a moving s/I interface in reference to a fixed s/l interface. This measurement technique is non-intrusive in nature and is based on the Seebeck principle. In this paper real time results obtained to characterize the dynamics of irregular eutectic spacing will be presented. As an example fiber to lamella or plate transition in the Al-Al2Au eutectic system will be discussed. Further, a resolution limit of 25 micron has permitted viewing in real time morphological instability and cellular growth in Al-Au and Al-Ag systems. Simultaneously, a systematic investigation has been carried out to measure interfacial undercooling for Pb-1 wt.% Sn at and near the marginal stability regime. In conjunction with the XTM observations this study attempts to validate existing relationships between undercooling and growth velocity during plane front growth, marginal stability regime, and stable cellular growth.

Sen, S.↗

Dry Deposition Methods Based on Turbulence Kinetic Energy: 1. Evaluation of Various Resistances and Sensitivity Studies Using a Single‐Point Model

Different functions are used to account for turbulence strength in the atmospheric boundary layer for different stability regimes. These functions are one of the sources for differences among different atmospheric models' predictions and associated biases. Also, turbulence strength is underrepresented in some of the resistance formulations. Here, to address these issues with dry deposition, firstly we take advantage of three‐dimensional (3‐D) turbulence information in estimating resistances by proposing and validating a 3‐D turbulence velocity scale that is relevant for different stability regimes of boundary layer. Secondly, we hypothesize and validate that friction velocity measured by 3‐D sonic anemometer can be effectively replaced by the new turbulence velocity scale multiplied by the von Karman constant. Finally, we (a) present a set of resistance formulations for ozone (O 3 ) based on the 3‐D turbulence velocity scale; (b) intercompare estimations of such resistances with those obtained using existing formulations; and, (c) evaluate simulated O 3 fluxes using a single‐point dry deposition model against long‐term observations of O 3 fluxes at the Harvard Forest (MA) site. Results indicate that the new resistance formulations work very well in simulating surface latent heat and O 3 fluxes when compared to respective existing formulations and measurements at a decadal time scale. Findings from this research may help to improve the capability of dry deposition schemes for better estimation of dry deposition fluxes and create opportunities for the development of a community dry deposition model for use in regional/global air quality models.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of idealized large-eddy simulations performed with the Weather Research and Forecasting model using turbulence measurements from a 250 m meteorological mast

Abstract. We investigate the ability of the Weather Research and Forecasting model to perform large-eddy simulation of canonical flows. This is achieved through comparison of the simulation outputs with measurements from sonic anemometers on a 250 m meteorological mast located at Østerild, in northern Denmark. Østerild is on a flat and rough area, and for the predominant wind directions, the atmospheric flow can be considered to be close to homogeneous. The idealized simulated flows aim at representing atmospheric boundary layer turbulence under unstable, neutral, and stable stability conditions at the surface, which are statistically significant conditions observed at Østerild. We found that the resolved fields from the simulations appear to have the characteristics of the three stability regimes. Vertical profiles of observed mean wind speeds and direction are well reproduced by the simulations, with the largest differences under near-neutral conditions, where the effect of the subgrid-scale model is evident on the vertical wind shear close to the surface. Vertical profiles of observed eddy fluxes are also well reproduced by the simulations, with the largest differences for the three velocity component variances under stable stability conditions, although nearly always within the observed variability. With regards to turbulent kinetic energy, we find good agreement between observations and simulations at all vertical levels. Simulated and observed velocity spectra match very well and show very similar behavior with height and with atmospheric stability within the low-frequency interval; at the effective resolution, the simulated spectra show the typical drop-off of finite differences. Our findings demonstrate that these idealized simulations reproduce the characteristics of atmospheric stability regimes often observed at a high turbulent and flat site within a direction sector, where the air flows over nearly homogeneous land.

17 WIND ENERGY↗

Dissecting the Roles of Stability and Dynamic Regimes in Tropical Cloud Feedbacks Simulated from an Upgraded Multiscale Modeling Framework (MMF)

In this study, we analyze the tropical cloud feedback from a 2 K sea surface temperature perturbation experiment performed with a multiscale modeling framework (MMF). We sort the monthly-mean cloud and radiative properties according to circulation and stability regimes. It is found that the dynamic component of cloud feedback for a given dynamic regime dominates the thermodynamic component in terms of the absolute magnitude.

Kuan-Man Xu↗

Dynamic Subgrid Turbulence Modeling for Shallow Cumulus Convection Simulations beyond LES Resolutions

A scale-dependent dynamic Smagorinsky model is implemented in the Met Office/NERC Cloud (MONC) model using two averaging flavors, along Lagrangian pathlines and local moving averages. The dynamic approaches were compared against the conventional Smagorinsky–Lilly scheme in simulating the diurnal cycle of shallow cumulus convection. The simulations spanned from the LES to the near-gray-zone and gray-zone resolutions and revealed the adaptability of the dynamic model across the scales and different stability regimes. The dynamic model can produce a scale- and stability-dependent profile of the subfilter turbulence length scale across the chosen resolution range. At gray-zone resolutions the adaptive length scales can better represent the early precloud boundary layer leading to temperature and moisture profiles closer to the LES compared to the standard Smagorinsky. As a result, the initialization and general representation of the cloud field in the dynamic model is in good agreement with the LES. In contrast, the standard Smagorinsky produces a less well-mixed boundary layer, which fails to ventilate moisture from the boundary layer, resulting in the delayed spinup of the cloud layer. Moreover, strong downgradient diffusion controls the turbulent transport of scalars in the cloud layer. However, the dynamic approaches rely on the resolved field to account for nonlocal transports, leading to overenergetic structures when the boundary layer is fully developed and the Lagrangian model is used. Introducing the local averaging version of the model or adopting a new Lagrangian time scale provides stronger dissipation without significantly affecting model behavior.

54 ENVIRONMENTAL SCIENCES↗

Formation of Martian Gullies by the Flow of Simultaneously Freezing and Boiling Liquid Water

Geomorphic evidence suggests that recent gullies on Mars were formed by fluvial activity. The Martian gully features are significant because their existence implies the presence of liquid water near the surface on Mars in geologically recent times. Irrespective of the ultimate source of the fluid carving the gullies, we seek to understand the behavior of this fluid after it reaches the Martian surface. We find that, contrary to popular belief, the fluvially-carved Martian gullies require formation conditions such as now occur on Mars, outside of the temperature-pressure stability regime of liquid water. Mars Global Surveyor observations of gully length and our modeling of water stability are consistent with gully formation from the action of pure liquid water that is simultaneously boiling and freezing.

Heldmann, Jennifer L.↗

ELM-free H-mode phase and decoupling of peeling–ballooning stability boundary in the MAST Upgrade tokamak

A linear magnetohydrodynamic (MHD) peeling–ballooning stability analysis of the edge-localized mode (ELM)-free phase of a MAST Upgrade (MAST-U) H-mode plasma is presented. In contrast to other similar discharges, #47018 is found to have a significantly higher and wider pedestal during its ELM-free H-mode phase that lasts for approximately 80 ms; this is made possible by the reduced core MHD mode activity on the q = 2 surface. During this period, there is sustained decoupling of peeling and ballooning branches of the stability boundary on J–α space, opening an access channel to the second stability regime with higher peaks in pedestal current density $J_{\mathrm{N,ped}}$ and pressure gradient (α). Such decoupling of the stability boundary has not previously been observed in MAST-U H-modes, and if such a condition can be readily reproduced, it opens a wide range of opportunities for MAST-U to explore low-collisionality peeling-limited pedestal regimes as well as advanced scenarios such as quiescent H-modes that are relevant to future reactors such as STEP and ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

36 Hz integral linewidth laser based on a photonic integrated 4.0 m coil resonator

Laser stabilization sits at the heart of many precision scientific experiments and applications, including quantum information science, metrology, and atomic timekeeping. Many of these systems narrow the laser linewidth and stabilize the carrier by use of Pound–Drever–Hall (PDH) locking to a table-scale, ultrahigh quality factor (Q), vacuum spaced Fabry–Perot reference cavity. Integrating these cavities to bring characteristics of PDH stabilization to the chip scale is critical to reducing their size, cost, and weight, and enabling a wide range of portable and system-on-chip applications. We report a significant advance in integrated laser linewidth narrowing, stabilization, and noise reduction by use of a photonic integrated 4.0 m long coil resonator to stabilize a semiconductor laser. We achieve a 36 Hz 1 / π <#comment/> -integral linewidth, Allan deviation of 1.8 × <#comment/> 10 − <#comment/> 13 at 10 ms measurement time, and a 2.3 kHz/s drift—to the best of our knowledge, the lowest integral linewidth and highest stability demonstrated for an integrated waveguide reference cavity. This performance represents over an order of magnitude improvement in integral linewidth and frequency noise over previous integrated waveguide PDH stabilized reference cavities and bulk-optic and integrated injection locked approaches, and over two orders of magnitude improvement in frequency and phase noise than integrated injection locked approaches. Two different wavelength coil designs are demonstrated, stabilizing lasers at 1550 nm and 1319 nm. The resonator is bus-coupled to a 4.0 m long coil, with a 49 MHz free spectral range, mode volume of 1.0 × <#comment/> 10 10 µ <#comment/> m 3 , and 142 million intrinsic Q , fabricated in a CMOS compatible, ultralow loss silicon nitride waveguide platform. Our measurements and simulations show that the thermorefractive noise floor for this particular cavity is reached for frequencies down to 20 Hz in an ambient environment with simple passive vibration isolation and without vacuum or thermal isolation. The thermorefractive noise limited performance is estimated to yield an 8 Hz 1 / π <#comment/> -integral linewidth and Allan deviation of 5 × <#comment/> 10 − <#comment/> 14 at 10 ms, opening a stability regime that heretofore has been available only in fundamentally non-integrated systems. These results demonstrate the potential to bring the characteristics of laboratory-scale stabilized lasers to the integrated, wafer-scale compatible chip scale, and are of interest for a number of applications in quantum technologies and atomic, molecular, and optical physics, and with further developments below 10 Hz linewidth, can be highly relevant to ultralow noise microwave generation.

Liu, Kaikai (ORCID:0000000213609646)↗

New methods to improve the vertical extrapolation of near-surface offshore wind speeds

Accurate characterization of the offshore wind resource has been hindered by a sparsity of wind speed observations that span offshore wind turbine rotor-swept heights. Although public availability of floating lidar data is increasing, most offshore wind speed observations continue to come from buoy-based and satellite-based near-surface measurements. The aim of this study is to develop and validate novel vertical extrapolation methods that can accurately estimate wind speed time series across rotor-swept heights using these near-surface measurements. We contrast the conventional logarithmic profile against three novel approaches: a logarithmic profile with a long-term stability correction, a single-column model, and a machine-learning model. These models are developed and validated using 1 year of observations from two floating lidars deployed in US Atlantic offshore wind energy areas. We find that the machine-learning model significantly outperforms all other models across all stability regimes, seasons, and times of day. Machine-learning model performance is considerably improved by including the air–sea temperature difference, which provides some accounting for offshore atmospheric stability. Finally, we find no degradation in machine-learning model performance when tested 83 km from its training location, suggesting promising future applications in extrapolating 10 m wind speeds from spatially resolved satellite-based wind atlases.

17 WIND ENERGY↗

Performance, Stability, and Plume Characterization of the HERMeS Thruster with Boron Nitride Silica Composite Discharge Channel

NASA's Hall Effect Rocket with Magnetic Shielding (HERMeS) 12.5kW Technology Demonstration Unit-3 (TDU-3) has been the subject of extensive technology maturation in preparation for flight system development. Detailed performance, stability, and plume characterization tests of the thruster were performed at NASA GRC's Vacuum Facility 5 (VF-5). The TDU-3 thruster implements a magnetic topology that is identical to TDU-1. The TDU-3 boron nitride silica composite discharge channel material is different than the TDU-1 heritage boron nitride discharge channel material. Performance and stability characterization of the TDU-3 thruster was performed at discharge voltages between 300V and 600V and at discharge currents between 5A and 21.8A. The thruster performance and stability were assessed for varying magnetic field strength, cathode flow fractions between 5% and 9%, varying harness inductance, and for reverse magnet polarity. Performance characterization test results indicate that the TDU-3 thruster performance is in family with the TDU-1 levels. TDU-3's thrust efficiency of 65% and specific impulse of 2,800sec at 600V and 12.5kW exceed performance levels of SOA Hall thrusters. Thruster stability regimes were characterized with respect to the thruster discharge current oscillations (discharge current peak-to-peak and root mean square magnitudes), discharge current waveform power spectral density analysis, and maps of the current-voltage-magnetic field. Stability characterization test results indicate a stability profile similar to TDU-1. Finally, comparison of the TDU-1 and TDU-3 plume profiles found that there were negligible differences in the plasma plume characteristics between the TDU with heritage boron nitride versus the boron nitride silica composite discharge channel.

ion engines↗