3D simulations of vertical displacement events in tokamaks: A benchmark of M3D-C1^1, NIMROD and JOREK

In recent years, the nonlinear 3D magnetohydrodynamic codes JOREK, M3D-C$^1$ and NIMROD developed the capability of modelling realistic 3D vertical displacement events (VDEs) including resistive walls. In this paper, a comprehensive 3D VDE benchmark is presented between these state of the art codes. The simulated case is based on an experimental NSTX plasma but with a simplified rectangular wall. In spite of pronounced differences between physics models and numerical methods, the comparison shows very good agreement in the relevant quantities used to characterize disruptions such as the 3D wall forces and energy decay. This benchmark does not only bring confidence regarding the use of the mentioned codes for disruption studies, but also shows differences with respect to the used models (e.g. reduced versus full MHD models). The simulations show important 3D features for a NSTX plasma such as the self-consistent evolution of the halo current and the origin of the wall forces. In contrast to other reduced MHD models based on an ordering in the aspect ratio, the ansatz based JOREK reduced MHD model allows capturing the 3D dynamics even in the spherical tokamak limit considered here

Axisymmetric simulations of vertical displacement events in tokamaks: A benchmark of M3D-C1, NIMROD and JOREK

A benchmark exercise for the modeling of vertical displacement events(VDEs) is presented and applied to the 3D nonlinear magneto-hydrodynamic codesM3D-C1, JOREK and NIMROD. The simulations are based on a vertically unstableNSTX equilibrium enclosed by an axisymmetric resistive wall with rectangular crosssection. A linear dependence of the linear VDE growth rates on the resistivity ofthe wall is recovered for sufficiently large wall conductivity and small temperatures inthe open field line region. The benchmark results show good agreement between theVDE growth rates obtained from linear NIMROD and M3D-C1simulations as wellas from the linear phase of axisymmetric nonlinear JOREK, NIMROD and M3D-C1simulations. Axisymmetric nonlinear simulations of a full VDE performed with thethree codes are compared and excellent agreement is found regarding plasma locationand plasma currents as well as eddy and halo currents in the wall.</p

Variación de 13C en el ambiente de Madrid (atmósfera, aguas superficiales y aguas subterráneas)

Stable isotopes of carbon are studied in atmosphere, freshwater (rivers) and groundwater of the region of Madrid. Sampling and analytical techniques are adapted for the isotopic study of carbon in the environment. The atmosphere samples show δ13C values from -13.92 to -15.67 ‰, closely related to the polution level of air. The water from wells, to the Tertiary detritic aquifer, with δ13C among -13.92 and -15.67 ‰. The dissolution of carbonatic rocks and the biologic activity could be the origin of the wide variations δ13C from -7.07 to -11.85 %0) found in the waters from rivers Jarama, Guadalix and Guadarrama.Se analizan el 13C de la atmósfera, ríos y pozos del área de Madrid, al objeto de evaluar su distribución y origen, mediante adaptación y desarrollo de técnicas de muestreo y análisis. Las muestras atmosféricas presentan valores de δ13C situados entre -9,81 y -14,64 ‰, que se corresponden gradualmente con los puntos de menor a mayor grado de polución por aportes principalmente de C02 procedente de la quema de combustibles fósiles. Las aguas de los pozos analizados, pertenecientes al acuífero terciario del detrítico, muestran un estrecho margen de variación, entre δ13C -13,92 y -15,67 ‰, relacionado con los procesos geoquímicos naturales. Las aguas de los ríos Jarama, Guadalix y Guadarrama, presentan una amplia variación en su composición isotópica (δ13C entre -7,07 y -11,85 ‰) debido a la disolución de carbonatos a su paso por las formaciones de calizas del Mesozoico y, en menor medida, a los aportes de efluentes contaminados

Simulations of COMPASS vertical displacement events with a self-consistent model for halo currents including neutrals and sheath boundary conditions

The understanding of the halo current properties during disruptions is key to design and operate large scale tokamaks in view of the large thermal and electromagnetic loads that they entail. For the first time, we present a fully self-consistent model for halo current simulations including neutral particles and sheath boundary conditions. The model is used to simulate vertical displacement events (VDEs) occurring in the COMPASS tokamak. Recent COMPASS experiments have shown that the parallel halo current density at the plasma-wall interface is limited by the ion saturation current during VDE-induced disruptions. We show that usual magneto-hydrodynamic boundary conditions can lead to the violation of this physical limit and we implement this current density limitation through a boundary condition for the electrostatic potential. Sheath boundary conditions for the density, the heat flux, the parallel velocity and a realistic parameter choice (e.g. Spitzer's resistivity and Spitzer-Harm parallel thermal conductivity) extend present VDE simulations beyond the state of the art. Experimental measurements of the current density, temperature and heat flux profiles at the COMPASS divertor are compared with the results obtained from axisymmetric simulations. Since the ion saturation current density (Jsat) is shown to be essential to determine the halo current profile, parametric scans are performed to study its dependence on different quantities such as the plasma resistivity and the particle and heat diffusion coefficients. In this respect, the plasma resistivity in the halo region broadens significantly the Jsat profile, increasing the halo width at a similar total halo current

Functionalized cyclophellitols are selective glucocerebrosidase inhibitors and induce a bona fide neuropathic Gaucher model in zebrafish

Gaucher disease is caused by inherited deficiency in glucocerebrosidase (GBA, a retaining β-glucosidase), and deficiency in GBA constitutes the largest known genetic risk factor for Parkinson's disease. In the past, animal models of Gaucher disease have been generated by treatment with the mechanism-based GBA inhibitors, conduritol B epoxide (CBE), and cyclophellitol. Both compounds, however, also target other retaining glycosidases, rendering generation and interpretation of such chemical knockout models complicated. Here we demonstrate that cyclophellitol derivatives carrying a bulky hydrophobic substituent at C8 are potent and selective GBA inhibitors and that an unambiguous Gaucher animal model can be readily generated by treatment of zebrafish with these

How Gibbs distributions may naturally arise from synaptic adaptation mechanisms. A model-based argumentation

This paper addresses two questions in the context of neuronal networks dynamics, using methods from dynamical systems theory and statistical physics: (i) How to characterize the statistical properties of sequences of action potentials ("spike trains") produced by neuronal networks ? and; (ii) what are the effects of synaptic plasticity on these statistics ? We introduce a framework in which spike trains are associated to a coding of membrane potential trajectories, and actually, constitute a symbolic coding in important explicit examples (the so-called gIF models). On this basis, we use the thermodynamic formalism from ergodic theory to show how Gibbs distributions are natural probability measures to describe the statistics of spike trains, given the empirical averages of prescribed quantities. As a second result, we show that Gibbs distributions naturally arise when considering "slow" synaptic plasticity rules where the characteristic time for synapse adaptation is quite longer than the characteristic time for neurons dynamics.Comment: 39 pages, 3 figure