-
ITER-IA 3D MHD Simulations of Shattered Pellet Injection(SPI)- D1.1 Optimization of the SPI model
Authors:
Charlson. C. Kim,
B. C. Lyons,
Y. Q. Liu,
J. T. McClenaghan,
P. B. Parks,
L. L. Lao
Abstract:
This report is in partial fulfillment of deliverable D1.1 Optimization of the SPI model and summarizes axisymmetric ITER SPI parameter scans performed by the NIMROD code for several ITER equilibria. These axisymmetric parameter scans are to assess the sensitivity of various injection parameters in preparation for 3D MHD SPI simulations. The scans are comprised of 5 scenarios: S1 - fragment size sc…
▽ More
This report is in partial fulfillment of deliverable D1.1 Optimization of the SPI model and summarizes axisymmetric ITER SPI parameter scans performed by the NIMROD code for several ITER equilibria. These axisymmetric parameter scans are to assess the sensitivity of various injection parameters in preparation for 3D MHD SPI simulations. The scans are comprised of 5 scenarios: S1 - fragment size scan : 3 uniform pencil beam, 1 distributed size pencil beam S2 - velocity scan : v = [250,500,750]m/s S3 - velocity dispersion scan : dv/v = [0.2,0.4] (linear distribution) S4 - poloidal extent of plume : [15' ,45' ] (linear distribution) (dv/v=0.2) S5 - poloidal injection angle : +/-[20' ,45' ] (dv/v=0.2) These scans are performed with several ITER equilibria representative of the operating range, from low current and thermal energy (H123 5MA, 29MJ Hydrogen H-mode) to high current and high thermal energy (DT24 15MA, 370MJ D-T H-mode).
△ Less
Submitted 19 October, 2023;
originally announced October 2023.
-
ITER-IA 3D MHD Simulations of Shattered Pellet Injection(SPI) -- D1.3 Code Validation (DIII-D)
Authors:
Charlson. C. Kim,
T. Bechtel,
J. L. Herfindal,
B. C. Lyons,
Y. Q. Liu,
P. B. Parks,
L. Lao
Abstract:
This report is in partial fulfillment of deliverable D1.3 Code Validation (DIII-D). These simulations focus on thermal quench phase of the SPI mitigation and are not typically carried beyond it to the current spike and subsequent current quench. NIMROD SPI simulations[1] are validated against DIII-D experiments. The target plasma for these simulations is DIII-D 160606@02990ms.
This report is in partial fulfillment of deliverable D1.3 Code Validation (DIII-D). These simulations focus on thermal quench phase of the SPI mitigation and are not typically carried beyond it to the current spike and subsequent current quench. NIMROD SPI simulations[1] are validated against DIII-D experiments. The target plasma for these simulations is DIII-D 160606@02990ms.
△ Less
Submitted 19 October, 2023;
originally announced October 2023.
-
Lagrangian particle model for 3D simulation of pellets and SPI fragments in tokamaks
Authors:
R. Samulyak,
S. Yuan,
N. Naitlho,
P. B. Parks
Abstract:
A 3D numerical model for the ablation of pellets and shattered pellet injection (SPI) fragments in tokamaks in the plasma disruption mitigation and fueling parameter space has been developed based on the Lagrangian particle code [R. Samulyak, X. Wang, H.-S. Chen, Lagrangian Particle Method for Compressible Fluid Dynamics, J. Comput. Phys., 362 (2018), 1-19]. The pellet code implements the low magn…
▽ More
A 3D numerical model for the ablation of pellets and shattered pellet injection (SPI) fragments in tokamaks in the plasma disruption mitigation and fueling parameter space has been developed based on the Lagrangian particle code [R. Samulyak, X. Wang, H.-S. Chen, Lagrangian Particle Method for Compressible Fluid Dynamics, J. Comput. Phys., 362 (2018), 1-19]. The pellet code implements the low magnetic Reynolds number MHD equations, kinetic models for the electronic heating, a pellet surface ablation model, an equation of state that supports multiple ionization states, radiation, and a model for grad-B drift of the ablated material across the magnetic field. The Lagrangian particle algorithm is highly adaptive, capable of simulating a large number of fragments in 3D while eliminating numerical difficulties of dealing with the tokamak background plasma. The code has achieved good agreement with theory for spherically symmetric ablation flows. Axisymmetric simulations of neon and deuterium pellets in magnetic fields ranging from 1 to 6 Tesla have been compared with previous simulations using the FronTier code, and very good agreement has also been obtained. The main physics contribution of the paper is a detailed study of the influence of 3D effects, in particular grad-B drift, on pellet ablation rates and properties of ablation clouds. Smaller reductions of ablation rates in magnetic fields compared to axially symmetric simulations have been demonstrated because the ablated material is not confined to narrowing channels in the presence of grad-B drift. Contribution of various factors in the grad-B drift model has also been quantified.
△ Less
Submitted 26 November, 2020; v1 submitted 13 November, 2020;
originally announced November 2020.
-
Near-field simulations of pellet ablation for disruptions mitigation in tokamaks
Authors:
Nicolas Bosviel,
Paul B. Parks,
Roman Samulyak
Abstract:
Detailed numerical studies of the ablation of a single neon pellet in the plasma disruption mitigation parameter space have been performed. Simulations were carried out using FronTier, a hydrodynamic and low magnetic Reynolds number MHD code with explicit tracking of material interfaces. FronTier's physics models resolve the pellet surface ablation and the formation of a dense, cold cloud of ablat…
▽ More
Detailed numerical studies of the ablation of a single neon pellet in the plasma disruption mitigation parameter space have been performed. Simulations were carried out using FronTier, a hydrodynamic and low magnetic Reynolds number MHD code with explicit tracking of material interfaces. FronTier's physics models resolve the pellet surface ablation and the formation of a dense, cold cloud of ablated material, the deposition of energy from hot plasma electrons passing through the ablation cloud, expansion of the ablation cloud along magnetic field lines and the radiation losses. A local thermodynamic equilibrium model based on Saha equations has been used to resolve atomic processes in the cloud and Redlich-Kwong corrections to the ideal gas equation of state for cold and dense gases have been used near the pellet surface. The FronTier pellet code is the next generation of the code described in [R. Samulyak, T. Lu, P. Parks, Nuclear Fusion, (47) 2007, 103--118]. It has been validated against the semi-analytic improved Neutral Gas Shielding model in the 1D spherically symmetric approximation. Main results include quantification of the influence of atomic processes and Redlich-Kwong corrections on the pellet ablation in spherically symmetric approximation and verification of analytic scaling laws in a broad range of pellet and plasma parameters. Using axially symmetric MHD simulations, properties of ablation channels and the reduction of pellet ablation rates in magnetic fields of increasing strength have been studied. While the main emphasis has been given to neon pellets for the plasma disruption mitigation, selected results on deuterium fueling pellets have also been presented.
△ Less
Submitted 15 September, 2020;
originally announced September 2020.
-
Why aren't more theories named after women?: Teaching women's history in physics
Authors:
Beth Parks
Abstract:
Barriers to women's education and employment in Europe and the U.S. in the nineteenth century made it unlikely that any women would be among the few physicists whose ideas are taught in high school and college courses. This paper explores the social settings in which three influential physicists worked--James Clerk Maxwell, Robert Millikan, and Albert Einstein--to better understand the limited opp…
▽ More
Barriers to women's education and employment in Europe and the U.S. in the nineteenth century made it unlikely that any women would be among the few physicists whose ideas are taught in high school and college courses. This paper explores the social settings in which three influential physicists worked--James Clerk Maxwell, Robert Millikan, and Albert Einstein--to better understand the limited opportunities available to women. By acknowledging and explaining why there weren't more women among these founding physicists, instructors may help students understand the barriers that still exist and feel more empowered to overcome them and pursue physics as a career.
△ Less
Submitted 1 July, 2020;
originally announced July 2020.
-
Studying the Underlying Event in Drell-Yan and High Transverse Momentum Jet Production at the Tevatron
Authors:
The CDF Collaboration,
T. Aaltonen,
J. Adelman,
B. Alvarez Gonzalez,
S. Amerio,
D. Amidei,
A. Anastassov,
A. Annovi,
J. Antos,
G. Apollinari,
A. Apresyan,
T. Arisawa,
A. Artikov,
J. Asaadi,
W. Ashmanskas,
A. Attal,
A. Aurisano,
F. Azfar,
W. Badgett,
A. Barbaro-Galtieri,
V. E. Barnes,
B. A. Barnett,
P. Barria,
P. Bartos,
G. Bauer
, et al. (554 additional authors not shown)
Abstract:
We study the underlying event in proton-antiproton collisions by examining the behavior of charged particles (transverse momentum pT > 0.5 GeV/c, pseudorapidity |η| < 1) produced in association with large transverse momentum jets (~2.2 fb-1) or with Drell-Yan lepton-pairs (~2.7 fb-1) in the Z-boson mass region (70 < M(pair) < 110 GeV/c2) as measured by CDF at 1.96 TeV center-of-mass energy. We u…
▽ More
We study the underlying event in proton-antiproton collisions by examining the behavior of charged particles (transverse momentum pT > 0.5 GeV/c, pseudorapidity |η| < 1) produced in association with large transverse momentum jets (~2.2 fb-1) or with Drell-Yan lepton-pairs (~2.7 fb-1) in the Z-boson mass region (70 < M(pair) < 110 GeV/c2) as measured by CDF at 1.96 TeV center-of-mass energy. We use the direction of the lepton-pair (in Drell-Yan production) or the leading jet (in high-pT jet production) in each event to define three regions of η-φspace; toward, away, and transverse, where φis the azimuthal scattering angle. For Drell-Yan production (excluding the leptons) both the toward and transverse regions are very sensitive to the underlying event. In high-pT jet production the transverse region is very sensitive to the underlying event and is separated into a MAX and MIN transverse region, which helps separate the hard component (initial and final-state radiation) from the beam-beam remnant and multiple parton interaction components of the scattering. The data are corrected to the particle level to remove detector effects and are then compared with several QCD Monte-Carlo models. The goal of this analysis is to provide data that can be used to test and improve the QCD Monte-Carlo models of the underlying event that are used to simulate hadron-hadron collisions.
△ Less
Submitted 16 March, 2010;
originally announced March 2010.
-
Linewidth of single photon transitions in Mn$_{12}$-acetate
Authors:
Beth Parks,
Joseph Loomis,
Evan Rumberger,
David N. Hendrickson,
George Christou
Abstract:
We use time-domain terahertz spectroscopy to measure the position and linewidth of single photon transitions in Mn$_{12}$-acetate. This linewidth is compared to the linewidth measured in tunneling experiments. We conclude that local magnetic fields (due to dipole or hyperfine interactions) cannot be responsible for the observed linewidth, and suggest that the linewidth is due to variations in th…
▽ More
We use time-domain terahertz spectroscopy to measure the position and linewidth of single photon transitions in Mn$_{12}$-acetate. This linewidth is compared to the linewidth measured in tunneling experiments. We conclude that local magnetic fields (due to dipole or hyperfine interactions) cannot be responsible for the observed linewidth, and suggest that the linewidth is due to variations in the anisotropy constants for different clusters. We also calculate a lower limit on the dipole field distribution that would be expected due to random orientations of clusters and find that collective effects must narrow this distribution in tunneling measurements.
△ Less
Submitted 31 October, 2001;
originally announced October 2001.