Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fixed grid”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

Making coal relevant for small scale applications: Modular gasification for syngas/engine CHP applications in challenging environments

Small-scale coal gasification technology, coupled to a reciprocating engine generator, has the potential for making coal a cost competitive resource for meeting the flexible energy needs and resiliency requirements of utilities across the United States. To maintain grid stability and reliability, electrical generation must be regulated to match the load at the proper voltage and frequency. With the expansion of intermittent sources into the grid, such as wind and solar, frequency and voltage regulation become increasingly important and challenging. This work presents the results of a Front End Engineering Design (FEED) effort to detail the engineering and preliminary economics of a small-scale, air blown, fixed-bed gasification process, operating at near-atmospheric pressure, with gas cleanup to provide syngas and pyrolysis liquid fuels for use in reciprocating engine generators for combined heat and power at the University of Alaska Fairbanks. A very detailed assessment of capital and operating costs allows the evaluation of a levelized cost of electricity of 208.06 $/MWh, which can be reduced through the selling of by-products (steam for heating purposes and pyrolysis liquids to be used for power generation in an existing diesel engine). Here, a combined sensitivity analysis, based on the Monte Carlo approach, has been carried out to evaluate the effects of the uncertainties (capital and operating costs and plant annual availability) on the LCOE.

01 COAL, LIGNITE, AND PEAT↗

Assessing CLUBB PDF Closure Assumptions for a Continental Shallow-to-Deep Convective Transition Case Over Multiple Spatial Scales

Assumed-PDF (probability density function) higher-order turbulence closures (APHOCs) are now widely used for parameterizing boundary layer turbulence and shallow convection in Earth system models (ESMs). A better understanding of the resolution-dependent behavior of APHOCs is essential for improving the performance of next-generation ESMs with intended horizontal resolutions finer than 10 km. In this study, we evaluate the PDF family of Analytic double-Gaussian 1 implemented in Cloud Layers Unified By Binormals (CLUBB) over a range of spatial scales (Dx) from 2 to 100 km. A 120-km-wide large eddy simulation (LES) for a continental convection case during 2016 Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign serves as benchmark to evaluate the PDF closure using an off-line approach. We find during the shallow convection period, the CLUBB PDF closure tends to produce positive biases of cloud properties and liquid water flux near cloud base for all scales of analysis. It produces negative biases for these variables near cloud top that are more severe for Dx larger than 25 km. Results show that replacing the CLUBB-parameterized moisture and temperature skewnesses with LES-derived ones can fix most of the biases if clipping of input moments is allowed to prevent the occurrence of unrealizable solutions. Overall, the performance of the PDF closure is better for smaller Dx = 2–5 km than for larger Dx = 50–100 km; for a given grid spacing, it is better when the convective clouds become deeper in the late afternoon. Likely causes for the resolution dependence and implications for improving the PDF closure are discussed.

54 ENVIRONMENTAL SCIENCES↗

Summary of the Third AIAA CFD Drag Prediction Workshop

The workshop focused on the prediction of both absolute and differential drag levels for wing-body and wing-al;one configurations of that are representative of transonic transport aircraft. The baseline DLR-F6 wing-body geometry, previously utilized in DPW-II, is also augmented with a side-body fairing to help reduce the complexity of the flow physics in the wing-body juncture region. In addition, two new wing-alone geometries have been developed for the DPW-II. Numerical calculations are performed using industry-relevant test cases that include lift-specific and fixed-alpha flight conditions, as well as full drag polars. Drag, lift, and pitching moment predictions from previous Reynolds-Averaged Navier-Stokes computational fluid Dynamics Methods are presented, focused on fully-turbulent flows. Solutions are performed on structured, unstructured, and hybrid grid systems. The structured grid sets include point-matched multi-block meshes and over-set grid systems. The unstructured and hybrid grid sets are comprised of tetrahedral, pyramid, and prismatic elements. Effort was made to provide a high-quality and parametrically consistent family of grids for each grid type about each configuration under study. The wing-body families are comprised of a coarse, medium, and fine grid, while the wing-alone families also include an extra-fine mesh. These mesh sequences are utilized to help determine how the provided flow solutions fair with respect to asymptotic grid convergence, and are used to estimate an absolute drag of each configuration.

Vassberg, John C.↗

NEXT Propellant Management System Integration With Multiple Ion Thrusters

As a critical part of the NEXT test validation process, a multiple-string integration test was performed on the NEXT propellant management system and ion thrusters. The objectives of this test were to verify that the PMS is capable of providing stable flow control to multiple thrusters operating over the NEXT system throttling range and to demonstrate to potential users that the NEXT PMS is ready for transition to flight. A test plan was developed for the sub-system integration test for verification of PMS and thruster system performance and functionality requirements. Propellant management system calibrations were checked during the single and multi-thruster testing. The low pressure assembly total flow rates to the thruster(s) were within 1.4 percent of the calibrated support equipment flow rates. The inlet pressures to the main, cathode, and neutralizer ports of Thruster PM1R were measured as the PMS operated in 1-thruster, 2-thruster, and 3-thruster configurations. It was found that the inlet pressures to Thruster PM1R for 2-thruster and 3-thruster operation as well as single thruster operation with the PMS compare very favorably indicating that flow rates to Thruster PM1R were similar in all cases. Characterizations of discharge losses, accelerator grid current, and neutralizer performance were performed as more operating thrusters were added to the PMS. There were no variations in these parameters as thrusters were throttled and single and multiple thruster operations were conducted. The propellant management system power consumption was at a fixed voltage to the DCIU and a fixed thermal throttle temperature of 75 C. The total power consumed by the PMS was 10.0, 17.9, and 25.2 W, respectively, for single, 2-thruster, and 3-thruster operation with the PMS. These sub-system integration tests of the PMS, the DCIU Simulator, and multiple thrusters addressed, in part, the NEXT PMS and propulsion system performance and functionality requirements.

Sovey, James S.↗

Interactive Visualization of High-Dimensional Petascale Ocean Data

We describe an application for interactive visualization of 5 petabytes of time-varying multivariate data from a high-resolution global ocean circulation model. The input data are 10311 hourly (ocean time) time steps of various 2D and 3D fields from a 22-billion point 1/48- degree “lat-lon cap” configuration of the MIT General Circulation Model (MITgcm). We map the global horizontal model domain onto our 128-screen (8x16) tiled display wall to produce a canonical tiling with approximately one MITgcm grid point per display pixel, and using this tiling we encode the entire time series for multiple native and computed scalar quantities at a collection of ocean depths. We reduce disk bandwidth requirements by converting the model’s floating point data to 16-bit fixed point values, and compressing those values with a lossless video encoder, which together allow synchronized playback at 24 time steps per second across all 128 displays. The application allows dynamic assignment of any two encoded tiles to any display, and has multiple interfaces for quickly specifying various orderly arrangements of tiles. All subsequent rendering is done on the fly, with run time control of colormaps, transfer functions, histogram equalization, and labeling. The two data streams on each screen can be rendered independently and combined in various ways, including blending, differencing, horizontal/ vertical wipes, and checkerboarding. The two data streams on any screen can optionally be displayed as a scatterplot in their joint attribute space. All scatterplots and map-view plots from the same x/y location and depth are linked so they all show the current brushable selection. Ocean scientists have used the system, and have found previously unidentified features in the data.

Interactive↗

Methods and systems for evaluating data transportability in distribution grids

This document describes systems and techniques for evaluating and improving distribution-grid data transportability. These systems and techniques allow engineers to quantify the data transportability of a communication system within or connected to a distribution grid, which represents an ability to transport in real-time telemetry from source locations (e.g., sensors in the distribution grid) to control mechanisms. Distribution engineers can use the sensor readings to perform grid analytics, control operating parameters, and operate protection systems. Distribution engineers can also use the transportability of the communication system to evaluate the observability of the distribution grid, which represents an ability to combine actual measurements and various types of computations (e.g., analytics, estimators, forecasters) from a system model. Distribution engineers can then generate a sensor allocation plan that indicates the number and location of sensors to maximize observability for a fixed sensor cost and/or minimize sensor cost for predetermined observability.

Taft, Jeffrey D.↗

On the numerical sensitivity of cellular automata grain structure predictions to large thermal gradients and cooling rates

Cellular automata (CA) models of as-solidified grain structure, originally developed and applied to casting, have become a common means of predicting grain structure resulting from Additive Manufacturing (AM) processes. The majority of these models are based on the decentered octahedron approach, which attempts to correct for the effect of grid anisotropy on the prediction of competitive solidification of dendritic grains. However, AM solidification occurs under cooling rates ($\dot{T}$) and thermal gradients (G) that are orders of magnitude larger than those encountered in casting, and no systematic investigation on the effect of the CA model cell size (Δx) and time step (Δt) on AM microstructure predictions has been performed. Here, in this study, such an investigation is first performed via simulation of individual grains of various crystallographic orientations with a fixed, unidirectional G, showing that CA prediction of the steady-state undercooling matched the expected values based on the interfacial response function at small G and deviated from the expected values at large G. Simulation of competitive growth of multiple grains showed a weakening of the predicted texture as G and Δx became large. Simulation of solidification under AM conditions, where G and $\dot{T}$ vary spatially across the melt pools, showed that not only does grain selection weaken and deviate from expectations at large Δx, but grains with crystallographic $\langle$100$\rangle$ aligned with the grid directions are more adversely affected by the temperature field discontinuities than grains with other crystallographic orientations. Despite the fact that the exact grain competition results depended on Δt, the overall texture development was notably less sensitive to Δt than Δx, provided that a reasonable value of Δt is selected based on the ratio of Δx to the maximum local solidification velocity in the simulation domain. Finally, from the directional solidification and AM simulation results, an analysis of computational cost compared to simulation resolution is performed based on an equation derived to quantify the relatively inaccuracy in grain selection based on the model and temperature field inputs. From this analysis, it is concluded that there is a need for algorithmic improvements to improve CA grain competition accuracy for large G processing conditions as sufficiently small Δx to resolve the necessary competition is intractable for many AM processing conditions.

36 MATERIALS SCIENCE↗

Microwave Power Combiners for Signals of Arbitrary Amplitude

Schemes for combining power from coherent microwave sources of arbitrary (unequal or equal) amplitude have been proposed. Most prior microwave-power-combining schemes are limited to sources of equal amplitude. The basic principle of the schemes now proposed is to use quasi-optical components to manipulate the polarizations and phases of two arbitrary-amplitude input signals in such a way as to combine them into one output signal having a specified, fixed polarization. To combine power from more than two sources, one could use multiple powercombining stages based on this principle, feeding the outputs of lower-power stages as inputs to higher-power stages. Quasi-optical components suitable for implementing these schemes include grids of parallel wires, vane polarizers, and a variety of waveguide structures. For the sake of brevity, the remainder of this article illustrates the basic principle by focusing on one scheme in which a wire grid and two vane polarizers would be used. Wire grids are the key quasi-optical elements in many prior equal-power combiners. In somewhat oversimplified terms, a wire grid reflects an incident beam having an electric field parallel to the wires and passes an incident beam having an electric field perpendicular to the wires. In a typical prior equal-power combining scheme, one provides for two properly phased, equal-amplitude signals having mutually perpendicular linear polarizations to impinge from two mutually perpendicular directions on a wire grid in a plane oriented at an angle of 45 with respect to both beam axes. The wires in the grid are oriented to pass one of the incident beams straight through onto the output path and to reflect the other incident beam onto the output path along with the first-mentioned beam.

Conroy, Bruce↗

An Updated Secondary Lunar Meteoroid Ejecta Model for Engineering Design

Introduction: The surface of the Moon is constantly being bombarded by a flux of meteoroids of various sizes. Impacts due to these meteoroids produce secondary ejecta material at much lower speeds but with a total mass larger than the original impactor. Details about the secondary ejecta are important for planning missions on the lunar surface. In this work, an updated ejecta model is presented called the Meteoroid Model of Secondary Ejecta (MeMoSeE), to re-place the Apollo-era ejecta model, NASA SP-8013 [1], in the SLS-SPEC 159 Design Specification for Natural Environments (DSNE) [2]. The model pro-duces secondary ejecta flux environments for a user-specified location on the lunar surface, and sorts the incoming secondary flux by angular direction and speed. Methods: MeMoSeE is separated into three parts: the inputs, the conversion step, and the integration of fluxes. Inputs to the model include the primary meteoroid fluxes and the primary near-Earth object (NEO) fluxes. Meteoroid fluxes, both asteroidal and cometary, are calculated using the Meteoroid Engineering Model (MEM3) [3] for different locations on the Moon. For each surface location, an ephemeris is generated using the JPL HORIZONS System [4] that feeds into MEM3. The NEO fluxes are approximated by the high-density population of MEM3 (i.e., only the directionality), where the speed distribution of the NEO fluxes is renormalized to match observations [5]. The regolith properties are used as defined in the DSNE [2]. The conversion step utilizes scaling laws given by Housen & Holsapple 2011 [6] to convert the primary impactor flux to the total mass of secondary ejecta. The ejecta distribution, at the point-of-impact (POI), is separated into a zenith angle and azimuthal angle distribution. The zenith angle distribution follows a beta distribution where the peak depends on the impact altitude angle and the impact azimuth [7]. We employ an ejecta azimuth distribution that is based on Rival & Mandeville 1999 [8] which focuses ejecta in the downstream direction for more oblique impacts. Finally, during the integration step, we sum secondary ejecta number fluxes at a particular region-of-interest (ROI) that originated from many POI locations over the entire surface of the Moon. We keep track of both altitude and azimuth angle bins as well as a range of speed bins, following the igloo gridding as done in MEM3 [9]. The ejecta particle size distribution and density is assumed to be the same as the lunar regolith [10, 2]. Results: The primary fluxes are computed for one Metonic cycle (19 years) for various locations over the lunar surface with a fixed orientation. Both the angular and speed distributions of the primary fluxes are dependent on the latitude and longitude. The over-all primary fluxes show a roughly 13% increase from the eastern limb to the western limb. In general, the speed distribution of the secondary fluxes span from a user-defined minimum speed to the escape speed of the Moon (2.38 km/s), roughly following a power-law relation [6]. Different parts of the speed distribution come from different primary impact locations on the Moon. The secondary ejecta is dominated by the slowest speeds, where these particles originate nearby the ROI. For speeds around 71% of the escape speed, the secondary ejecta originates from locations near the antipodal point. On the other hand, for secondary ejecta speeds that exceed roughly 90% the escape speed, the ejecta particles come from all over the lunar surface to the ROI. Comparing the secondary ejecta fluxes from MeMoSeE with NASA SP-8013 [1], there is a reduction by about 2-3 orders of magnitude for secondary ejecta particles greater than 1 µg. These estimates agree with recent findings from Bjorkman & Christiansen 2019 [11]. The secondary ejecta fluxes are also compared with the primary fluxes, where the ejecta fluxes are roughly an order of magnitude greater than the primary fluxes. References: [1] Cour-Palais, B. G., (1969) NASA SP-8013. [2] NASA SLS-SPEC-159 Rev. H (2020). [3] Moorhead, A. V., et al. (2019) JS&R, 1-17. [4] Giorgini, J. D., (2015) IAUGA, 29, 2256293. [5] Moorhead, A. V., (2020) Memo OSMA/MEO/Lunar-001. [6] Housen, K. R., and Holsapple, K. A., (2011) Icarus, 211(1), 856-875. [7] Gault, D. E., and Wedekind, J. A., (1978) L&PSCP, 9, 3843-3875. [8] Rival, M., and Mandeville, J., (1999) Space Debris, 1(1), 45-57. [9] Moorhead, A. V., (2019) MEM3 User Guide. [10] Carrier, W. D., (1973) The Moon, 6(3-4), 250-263. [11] Bjorkman, M. D. & Christiansen, E. L., (2019) ODC

impacts↗

Application of adaptive grids to fluid-flow problems with asymptotic solutions

Coordinate system selection is an important consideration in the asymptotic numerical solution of any fluid-flow or heat transfer problem. This paper uses a new technique that provides a simple way of moving the mesh points in physical space in order to reduce the error in the computed asymptotic solution relative to that obtained using a fixed mesh. Applications to fluid-flow problems are presented, including boundary layer flow and inviscid supersonic flow over cylinders, and wedges with associated detached shocks. The treatment of curved boundaries, stationary and nonstationary boundaries, and systems of PDE's is discussed. Significant error reductions are demonstrated.

Rai, M. M.↗

Methane monitoring and detection apparatus and methods

A low cost, low power, passive optical methane monitoring system for fixed-position installation at oil and gas production well pads and gathering centers is disclosed. The optical methane monitoring system disclosed can be a scannable field of view Near Infrared (NIR) filter photometer to detect and quantify methane concentration in a two dimensional or a three dimensional grid above and around a facility. A randomized fiber optic bundle is disclosed that can be used to direct the total optical power from a collection lens to two or more isolated optical channels. Band pass filters isolate a desired wavelength range for transmission measurements for the two or more channels. Also disclosed is an absorption algorithm which accounts for variable background spectral intensity as well as correcting for water vapor and overall scattering effects to measure methane concentration for a given field of view.

Lievois, John↗

Techno-Economic and Deployment Analysis of Fossil Fuel-Based Power Generation with Integrated Energy Storage

Most existing coal-fired power plants were designed for sustained operation at full load to maximize efficiency, reliability, and revenues. Depending on plant type and design, these plants can adjust output within a fixed range in response to plant or market conditions. The need for flexibility driven by increased penetration of variable and non-dispatchable power generation such as wind and solar is shifting traditional mission profile of thermoelectric power plants in three ways: more frequent shutdowns when market or grid conditions warrant, more aggressive load ramp rates (rate of output change), and lower minimum sustainable load, which provides a wider operating range and helps avoid costly plant shutdowns. The recent studies have shown that flexibility of a coal-fired power plant can be improved by the energy storage. The objective of this project was to analyze a set of energy storage options (technologies) and determine their impact on flexibility and economics of a representative coal-fired power plant. The effect of five Energy Storage Systems (ESSs) integrated with a coal power plant on plant flexibility and economics was investigated in this study. The results obtained in this project and presented in this report showed that ESS integrated with a thermal power plant improves plant flexibility and participation in the Energy and Ancillary Services markets and also improves plant financial performance.

20 FOSSIL-FUELED POWER PLANTS↗

Energy Storage Improves Power Plant Flexibility and Economic Performance

Most existing coal-fired power plants were designed for sustained operation at full load to maximize efficiency, reliability, and revenue, as well as to operate air pollution control devices at design conditions. Depending on plant type and design, these plants can adjust output within a fixed range in response to plant operating or market conditions. The need for flexibility driven by increased penetration of variable and non-dispatchable power generation, such as wind and solar, is shifting the traditional mission profile of thermoelectric power plants in three ways: more frequent shutdowns when market or grid conditions warrant, more aggressive load ramp rates (rate of output change), and a lower minimum sustainable load, which provides a wider operating range and helps avoid costly plant shutdowns. Recent studies have shown that the flexibility of a coal-fired power plant can be improved by energy storage. The objective of this work was to analyze a set of energy storage options and determine their impact on the flexibility and economics of a representative coal-fired power plant. The effect of three energy storage systems integrated with a coal power plant on plant flexibility and economics was investigated. The results obtained in this project show that energy storage systems integrated with a thermal power plant improve plant flexibility and participation in the energy and ancillary services markets, which improves plant financial performance.

20 FOSSIL-FUELED POWER PLANTS↗

Adaptive bandwidth management for entanglement distribution in quantum networks

Flexible grid wavelength division multiplexing is a powerful tool in lightwave communications to maximize spectral efficiency. In the emerging field of quantum networking, the need for effective resource provisioning is particularly acute, given the generally lower power levels, higher sensitivity to loss, and inapplicability of optical detection and retransmission. In this letter, we leverage flex grid technology to demonstrate reconfigurable distribution of quantum entanglement in a four-user tabletop network. By adaptively partitioning bandwidth with a single wavelength-selective switch, we successfully equalize two-party coincidence rates that initially differ by over two orders of magnitude. Our scalable approach introduces loss that is fixed with the number of users, offering a practical path for the establishment and management of quality-of-service guarantees in large quantum networks.

42 ENGINEERING↗

PopGNN: Graph Neural Network-Based Flexible Future Population Forecasting Model

Accurate population forecasts is important to plan critical infrastructure and services, from housing and education to healthcare and transport. However, traditional population prediction studies have only employed traditional machine learning models limited to capture complex spatial interdependencies and patterns. Althogh recently computer vision-based framework was introduced with with promising accuracy, it has critical limitations for real-world planning applications: it function only at fixed spatial resolutions, restricting their use in diverse boundaries such as census tracts, neighborhoods, or administrative zones. Therefore, this study suggests a Graph Neural Network (GNN)-based population prediction framework, called PopGNN. This model recorded remarkable performance compared with state-of-the-art models and traditional baseline models in the grid and administrative boundaries. Furthermore, our framework achieved comparable predictive accuracy to a computer vision-based model in both the South Korea and Tennessee case studies. Consequently, this study is valuable in that a single model can provide accurate population forecasts that address diverse planning demands, ranging from granular grid-level estimates for precise service allocation and facility location planning to aggregate administrative-level forecasts for macro-scale regional policy and resource distribution.

97 MATHEMATICS AND COMPUTING↗

An upwind-biased, point-implicit relaxation algorithm for viscous, compressible perfect-gas flows

An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for three-dimensional, viscous, compressible, perfect-gas flows is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. Viscous terms are discretized using central differences. The relaxation strategy is well suited for computers employing either vector or parallel architectures. It is also well suited to the numerical solution of the governing equations on unstructured grids. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large, block tri-diagonal systems. Convergence rates and grid refinement studies are conducted for Mach 5 flow through an inlet with a 10 deg compression ramp and Mach 14 flow over a 15 deg ramp. Predictions for pressure distributions, surface heating, and aerodynamics coefficients compare well with experiment data for Mach 10 flow over a blunt body.

Gnoffo, Peter A.↗

An Accelerated Clip Algorithm for Unstructured Meshes: A Batch-Driven Approach

The clip technique is a popular method for visualizing complex structures and phenomena within 3D unstructured meshes. Meshes can be clipped by specifying a scalar isovalue to produce an output unstructured mesh with its external surface as the isovalue. Similar to isocontouring, the clipping process relies on scalar data associated with the mesh points, including scalar data generated by implicit functions such as planes, boxes, and spheres, which facilitates the visualization of results interior to the grid. In this paper, we introduce a novel batch-driven parallel algorithm based on a sequential clip algorithm designed for high-quality results in partial volume extraction. Our algorithm comprises five passes, each progressively processing data to generate the resulting clipped unstructured mesh. The novelty lies in the use of fixed-size batches of points and cells, which enable rapid workload trimming and parallel processing, leading to a significantly improved memory footprint and run-time performance compared to the original version. On a 32-core CPU, the proposed batch-driven parallel algorithm demonstrates a run-time speed-up of up to 32.6x and a memory footprint reduction of up to 4.37x compared to the existing sequential algorithm. The software is currently available under an open-source license in the VTK visualization system.

Tsalikis, Spiros↗

Energy and particle balance during plasma detachment in a long-leg divertor configuration

Abstract Comprehensive studies of energy and particle balances in the transition to plasma detachment in an alternative divertor configuration with long outer legs are shown. Numerical simulations are performed with the 2D code suite SOLPS 4.3, using a disconnected double null grid with narrow, tightly baffled long poloidal leg divertors at the outer lower target and outer upper target. A particle count scan is performed using the ‘closed gas box’ model, where the tunable parameter in the simulations is the total number of deuterium particles in the simulation space and all other parameters are held fixed, including a constant input power and trace neon impurity radiation, to assess the physics of the transition to detachment in the system as the particle count increases. Three main aspects of the physics of divertor detachment are addressed: the criteria for the local onset of divertor detachment in each of the divertors, the distribution of heat flux and other plasma parameters between the four divertors as each divertor transitions to detachment, and the role of perpendicular transport in the transition to the detached regime. A synergistic mechanism by which the cross-field transport is reduced by factors associated with the onset of plasma recombination effects is identified. These results are compared to the existing understanding of the physics of the transition to plasma detachment in standard divertors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗