Plasma Edge Kinetic-MHD Modeling in Tokamaks Using Kepler Workflow for Code Coupling, Data Management and Visualization

A new predictive computer simulation tool targeting the development of the H-mode pedestal at the plasma edge in tokamaks and the triggering and dynamics of edge localized modes (ELMs) is presented in this report. This tool brings together, in a coordinated and effective manner, several first-principles physics simulation codes, stability analysis packages, and data processing and visualization tools. A Kepler workflow is used in order to carry out an edge plasma simulation that loosely couples the kinetic code, XGC0, with an ideal MHD linear stability analysis code, ELITE, and an extended MHD initial value code such as M3D or NIMROD. XGC0 includes the neoclassical ion-electron-neutral dynamics needed to simulate pedestal growth near the separatrix. The Kepler workflow processes the XGC0 simulation results into simple images that can be selected and displayed via the Dashboard, a monitoring tool implemented in AJAX allowing the scientist to track computational resources, examine running and archived jobs, and view key physics data, all within a standard Web browser. The XGC0 simulation is monitored for the conditions needed to trigger an ELM crash by periodically assessing the edge plasma pressure and current density profiles using the ELITE code. If an ELM crash is triggered, the Kepler workflow launches the M3D code on a moderate-size Opteron cluster to simulate the nonlinear ELM crash and to compute the relaxation of plasma profiles after the crash. This process is monitored through periodic outputs of plasma fluid quantities that are automatically visualized with AVS/Express and may be displayed on the Dashboard. Finally, the Kepler workflow archives all data outputs and processed images using HPSS, as well as provenance information about the software and hardware used to create the simulation. The complete process of preparing, executing and monitoring a coupled-code simulation of the edge pressure pedestal buildup and the ELM cycle using the Kepler scientific workflow system is described in this paper

Time dependence of triplet-singlet excitation transfer from compact poly rA to bound dye at 77 K

The nonexponential phosphorescence decay of a highly folded form of poly-riboadenylic acid (poly rA) with noncovalently bound dye is explained by a novel application of a well-known theory of electronic excitation transfer based on the Förster mechanism. This theory, originally used to describe singlet-singlet energy transfer from donor molecules to an acceptor in a solution, is here applied to the transfer of triplet excitation from the adenine (in poly rA) to the singlet manifold of either of the bound dyes, ethidium bromide or proflavine. New experimental data are presented that allow straight-forward theoretical interpretation. These data fit the form predicted by the theory, U(t) exp (-Bt[superscript 1/2]), where U(t) is the decay of the poly rA phosphorescence in the absence of dye, for a range of relative concentrations of either dye. The self-consistency of these theoretical fits is demonstrated by the proportionality of B to the square root of the Förster triplet-singlet overlap integrals for transfer from poly rA to each of the dyes, as demanded by the theory. From these self-consistent values of B, the theory enables one to deduce the mean packing density of nucleotides in this folded poly rA, which we estimate to be approximately ~ 1 nm⁻³. We conclude that some variations of the method described here may be useful for deducing packing densities of nucleotides in other compact nucleic acid structures.This is the publisher’s final pdf. The article is copyrighted by The Biophysical Society and published by Cell Press/Elsevier. It can be found at: http://www.cell.com/biophysj

Analysis and Modeling of DIII-D Hybrid Discharges and their Extrapolation to ITER

Recent experiments on tokamaks around the world [1-5] have demonstrated discharges with moderately high performance in which the q-profile remains stationary, as measured by the motional Stark effect diagnostic, for periods up to several {tau}{sub R}. Hybrid discharges are characterize by q{sub min} {approx} 1, high {beta}{sub N}, and good confinement. These discharges have been termed hybrid because of their intermediate nature between that of an ordinary H-mode and advanced tokamak discharges. They form an attractive scenario for ITER as the normalized fusion performance ({beta}{sub N}H{sub 89P}/q{sub 95}{sup 2}) is at or above that for the ITER baseline Q{sub fus} = 10 scenario, even for q{sub 95} as high as 4.6. The startup phase is thought to be crucial to the ultimate evolution of the hybrid discharge. An open question is how hybrid discharges achieve and maintain their stationary state during the initial startup phase. To investigate this aspect of hybrid discharges, we have used the CORSICA code to model the early stages of a discharge. Results clearly indicate that neoclassical current evolution alone is insufficient to account for the time evolution of the q-profile and that an addition of non-inductive current source must be incorporated into the model to reproduce the experimental time history. We include non-inductive neutral beam and bootstrap current sources in the model, and investigate the difference between simulations with these sources and the experimentally inferred q-profile. Further, we have made preliminary estimates of the spatial structure of the current needed to bring the simulation and experiment into agreement. This additional non-inductive source has not been tied to any physical mechanism as yet. We present these results and discuss the implications for hybrid startup on ITER

Model dependence of single-energy fits to pion photoproduction data

Model dependence of multipole analysis has been explored through energy-dependent and single-energy fits to pion photoproduction data. The MAID energy-dependent solution has been used as input for an event generator producing realistic pseudo data. These were fitted using the SAID parametrization approach to determine single-energy and energy-dependent solutions over a range of lab photon energies from 200 to 1200 MeV. The resulting solutions were found to be consistent with the input amplitudes from MAID. Fits with a $\chi$-squared per datum of unity or less were generally achieved. We discuss energy regions where consistent results are expected, and explore the sensitivity of fits to the number of included single- and double-polarization observables. The influence of Watson's theorem is examined in detail.Comment: 12 pages, 8 figure

Stable bootstrap-current driven equilibria for low aspect ratio tokamaks

Low aspect ratio tokamaks can potentially provide a high ratio of plasma pressure to magnetic pressure {beta} and high plasma current I at a modest size, ultimately leading to a high power density compact fusion power plant. For the concept to be economically feasible, bootstrap current must be a major component of the plasma current. A high value of the Troyon factor {beta}{sub N} and strong shaping are required to allow simultaneous operation at high {beta} and high bootstrap current fraction. Ideal magnetohydrodynamic stability of a range of equilibria at aspect ratio 1.4 is systematically explored by varying the pressure profile and shape. The pressure and current profiles are constrained in such a way as to assure complete bootstrap current alignment. Both {beta}{sub N} and {beta} are defined in terms of the vacuum toroidal field. Equilibria with {beta}{sub N} {ge} 8 and {beta} - 35% to 55% exist which are stable to n = {infinity} ballooning modes, and stable to n = 0, 1,2,3 kink modes with a conducting wall. The dependence of {beta} and {beta}{sub N} with respect to aspect ratio is also considered

Studies on Neutral Beam Injection into the SSPX Spheromak Plasma

In the Sustained Spheromak Physics Experiment, (SSPX) ['Improved operation of the SSPX spheromak', R.D. Wood, D.N. Hill, E.B. Hooper, S. Woodruff1, H.S. McLean and B.W. Stallard, Nucl. Fusion 45 1582-1588 (2005)], plasmas with core electron temperatures reaching up to 500 eV at densities of 10{sup 20}/m{sup 3} have been sustained for several milliseconds, making them suitable as targets for neutral beam injection. High performance and further progress in understanding Spheromak plasma physics are expected if neutral beams are injected into the plasma. This paper presents the results of numerical 1.5 D modeling of the plasma to calculate neutral beam current drive and ion and electron heating. The results are presented for varying initial conditions of density, temperatures and profiles and beam energy, injection angle and power. Current drive efficiency (Ampere/Watt of absorbed power) of up to 0.08 can be achieved with best performance SSPX shots as target. Analyses of neutral beam heating indicate that ion temperatures of up to 1.5 keV and electron temperatures of up to 750 eV can be obtained with injection of about 1 MW of neutral beam for 5-10 ms and with diffusivities typically observed in SSPX. Injection targeting near the magnetic axis appears to be the best for heating and current drive. Effect of the current drive and evolution of SSPX equilibrium are discussed

Evolution of breeding plumages in birds: A multiple-step pathway to seasonal dichromatism in New World warblers (Aves: Parulidae)

Ecology and Evolution published by John Wiley & Sons Ltd Many species of birds show distinctive seasonal breeding and nonbreeding plumages. A number of hypotheses have been proposed for the evolution of this seasonal dichromatism, specifically related to the idea that birds may experience variable levels of sexual selection relative to natural selection throughout the year. However, these hypotheses have not addressed the selective forces that have shaped molt, the underlying mechanism of plumage change. Here, we examined relationships between life-history variation, the evolution of a seasonal molt, and seasonal plumage dichromatism in the New World warblers (Aves: Parulidae), a family with a remarkable diversity of plumage, molt, and life-history strategies. We used phylogenetic comparative methods and path analysis to understand how and why distinctive breeding and nonbreeding plumages evolve in this family. We found that color change alone poorly explains the evolution of patterns of biannual molt evolution in warblers. Instead, molt evolution is better explained by a combination of other life-history factors, especially migration distance and foraging stratum. We found that the evolution of biannual molt and seasonal dichromatism is decoupled, with a biannual molt appearing earlier on the tree, more dispersed across taxa and body regions, and correlating with separate life-history factors than seasonal dichromatism. This result helps explain the apparent paradox of birds that molt biannually but show breeding plumages that are identical to the nonbreeding plumage. We find support for a two-step process for the evolution of distinctive breeding and nonbreeding plumages: That prealternate molt evolves primarily under selection for feather renewal, with seasonal color change sometimes following later. These results reveal how life-history strategies and a birds\u27 environment act upon multiple and separate feather functions to drive the evolution of feather replacement patterns and bird coloration

Confinement Studies in High Temperature Spheromak Plasmas

Recent results from the SSPX spheromak experiment demonstrate the potential for obtaining good energy confinement (Te > 350eV and radial electron thermal diffusivity comparable to tokamak L-mode values) in a completely self-organized toroidal plasma. A strong decrease in thermal conductivity with temperature is observed and at the highest temperatures, transport is well below that expected from the Rechester-Rosenbluth model. Addition of a new capacitor bank has produced 60% higher magnetic fields and almost tripled the pulse length to 11ms. For plasmas with T{sub e} > 300eV, it becomes feasible to use modest (1.8MW) neutral beam injection (NBI) heating to significantly change the power balance in the core plasma, making it an effective tool for improving transport analysis. We are now developing detailed designs for adding NBI to SSPX and have developed a new module for the CORSICA transport code to compute the correct fast-ion orbits in SSPX so that we can simulate the effect of adding NBI; initial results predict that such heating can raise the electron temperature and total plasma pressure in the core by a factor of two

Listening to the voices of women suffering perinatal psychological distress

This article suggests that transactional analysis can be an effective treatment approach for women suffering from mental health conditions and the emotional and life disturbances that may occur during the perinatal period. It offers a brief introduction to perinatal psychological distress followed by a description of the use of transactional analysis psychotherapy for this condition. The article outlines a new model for a research project that aims to ascertain women’s views on the helpfulness of the treatment and to gain a better understanding of the stigma often associated with perinatal mental health issues. The author argues for the necessity of qualitative research to assess the efficacy of transactional-analysis-based treatment and to increase our knowledge about the change process in transactional analysis psychotherapy with this client population as well as to inform future transactional analysis treatment protocols

High In Situ Repeatability of Behaviour Indicates Animal Personality in the Beadlet Anemone Actinia equina (Cnidaria)

‘Animal personality’ means that individuals differ from one another in either single behaviours or suites of related behaviours in a way that is consistent over time. It is usually assumed that such consistent individual differences in behaviour are driven by variation in how individuals respond to information about their environment, rather than by differences in external factors such as variation in microhabitat. Since behavioural variation is ubiquitous in nature we might expect ‘animal personality’ to be present in diverse taxa, including animals with relatively simple nervous systems. We investigated in situ startle responses in a sea anemone, Actinia equina, to determine whether personalities might be present in this example of an animal with a simple nervous system. We found very high levels of repeatability among individuals that were re-identified in the same locations over a three week sampling period. In a subset of the data, where we used tide-pool temperature measurements to control for a key element of variation in microhabitat, these high levels of repeatability remained. Although a range of other consistent differences in micro-habitat features could have contributed to consistent differences between the behaviour of individuals, these data suggest the presence of animal personality in A. equina. Rather than being restricted to certain groups, personality may be a general feature of animals and may be particularly pronounced in species with simple nervous systems