Engineering topics

Leonard, A.

Publications and source records attributed to Leonard, A..

At least 19 records

DOE OSTI2024

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↗

DOE OSTI2022

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↗

DOE OSTI2021

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↗

DOE OSTI2020

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↗

NASA NTRS1989

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.↗

NASA NTRS1988

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.↗

NASA NTRS1988

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.↗

NASA NTRS1986

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.↗

NASA NTRS1986

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.↗

NASA NTRS1986

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.↗

NASA NTRS1985

Evolution of a curved vortex filament into a vortex ring

Moin, P.↗

NASA NTRS1985

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.↗

NASA NTRS1985

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.↗

NASA NTRS1985

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.↗

NASA NTRS1984

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.↗

NASA NTRS1983

A spectral numerical method for the Navier-Stokes equations with applications to Taylor-Couette flow

A new spectral method for solving the incompressible Navier-Stokes equations in a plane channel and between concentric cylinders is presented. The method uses spectral expansions which inherently satisfy the boundary conditions and the continuity equation and yield banded matrices which are efficiently solved at each time step. In addition, the number of dependent variables is reduced, resulting in a reduction in computer memory requirements. Several test problems have been computed for the channel flow and for flow between concentric cylinders, including Taylor-Couette flow with axisymmetric Taylor vortices and wavy vortices. In all cases, agreement with available experimental and theoretical results is very good.

Moser, R. D.↗

NASA NTRS1983

Numerical simulation of vortex breakdown by the vortex-filament method

The vortex filament method was applied to the simulation of vortex breakdown. The principal vortex region was represented by multiple filaments, and an axial velocity component was induced by a spiral winding of the filaments. First, an accuracy check was performed for a cylindrical swirling flow with simple analytical expressions for the axial and theta velocities. The result suggests that the flow field is simulated to any accuracy by increasing the number of filaments. Second, an axisymmetric type vortex breakdown was simulated, with experimental data serving as upstream conditions. The calculated axial and theta velocity contours show the breakdown of the vortex, including a rapid change in the vortex core, followed axially by a recovery zone and then a second breakdown. When three dimensional initial data are used the second breakdown appears to be of the spiral type in correspondence with experimental observations. The present method is easily used to simulate other types of vortex breakdown or other vortex flows with axial velocity.

Nakamura, Y.↗

NASA NTRS1983

Numerical simulation of vortex breakdown by the vortex-filament method

The vortex-filament method was applied to the simulation of vortex breakdown. The principal vortex region was represented by multiple filaments, and an axial velocity component was induced by a spiral winding of the filaments. An accuracy check was performed for a cylindrical swirling flow with simple analytical expressions for the axial and theta velocities. The result suggests that the flow field can be simulated to any accuracy by increasing the number of filaments. An axisymmetric-type vortex breakdown was simulated, with experimental data serving as upstream conditions. The calculated axial- and theta-velocity contours show the breakdown of the vortex, including a rapid change in the vortex core, followed axially by a recovery zone and then a second breakdown. When three dimensional initial data are used the second breakdown appears to be of the spiral type in correspondence with experimental observations. The present method can easily be used to simulate other types of vortex breakdown or other vortex flows with axial velocity.

Nakamura, Y.↗