Engineering Papers⌕ Search

Engineering topics

Leonard, A.

Publications and source records attributed to Leonard, A..

At least 19 records

Real-time plasma equilibrium reconstruction and shape control for the MAST Upgrade tokamak

Real-time magnetic control has been developed to deliver precise control of multiple plasma shape parameters for advanced divertor configurations, including double-null, Super-X, X-point target and X-divertor for the first time on the MAST Upgrade (MAST-U) spherical tokamak. Successful real-time magnetic equilibrium control of different plasma shape variables has been accomplished in the 2022–2023 MAST-U experimental campaign for the advanced MAST-U divertor configurations. Application of the MAST-U boundary reconstruction algorithm, LEMUR, is described and compared with off-line equilibrium reconstruction and diagnostic measurements. The process of design and verification of the axisymmetric magnetic control schemes using a suite of control analysis tools (known collectively as TokSys) is also described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impurity leakage and radiative cooling in the first nitrogen and neon seeding study in the closed DIII-D SAS configuration

A comparative study of nitrogen versus neon has been carried out to analyze the impact of the two radiative species on power dissipation, SOL impurity distribution, divertor and pedestal characteristics. The experimental results show that N remains compressed in the divertor, thereby providing high radiative losses without affecting the pedestal profiles and displacing carbon as dominant radiator. Neon, instead, radiates more upstream than N thus reducing the power flux through the separatrix leading to a reduced ELM frequency and compression in the divertor. A significant amount of neon is measured in the plasma core leading to a steeper density gradient. The different behavior between the two impurities is confirmed by SOLPS-ITER modeling which for the first time at DIII-D includes multiple impurity species and a treatment of full drifts, currents and neutral–neutral collisions. The impurity transport in the SOL is studied in terms of the parallel momentum balance showing that N is mostly retained in the divertor whereas Ne leaks out consistent with its higher ionization potential and longer mean free path. This is also in agreement with the enrichment factor calculations which indicate lower divertor enrichment for neon. The strong ionization source characterizing the SAS divertor causes a reversal of the main ions and impurity flows. The flow reversal together with plasma drifts and the effect of the thermal force contribute significantly in the shift of the impurity stagnation point affecting impurity leakage. This work provides a demonstration of the impurity leakage mechanism in a closed divertor structure and the consequent impact on pedestal. Since carbon is an intrinsic radiator at DIII-D, in this paper we have also demonstrated the different role of carbon in the N vs Ne seeded cases both in the experiments and in the numerical modeling. Here, carbon contributes more when neon seeding is injected compared to when nitrogen is used. Finally, the results highlight the importance of accompanying experimental studies with numerical modeling of plasma flows, drifts and ionization profile to determine the details of the SOL impurity transport as the latter may vary with changes in divertor regime and geometry. In the cases presented here, plasma drifts and flow reversal caused by high level of closure in the slot upper divertor at DIII-D play an important role in the underlined mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical assessment of the new V-shape small-angle slot divertor on DIII-D

The small-angle slot (SAS) divertor of the DIII-D tokamak, and its upcoming upgrade, the V-shape small-angle slot (SAS-V) divertor, are numerically investigated using the SOLPS-ITER code package, including the effect of particle drifts, for a range of plasma density, heating power, strike point position in the slot, and for both magnetic field directions. The simulations show that the electron temperature near the strike point is reduced in SAS-V compared to SAS, for both magnetic field directions, such that SAS-V achieves divertor detachment at a lower value of the outboard mid-plane separatrix electron density. The detachment threshold is lower because the V-shape focuses recycling neutrals on the V-end, densifying and cooling the plasma in the slot. At sufficiently high density, the V-shape also reduces the radial gradient of the temperature profile at the target, which in turns reduces the radial electric field and the E x B drift velocities, further densifying and cooling the plasma in the slot and leading to detachment. The V-shape effect, however, is reduced for higher heating power. With more heating power, the detachment density increases, reducing the ionization mean free path of recycled neutrals, which therefore become less sensitive to target shape changes. Furthermore, this suggests that in a fusion reactor, where the heating power is high, optimization of the divertor target shape needs to be combined with other strategies to lower the detachment density, such as in-slot injection of low-Z impurities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Plasma flux expansion control on the DIII-D tokamak

A new controller has been developed with help of the flexible divertor poloidal-field coil set of the DIII-D tokamak, to aid in the precise control of the flux expansion in the scrape-off layer (SOL). Here, the single-input multiple-output (SIMO) architecture ensures flexibility through a complementary set of orthogonal actuator direction to guarantee minimum effect on existing controlled variables, e.g. radial and vertical position control of the X-point. A non-linear free-boundary simulation code (GSevolve) is used for simulating the closed-loop response and for verifying the implementation of the control algorithm on the DIII-D plasma control system (PCS). First results of the experimental commissioning of the new controller during 2020 DIII-D campaign are also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Layout and assembly technique of the GEM chambers for the upgrade of the CMS first muon endcap station

Triple-GEM detector technology was recently selected by CMS for a part of the upgrade of its forward muon detector system as GEM detectors provide a stable operation in the high radiation environment expected during the future High-Luminosity phase of the Large Hadron Collider (HL-LHC). In a first step, GEM chambers (detectors) will be installed in the innermost muon endcap station in the 1.6<|η|<2.2 pseudo-rapidity region, mainly to control level-1 muon trigger rates after the second LHC Long Shutdown. These new chambers will add redundancy to the muon system in the η-region where the background rates are high, and the bending of the muon trajectories due to the CMS magnetic field is small. A novel construction technique for such chambers has been developed in such a way where foils are mounted onto a single stack and then uniformly stretched mechanically, avoiding the use of spacers and glue inside the active gas volume. We describe the layout, the stretching mechanism and the overall assembly technique of such GEM chambers.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Operational Experience With the GEM Detector Assembly Lines for the CMS Forward Muon Upgrade

The CMS Collaboration has been creating large-area triple-gas electron multiplier (GEM) detectors to be installed in the muon Endcap regions of the CMS experiment in 2019 to maintain forward muon trigger and tracking performance at the High-Luminosity upgrade of the Large Hadron Collider (LHC); 10 preproduction detectors were built at CERN to commission the first assembly line and the quality controls (QCs). These were installed in the CMS detector in early 2017 and participated in the 2017 LHC run. The collaboration has prepared several additional assembly and QC lines for distributed mass production of 160 GEM detectors at various sites worldwide. In 2017, these additional production sites have optimized construction techniques and QC procedures and validated them against common specifications by constructing additional preproduction detectors. Using the specific experience from one production site as an example, we report on how the QCs make use of independent hardware and trained personnel to ensure fast and reliable production. Preliminary results on the construction status of CMS GEM detectors are presented with details of the assembly sites involvement.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Improving Weak Lensing Mass Map Reconstructions using Gaussian and Sparsity Priors: Application to DES SV

Mapping the underlying density field, including non-visible dark matter, using weak gravitational lensing measurements is now a standard tool in cosmology. Due to its importance to the science results of current and upcoming surveys, the quality of the convergence reconstruction methods should be well understood. We compare three methods: Kaiser–Squires (KS), Wiener filter, and Glimpse. Kaiser–Squires is a direct inversion, not accounting for survey masks or noise. The Wiener filter is well-motivated for Gaussian density fields in a Bayesian framework. Glimpse uses sparsity, aiming to reconstruct non-linearities in the density field. We compare these methods with several tests using public Dark Energy Survey (DES) Science Verification (SV) data and realistic DES simulations. The Wiener filter and Glimpse offer substantial improvements over smoothed Kaiser–Squires with a range of metrics. Both the Wiener filter and Glimpse convergence reconstructions show a 12 per cent improvement in Pearson correlation with the underlying truth from simulations. To compare the mapping methods’ abilities to find mass peaks, we measure the difference between peak counts from simulated ΛCDM shear catalogues and catalogues with no mass fluctuations (a standard data vector when inferring cosmology from peak statistics); the maximum signal-to-noise of these peak statistics is increased by a factor of 3.5 for the Wiener filter and 9 for Glimpse. With simulations, we measure the reconstruction of the harmonic phases; the phase residuals’ concentration is improved 17 per cent by Glimpse and 18 per cent by the Wiener filter. The correlationbetween reconstructions from data and foreground redMaPPer clusters is increased 18 per cent by the Wiener filter and 32 per cent by Glimpse.

79 ASTRONOMY AND ASTROPHYSICS↗

Improved vortex methods for three-dimensional flows

Robust numerical methods are developed for three-dimensional incompressible vortical flows, using Lagrangian vortex elements. A successful scheme must be able to handle regions of intense vortex stretching and vortex reconnection with reasonable accuracy (without diverging). Here, consideration is given to vortex particles, also commonly called vortons or vortex sticks. The following issues are discussed: (1) use of delta-function elements and weak solutions of the vorticity equation; (2) use of smoothed elements and the choice of the smoothing function; (3) representation of viscous effects and the redistribution of element strength; and (4) conservation laws (are they satisfied?). The various proposed schemes have been tested on flows involving a strong interaction between two vortex rings.

Winckelmans, G.↗

Weak solutions of the three-dimensional vorticity equation with vortex singularities

The extension of the concept of vortex singularities, developed by Saffman and Meiron (1986) for the case of two-dimensional point vortices in an incompressible vortical flow, to the three-dimensional case of vortex sticks (vortons) is investigated analytically. The derivation of the governing equations is explained, and it is demonstrated that the formulation obtained conserves total vorticity and is a weak solution of the vorticity equation, making it an appropriate means for representing three-dimensional vortical flows with limited numbers of vortex singularities.

Winckelmans, G.↗

Robust vortex methods for three-dimensional incompressible flows

Vortex methods for the numerical simulation of incompressible three-dimensional flows at high Reynolds number are discussed. Particular emphasis is placed on schemes that prevent the excessive generation of mesh elements, include viscous effects, and are relatively inexpensive even as the number of elements becomes large. Several simulations are presented involving the interaction of vortex rings.

Chua, K.↗

Vortex simulation of forced mixing layers

Two-dimensional, spatially growing, turbulent mixing layers are simulated numerically by a vortex method and the results are compared with those determined experimentally. The effects of artificial forcing on flow development are also studied. Many of the flow features which have been observed experimentally are reproduced, and good quantitative agreements between experiments and computations are obtained.

Inoue, O.↗

Internal structure of a vortex breakdown

An axisymmetric vortex breakdown was well simulated by the vortex filament method. The agreement with the experiment was qualitatively good. In particular, the structure in the interior of the vortex breakdown was ensured to a great degree by the present simulation. The second breakdown, or spiral type, which occurs downstream of the first axisymmetric breakdown, was simulated more similarly to the experiment than before. It shows a kink of the vortex filaments and strong three-dimensionality. Furthermore, a relatively low velocity region was observed near the second breakdown. It was also found that it takes some time for this physical phenomenon to attain its final stage. The comparison with the experiment is getting better as time goes on. In this paper, emphasis is placed on the comparison of the simulated results with the experiment. The present results help to make clear the mechanism of a vortex breakdown.

Nakamura, Y.↗

Evolution of a curved vortex filament into a vortex ring

The deformation of a hairpin-shaped vortex filament under self-induction and in the presence of shear is studied numerically using the Biot-Savart law. It is shown that the tip region of an elongated hairpin vortex evolves into a vortex ring and that the presence of mean shear impedes the process. Evolution of a finite-thickness vortex sheet under self-induction is also investigated using the Navier-Stokes equations. The layer evolves into a hairpin vortex which in turn produces a vortex ring of high Reynolds stress content. These results indicate a mechanism for the generation of ring vortices in turbulent shear flows, and a link between the experimental and numerical observation of hairpin vortices and the observation of ring vortices in the outer regions of turbulent boundary layers.

Moin, P.↗

Evolution of a curved vortex filament into a vortex ring

The deformation of a hairpin-shaped vortex filament under self-induction and in the presence of shear is studied numerically using the Biot-Savart law. It is shown that the tip region of an elongated hairpin vortex evolves into a vortex ring and that the presence of mean shear impedes the process. Evolution of a finite-thickness vortex sheet under self-induction is also investigated using the Navier-Stokes equations. The layer evolves into a hairpin vortex which in turn produces a vortex ring of high Reynolds stress content. These results indicate a mechanism for the generation of ring vortices in turbulent shear flows, and a link between the experimental and numerical observation of hairpin vortices and the observation of ring vortices in the outer regions of turbulent boundary layers.

Moin, P.↗

Vortex breakdown simulation

A vortex breakdown was simulated by the vortex filament method, and detailed figures are presented based on the results. Deformations of the vortex filaments showed clear and large swelling at a particular axial station which implied the presence of a recirculation bubble at that station. The tendency for two breakdowns to occur experimentally was confirmed by the simulation, and the jet flow inside the bubble was well simulated. The particle paths spiralled with expansion, and the streamlines took spiral forms at the breakdown with expansion.

Nakamura, Y.↗

Computing three-dimensional incompressible flows with vortex elements

The techniques, capabilities and applicability of numerical models of three-dimensional, unsteady vortical flows with high Re are assessed. Vorticity is calculated only in appropriate regions and the velocity field is derived from the boundary conditions. Vorticity is assumed to take the shape of tubes with uniform core structures in the case of turbulence. The efforts being made to simplify equations for dense collections of vortex filaments in order to make them tractable to computer simulations are described. The effectiveness of vorticity arrow representations for accurately describing vorticity fields near surfaces is discussed, along with Lagrangian vortex elements, which may be of use in modelling the rotational part of flows around bluff bodies, nonuniform density flows and chemically reacting flows.

Leonard, A.↗

Direct numerical simulation of equilibrium turbulent boundary layers

This paper describes the simulation of turbulent boundary layers by direct numerical solution of the three-dimensional, time-dependent Navier-Stokes equations, using a spectral method. The flow is incompressible, with Re sub delta = 1000. The equations are written in the self-similar coordinate system and periodic streamwise and spanwise boundary conditions are imposed. A family of nine 'equilibrium' boundary layers, from the strongly accelerated 'sink' flow to Stratford's separating boundary layer is treated. Good general agreement with experiments is observed. The effects of pressure gradients on the structures and statistics, both in the wall and wake regions are discussed.

Spalart, P. R.↗

Simulating a Three-Dimensional Flow in Pipes

Numerical method simulates in unsteady, incompressible flow in axisymmetric pipe in three dimensions. Method relies on velocity-vector expansion technique in which each vector function in expansion set is divergence-free and satisfies boundary conditions for viscous flow.

Leonard, A.↗