Search NASA⌕ Search

SEARCH · Search NASA

Results for “core optimization”

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 109 records · Page 6

Fermilab Automation of Coil Winding System

The automated coil winder project creates a repeatable, scalable method for producing superconducting magnet coils with increased accuracy and efficiency. The design uses an ortho cyclic winding method to achieve high fill factors by combining computer numerical control (CNC), microcontroller-based controls, and sophisticated optimization techniques. The method optimizes coil turns, core geometry, and current consumption by utilizing analytical models based on magnetic moment, Ohm’s law, and power constraints. For particle accelerator applications, prototype testing using a 3D-printed coil winder has shown enhanced winding repeatability, thermal management, and manufacturability.

Dzida, Rafal [Northern Illinois U.]↗

ELM control optimization for various ITER scenarios based on linear and quasi-linear figures of merit

For the purpose of better understanding type-I edge localized mode (ELM) control in ITER with resonant magnetic perturbation (RMP) fields, the plasma response to RMP is computed by a resistive full magneto-hydrodynamic model in toroidal geometry. Five scenarios designed for ITER are considered, ranging from the pre-nuclear to nuclear phases. The plasma response to RMP is quantified by the plasma surface displacement near the X-point of the divertor plasma and at the outboard mid-plane. The optimal coil configurations between two high- Q deuterium-tritium (DT) scenarios (at the same plasma current of 15 MA and the same magnetic field of 5.3 T but different fusion gains, Q = 5 and 10) are predicted to be similar. For the other ITER scenarios with similar edge safety factor q 95 ~ 3 to that of the baseline scenario, the optimal coil phasing is also similar. The optimization results are different for a half-current full-field (7.5 MA/5.3 T) scenario, largely due to the difference in q 95 . The RMP coil currents are also optimized to tailor the core vs edge toroidal torques exerted by the 3D RMP fields on the plasma column. Torque optimization, with various objective functions proposed in the study, is useful for minimizing the side effects of RMP on the plasma core flow in ITER, while still maintaining the ELM control capability. Full utilization of three rows of ELM control coils in ITER is found to be essential to ensure both flexibility and robustness of ELM control, in terms of both linear and quasilinear plasma responses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bayesian optimization to design a novel x-ray shaping device

In radiation therapy, x-ray dose must be precisely sculpted to the tumor, while simultaneously avoiding surrounding organs at risk. This requires modulation of x-ray intensity in space and/or time. Typically, this is achieved using a multi leaf collimator (MLC) - a complex mechatronic device comprising over one hundred individually powered tungsten ‘leaves’ that move in or out of the radiation field as required. Here, an all-electronic x-ray collimation concept with no moving parts is presented, termed “SPHINX”: Scanning Pencil-beam High-speed Intensity-modulated X-ray source. SPHINX utilizes a spatially distributed bremsstrahlung target and collimator array in conjunction with magnetic scanning of a high energy electron beam to generate a plurality of small x-ray “beamlets.” A simulation framework was developed in Topas Monte Carlo incorporating a phase space electron source, transport through user defined magnetic fields, bremsstrahlung x-ray production, transport through a SPHINX collimator, and dose in water. This framework was completely parametric, meaning a simulation could be built and run for any supplied geometric parameters. This functionality was coupled with Bayesian optimization to find the best parameter set based on an objective function which included terms to maximize dose rate for a user defined beamlet width while constraining inter-channel cross talk and electron contamination. Designs for beamlet widths of 5, 7, and 10 mm 2 were generated. Each optimization was run for 300 iterations and took approximately 40 h on a 24-core computer. For the optimized 7-mm model, a simulation of all beamlets in water was carried out including a linear scanning magnet calibration simulation. Finally, a back-of-envelope dose rate formalism was developed and used to estimate dose rate under various conditions. The optimized 5–, 7–, and 10-mm models had beamlet widths of 5.1 , 7.2 , and 10.1 mm 2 and dose rates of 3574, 6351, and 10 015 Gy/C, respectively. The reduction in dose rate for smaller beamlet widths is a result of both increased collimation and source occlusion. For the simulation of all beamlets in water, the scanning magnet calibration reduced the offset between the collimator channels and beam centroids from 2.9 ±1.9 mm to 0.01 ±0.03 mm. A slight reduction in dose rate of approximately 2% per degree of scanning angle was observed. Based on a back-of-envelope dose rate formalism, SPHINX in conjunction with next-generation linear accelerators has the potential to achieve substantially higher dose rates than conventional MLC-based delivery, with delivery of an intensity modulated 100 x 100 mm 2 field achievable in 0.9 to 10.6 s depending on the beamlet widths used. Bayesian optimization was coupled with Monte Carlo modeling to generate SPHINX geometries for various beamlet widths. A complete Monte Carlo simulation for one of these designs was developed, including electron beam transport of all beamlets through scanning magnets, x-ray production and collimation, and dose in water. These results demonstrate that SPHINX is a promising candidate for sculpting radiation dose with no moving parts, and has the potential to vastly improve both the speed and robustness of radiotherapy delivery. A multi-beam SPHINX system may be a candidate for delivering magavoltage FLASH RT in humans.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Optimization on fixed low latency implementation of the GBT core in FPGA

We present that in the upgrade of ATLAS experiment, the front-end electronics components are subjected to a large radiation background. Meanwhile high speed optical links are required for the data transmission between the on-detector and off-detector electronics. The GBT architecture and the Versatile Link (VL) project are designed by CERN to support the 4.8 Gbps line rate bidirectional high-speed data transmission which is called GBT link. In the ATLAS upgrade, besides the link with on-detector, the GBT link is also used between different off-detector systems. The GBTX ASIC is designed for the on-detector front-end, correspondingly for the off-detector electronics, the GBT architecture is implemented in Field Programmable Gate Arrays (FPGA). CERN launches the GBT-FPGA project to provide examples in different types of FPGA. In the ATLAS upgrade framework, the Front-End LInk eXchange (FELIX) system is used to interface the front end electronics of several ATLAS subsystems. The GBT link is used between them, to transfer the detector data and the timing, trigger, control and monitoring information. The trigger signal distributed in the down-link from FELIX to the front-end requires a fixed and low latency. In this paper, several optimizations on the GBT-FPGA IP core are introduced, to achieve a lower fixed latency. For FELIX, a common firmware will be used to interface different front-ends with support of both GBT modes: the forward error correction mode and the wide mode. The modified GBT-FPGA core has the ability to switch between the GBT modes without FPGA reprogramming. Finally, the system clock distribution of the multi-channel FELIX firmware is also discussed in this paper.

47 OTHER INSTRUMENTATION↗

Diii-D Research Advancing The Physics Basis For Optimizing The Tokamak Approach To Fusion Energy

DIII-D physics research addresses critical challenges for operation of ITER and the next generation of fusion energy devices. This is done through a focus on innovations to provide solutions for high performance long pulse operation, coupled with fundamental plasma physics understanding and model validation, to drive scenario development by integrating high performance core and boundary plasmas. Substantial increases in off-axis current drive efficiency from an innovative top launch system for EC power, and in pressure broadening for Alfven eigenmode control from a co-/counter-Ip steerable off-axis neutral beam, all improve the prospects for optimization of future long pulse/steady state high performance tokamak operation. Fundamental studies into the modes that drive the evolution of the pedestal pressure profile and electron vs. ion heat flux validate predictive models of pedestal recovery after ELMs. Understanding the physics mechanisms of ELM control and density pumpout by 3D magnetic perturbation fields leads to confident predictions for ITER and future devices. Validated modeling of high-Z shattered pellet injection for disruption mitigation, runaway electron dissipation, and techniques for disruption prediction and avoidance including machine learning, give confidence in handling disruptivity for future devices. With this physics understanding and suite of actuators, a high poloidal beta optimized-core scenario with an internal transport barrier that projects nearly to Q=10 in ITER at 9 MA was coupled to a detached divertor, and a near Super H-mode optimized-pedestal scenario with co-Ip beam injection, was coupled to a radiative divertor. The Hybrid core scenario was achieved directly, without the need for internal flux pumping, using off-axis current drive actuators, and a controller to avoid mode locking and disruptions in the ITER baseline scenario, based on plasma response to probing 3D fields, was demonstrated. Finally innovative tokamak operation using a negative triangularity shape showed many attractive features for future pilot plant operation.

Fenstermacher, M.↗

Graded Infill Structure of Wind Turbine Blade Accounting for Internal Stress in Big Area Additive Manufacturing

In additively manufactured (AM) components, infill structure significantly affects the mechanical performance of the final printed part. However, mechanical stress induced by operation loads has been so far neglected for infill patterning. Most slicers currently available in the market provide infill patterns that are uniform in shape and size regardless of the operational loading. We develop a design approach for infill patterns that accounts for the induced stress. This approach differs from topology optimization as it focuses on the porous infill, which allows the external shape of the printed part to remain intact. The proposed approach uses a computational stress analysis to control the distribution of the local density of the infill pattern. We have applied the approach to a wind turbine blade core with infill densities optimized based on the structural loads. The blade core is fabricated in our big area additive manufacturing (BAAM) system. To ensure less warpage and better inter-layer bonding, fast layer deposition is critical in BAAM system. We have optimized the tool path sequences to minimize the deposition time via the solution to the Chinese Postman Problem (CPP). For the application of wind turbine infill, the deposition from the CPP method is twice faster than the deposition from conventional slicing.

Kim, Pum↗

DIII-D research advancing the physics basis for optimizing the tokamak approach to fusion energy

DIII-D physics research addresses critical challenges for the operation of ITER and the next generation of fusion energy devices. This is done through a focus on innovations to provide solutions for high performance long pulse operation, coupled with fundamental plasma physics understanding and model validation, to drive scenario development by integrating high performance core and boundary plasmas. Substantial increases in off-axis current drive efficiency from an innovative top launch system for EC power, and in pressure broadening for Alfven eigenmode control from a co-/counter-I p steerable off-axis neutral beam, all improve the prospects for optimization of future long pulse/steady state high performance tokamak operation. Fundamental studies into the modes that drive the evolution of the pedestal pressure profile and electron vs ion heat flux validate predictive models of pedestal recovery after ELMs. Understanding the physics mechanisms of ELM control and density pumpout by 3D magnetic perturbation fields leads to confident predictions for ITER and future devices. Validated modeling of high-Z shattered pellet injection for disruption mitigation, runaway electron dissipation, and techniques for disruption prediction and avoidance including machine learning, give confidence in handling disruptivity for future devices. For the non-nuclear phase of ITER, two actuators are identified to lower the L–H threshold power in hydrogen plasmas. With this physics understanding and suite of capabilities, a high poloidal beta optimized-core scenario with an internal transport barrier that projects nearly to Q = 10 in ITER at ~8 MA was coupled to a detached divertor, and a near super H-mode optimized-pedestal scenario with co-I p beam injection was coupled to a radiative divertor. The hybrid core scenario was achieved directly, without the need for anomalous current diffusion, using off-axis current drive actuators. Also, a controller to assess proximity to stability limits and regulate β N in the ITER baseline scenario, based on plasma response to probing 3D fields, was demonstrated. Finally, innovative tokamak operation using a negative triangularity shape showed many attractive features for future pilot plant operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of Genetic Algorithm Based Multi-Objective Plant Reload Optimization Platform

The U.S. nuclear industry is facing a challenge in maintaining required levels of safety while ensuring economic competitiveness to stay in business. Safety remains a key parameter for all aspects of light-water reactor nuclear power plant operations. Safety can become more economical by using a risk-informed ecosystem, such as the one being developed in the Risk-Informed Systems Analysis Pathway under the U.S. Department of Energy Light Water Reactor Sustainability Program. The Light Water Reactor Sustainability Program promotes a wide range of research and development activities to maximize both the safety and economic efficiency of nuclear power plants through improved scientific understanding, especially given that many plants are now considering second license renewals. The Risk-Informed Systems Analysis Pathway has two main goals: Deploy methodologies and technologies that better represent safety margins and cost and safety factors; Develop advanced applications that enable cost-effective plant operations. The Plant Reload Optimization Platform development project aims to build a reactor core design tool that includes reactor safety and fuel performance analyses and uses artificial intelligence to support the optimization of core design solutions. This report summarizes genetic-algorithm-based multi-objective fuel reload optimization activities, specifically: Developing the non-dominated sorting genetic algorithm II optimizer in the Risk Analysis and Virtual ENviroment (RAVEN); Demonstrating and validating the developed non-dominated sorting genetic algorithm II optimizer using benchmark optimization problems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Advanced Manufacturing for Nuclear Core Design

Advanced manufacturing technologies are fundamentally altering the way components are designed and manufactured. The potential application space for these technologies within the nuclear industry is broad, with the most benefit realized from the ability to approach the nuclear core design problem from an advanced manufacturing perspective. One of the objectives of the US Department of Energy Office of Nuclear Energy (DOE-NE) Transformational Challenge Reactor (TCR) design and analysis activities is to demonstrate this approach. A mitigation of conventional constraints realized through advanced manufacturing methods opens up the design space, enabling a departure from typical axially uniform pin- and plate-type geometries common to operating nuclear reactors worldwide. The ability to use more heterogenous spatial and material flexibilities increases the demand on the predictability of higher fidelity modeling and simulation tools to be used in the design, analysis, and optimization processes. Leading candidate advanced manufacturing technologies for nuclear applications are additive methods that provide for high dimensional control and incorporation of small design features such as powder bed fusion, laser directed energy deposition, and binder jetting. Potential applications to core design identified during DOE-NE TCR design and analysis activities fall into two categories: those that are uniquely enabled by additive manufacturing technologies and those that are enabled by the characteristics (e.g., repeatability and instrumentation) inherent to these technologies. Identified applications include reducing in-core parasitic materials, improving heat conduction across fuel-cladding gaps, optimizing thermal hydraulic design with flexible coolant flow geometries, reducing uncertainties in manufactured components, optimizing core conditions with design heterogeneities, designing feedback mechanisms using material characteristics, and embedding burnable poisons. These potential applications join a host of nuclear industry advanced manufacturing activities to recapture manufacturing capabilities, reduce manufacturing costs and complexities, reduce deployment time, and enable performance and operational enhancements.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Optimization of the Post-Operational Phase on Two Belgian Multi-Unit Nuclear Power Plants: the Case of the Non-Fissile Irradiated Core Items - 20156

The current legal framework in Belgium foresees the progressive phase out of nuclear power between October 2022 (Doel 3) and December 2025 (Doel 2). Upon its definitive shutdown, each unit of the Tihange and Doel sites will enter a Post-Operational Phase (POP) and be prepared for its Decontamination and Decommissioning (D and D). Prior to obtaining the D and D license, the Operator Electrabel is legally required to remove any non-fissile irradiated core items stored in the deactivation pools. The non-fissile irradiated core items consist essentially of control rods, poison rods and source thimbles as well as thimble plugs and foreign materials irradiated during operation: - Their significant content in highly radiant radionuclides (up to 6 TBq of Co-60 per kg of irradiated material) renders all existing operational waste management processes inadequate due to insufficient biological shielding; - Their high concentrations in long-lived radionuclides call for their disposal in a geological repository for which no final design nor waste acceptance criteria are expected prior to 2050. Uncertainties in the Belgian energy supply and security, however, require the Operator to be prepared for a partial nuclear phase out, where one or more units would benefit from lifetime extension while the remaining units would undergo decommissioning. The present paper aims at presenting how Electrabel, in partnership with Tractebel, addressed this challenge by maximizing the use of synergies within the respective sites as well as between both sites themselves, all the while accounting for site specificities. The most recent results and state of progress of the project will be detailed and the first lessons learned will be shared. The project has been split in multiple tasks and phased as follows: - An inventory phase aimed at mapping the contents, origin, composition and history of the non-fissile irradiated core items; - A pre-characterization phase based on neutron activation models; - A waste sorting phase aimed at separating waste forms for which an evacuation route exists from those for which such route does not exist; - A feasibility phase aimed at exploring all possible scenarios for the management of non-fissile irradiated core items and identifying the optimal feasible solution for each site; - A preparation phase (currently ongoing), developing further the optimal solution and ensuring that back-up solutions are available for any foreseeable change of context (licensing issue, modification in the nuclear phase-out program, etc.) and initiating early contacts with potential subcontractors for segmentation works and cask manufacturers, as well as the Belgian regulatory body and waste management agency. This phase also foresees the investigation of destructive and non-destructive radiological measurements to support the detailed characterization of the waste forms; - A realization phase (future work). (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

AI-based design of a nuclear reactor core

The authors developed an artificial intelligence (AI)-based algorithm for the design and optimization of a nuclear reactor core based on a flexible geometry and demonstrated a 3× improvement in the selected performance metric: temperature peaking factor. The rapid development of advanced, and specifically, additive manufacturing (3-D printing) and its introduction into advanced nuclear core design through the Transformational Challenge Reactor program have presented the opportunity to explore the arbitrary geometry design of nuclear-heated structures. The primary challenge is that the arbitrary geometry design space is vast and requires the computational evaluation of many candidate designs, and the multiphysics simulation of nuclear systems is very time-intensive. Therefore, the authors developed a machine learning-based multiphysics emulator and evaluated thousands of candidate geometries on Summit, Oak Ridge National Laboratory’s leadership class supercomputer. The results presented in this work demonstrate temperature distribution smoothing in a nuclear reactor core through the manipulation of the geometry, which is traditionally achieved in light water reactors through variable assembly loading in the axial direction and fuel shuffling during refueling in the radial direction. The conclusions discuss the future implications for nuclear systems design with arbitrary geometry and the potential for AI-based autonomous design algorithms.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Comparing the Performance of Julia on CPUs versus GPUs and Julia-MPI versus Fortran-MPI: a case study with MPAS-Ocean (Version 7.1)

Abstract. Some programming languages are easy to develop at the cost of slow execution, while others are fast at runtime but much more difficult to write. Julia is a programming language that aims to be the best of both worlds – a development and production language at the same time. To test Julia's utility in scientific high-performance computing (HPC), we built an unstructured-mesh shallow water model in Julia and compared it against an established Fortran-MPI ocean model, the Model for Prediction Across Scales–Ocean (MPAS-Ocean), as well as a Python shallow water code. Three versions of the Julia shallow water code were created: for single-core CPU, graphics processing unit (GPU), and Message Passing Interface (MPI) CPU clusters. Comparing identical simulations revealed that our first version of the Julia model was 13 times faster than Python using NumPy, where both used an unthreaded single-core CPU. Further Julia optimizations, including static typing and removing implicit memory allocations, provided an additional 10–20× speed-up of the single-core CPU Julia model. The GPU-accelerated Julia code was almost identical in terms of performance to the MPI parallelized code on 64 processes, an unexpected result for such different architectures. Parallelized Julia-MPI performance was identical to Fortran-MPI MPAS-Ocean for low processor counts and ranges from 2× faster to 2× slower for higher processor counts. Our experience is that Julia development is fast and convenient for prototyping but that Julia requires further investment and expertise to be competitive with compiled codes. We provide advice on Julia code optimization for HPC systems.

54 ENVIRONMENTAL SCIENCES↗

Phase-based design of CO 2 capture, transport, and storage infrastructure via SimCCS 3.0

The design of optimal infrastructure is essential for the deployment of commercial and large-scale carbon capture and storage (CCS) technology. During the design process, it is important to consider CO 2 capture and storage locations and CO 2 transportation pipelines to minimize the total project cost. SimCCS, first introduced in 2009, is an integrated open-source tool to optimize CCS infrastructure. The core CCS infrastructure design problem in SimCCS is structured as a mixed-integer linear programming problem by selecting the optimal pipeline routes, searching CO 2 source capture and storage locations, and determining the corresponding CO 2 amounts to meet desired capture targets. Multiple important and practical features have been developed to the latest version of SimCCS, SimCCS 3.0 . One of these features is phase-based modeling which enables users to dynamically design the CCS infrastructure. We demonstrate the phased-based modeling capability using two CCS infrastructure optimization case studies. The results from these case studies reveal that the phase-based modeling capability in SimCCS is particularly useful to optimize the dynamic deployment of CCS projects.

54 ENVIRONMENTAL SCIENCES↗

Investigation of the neoclassical ambipolar electric field in ion-root plasmas on W7-X

The role of the radial electric field in high-performance ion-root plasmas on Wendelstein 7-X (W7-X) is examined and compared with neoclassical predictions herein. The W7-X stellarator is the world’s first large-scale optimized stellarator. One of the important targets chosen for optimization during the W7-X design process was the reduction of core neoclassical heat transport. This optimization was targeted for reactor-relevant high-density plasmas with $T_e ≈ T_i$ in which the neoclassical ambipolar radial electric field is expected to be negative throughout the plasmas core. Measurements of the core radial electric field $(E_r)$ have confirmed that ion-root conditions (negative $(E_r)$ in the plasma core) have been achieved in W7-X with high-density plasmas and central ECRH. These measured $(E_r)$ profiles agree well with the neoclassical ambipolar $(E_r)$ predicted by the code SFINCS . This good agreement provides confidence in the validity of neoclassical calculations in high-density ion-root conditions, and enables initial studies on the role of neoclassical transport in the optimized high-density regime of W7-X. Profile measurements of electron temperature $(T_e)$, ion temperature $(T_i)$ and electron density $(n_e)$ along with approximations for the average value of $Z$ eff have been used as inputs to the sfincs code to calculate the ambipolar $E_r$ profile along with neoclassical ion and electron energy flux profiles $(Q$ NCi , $Q$ NCe $)$. Finally the total experimental energy input to the plasma from ECRH heating has been compared to the neoclassical energy fluxes to provide a first estimate for the fraction of transport that can be attributed to neoclassical processes in reactorrelevant high-density ion-root conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CMS readiness for multi-core workload scheduling

In the present run of the LHC, CMS data reconstruction and simulation algorithms benefit greatly from being executed as multiple threads running on several processor cores. The complexity of the Run 2 events requires parallelization of the code to reduce the memory-per- core footprint constraining serial execution programs, thus optimizing the exploitation of present multi-core processor architectures. The allocation of computing resources for multi-core tasks, however, becomes a complex problem in itself. The CMS workload submission infrastructure employs multi-slot partitionable pilots, built on HTCondor and GlideinWMS native features, to enable scheduling of single and multi-core jobs simultaneously. This provides a solution for the scheduling problem in a uniform way across grid sites running a diversity of gateways to compute resources and batch system technologies. This paper presents this strategy and the tools on which it has been implemented. The experience of managing multi-core resources at the Tier-0 and Tier-1 sites during 2015, along with the deployment phase to Tier-2 sites during early 2016 is reported. The process of performance monitoring and optimization to achieve efficient and flexible use of the resources is also described.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Merging of the superbanana plateau and $\sqrt{v}$ transport regimes in nearly quasisymmetric stellarators

We report alpha particle confinement is one of the most demanding issues for stellarators. It now seems clear that it is possible to design optimized stellarators that confine the background plasma at near tokamak radial transport levels. Moreover, adequate collisionless alpha particle confinement is possible in the core of a highly optimized stellarator. Here, the collisional confinement of barely trapped alphas in an optimized stellarator is considered by accounting for the resonance due to the reversal in direction of the drift within a flux surface and investigating the sensitive role of magnetic shear in keeping this resonance close to the passing boundary in some nearly quasisymmetric stellarator configurations. The treatment relies on a narrow collisional boundary layer formulation that combines the responses of both these resonant pitch angle alphas and the remaining barely trapped alphas. A novel merged regime treatment leads to explicit expressions for the energy diffusivity for both superbanana plateau (or resonant plateau) and $\sqrt{v}$ transport in the large aspect ratio limit for a slowing down tail alpha distribution function, where v is the effective pitch angle scattering collision frequency of the trapped alphas off the background ions. Depending on the details of the optimization scheme and the sign of the magnetic shear, modest magnetic shear can be used to reduce superbanana (or resonant) plateau transport to below the $\sqrt{v}$ transport level. In addition, a quasilinear equation retaining spatial diffusion is derived for a general alpha distribution function that allows the radial alpha transport to modify the distribution so it is no longer isotropic in velocity space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗