Search NASA⌕ Search

SEARCH · Search NASA

Results for “P CODES”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Vadose Zone Model for C-9 Pond Area for Composite Analysis

The objectives of the vadose modeling for the updated Hanford Site composite analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates for the plateau to river (P2R) model, version 8.3 (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The following 16 radionuclides were selected for this modeling effort; carbon-14 (C-14), chlorine-36 (Cl-36), tritium (H-3), iodine-129 (I-129), neptunium-237 (Np-237), rhenium-187 (Re-187), strontium-90 (Sr-90), technetium-99 (Tc-99), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), uranium-238 (U-238), radium-226 (Ra-226), and thorium-230 (Th-230). The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The parallel version of the Subsurface Transport Over Multiple Phases (STOMP1) simulator, officially named the exascale Subsurface Transport Over Multiple Phases (eSTOMP), is used to simulate flow and transport for the vadose models. The documentation for the STOMP code is comprehensive. The theoretical and numerical approaches applied in the STOMP code are documented in a published theory guide (PNNL-12030, STOMP Subsurface Transport Over Multiple Phases Version 2.0 Theory Guide). The code has undergone a rigorous verification procedure against analytical solutions, laboratory-scale experiments, and field-scale demonstrations. The application guide (PNNL-11216, STOMP Subsurface Transport Over Multiple Phases Application Guide) provides instructive examples in the application of the code to classical groundwater problems. The user’s guide (PNNL-15782, STOMP: Subsurface Transport Over Multiple Phases Version 4.0: User’s Guide) describes the general use, input file formatting, compilation, and execution of the code. The primary output of the vadose zone modeling is radionuclide transfer rates to the groundwater for input into the P2R model. The rates will be summed over the 100 by 100 m P2R grid cells that fall within the vadose zone model source domain.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Runaway electron transport in stochastic toroidal magnetic fields

Here, we study the transport and confinement properties of runaway electrons (RE) in the presence of magnetic fields with perturbations producing different levels of stochasticity. We use Kinetic Orbit Runaway Electron Code (KORC) [Carbajal et al. , Phys. Plasmas 24, 042512 (2017) and del-Castillo-Negrete et al. , Phys. Plasmas 25, 056104 (2018)] for simulating the full-orbit (FO) and guiding-center (GC) dynamics of RE in perturbed magnetic fields that exhibit magnetic islands. We extend previous works on this problem [Wingen et al. , Nucl. Fusion 46, 941 (2006); Izzo et al. , Nucl. Fusion 51, 063032 (2011); Papp et al. , Nucl. Fusion 51, 043004 (2011); V. Izzo and P. Parks, Phys. Plasmas 24, 060705 (2017); and Sommariva et al. , Nucl. Fusion 58, 016043 (2018)] by studying in detail full-orbit effects on the RE dynamics. We quantify FO effects on RE transport by performing one-to-one comparisons between FO and GC simulations. It is found that, for the magnetic field configurations considered, GC simulations predict twice the RE losses of FO simulations for 1 MeV and four times the RE losses of FO simulations for 25 MeV. Similarly, we show how different GC and FO dynamics of RE moving around magnetic islands can be, especially in the scenario where the RE Larmor radius is on the order of the size of the magnetic island. We also study the role of rotation of the magnetic islands on RE confinement, and we find that low-frequency toroidal rotation has no observable effect on RE transport in the cases considered. These results shed some light into the potential of avoidance or mitigation mechanisms based on magnetic perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Dynamic load balancing with enhanced shared-memory parallelism for particle-in-cell codes

Furthering our understanding of many of today’s interesting problems in plasma physics – including plasma based acceleration and magnetic reconnection with pair production due to quantum electrodynamic effects – requires large-scale kinetic simulations using particle-in-cell (PIC) codes. However, these simulations are extremely demanding, requiring that contemporary PIC codes be designed to efficiently use a new fleet of exascale computing architectures. To this end, the key issue of parallel load balance across computational nodes must be addressed. We discuss the implementation of dynamic load balancing by dividing the simulation space into many small, self-contained regions or ‘‘tiles,’’ along with shared-memory (e.g., OpenMP) parallelism both over many tiles and within single tiles. The load balancing algorithm can be used with three different topologies, including two space-filling curves. Here, we tested this implementation in the code Osiris and show low overhead and improved scalability with OpenMP thread number on simulations with both uniform load and severe load imbalance. Compared to other load-balancing techniques, our algorithm gives order-of-magnitude improvement in parallel scalability for simulations with severe load imbalance issues.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Magnetic state switching in FeGa microstructures

This work demonstrates that magnetoelectric composite heterostructures can be designed at the length scale of 10µms that can be switched from a magnetized state to a vortex state, effectively switching the magnetization off, using electric field induced strain. This was accomplished using thin film magnetoelectric heterostructures of Fe 81.4 Ga 18.6 on a single crystal (011) [Pb(Mg 1/3 Nb 2/3 )O 3 ] 0.68 -[PbTiO 3 ] 0.32 (PMN-32PT) ferroelectric substrate. Here, the heterostructures were tripped from a multi-domain magnetized state to a flux closure vortex state using voltage induced strain in a piezoelectric substrate. FeGa heterostructures were deposited on a Si-substrate for superconducting quantum interference device magnetometry characterization of the magnetic properties. The magnetoelectric coupling of a FeGa continuous film on PMN-32PT was characterized using a magneto optical Kerr effect magnetometer with bi-axial strain gauges, and magnetic multi-domain heterostructures were imaged using x-ray magnetic circular dichroism—photoemission electron microscopy during the transition to the vortex state. The domain structures were modelled using MuMax 3 , a micromagnetics code, and compared with observations. The results provide considerable insight into designing magnetoelectric heterostructures that can be switched from an ‘on’ state to an ‘off’ state using electric field induced strain.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Outcomes of OpenMP Hackathon: OpenMP Application Experiences with the Offloading Model (Part I)

This paper reports on experiences gained and practices adopted when using the latest features of OpenMP to port a variety of HPC applications and mini-apps based on different computational motifs (BerkeleyGW, WDMApp/XGC, GAMESS, GESTS, and GridMini) to accelerator-based, leadership-class, high-performance supercomputer systems at the Department of Energy. As recent enhancements to OpenMP become available in implementations, there is a need to share the results of experimentation with them in order to better understand their behavior in practice, to identify pitfalls, and to learn how they can be effectively deployed in scientific codes. Additionally, we identify best practices from these experiences that we can share with the rest of the OpenMP community.

Chapman, Barbara↗
DOE OSTI2020

BISON Capability, Validation and Demonstration for Reactivity-Initiated Accidents

Reactivity Initiated Accidents (RIA) are design basis accidents that could have adverse impacts on the core coolability. In the unlikely event that sufficient reactivity is inserted into the reactor core by the ejected/dropped control rod, prompt energy deposition into the fuel can occur, which when sufficiently high can lead to fuel rod failure or, at large energy deposition levels, expulsion of UO2 fragments or molten UO2 material from the fuel rod. This results in a release of fission product and fuel into the coolant and has the potential to compromise core coolability and threatening the pressure boundary of the primary coolant system. The design basis RIA is one of the industry challenging problems that the CASL aimed to address; the CASL RIA Challenge Problem Charter [1] states, ?The Pressurized Water Reactor (PWR) Rod Ejection Accident (REA) and Boiling Water Reactor (BWR) Control Rod Drop Accident (CRDA) are postulated accidents with consequences that are important to nuclear safety (fuel rod integrity and core coolability). Currently each reload core design must be analyzed to meet regulatory acceptance criteria. The goal of CASLs ModSim capability for RIA is to model the event at a higher fidelity, with validation to existing tests, to better model the transient neutronics and the progression of the fuel and cladding thermal-mechanical behavior. These improved analytical capabilities can be used to better inform reload core design, limits on fuel assembly discharge burnup, restrictions on placement of fuel in the reactor, control rod insertion limits, operating margin, and performance sensitivities.? In support of that charter, BISON, the fuel performance code has been used to demonstrate the simulation of thermal-mechanical behavior of LWR fuels during a reactivity-initiated accident. The combination of mechanical, thermal, and thermal-hydraulic phenomena all present during an RIA makes a multi-physics code such as BISON a valuable tool for modeling these scenarios. This paper will highlight many of the activities associated with BISON relevant to reactivity-initiated accident capability development and validation efforts. These efforts have been performed under both CASL and the Nuclear Energy Advanced Modeling and Simulation (NEAMS) programs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

High speed two-dimensional event detection and imaging using an analog interface and a massively parallel processor

A quantitative pulse count (event detection) algorithm with linearity to high count rates is accomplished by combining a high-speed, high frame rate camera with simple logic code run on a massively parallel processor such as a GPU or a FPGA. The parallel processor elements examine frames from the camera pixel by pixel to find and tag events or count pulses. The tagged events are combined to form a combined quantitative event image.

Waugh, Justin↗
DOE OSTI2023

Using the Stix finite element RF code to investigate operation optimization of the ICRF antenna on Alcator C-Mod

Abstract As the Ion Cyclotron Radio Frequency range (ICRF) heating becomes more favorable in fusion devices, the urgency of predicting and mitigating impurity generation that arises from it becomes more pressing. In the ICRF regime, rectified Radio Frequency (RF) sheaths are known to form at antenna and material edges that influence negative effects like sputtering and a decrease in heating efficiency. Methods to mitigate the formation of these RF sheaths through RF image currents cancellation have been experimentally studied. A power-phasing scan done on Alcator C-Mod in which the amount of power on the two inner straps ( P in ) versus the total 4 straps ( P tot ) was varied showed a minimization of enhanced potentials between P in / P tot ∼ 0.7–0.9 while impurities were minimized for P in / P tot ∼ 0.5–0.8. New capabilities in the realm of representing the RF sheath numerically now allow for these experiments to be simulated. Given the size of the sheath relative to the scale of the device, it can be approximated as a Boundary Condition (BC). A new parallelized cold-plasma wave equation solver called Stix implements a non-linear sheath impedance model BC formulated by Myra et al (2015 Phys. Plasmas 22 062507) through the method of finite elements using the MFEM library [ http://mfem.org ]. It is seen that Stix shows qualitative agreement with the measured C-Mod enhanced potentials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Validation of prediction capability of operating space for plasma initiation in MAST-U

DYON is a plasma initiation modelling code that solves the differential equation system of the full circuit equations (plasma current, active coil currents and eddy currents in full passive structures) and 0D global energy and particle balance equations (Kim 2022 Nucl. Fusion 62 126012). In order to test the capability of the full electromagnetic plasma initiation model to predict individual discharges in experiments and thus the operating space in the device, a dedicated experimental database was built in MAST-U by scanning the prefilled gas pressure p 0 and the induced loop voltage V loop . In the experimental operating space of p 0 and V loop the lower and the upper limits of p 0 are determined by the plasma breakdown failure and the plasma burn-through failure, respectively. The lower limit of V loop is determined by the plasma burn-through failure. By directly reading the control room data used in each discharge (i.e. currents in the solenoid, poloidal field coils, and toroidal field coils, p 0 , and gas puffing rate), the full electromagnetic DYON consistently predicted the failed breakdown, failed burn-through, and successful plasma initiation discharges in the experimental database, demonstrating its capability to predict the operating space for inductive plasma initiation. The Paschen curve calculated with the effective connection length in MAST-U indicates a much higher p 0 required for plasma breakdown than the experimental data, indicating that individual field line evaluation is necessary to calculate the quantitative requirements for Townsend breakdown. The demonstration in this paper shows that the full electromagnetic DYON could be a useful simulation tool to assess the feasibility of inductive plasma initiation and to optimise operating scenarios in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

DSPS: Differentiable stellar population synthesis

ABSTRACT Models of stellar population synthesis (SPS) are the fundamental tool that relates the physical properties of a galaxy to its spectral energy distribution (SED). In this paper, we present DSPS: a python package for SPS. All of the functionality in DSPS is implemented natively in the JAX library for automatic differentiation, and so our predictions for galaxy photometry are fully differentiable, and directly inherit the performance benefits of JAX, including portability onto GPUs. DSPS also implements several novel features, such as i) a flexible empirical model for stellar metallicity that incorporates correlations with stellar age, ii) support for the Diffstar model that provides a physically-motivated connection between the star formation history of a galaxy (SFH) and the mass assembly of its underlying dark matter halo. We detail a set of theoretical techniques for using autodiff to calculate gradients of predictions for galaxy SEDs with respect to SPS parameters that control a range of physical effects, including SFH, stellar metallicity, nebular emission, and dust attenuation. When forward modelling the colours of a synthetic galaxy population, we find that DSPS can provide a factor of 5 speed-up over standard SPS codes on a CPU, and a factor of 300-400 on a modern GPU. When coupled with gradient-based techniques for optimization and inference, DSPS makes it practical to conduct expansive likelihood analyses of simulation-based models of the galaxy–halo connection that fully forward model galaxy spectra and photometry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Multi-Physics Investigation of a Natural Circulation Molten Salt Micro-Reactor that Utilizes an Experimental In-Pile Device to Improve Core Physics and System Thermal-Hydraulic Performance.

INL employee PhD Dissertation - The Molten Salt Reactor (MSR) concept is a rapidly evolving Generation IV design that has recently attracted favorable attention due to the potential for reducing waste generation, realizing passive safety features, and seizing on the opportunity for cost effective economics. This thesis investigates the performance benefit of a new device invented by the doctoral candidate. The device is referred to as a Wrapped Helix around an Inclined Plane or WHIP. The WHIP is protected under a provisional patent filed with the USA Patent Office on 28 September 2021 under application number 63/261,776, BEA docket number BA-1254. The WHIP device can be located in-core or near-core to promoted enhanced thermal-hydraulics and neutronics performance. While the WHIP can be employed in a variety of solid or liquid fuel designs, this thesis investigates the device’s benefit in the application of a natural-circulation, micro-molten salt nuclear battery concept (MsNB). This thesis will specifically investigate the temperature coefficient of the MsNB fuel (FLiNaK) using novel temperature sensitivity techniques unique to Serpent particle transport code, evaluate the WHIP’s thermal-hydraulic and neutronic performance effects using both established (STAR-CCM+) and novel (Python code developed to estimate the circulation effective delayed neutron fraction, ßef f ) analytical and numerical methods, evaluate the neutron noise behavior of the MsNB and how the WHIP may alter the character of the MsNB’s transfer function, and how the WHIP affects the autonomous, load-following performance under transient power conditions using Python code developed by the candidate. The sum body of this work, in part, has been published in three journal articles as the timing of the provisional patent process has allowed. Results show that creative utilization of WHIP engineering design and function reduces reactor system volume, fuel loading, control/stability in the buoyant regime

buoyant flow↗
DOE OSTI2025

Understanding the formation of a low-pressure pedestal in the presence of a strong internal transport barrier in DIII-D high β p plasmas

As a promising scenario for fusion reactors, the high poloidal-beta (β P ) scenario is characterized by a strong large radius internal transport barrier (ITB), which significantly enhances the overall confinement quality and the bootstrap current fraction for fully non-inductive operation. It is frequently observed that in the presence of a strong ITB, the pedestal height is lower and is accompanied by small edge localized modes (ELMs), which further improves the compatibility of a high performance core with an edge solution. A mechanism for the formation of the low pedestal is proposed in this paper. It is found that the strong ITB creates an off-axis bootstrap current to clamp the local safety factor q, and thus the magnetic shear in the outer core/pedestal region is increased. Gyrokinetic simulations with the CGYRO code show that the higher magnetic shear brings the experimental profiles into the range where the growth rate of drift-wave instabilities and thus transport is higher, and therefore a lower pedestal gradient is expected. Here, the combination of low pedestal and high magnetic shear further enhances the turbulent transport across the whole pedestal, consistent with power balance analysis. Such a positive feedback mechanism ultimately results in a lower pressure pedestal as observed in experiments. Under such a low pedestal, linear simulations with BOUT++ predict the growth rates of peeling–ballooning modes to be lower across the whole toroidal mode number spectra, and the nonlinear BOUT++ simulation exhibits lower saturated fluctuation intensity as well, consistent with the experimentally observed lower ELM size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Tracking Active Phase Behavior on Boron Nitride during the Oxidative Dehydrogenation of Propane Using Operando X-ray Raman Spectroscopy

Hexagonal boron nitride (hBN) is a highly selective catalyst for the oxidative dehydrogenation of propane (ODHP) to propylene. Using a variety of ex situ characterization techniques, the activity of the catalyst has been attributed to the formation of an amorphous boron oxyhydroxide surface layer. The ODHP reaction mechanism proceeds via a combination of surface mediated and gas phase propagated radical reactions with the relative importance of both depending on the surface-to-void-volume ratio. Here we demonstrate the unique capability of operando X-ray Raman spectroscopy (XRS) to investigate the oxyfunctionalization of the catalyst under reaction conditions (1 mm outer diameter reactor, 500 to 550 °C, P = 30 kPa C 3 H 8 , 15 kPa O 2 , 56 kPa He). We probe the effect of a water cofeed on the surface of the activated catalyst and find that water removes boron oxyhydroxide from the surface, resulting in a lower reaction rate when the surface reaction dominates and an enhanced reaction rate when the gas phase contribution dominates. Computational description of the surface transformations at an atomic-level combined with high precision XRS spectra simulations with the OCEAN code rationalize the experimental observations. Finally, this work establishes XRS as a powerful technique for the investigation of light element-containing catalysts under working conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Demonstration of low-mode shape control in indirect-drive double shell implosions at the NIF

Double shell inertial confinement fusion is a concept for achieving robust thermonuclear burn that uses dense metal shells to compress deuterium-tritium (DT) fuel to fusion conditions. Double shell implosions are typically indirectly driven and involve a target that consists of a low-Z ablator, a foam layer, and a high-Z pusher surrounding the DT fuel. The goal of the campaign is to achieve a volumetric burn as radiation losses from the DT fuel are trapped by the opaque high-Z shell. The overall performance of double shell implosions relies on the efficient collisional transfer of kinetic energy between layers. The efficiency of this transfer (and therefore the overall performance of a given implosion) is degraded by the presence of low-mode asymmetries. P2 asymmetries are often observed in spatially resolved 2D radiographs of nominal double shell implosions. This work discusses three such experiments: one with an oblate P2 asymmetry, one with a prolate P2 asymmetry, and one with an approximate spherical symmetry. After performing a shape analysis of the oblate and prolate implosions to quantify asymmetries, these experimental results were compared with the results of hydrodynamic simulations for the two experiments. Differences between the experiment and simulation were then used to design an approximately spherical implosion by altering the incident laser cone fraction. Radiographs from the experiment that implemented the modified cone fraction show evidence of an implosion that is approximately spherical until bang time. This design is intended to serve as a point design for future studies that will seek to optimize various aspects of the double shell target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A Synoptic VLBI Technique for Localizing Nonrepeating Fast Radio Bursts with CHIME/FRB

We demonstrate the blind interferometric detection and localization of two fast radio bursts (FRBs) with subarcminute precision on the 400 m baseline between the Canadian Hydrogen Intensity Mapping Experiment (CHIME) and the CHIME Pathfinder. In the same spirit as Very Long Baseline Interferometry (VLBI), the telescopes were synchronized to separate clocks, and the channelized voltage (herein referred to as baseband) data were saved to a disk with correlation performed offline. The simultaneous wide field of view and high sensitivity required for blind FRB searches implies a high data rate—6.5 terabits per second (Tb/s) for CHIME and 0.8 Tb s{sup −1} for the Pathfinder. Since such high data rates cannot be continuously saved, we buffer data from both telescopes locally in memory for ≈40 s, and write to the disk upon receipt of a low-latency trigger from the CHIME Fast Radio Burst Instrument (CHIME/FRB). The ≈200 deg{sup 2} field of view of the two telescopes allows us to use in-field calibrators to synchronize the two telescopes without needing either separate calibrator observations or an atomic timing standard. In addition to our FRB observations, we analyze bright single pulses from the pulsars B0329+54 and B0355+54 to characterize systematic localization errors. Our results demonstrate the successful implementation of key software, triggering, and calibration challenges for CHIME/FRB Outriggers: cylindrical VLBI outrigger telescopes which, along with the CHIME telescope, will localize thousands of single FRB events with sufficient precision to unambiguously associate a host galaxy with each burst.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Global Distribution of EMIC Waves and Its Association to Subauroral Proton Precipitation During the 27 May 2017 Storm: Modeling and Multipoint Observations

Recent simulation studies using the RAM-SCB model showed that proton precipitation contributes significantly to the total energy flux deposited into the subauroral ionosphere thereby affecting the magnetosphere-ionosphere coupling. Here, in this study, we use the BATS-R-US + RAM-SCB model to understand the evolution of ElectroMagnetic Ion Cyclotron (EMIC) waves in the inner magnetosphere, their correspondence to the proton precipitation into the subauroral ionosphere, and to assess the performance of the model in reproducing the EMIC wave-particle interactions. During the 27 May 2017 storm, Arase and RBSP-A satellites observed typical signatures of EMIC waves in the inner magnetosphere. Within this interval, Defense Meteorological Satellite Program (DMSP) and National Oceanic and Atmospheric Administration (NOAA)/MetOp satellites observed significant proton precipitation in the dusk-midnight sector. Simulation results show that H- and He-band EMIC waves are excited within regions of strong temperature anisotropy near the plasmapause. The simulated growth rates of EMIC waves show a similar trend to that of the EMIC wave power observed by the Arase and RBSP-A satellites, suggesting that the model can reproduce the EMIC wave activity qualitatively. The simulated H-band waves in the dusk sector are stronger than He-band waves possibly due to the presence of excess protons in the boundary conditions obtained from the BATS-R-US code. The precipitating proton fluxes reproduced by the simulation with EMIC waves are found to agree reasonably well with the DMSP and NOAA/MetOp satellite observations. It is suggested that EMIC wave scattering of ring current ions can account for proton precipitation observed by the DMSP and MetOp satellites during the 27 May 2017 storm.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Experimental observation of the geodesic acoustic frequency limit for the NBI-driven Alfvén eigenmodes in TJ-II

We study Alfvén Eigenmodes (AEs) in the TJ-II heliac in hydrogen plasmas heated by hydrogen cofield NBI. Taking advantage of the unique TJ-II flexibility in varying plasma current, we have observed strong variation of the AE frequency from f AE ~ 30 to ~ 220 kHz for selected modes. An advanced heavy ion beam probe diagnostic determines the spatial location and internal amplitudes of the modes. The modes satisfy a local AE dispersion relation including the geodesic acoustic frequency that represents the lowest frequency of the mode. Linear MHD-modeling with STELLGAP and FAR3D codes shows that the calculated temporal evolution of the mode frequency reproduces the observed maxima and minima at the same time intervals with a similar frequency range and the radial profile peaks near the outer edge of the observed one.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗