Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative limiters”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Numerical error analysis of SOLPS-ITER simulations of EAST

Abstract Plasma edge simulations with codes like SOLPS-ITER are widely employed to interpret fusion experiments. However, numerical errors appearing in such simulations are rarely investigated, despite their potential large impact on simulation results. These errors consist of the statistical error and the bias, both resulting from the finite number of employed EIRENE Monte Carlo particles and incomplete convergence, and the discretization error due to the finite resolution of the computational grids. In this contribution, the resulting numerical errors on simulations of pure deuterium and neon seeded H-mode EAST discharges are examined. The statistical error can be kept small compared to other numerical error contributions by averaging the plasma profiles. This allows investigating the bias and discretization errors using Richardson extrapolation. It is shown that grid refinement and the number of employed Monte Carlo particles have the largest influence on the result, in agreement with similar studies of an ITER deuterium case. For the first time, numerical error bars on the entire simulated target profiles are determined showing that the largest numerical error is 17.9%, mainly due to the plasma grid discretization. On top, also numerical errors on simulated neutral pressures are investigated in detail, for which the statistical error is dominant. The analysis demonstrates which setup is needed to keep numerical errors limited: the SOLPS-ITER averaging procedure should be employed including enough EIRENE particles, and the involved grids should be sufficiently refined to reduce discretization errors.

Boeyaert, Dieter (ORCID:0000000309208660)↗
DOE OSTI2025

Feasibility of fusion plasma burn control via real-time, sub-divertor neutral gas isotopic and compositional analysis

The ability to provide fusion burn control without requiring physical access through the first wall and fuel breeding blankets, would be vital for any future, magnetically confined fusion power reactor. A multi-sensor, fusion fuel cycle exhaust, neutral gas analysis system on JET, capable of delivering real time data, and accessing only the sub-divertor region, provides an excellent example of such capability. Optimized for and operated during the deuterium–tritium experimental campaigns 2 and 3 (DTE2, DTE3), it is proving valuable for planning to explore fusion reactor burn control in ITER with a comparable diagnostic system called the Diagnostic Residual Gas Analyzer (DRGA). This paper aims to show feasibility of developing model-based controllers for ITER and next generation, reactor-relevant devices, by building both on the empirical experience in JET-DTE2, and on the already emerging experience on developing such models specifically for ITER. The paper begins with a specific use-case from JET-DTE2, pertaining to the observed sensitivity of the fusion neutron yield on the concentration of isotopic helium-3 ( 3 He), with data from one of the high-performance DT shots exhibited with emphasis on the 3 He measurement via the sub-divertor. Then, a first model is developed and then explored with simulations that aim to discover how well the controllers in the model react to either insufficient levels of 3 He or excessive levels of 3 He. The simulations then explore potential impact from a delay in the measurement (or the response) that would be comparable to the ∼1 s, conductance limited response for the ITER DRGA system, currently in its final design. The simulations show that control is feasible, and that its effectiveness is not significantly impacted by such delay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Progress towards edge-localized mode suppression via magnetic perturbations in hydrogen plasmas

Abstract The suppression of edge-localized modes (ELMs) by applying resonant magnetic perturbations (RMPs) is well studied in low collisionality deuterium plasmas as a measure to reduce transient divertor heat loads. However, ELM suppression has yet to be demonstrated in non-nuclear fuels such as hydrogen and hydrogen + helium mixtures which are the main ion species to be used in the ITER pre-fusion power operation (PFPO) phase. For the first time, attempts have been made to access ELM suppression with RMPs in ITER-like low collisionality hydrogen plasmas at DIII-D and ASDEX Upgrade. The DIII-D experiments focused on operation with injected power slightly above the L–H power threshold similar to the expected conditions in the ITER PFPO phase with limited external heating power. The RMPs were found to trigger H–L backtransitions, which is shown to be avoided by reducing the L–H power threshold by diluting the plasma with helium. The additional helium combined with a larger measured neutral density of hydrogen inside the separatrix compared to ELM suppressed deuterium plasmas precluded access to a pedestal top density below the known RMP-ELM suppression threshold. At ASDEX Upgrade, RMP-ELM suppression has been achieved when the concentration of 1 H in the hydrogen isotope mix is below 40 % . While all known access criteria for RMP-ELM suppression were met above this threshold, full ELM suppression was replaced by strong mitigation. The most prominent difference between the hydrogen and deuterium plasmas was a change of turbulence characteristics in the pedestal where Doppler reflectometry measurements suggest a significant reduction of turbulence even at small hydrogen concentrations. In conclusion, these experiments not only identify issues that may prevent access to RMP-ELM suppression in the ITER PFPO phase, but also highlight missing physics in our current understanding of RMP-ELM suppression such as potentially the role of turbulence in the pedestal gradient region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2013

A Framework for Robust Multivariable Optimization of Integrated Circuits in Space Applications

Application Specific Integrated Circuit (ASIC) design for space applications involves multiple challenges of maximizing performance, minimizing power and ensuring reliable operation in extreme environments. This is a complex multidimensional optimization problem which must be solved early in the development cycle of a system due to the time required for testing and qualification severely limiting opportunities to modify and iterate. Manual design techniques which generally involve simulation at one or a small number of corners with a very limited set of simultaneously variable parameters in order to make the problem tractable are inefficient and not guaranteed to achieve the best possible results within the performance envelope defined by the process and environmental requirements. What is required is a means to automate design parameter variation, allow the designer to specify operational constraints and performance goals, and to analyze the results in a way which facilitates identifying the tradeoffs defining the performance envelope over the full set of process and environmental corner cases. The system developed by the Mixed Signal ASIC Group (MSAG) at the Goddard Space Flight Center is implemented as framework of software modules, templates and function libraries. It integrates CAD tools and a mathematical computing environment, and can be customized for new circuit designs with only a modest amount of effort as most common tasks are already encapsulated. Customization is required for simulation test benches to determine performance metrics and for cost function computation. Templates provide a starting point for both while toolbox functions minimize the code required. Once a test bench has been coded to optimize a particular circuit, it is also used to verify the final design. The combination of test bench and cost function can then serve as a template for similar circuits or be re-used to migrate the design to different processes by re-running it with the new process specific device models. The system has been used in the design of time to digital converters for laser ranging and time-of-flight mass spectrometry to optimize analog, mixed signal and digital circuits such as charge sensitive amplifiers, comparators, delay elements, radiation tolerant dual interlocked (DICE) flip-flops and two of three voter gates.

DuMonthier, Jeffrey↗
NASA NTRS2013

Robust Multivariable Optimization and Performance Simulation for ASIC Design

Application-specific-integrated-circuit (ASIC) design for space applications involves multiple challenges of maximizing performance, minimizing power, and ensuring reliable operation in extreme environments. This is a complex multidimensional optimization problem, which must be solved early in the development cycle of a system due to the time required for testing and qualification severely limiting opportunities to modify and iterate. Manual design techniques, which generally involve simulation at one or a small number of corners with a very limited set of simultaneously variable parameters in order to make the problem tractable, are inefficient and not guaranteed to achieve the best possible results within the performance envelope defined by the process and environmental requirements. What is required is a means to automate design parameter variation, allow the designer to specify operational constraints and performance goals, and to analyze the results in a way that facilitates identifying the tradeoffs defining the performance envelope over the full set of process and environmental corner cases. The system developed by the Mixed Signal ASIC Group (MSAG) at the Goddard Space Flight Center is implemented as a framework of software modules, templates, and function libraries. It integrates CAD tools and a mathematical computing environment, and can be customized for new circuit designs with only a modest amount of effort as most common tasks are already encapsulated. Customization is required for simulation test benches to determine performance metrics and for cost function computation.

DuMonthier, Jeffrey↗
DOE OSTI2021

Simulation of pellet ELM triggering in low-collisionality, ITER-like discharges

3D nonlinear, as well as 2D linear M3D-C1 simulations are used to model ELM triggering by small pellets in DIII-D discharges in the ITER relevant, peeling-limited pedestal stability regime. A critical pellet size threshold is found in both experiment and modeling depending on pedestal conditions, pellet velocity and injection direction. Using radial injection at the outboard midplane, the threshold is determined by M3D-C1 for multiple time slices of a DIII-D low-collisionality discharge that has pellet ELM triggering. Experimental observations show that a larger pellet size than the standard 1.3 mm diameter is necessary for ELM triggering; 1.8 mm pellets triggered several ELMs in cases where a smaller pellet failed. The M3D-C1 simulations are in good agreement with these observations. While the 2D linear simulations give insight into the change of growth rates for various toroidal modes with pellet size, the 3D nonlinear simulations apply a pellet ablation model that mimics the actual injection with good match to the experiment. The 3D nonlinear simulation confirms the pellet ELM triggering for a pellet size larger than the threshold found by the linear simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1974

Computation of steady axisymmetric flow using a one-dimensional time-dependent method

An iterative numerical method for computing steady, three dimensional, viscous, compressible flow fields, about aerodynamic bodies was studied. In order to develop the iterative method economically, the primary emphasis was directed towards supersonic, axisymmetric flow. However, the technique readily extends to three spatial dimensions. The viscous flow field about a cone-cylinder-flare body was calculated and compared to existing experimental data. Numerical predictions of the cone boundary layer and the flow field shock wave structure agreed with corresponding measurements. A separation was calculated at the cylinder-flare junction in six iterations; however, the size of the vortex did not correspond to the measured size. It was estimated that fifty iterations would be required to properly define the vortex. It was concluded that the iteration technique is of limited value for plane two dimensional and axisymmetrix flows, but of great value for three-dimensional flows.

Walitt, L.↗
NASA NTRS1976

Advances in the application of fast semidirect computational methods in transonic flow

The paper uses finite-difference algorithms called 'fast direct elliptic solvers' within an iteration scheme for the rapid solution of the equations of inviscid transonic aerodynamics. The methods are called 'direct' (or semidirect) because the entire computational field is solved at once rather than in successive traverses over the field. These semidirect iterative methods have been limited here to the investigation of two-dimensional steady inviscid flow over airfoils in subsonic free stream.

Martin, E. D.↗
NASA NTRS1992

Schroedinger's radial equation - Solution by extrapolation

A high-accuracy numerical method for the solution of a 1D Schroedinger equation that is suitable for a diatomic molecule, obtained by combining a finite-difference method with iterative extrapolation to the limit, is presently shown to have several advantages over more conventional methods. Initial guesses for the term values are obviated, and implementation of the algorithm is straightforward. The method is both less sensitive to round-off error, and faster than conventional methods for equivalent accuracy. These advantages are illustrated through the solution of Schroedinger's equation for a Morse potential function suited for HCl and a numerically derived Rydberg-Klein-Rees potential function for the X 1Sigma(+) state of CO.

Goorvitch, D.↗
NASA NTRS2003

Factors Affecting the Evolution of Hurricane Erin and the Distributions of Hydrometeors: Role of Microphysical Processes

Fine-resolution simulations of Hurricane Erin 2001 are conducted using the Penn State University/National Center for Atmospheric Research mesoscale model version 3.5 to investigate the role of thermodynamic, boundary layer and microphysical processes in Erin's growth and maintenance, and their effects on the horizontal and vertical distributions of hydrometeors. Through comparison against radar, radiometer, and dropsonde data collected during the Convection and Moisture Experiment 4, it is seen that realistic simulations of Erin are obtained provided that fine resolution simulations with detailed representations of physical processes are conducted. The principle findings of the study are as follows: 1) a new iterative condensation scheme, which limits the unphysical increase of equivalent potential temperature associated with most condensation schemes, increases the horizontal size of the hurricane, decreases its maximum rainfall rate, reduces its intensity, and makes its eye more moist; 2) in general, microphysical parameterization schemes with more categories of hydrometeors produce more intense hurricanes, larger hydrometeor mixing ratios, and more intense updrafts and downdrafts; 3) the choice of coefficients describing hydrometeor fall velocities has as big of an impact on the hurricane simulations as does choice of microphysical parameterization scheme with no clear relationship between fall velocity and hurricane intensity; and 4) in order for a tropical cyclone to adequately intensify, an advanced boundary layer scheme (e.g., Burk-Thompson scheme) must be used to represent boundary layer processes. The impacts of varying simulations on the horizontal and vertical distributions of different categories of hydrometeor species, on equivalent potential temperature, and on storm updrafts and downdrafts are examined to determine how the release of latent heat feedbacks upon the structure of Erin. In general, all simulations tend to overpredict precipitation rate and hydrometeor mixing ratios. The ramifications of these findings for quantitative precipitation forecasts (QPFs) of tropical cyclones are discussed.

McFarquhar, Greg M.↗
NASA NTRS2011

AGN Obscuration Through Dusty Infrared Dominated Flows: Radiation-Hydrodynamics Solution for the Wind - 1

We construct a radiation-hydrodynamics model for the obscuring toroidal structure in active galactic nuclei. In this model the obscuration is produced at parsec scale by a dense, dusty wind which is supported by infrared radiation pressure on dust grains. To find the distribution of radiation pressure, we numerically solve the 2D radiation transfer problem in a flux limited diffusion approximation. We iteratively couple the solution with calculations of stationary 1D models for the wind, and obtain the z-component of the velocity. Our results demonstrate that for AGN luminosities greater than 0.1 L(sub edd) external illumination can support a geometrically thick obscuration via outflows driven by infrared radiation pressure. The terminal velocity of marginally Compton-thin models (0.2 < tau(sub T) < 0.6), is comparable to or greater than the escape velocity. In Compton thick models the maximum value of the vertical component of the velocity is lower than the escape velocity, suggesting that a significant part of our torus is in the form of failed wind. The results demonstrate that obscuration via normal or failed infrared-driven winds is a viable option for the AGN torus problem and AGN unification models. Such winds can also provide an important channel for AGN feedback.

Dorodnitsyn, A.↗
DOE OSTI2026

wa-hls4ml: A Benchmark and Surrogate Models for hls4ml Resource and Latency Estimation

As machine learning (ML) is increasingly implemented in hardware to address real-time challenges in scientific applications, the development of advanced toolchains has significantly reduced the time required to iterate on various designs. These advancements have solved major obstacles, but also exposed new challenges. For example, processes that were not previously considered bottlenecks, such as hardware synthesis, are becoming limiting factors in the rapid iteration of designs. To mitigate these emerging constraints, multiple efforts have been undertaken to develop an ML-based surrogate model that estimates resource usage of ML accelerator architectures. We introduce wa-hls4ml, a benchmark for ML accelerator resource and latency estimation, and its corresponding initial dataset of over 680,000 fully connected and convolutional neural networks, all synthesized using hls4ml and targeting Xilinx FPGAs. The benchmark evaluates the performance of resource and latency predictors against several common ML model architectures, primarily originating from scientific domains, as exemplar models, and the average performance across a subset of the dataset. Additionally, we introduce GNN- and transformer-based surrogate models that predict latency and resources for ML accelerators. We present the architecture and performance of the models and find that the models generally predict latency and resources for the 75% percentile within several percent of the synthesized resources on the synthetic test dataset.

Hawks, Benjamin [Fermilab] (ORCID:0000000157000288↗
DOE OSTI2025

wa-hls4ml: A Benchmark and Surrogate Models for hls4ml Resource and Latency Estimation

As machine learning (ML) is increasingly implemented in hardware to address real-time challenges in scientific applications, the development of advanced toolchains has significantly reduced the time required to iterate on various designs. These advancements have solved major obstacles, but also exposed new challenges. For example, processes that were not previously considered bottlenecks, such as hardware synthesis, are becoming limiting factors in the rapid iteration of designs. To mitigate these emerging constraints, multipleefforts have been undertaken to develop an ML-based surrogate model that estimates resource usage of synthesized ML accelerator architectures. We introduce wa-hls4ml, a benchmark for ML accelerator resource and latency estimation, and its corresponding initial dataset of over 680 000 fully connected and convolutional neural networks, all synthesized using hls4ml and targeting Xilinx FPGAs. The benchmark evaluates the performance of resource and latency predictors against several common ML model architectures, primarily originating from scientific domains, as exemplar models, and the average performance across a subset of the dataset. Additionally, we introduce GNN- and transformer-based surrogate models that predict latency and resources for ML accelerators. We present the architecture and performance of the models and find that the models generally predict latency and resources for the 75% percentile within several percent of the synthesized resources on the synthetic test dataset.

Hawks, Benjamin G. [Fermilab]↗
DOE OSTI2025

wa-hls4ml and lui-gnn: A benchmark and GNN-based surrogate model for hls4ml resource and latency estimation

As machine learning (ML) increasingly serves as a tool for addressing real-time challenges in scientific applications, the development of advanced tooling has significantly reduced the time required to iterate on various designs. These advancements have solved major obstacles, but also exposed new challenges. For example, processes that were not previously considered bottlenecks, such as model synthesis, are now becoming limiting factors in the rapid iteration of designs. To reduce these emerging constraints, multiple efforts are being launched toward designing an ML-based surrogate model that estimates resource usage of synthesized accelerator architectures. This model would reduce the design iteration time, especially when designing within a set of given hardware constraints. This approach shows considerable potential, but as it stands, the effort is early and would benefit from coordination and standardization to assist future work as it emerges. We introduce wa-hls4ml, a benchmark for ML accelerator resource and latency estimation, and its corresponding initial dataset of more than 100,000 fully connected neural networks, all synthesized using hls4ml and targeting Xilinx FPGAs. In addition to the resource utilization and latency data provided, the dataset includes generated artifacts and log files for many of the synthesized neural networks, in order to support future research in ML-based code generation. The benchmark evaluates the performance of resource and latency predictors against several common ML model architectures, primarily originating from scientific domains, as exemplar models, as well as the average performance across a subset of the dataset. We measure the performance of a given predictor model through multiple metrics, including $R^2$ score and SMAPE on regression tasks, as well as inference time to further characterize the estimator under test. Additionally, we introduce the latency/utilization inference graph neural network (lui-gnn), a surrogate model that uses a graph neural network to represent input architectures in the form of a directed graph. This graph representation allows for a diverse set of model architectures to all be effectively handled by a surrogate model. We present the architecture and performance of the model, as evaluated by the new proposed benchmark, including SMAPE, $R^2$ score, and inference times, and find that lui-gnn generally predicts latency and utilization for the 75\% quantile within several percent of the synthesized resources on the synthetic test dataset, indicating that this approach of estimating resource and latency via a surrogate models has promise and warrants further research.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Assessment of the burning-plasma operational space in ITER by using a control-oriented core-SOL-divertor model

In future tokamaks, the control of burning plasmas will require careful regulation of the plasma density and temperature. Along with the design of effective burn-control systems, understanding how the fusion power varies in the density-temperature space is vital for the operation of fusion power plants. Here in this work, the steady-state operational space of ITER is studied using a control-oriented core-plasma model coupled to a two-point model of the scrape-off-layer (SOL) and divertor regions. The two models are coupled through the exchange of input-output parameters. The deuterium and tritium recycling from the wall are output parameters of the SOL-divertor model that are used as input parameters in the core-plasma density balance. Furthermore, the separatrix temperature, which is an output parameter of the SOL-divertor model, is incorporated into the radial core-plasma temperature profiles. Therefore, the temperature-dependent power balance of the plasma core is intimately linked to the SOL-divertor model. Both the power entering the SOL from the core, as determined by the core-plasma power balance, and the separatrix density, as dictated by the core-plasma density balance, are input parameters to the SOL-divertor model. They are control knobs in the SOL-divertor model that can be regulated using the core-plasma actuators: auxiliary power and pellet injection. There are various operational limitations, such as the saturation of the aforementioned actuators, that will prevent ITER from accessing certain high-fusion plasma regimes. The achievable tritium concentration in the fueling lines and the maximum sustainable heat load on the divertor will impose further restrictions. By accounting for these limitations, the ITER operational space is computed based on the coupled core-SOL-divertor model and visualized using Plasma Operation Contour (POPCON) plots that map performance metrics, such as the fusion to auxiliary power ratio, over the density-temperature space. Comparisons are drawn between plasmas with different recycling, confinement, and SOL-divertor conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1999

Ames Optimized TCA Configuration

Configuration design at Ames was carried out with the SYN87-SB (single block) Euler code using a 193 x 49 x 65 C-H grid. The Euler solver is coupled to the constrained (NPSOL) and the unconstrained (QNMDIF) optimization packages. Since the single block grid is able to model only wing-body configurations, the nacelle/diverter effects were included in the optimization process by SYN87's option to superimpose the nacelle/diverter interference pressures on the wing. These interference pressures were calculated using the AIRPLANE code. AIRPLANE is an Euler solver that uses a unstructured tetrahedral mesh and is capable of computations about arbitrary complete configurations. In addition, the buoyancy effects of the nacelle/diverters were also included in the design process by imposing the pressure field obtained during the design process onto the triangulated surfaces of the nacelle/diverter mesh generated by AIRPLANE. The interference pressures and nacelle buoyancy effects are added to the final forces after each flow field calculation. Full details of the (recently enhanced) ghost nacelle capability are given in a related talk. The pseudo nacelle corrections were greatly improved during this design cycle. During the Ref H and Cycle 1 design activities, the nacelles were only translated and pitched. In the cycle 2 design effort the nacelles can translate vertically, and pitch to accommodate the changes in the lower surface geometry. The diverter heights (between their leading and trailing edges) were modified during design as the shape of the lower wing changed, with the drag of the diverter changing accordingly. Both adjoint and finite difference gradients were used during optimization. The adjoint-based gradients were found to give good direction in the design space for configurations near the starting point, but as the design approached a minimum, the finite difference gradients were found to be more accurate. Use of finite difference gradients was limited by the CPU time limit available on the Cray machines. A typical optimization run using finite difference gradients can use only 30 to 40 design variables and one optimization iteration within the 8 hour queue limit for the chosen grid size and convergence level. The efficiency afforded by the adjoint method allowed for 50-120 design variables and 5-10 optimization iterations in the 8 hour queue. Geometric perturbations to the wing and fuselage were made using the Hicks/Henne (HH) shape functions. The HH functions were distributed uniformly along the chords of the wing defining sections and lofted linearly. During single-surface design, constraints on thickness and volume at selected wing stations were imposed. Both fuselage camber and cross-sectional area distributions were permitted to change during design. The major disadvantage to the use of these functions is the inherent surface waviness produced by repeated use of such functions. Many smoothing operations were required following optimization runs to produce a configuration with reasonable smoothness. Wagner functions were also used on the wing sections but were never used on the fuselage. The Wagner functions are a family of increasingly oscillatory functions that have also been used extensively in airfoil design. The leading and trailing edge regions of the wing were designed by use of polynomial and monomial functions respectively. Twist was attempted but was abandoned because of little performance improvement available from changing the baseline twist.

Cliff, Susan E.↗
DOE OSTI2021

Iterative hybrid manufacture of a titanium alloy component

Here, this paper describes an iterative hybrid (additive + subtractive) manufacturing approach for a titanium alloy (Ti6Al4V) part using a laser hotwire directed energy deposition system (LHWDED) and a traditional four-axis milling machine tool. The term iterative hybrid manufacturing is used to described hybrid manufacturing where the additive and subtractive operations occur in multiple stages rather than sequentially. It is currently common to produce an entire part by sequential hybrid manufacturing by additively manufacturing (AM) an entire preform geometry that then requires post processing by another machine tool to create final part features. By contrast, a part produced by iterative hybrid manufacturing (IHM) does not produce the entire preform geometry in a single AM process. Instead, a portion of the entire preform geometry is manufactured by an AM process, then that portion is transferred to another machine tool which creates features in that portion, and then that machined portion is transferred back to the AM machine to complete another AM process. There is no limit to the number of iterations that an IHM process can have. IHM offers several advantages over sequential hybrid manufacturing such as the use of shorter and stiffer subtractive tooling, better access to part geometries that require subtractive processes, and the separation of the AM heat source from the subtractive machine tool. A titanium alloy demonstration part was successfully manufactured by IHM with three iterations using a shared pallet system between the AM machine tool and the subtractive machine tool.

Hybrid manufacturing↗
DOE OSTI2021

The role of edge plasma parameters in H-mode density limit on the JET-ILW

A study of a dataset of JET H-mode plasma with the Be/W ITER-like wall (JET-ILW) shows that reaching the edge MHD ballooning limit leads to confinement degradation. However, unlike JET plasmas with a carbon wall (JET-C), the JET-ILW plasmas stay in a marginal dithering phase for a relatively long period, associated with a higher (≈20%) H-mode density limit (HDL) than JET-C equivalents. This suggests that ITER could be operated in H-mode with higher density than the scaling based on carbon wall devices, but likely with a dithering phase plasma with lower confinement. A new, reliable estimator for JET E r, min has been derived by combining HRTS measurements of pedestal gradient and edge-SOL decay lengths. JET radial E r ETB wells are observed in the range of -15 to -60 kV m -1 in high performance H-modes, consistent with previous CXRS results in ASDEX Upgrade. The results imply that a higher positive E &#x00d7; B shear in the near SOL plays a role in sustaining a marginal phase in JET-ILW which leads to a higher HDL than that in JET-C. The results of the JET-ILW dataset show agreement with the Goldston finite collisionality HD model for SOL broadening at high collisionality. A hypothesis for the dithering H-mode phase is proposed: as n e,SOL increases, ν *,SOL increases, SOL broadens, E r shear decreases, triggers L-mode; n e drops, ν *,SOL decreases, SOL becomes narrower, and E r shear increases, triggering H-mode, resulting in a cycle of H–L–H- oscillations. For burning plasma devices, such as ITER, operating just below the MHD limit for the dithering phase could be a promising regime for maximising core density, and fusion performance while minimising plasma-material interaction. The oscillatory signal during the dithering phase could be used as a precursor of undesirable plasma performance for control purposes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗