Source author record

D. Del-Castillo-Negrete

D. Del-Castillo-Negrete appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2022arXiv

Off-axis runaway-electron seed formation, growth and suppression

Novel x-ray detection technology enabled the first profile measurements of the birth and growth dynamics of runaway electrons (REs) at the edge of tokamaks during quiescent RE studies at the Madison Symmetric Torus. The formation of an off-axis RE seed with linear growth rates has been resolved for low energies, a hollow streaming parameter and large electric fields ($E_{\parallel}/E_{D}$) in agreement with theory and simulations. Secondary exponential growth rates have also been spatially resolved for the first time and are consistent with a convective transport of the order of the Ware pinch and energies up to $10^3\times T_{e,0}$. Numerical simulations are shown to reproduce the experimental observations including the off-axis runaway electron generation, radial transport and exponential growth at the core, as well as suppression due to $m=3$ resonant magnetic perturbations.

preprint2020arXiv

Spatially-dependent modeling and simulation of runaway electron mitigation in DIII-D

New simulations with the Kinetic Orbit Runaway electron (RE) Code KORC show RE deconfinement losses to the wall during plasma scrape off are the primary current dissipation mechanism in DIII-D experiments with high-Z impurity injection, and not collisional slowing down. The majority of simulations also exhibit an increase in the RE beam energy due to acceleration by the induced toroidal electric field, even while the RE beam current is decreasing. In this study, KORC integrates RE orbits using the relativistic guiding center equations of motion, and incorporates time-sequenced, experimental reconstructions of the magnetic and electric fields and line integrated electron density to construct spatiotemporal models of electron and partially-ionized impurity transport in the companion plasma. Comparisons of experimental current evolution and KORC results demonstrate the importance of including Coulomb collisions with partially-ionized impurity physics, initial RE energy, pitch angle, and spatial distributions, and spatiotemporal electron and partially-ionized impurity transport. This research provides an initial quantification of the efficacy of RE mitigation via injected impurities, and identification of the critical role played by loss of confinement due to plasma scrape off on the inner wall as compared to the relatively slow collisional damping.

preprint2016arXiv

A statistical study of gyro-averaging effects in a reduced model of drift-wave transport

A statistical study of finite Larmor radius (FLR) effects on transport driven by electrostatic drift-waves is presented. The study is based on a reduced discrete Hamiltonian dynamical system known as the gyro-averaged standard map (GSM). In this system, FLR effects are incorporated through the gyro-averaging of a simplified weak-turbulence model of electrostatic fluctuations. Formally, the GSM is a modified version of the standard map in which the perturbation amplitude, $K_0$, becomes $K_0 J_0(\hatρ)$, where $J_0$ is the zeroth-order Bessel function and $\hatρ$ is the Larmor radius. Assuming a Maxwellian probability density function (pdf) for $\hatρ$, we compute analytically and numerically the pdf and the cumulative distribution function of the effective drift-wave perturbation amplitude $K_0 J_0(\hatρ)$. Using these results we compute the probability of loss of confinement (i.e., global chaos), $P_{c}$, and the probability of trapping in the main drift-wave resonance, $P_{t}$. It is shown that $P_{c}$ provides an upper bound for the escape rate, and that $P_{t}$ provides a good estimate of the particle trapping rate. The analytical results are compared with direct numerical Monte-Carlo simulations of particle transport.

preprint2016arXiv

Modulated heat pulse propagation and partial transport barriers in chaotic magnetic fields

Direct numerical simulations of the time dependent parallel heat transport equation modeling heat pulses driven by power modulation in 3-dimensional chaotic magnetic fields are presented. The numerical method is based on the Fourier formulation of a Lagrangian-Green's function method that provides an accurate and efficient technique for the solution of the parallel heat transport equation in the presence of harmonic power modulation. The numerical results presented provide conclusive evidence that even in the absence of magnetic flux surfaces, chaotic magnetic field configurations with intermediate levels of stochasticity exhibit transport barriers to modulated heat pulse propagation. In particular, high-order islands and remnants of destroyed flux surfaces (Cantori) act as partial barriers that slow down or even stop the propagation of heat waves at places where the magnetic field connection length exhibits a strong gradient. Results on modulated heat pulse propagation in fully stochastic fields and across magnetic islands are also presented. In qualitative agreement with recent experiments in LHD and DIII-D, it is shown that the elliptic (O) and hyperbolic (X) points of magnetic islands have a direct impact on the spatio-temporal dependence of the amplitude of modulated heat pulses.

preprint2016arXiv

Self-consistent chaotic transport in a high-dimensional mean-field Hamiltonian map model

Self-consistent chaotic transport is studied in a Hamiltonian mean-field model. The model provides a simplified description of transport in marginally stable systems including vorticity mixing in strong shear flows and electron dynamics in plasmas. Self-consistency is incorporated through a mean-field that couples all the degrees-of-freedom. The model is formulated as a large set of $N$ coupled standard-like area-preserving twist maps in which the amplitude and phase of the perturbation, rather than being constant like in the standard map, are dynamical variables. Of particular interest is the study of the impact of periodic orbits on the chaotic transport and coherent structures. Numerical simulations show that self-consistency leads to the formation of a coherent macro-particle trapped around the elliptic fixed point of the system that appears together with an asymptotic periodic behavior of the mean field. To model this asymptotic state, we introduced a non-autonomous map that allows a detailed study of the onset of global transport. A turnstile-type transport mechanism that allows transport across instantaneous KAM invariant circles in non-autonomous systems is discussed. As a first step to understand transport, we study a special type of orbits referred to as sequential periodic orbits. Using symmetry properties we show that, through replication, high-dimensional sequential periodic orbits can be generated starting from low-dimensional periodic orbits. We show that sequential periodic orbits in the self-consistent map can be continued from trivial (uncoupled) periodic orbits of standard-like maps using numerical and asymptotic methods. Normal forms are used to describe these orbits and to find the values of the map parameters that guarantee their existence. Numerical simulations are used to verify the prediction from the asymptotic methods.

preprint2014arXiv

Area-preserving maps models of gyro-averaged ${\bf E} \times {\bf B}$ chaotic transport

Discrete maps have been extensively used to model 2-dimensional chaotic transport in plasmas and fluids. Here we focus on area-preserving maps describing finite Larmor radius (FLR) effects on ${\bf E} \times {\bf B}$ chaotic transport in magnetized plasmas with zonal flows perturbed by electrostatic drift waves. FLR effects are included by gyro-averaging the Hamiltonians of the maps which, depending on the zonal flow profile, can have monotonic or non-monotonic frequencies. In the limit of zero Larmor radius, the monotonic frequency map reduces to the standard Chirikov-Taylor map, and, in the case of non-monotonic frequency, the map reduces to the standard nontwist map. We show that in both cases FLR leads to chaos suppression, changes in the stability of fixed points, and robustness of transport barriers. FLR effects are also responsible for changes in the phase space topology and zonal flow bifurcations. Dynamical systems methods based on recurrence time statistics are used to quantify the dependence on the Larmor radius of the threshold for the destruction of transport barriers.

preprint2014arXiv

Front propagation in reaction-diffusion systems with anomalous diffusion

A numerical study of the role of anomalous diffusion in front propagation in reaction-diffusion systems is presented. Three models of anomalous diffusion are considered: fractional diffusion, tempered fractional diffusion, and a model that combines fractional diffusion and regular diffusion. The reaction kinetics corresponds to a Fisher-Kolmogorov nonlinearity. The numerical method is based on a finite-difference operator splitting algorithm with an explicit Euler step for the time advance of the reaction kinetics, and a Crank-Nicholson semi-implicit time step for the transport operator. The anomalous diffusion operators are discretized using an upwind, flux-conserving, Grunwald-Letnikov finite-difference scheme applied to the regularized fractional derivatives. With fractional diffusion of order $α$, fronts exhibit exponential acceleration, $a_L(t) \sim e^{γt/α}$, and develop algebraic decaying tails, $ϕ\sim 1/x^α$. In the case of tempered fractional diffusion, this phenomenology prevails in the intermediate asymptotic regime $\left(χt \right)^{1/α} \ll x \ll 1/λ$, where $1/λ$ is the scale of the tempering. Outside this regime, i.e. for $x > 1/λ$, the tail exhibits the tempered decay $ϕ\sim e^{-λx}/x^{α+1}$, and the front velocity approaches the terminal speed $v_*= \left(γ-λ^αχ\right)/ λ$. Of particular interest is the study of the interplay of regular and fractional diffusion. It is shown that the main role of regular diffusion is to delay the onset of front acceleration. In particular, the crossover time, $t_c$, to transition to the accelerated fractional regime exhibits a logarithmic scaling of the form $t_c \sim \log \left(χ_d/χ_f\right)$ where $χ_d$ and $χ_f$ are the regular and fractional diffusivities.

preprint2014arXiv

Heat pulse propagation in chaotic 3-dimensional magnetic fields

Heat pulse propagation in $3$-D chaotic magnetic fields is studied by solving the parallel heat transport equation using a Lagrangian-Green's function (LG) method. The LG method provides an efficient and accurate technique that circumvents limitations of finite elements and finite difference methods. The main two problems addressed are: (i) The dependence of the radial transport on the magnetic field stochasticity (controlled by the amplitude of the perturbation, $ε$); and (ii) The role of reversed shear configurations on pulse propagation. In all the cases considered there are no magnetic flux surfaces. However, radial transport is observed to depend strongly on $ε$ due to the presence of high-order magnetic islands and Cantori that act as quasi-transport barriers that preclude the radial penetration of heat pulses within physically relevant time scale. The dependence of the magnetic field connection length, $\ell_B$, on $ε$ is studied in detail. The decay rate of the temperature maximum, $\langle T \rangle_{max}(t)$, the time delay of the temperature response as function of the radius, $τ$, and the radial heat flux $\langle {{\bf q}\cdot {\hat e}_ψ} \rangle$, are also studied as functions of the magnetic field stochasticity and $\ell_B$. In all cases, the scaling of $\langle T \rangle_{max}$ with $t$ transitions from sub-diffusive, $\langle T \rangle_{max} \sim t^{-1/4}$, at short times ($χ_\parallel t< 10^5$) to a significantly slower scaling at longer times ($χ_\parallel t > 10^5$). A strong dependence on $ε$ is also observed on $τ$ and $\langle {{\bf q}\cdot {\hat e}_ψ} \rangle$. The radial propagation of pulses in fully chaotic fields considerably slows down in the shear reversal region and, as a result, $τ$, in reversed shear configurations is an order of magnitude longer than the one in monotonic $q$-profiles.

preprint2010arXiv

Lagrangian statistics and flow topology in forced two-dimensional turbulence

A study of the relationship between Lagrangian statistics and flow topology in fluid turbulence is presented. The topology is characterized using the Weiss criterion that provides a simplified tool to partition the flow into topologically different regions: elliptic (vortex dominated), hyperbolic (deformation dominated), and intermediate (turbulent background). The flow corresponds to forced two-dimensional Navier-Stokes turbulence in doubly periodic and circular bounded domains with non-slip boundary conditions. In the double periodic domain, the probability density function (pdf) of the Weiss field exhibits a negative skewness consistent with the fact that in periodic domains the flow is dominated by coherent vortex structures. On the other hand, in the circular domain, the elliptic and hyperbolic regions seem to be statistically similar. We follow a Lagrangian approach and obtain the statistics by tracking large ensembles of passively advected tracers. The pdfs of residence time in the topologically different regions are computed using the Lagrangian Weiss field, i.e., the Weiss field computed along the particles' trajectories. In elliptic and hyperbolic regions, the pdfs of the residence time have self-similar algebraic decaying tails. On the other hand, in the intermediate regions the pdf has exponential decaying tails. The conditional (with respect to the flow topology) pdfs of the Lagrangian velocity exhibit Gaussian behavior in the periodic and in the bounded domains. In contrast to the freely decaying turbulence case, the conditional pdfs of the Lagrangian acceleration in forced turbulence show a comparable level of intermittency in the periodic and the bounded domains. The conditional pdfs of the Lagrangian curvature are characterized, in all cases, by self-similar power law behavior with a decay exponent of order -2.

preprint2010arXiv

Levy ratchets in the spatially tempered fractional Fokker-Planck equation

Lévy ratchets are minimal models of fluctuation-driven transport in the presence of Lévy noise and periodic external potentials with broken spatial symmetry. In these systems, a net ratchet current can appear even in the absence of time dependent perturbations, external tilting forces, or a bias in the noise. The majority of studies on the interaction of Lévy noise with external potentials have assumed $α$-stable Lévy statistics in the Langevin description, which in the continuum limit corresponds to the fractional Fokker-Planck equation. However, the divergence of the low order moments is a potential drawback of $α$-stable distributions because, in applications, the moments represent physical quantities. For example, for $α<1$, the current $J$, in $α$-stable Lévy ratchets is unbounded. To overcome this limitation, we study ratchet transport using truncated Lévy distributions which in the continuum limit correspond to the spatially tempered fractional Fokker-Planck equation. The main object of study is the dependence of the ratchet current on the level of tempering, $λ$. For $λ\neq 0$, the statistics ultimately converges (although very slowly) to Gaussian diffusion in the absence of a potential. However, it is shown here that in the presence of a ratchet potential a finite current persists asymptotically for any finite value of $λ$. The current converges exponentially in time to the steady state value. The steady state current exhibits algebraically decay, $J\sim λ^{-ζ}$, for $α\geq 1.75$. However, for $α\leq 1.5$, the decay is exponential, $J \sim e^{-ξλ}$. In the presence of a bias in the Lévy noise, it is shown that the tempering can lead to a current reversal. A detailed numerical study is presented on the dependence of the current on $λ$ and the physical parameters of the system.

preprint2010arXiv

Local and nonlocal parallel heat transport in general magnetic fields

A novel approach that enables the study of parallel transport in magnetized plasmas is presented. The method applies to general magnetic fields with local or nonlocal parallel closures. Temperature flattening in magnetic islands is accurately computed. For a wave number $k$, the fattening time scales as $χ_{\parallel} τ\sim k^{-α}$ where $χ$ is the parallel diffusivity, and $α=1$ ($α=2$) for non-local (local) transport. The fractal structure of the devil staircase temperature radial profile in weakly chaotic fields is resolved. In fully chaotic fields, the temperature exhibits self-similar evolution of the form $T=(χ_{\parallel} t)^{-γ/2} L \left[ (χ_{\parallel} t)^{-γ/2} δψ\right]$, where $δψ$ is a radial coordinate. In the local case, $f$ is Gaussian and the scaling is sub-diffusive, $γ=1/2$. In the non-local case, $f$ decays algebraically, $L (η) \sim η^{-3}$, and the scaling is diffusive, $γ=1$.

preprint2010arXiv

Non-diffusive, non-local transport in fluids and plasmas

A review of non-diffusive transport in fluids and plasmas is presented. In the fluid context, non-diffusive chaotic transport by Rossby waves in zonal flows is studied following a Lagrangian approach. In the plasma physics context the problem of interest is test particle transport in pressure-gradient-driven plasma turbulence. In both systems the probability density function (PDF) of particle displacements is strongly non-Gaussian and the statistical moments exhibit super-diffusive anomalous scaling. Fractional diffusion models are proposed and tested in the quantitative description of the non-diffusive Lagrangian statistics of the fluid and plasma problems. Also, fractional diffusion operators are used to construct non-local transport models exhibiting up-hill transport, multivalued flux-gradient relations, fast pulse propagation phenomena, and "tunneling" of perturbations across transport barriers.

preprint2009arXiv

Universal Probability Distribution Function for Bursty Transport in Plasma Turbulence

Bursty transport phenomena associated with convective motion present universal statistical characteristics among different physical systems. In this letter, a stochastic univariate model and the associated probability distribution function for the description of bursty transport in plasma turbulence is presented. The proposed stochastic process recovers the universal distribution of density fluctuations observed in plasma edge of several magnetic confinement devices and the remarkable scaling between their skewness $S$ and kurtosis $K$. Similar statistical characteristics of variabilities have been also observed in other physical systems that are characterized by convection such as the X-ray fluctuations emitted by the Cygnus X-1 accretion disc plasmas and the sea surface temperature fluctuations.