Search

The STEP (Stability, Transport, Equilibrium, and Pedestal) integrated-modeling tool has been developed in OMFIT to predict stable, tokamak equilibria self-consistently with core-transport and pedestal calculations. STEP couples theory-based codes to integrate a variety of physics, including MHD stability, transport, equilibrium, pedestal formation, and current-drive, heating, and fueling. The input/output of each code is interfaced with a centralized ITER-IMAS data structure, allowing codes to be run in any order and enabling open-loop, feedback, and optimization workflows. This paradigm simplifies the integration of new codes, making STEP highly extensible. STEP has been verified against a published benchmark of six different integrated models. Core-pedestal calculations with STEP have been successfully validated against individual DIII-D H-mode discharges and across more than 500 discharges of the $H_{98,y2}$ database, with a mean error in confinement time from experiment less than 19%. STEP has also reproduced results in less conventional DIII-D scenarios, including negative-central-shear and negative-triangularity plasmas. Predictive STEP modeling has been used to assess performance in several tokamak reactors. Simulations of a high-field, large-aspect-ratio reactor show significantly lower fusion power than predicted by a zero-dimensional study, demonstrating the limitations of scaling-law extrapolations. STEP predictions have found promising EXCITE scenarios, including a high-pressure, 80%-bootstrap-fraction plasma. ITER modeling with STEP has shown that pellet fueling enhances fusion gain in both the baseline and advanced-inductive scenarios. Finally, STEP predictions for the SPARC baseline scenario are in good agreement with published results from the physics basis., Comment: 15 pages, 11 figures Associated with invited talk at 63nd Annual Meeting of the APS Division of Plasma Physics: https://meetings.aps.org/Meeting/DPP21/Session/NI02.1 . The following article has been submitted to Physics of Plasmas. After it is published, it will be found at https://publishing.aip.org/resources/librarians/products/journals/

This paper develops two non-inductive steady state scenarios for larger size configuration of China Fusion Engineering Test Reactor (CFETR) with integrated modeling simulations. A self-consistent core-pedestal coupled workflow for CFETR is developed under integrated modeling framework OMFIT, which allows more accurate evaluation of CFETR performance. The workflow integrates equilibrium code EFIT, transport codes ONETWO and TGYRO, and pedestal code EPED. A fully non-inductive baseline phase I scenario is developed with the workflow, which satisfies the minimum goal of Fusion Nuclear Science Facility. Compared with previous work, which proves the larger size and higher toroidal field CFETR configuration than has the advantages of reducing heating and current drive requirements, lowering divertor and wall power loads, allowing higher bootstrap current fraction and better confinement. A fully non-inductive high-performance phase II scenario is developed, which explores the alpha-particle dominated self-heating regime. Phase II scenario achieves the target of fusion power Pfus>1GW and fusion gain Qfus>20, and it largely reduces auxiliary heating and current drive power. Moreover, the large neutron production of phase II increases the energy generation power and tritium breeding rate.

Extended-MHD modeling of DIII-D tokamak [J. L. Luxon, Nucl. Fusion 42, 614 (2002)] quiescent H-mode (QH-mode) discharges with nonlinear NIMROD [C. R. Sovinec et al., J. Comput. Phys. 195, 355 (2004)] simulations saturates into a turbulent state but does not saturate when the steady-state flow inferred from measurements is not included. This is consistent with the experimental observations of the quiescent regime on DIII-D. The simulation with flow develops into a saturated turbulent state where the n=1 and 2 toroidal modes become dominant through an inverse cascade. Each mode in the range of n=1-5 is dominant at a different time. Consistent with experimental observations during QH-mode, the simulated state leads to large particle transport relative to the thermal transport. Analysis shows that the amplitude and phase of the density and temperature perturbations differ resulting in greater fluctuation-induced convective particle transport relative to the convective thermal transport. Comparison to magnetic-coil measurements shows that rotation frequencies differ between the simulation and experiment, which indicates that more sophisticated extended-MHD two-fluid modeling is required., Comment: 11 pages, 10 figures, manuscript associated with invited talk at APS-DPP Annual Meeting, 2016, in San Jose

The extended-MHD NIMROD code [C.R. Sovinec and J.R. King, J. Comput. Phys. 229, 5803 (2010)] models broadband-MHD activity from a reconstruction of a quiescent H-mode shot on the DIII-D tokamak [J. L. Luxon, Nucl. Fusion 42, 614 (2002)]. Computations with the reconstructed toroidal and poloidal ion flows exhibit low-n perturbations (n=1-5) that grow and saturate into a turbulent-like MHD state. The workflow used to project the reconstructed state onto the NIMROD basis functions re-solves the Grad-Shafranov equation and extrapolates profiles to include scrape-off-layer currents. Evaluation of the transport from the turbulent-like MHD state leads to a relaxation of the density and temperature profiles., Comment: 8 pages

The extended-MHD NIMROD code [C.R. Sovinec and J.R. King, J. Comput. Phys. 229, 5803 (2010)] is verified against the ideal-MHD ELITE code [H.R. Wilson et al. Phys. Plasmas 9, 1277 (2002)] on a diverted tokamak discharge. When the NIMROD model complexity is increased incrementally, resistive and first-order finite-Larmour radius effects are destabilizing and stabilizing, respectively. The full result is compared to local analytic calculations which are found to overpredict both the resistive destabilization and drift stabilization in comparison to the NIMROD computations., Comment: 9 pages

The transition quadrupole moments, $Q_{t}$, of rotational bands in the neutron-rich, even-mass $^{102-108}$Mo and $^{108-112}$Ru nuclei were measured in the 8 to 16 $\hbar $ spin range with the Doppler-shift attenuation method. The nuclei were populated as fission fragments from $^{252}$Cf fission. The detector setup consisted of the Gammasphere spectrometer and the HERCULES fast-plastic array. At moderate spin, the $Q_{t}$ moments are found to be reduced with respect to the values near the ground states. Attempts to describe the observations in mean-field-based models, specifically cranked relativistic Hartree-Bogoliubov theory, illustrate the challenge theory faces and the difficulty to infer information on $\gamma $ softness and triaxiality from the data., Comment: 17 pages, 4 figures

Global electromagnetic gyrokinetic simulations show the existence of near threshold conditions for both a high-$n$ kinetic ballooning mode (KBM) and an intermediate-$n$ kinetic version of peeling-ballooning mode (KPBM) in the edge pedestal of two DIII-D H-mode discharges. When the magnetic shear is reduced in a narrow region of steep pressure gradient, the KPBM is significantly stabilized, while the KBM is weakly destabilized and hence becomes the most-unstable mode. Collisions decrease the KBM's critical $\beta$ and increase the growth rate., Comment: 4 pages, 6 figures

The BOUT++ code is used to simulate ELMs in a shifted circle equilibrium. Reduced ideal MHD simulations are first benchmarked against the linear ideal MHD code ELITE, showing good agreement. Diamagnetic drift effects are included finding the expected suppression of high toroidal mode number modes. Nonlinear simulations are performed, making the assumption that the anomalous kinematic electron viscosity is comparable to the anomalous electron thermal diffusivity. This allows simulations with realistically high Lundquist numbers S = 1e8, finding ELM sizes of 5-10% of the pedestal stored thermal energy. Scans show a strong dependence of the ELM size resistivity at low Lundquist numbers, with higher resistivity leading to more violent eruptions. At high Lundquist numbers relevant to high-performance discharges, ELM size is independent of resistivity as hyper-resistivity becomes the dominant dissipative effect., Comment: 14 pages, 14 figures. Theory of Fusion Plasmas workshop, Varenna 2010. Submitted to Plas. Phys. Control. Fusion

A new modular code called BOUT++ is presented, which simulates 3D fluid equations in curvilinear coordinates. Although aimed at simulating Edge Localised Modes (ELMs) in tokamak X-point geometry, the code is able to simulate a wide range of fluid models (magnetised and unmagnetised) involving an arbitrary number of scalar and vector fields, in a wide range of geometries. Time evolution is fully implicit, and 3rd-order WENO schemes are implemented. Benchmarks are presented for linear and non-linear problems (the Orszag-Tang vortex) showing good agreement. Performance of the code is tested by scaling with problem size and processor number, showing efficient scaling to thousands of processors. Linear initial-value simulations of ELMs using reduced ideal MHD are presented, and the results compared to the ELITE linear MHD eigenvalue code. The resulting mode-structures and growth-rate are found to be in good agreement (BOUT++ = 0.245, ELITE = 0.239). To our knowledge, this is the first time dissipationless, initial-value simulations of ELMs have been successfully demonstrated., Comment: Submitted to Computer Physics Communications. Revised to reduce page count. 18 pages, 16 figures

Recent advances in the integrated modeling of ELMy H-mode plasmas are presented. A model for the H-mode pedestal and for the triggering of ELMs predicts the height, width, and shape of the H-mode pedestal and the frequency and width of ELMs. Formation of the pedestal and the L-H transition is the direct result of ExB flow shear suppression of anomalous transport. The periodic ELM crashes are triggered by either the ballooning or peeling MHD instabilities. The BALOO, DCON, and ELITE ideal MHD stability codes are used to derive a new parametric expression for the peeling-ballooning threshold. The new dependence for the peeling-ballooning threshold is implemented in the ASTRA transport code. Results of integrated modeling of DIII-D like discharges are presented and compared with experimental observations. The results from the ideal MHD stability codes are compared with results from the resistive MHD stability code NIMROD., Comment: 12th International Congress on Plasma Physics, 25-29 October 2004, Nice (France)

A study investigated the effects of the newly introduced Surgeon General's alcohol warnings and advertisements on college students. One hundred fifty-nine undergraduates in communication sciences at the University of Connecticut viewed slides of alcohol products, with or without advertisements and warnings. Following the viewings, subjects filled out questionnaires on product evaluation and on alcohol consumption and risk. Results indicated that the alcohol warnings boomeranged, causing drinkers to perceive greater benefits from the alcoholic beverages. Results also indicated that the advertisements had powerful effects on both drinkers and non-drinkers, increasing perceptions of benefits and decreasing perceptions of risks. Results suggest that, for young people, the warning labels in their present form may be counterproductive, and that advertisements make drinking more attractive. (Two tables of data are included; 46 references are attached.) (Author/PRA)