scispace - formally typeset
Search or ask a question

Showing papers in "arXiv: Plasma Physics in 2014"


Journal ArticleDOI
TL;DR: In this paper, the curvature drift and the parallel electric field dominate the dynamics and drive parallel heating in collisionless magnetic reconnection, and an upper limit on electron energy gain is obtained by balancing reconnection drive with radiative loss.
Abstract: The heating of electrons in collisionless magnetic reconnection is explored in particle-in-cell (PIC) simulations with non-zero guide fields so that electrons remain magnetized. In this regime electric fields parallel to B accelerate particles directly while those perpendicular to B do so through gradient-B and curvature drifts. The curvature drift drives parallel heating through Fermi reflection while the gradient B drift changes the perpendicular energy through betatron acceleration. We present simulations in which we evaluate each of these mechanisms in space and time in order to quantify their role in electron heating. For a case with a small guide field (20 % of the magnitude of the reconnecting component) the curvature drift is the dominant source of electron heating. However, for a larger guide field (equal to the magnitude of the reconnecting component) electron acceleration by the curvature drift is comparable to that of the parallel electric field. In both cases the heating by the gradient B drift is negligible in magnitude. Heating by the curvature-drift dominates in the outflow exhausts where bent field lines expand to relax their tension and is therefore distributed over a large area. In contrast, the parallel electric field is localized near X-lines. Acceleration by parallel electric fields may play a smaller role in large systems where the X-line occupies a vanishing fraction of the system. The curvature drift and the parallel electric field dominate the dynamics and drive parallel heating. A consequence is that the electron energy spectrum becomes extremely anisotropic at late time, which has important implications for quantifying the limits of electron acceleration due to synchrotron emission. An upper limit on electron energy gain that is substantially higher than earlier estimates is obtained by balancing reconnection drive with radiative loss.

131 citations


Journal ArticleDOI
TL;DR: In this article, a multidiagnostic approach, utilizing Langmuir probes in the midplane, X-point and divertor walls, along with Lithium beam and infrared measurements is employed to evaluate the evolution of the scrape-off layer (SOL) of ASDEX Upgrade across the L-mode density transition leading to the formation of a density shoulder.
Abstract: A multidiagnostic approach, utilizing Langmuir probes in the midplane, X-point and divertor walls, along with Lithium beam and infrared measurements is employed to evaluate the evolution of the Scrape-off Layer (SOL) of ASDEX Upgrade across the L-mode density transition leading to the formation of a density shoulder. The flattening of the SOL density profiles is linked to a regime change of filaments, which become faster and larger, and to a similar flattening of the $q_{\parallel}$ profile. This transition is related to the beginning of outer divertor detachment and leads to the onset of a velocity shear layer in the SOL. Experimental measurements are in good agreement with several filament models which describe the process as a transition from conduction to convection-dominated SOL perpendicular transport caused by an increase of parallel collisionality. These results could be of great relevance since both ITER and DEMO will feature detached divertors and densities largely over the transition values, and might therefore exhibit convective transport levels different to those observed typically in present-day devices.

127 citations


Journal ArticleDOI
TL;DR: In this article, a theory-based modelling of core heavy impurity transport is presented, and shown to be necessary for quantitative description of present experiments in JET and ASDEX Upgrade.
Abstract: Recent developments in theory-based modelling of core heavy impurity transport are presented, and shown to be necessary for quantitative description of present experiments in JET and ASDEX Upgrade. The treatment of heavy impurities is complicated by their large mass and charge, which result in a strong response to plasma rotation or any small background electrostatic field in the plasma, such as that generated by anisotropic external heating. These forces lead to strong poloidal asymmetries of impurity density, which have recently been added to numerical tools describing both neoclassical and turbulent transport. Modelling predictions of the steady-state two-dimensional tungsten impurity distribution are compared with experimental densities interpreted from soft X-ray diagnostics. The modelling identifies neoclassical transport enhanced by poloidal asymmetries as the dominant mechanism responsible for tungsten accumulation in the central core of the plasma. Depending on the bulk plasma profiles, neoclassical temperature screening can prevent accumulation, and can be enhanced by externally heated species, demonstrated here in ICRH plasmas.

106 citations


Journal ArticleDOI
TL;DR: In this paper, the effects of including fast neutrals and using realistic energy-dependent secondary electron emission coefficients for ions and neutrals in simulations of CCPs operated in argon at 13.56 MHz and at neutral gas pressures between 3 Pa and 100 Pa.
Abstract: In most PIC/MCC simulations of radio frequency capacitively coupled plasmas (CCPs) several simplifications are made: (i) fast neutrals are not traced, (ii) heavy particle induced excitation and ionization are neglected, (iii) secondary electron emission from boundary surfaces due to neutral particle impact is not taken into account, and (iv) the secondary electron emission coefficient is assumed to be constant, i.e. independent of the incident particle energy and the surface conditions. Here we question the validity of these simplifications under conditions typical for plasma processing applications. We study the effects of including fast neutrals and using realistic energy-dependent secondary electron emission coefficients for ions and fast neutrals in simulations of CCPs operated in argon at 13.56 MHz and at neutral gas pressures between 3 Pa and 100 Pa. We find a strong increase of the plasma density and the ion flux to the electrodes under most conditions, if these processes are included realistically in the simulation. The sheath widths are found to be significantly smaller and the simulation is found to diverge at high pressures for high voltage amplitudes in qualitative agreement with experimental findings. By switching individual processes on and off in the simulation we identify their individual effects on the ionization dynamics and plasma parameters. We conclude that fast neutrals and energy-dependent secondary electron emission coefficients must be included in simulations of CCPs in order to yield realistic results.

90 citations


Journal ArticleDOI
TL;DR: In this paper, nanosecond pulsed discharges in N$2}$ and N$ 2}$/0.9% H$ 2 O at atmospheric pressure (at 300 K) were studied with time-resolved imaging, optical emission spectroscopy and Rayleigh scattering.
Abstract: In this contribution, nanosecond pulsed discharges in N$_{2}$ and N$_{2}$/0.9% H$_{2}$O at atmospheric pressure (at 300 K) are studied with time-resolved imaging, optical emission spectroscopy and Rayleigh scattering. A 170 ns high voltage pulse is applied across two pin-shaped electrodes at a frequency of 1 kHz. The discharge consists of three phases: an ignition phase, a spark phase and a recombination phase. During the ignition phase the emission is mainly caused by molecular nitrogen (N$_{2}$(C-B)). In the spark and recombination phase mainly atomic nitrogen emission is observed. The emission when H$_{2}$O is added is very similar, except the small contribution of H$_{\alpha}$ and the intensity of the molecular N$_{2}$(C-B) emission is less.

87 citations


Journal ArticleDOI
TL;DR: In this article, the authors show that the increase in ELM frequency is correlated with the size of the edge peeling-tearing like response of the plasma and the distortions of the X-point region.
Abstract: Sustained ELM mitigation has been achieved on MAST and AUG using RMPs with a range of toroidal mode numbers over a wide region of low to medium collisionality discharges. The ELM energy loss and peak heat loads at the divertor targets have been reduced. The ELM mitigation phase is typically associated with a drop in plasma density and overall stored energy. In one particular scenario on MAST, by carefully adjusting the fuelling it has been possible to counteract the drop in density and to produce plasmas with mitigated ELMs, reduced peak divertor heat flux and with minimal degradation in pedestal height and confined energy. While the applied resonant magnetic perturbation field can be a good indicator for the onset of ELM mitigation on MAST and AUG there are some cases where this is not the case and which clearly emphasise the need to take into account the plasma response to the applied perturbations. The plasma response calculations show that the increase in ELM frequency is correlated with the size of the edge peeling-tearing like response of the plasma and the distortions of the plasma boundary in the X-point region.

81 citations


Journal ArticleDOI
TL;DR: In this paper, simulations of isolated 3D filaments in a slab geometry obtained using a 3D reduced fluid code were performed to investigate how the dynamics of a filament are affected by its amplitude, perpendicular size and parallel extent.
Abstract: This paper presents simulations of isolated 3D filaments in a slab geometry obtained using a 3D reduced fluid code. First, systematic scans were performed to investigate how the dynamics of a filament are affected by its amplitude, perpendicular size and parallel extent. The perpendicular size of the filament was found to have a strong influence on its motions, as it determined the relative importance of parallel currents to polarisation and viscous currents, whilst drift-wave instabilities were observed if the initial amplitude of the blob was increased sufficiently. Next, the 3D simulations were compared to 2D simulations using different parallel closures; namely, the sheath dissipation closure, which neglects parallel gradients, and the vorticity advection closure, which neglects the influence of parallel currents. The vorticity advection closure was found to not replicate the 3D perpendicular dynamics and overestimated the initial radial acceleration of all the filaments studied. In contrast, a more satisfactory comparison with the sheath dissipation closure was obtained, even in the presence of significant parallel gradients, where the closure is no longer valid. Specifically it captured the contrasting dynamics of filaments with different perpen- dicular sizes that were observed in the 3D simulations which the vorticity advection closure failed to replicate. However, neither closure successfully replicated the Boltzmann spinning effects and associated poloidal drift of the blob that was observed in the 3D simulations. Although the sheath dissipation closure was concluded to be more successful in replicating the 3D dynamics, it is emphasised that the vorticity closure may still be relevant for situations where the parallel current is inhibited from closing through the sheath due to effects such as increased resistivity near the targets.

76 citations


Journal ArticleDOI
TL;DR: In this paper, the impact of electromagnetic stabilization and flow shear stabilization on ITG turbulence is investigated, and it is found that at the inner radius, nonlinear electromagnetic stabilization is dominant, and is critical for achieving simulated heat fluxes in agreement with the experiment.
Abstract: The impact of electromagnetic stabilization and flow shear stabilization on ITG turbulence is investigated. Analysis of a low-$\beta$ JET L-mode discharge illustrates the relation between ITG stabilization, and proximity to the electromagnetic instability threshold. This threshold is reduced by suprathermal pressure gradients, highlighting the effectiveness of fast ions in ITG stabilization. Extensive linear and nonlinear gyrokinetic simulations are then carried out for the high-$\beta$ JET hybrid discharge 75225, at two separate locations at inner and outer radii. It is found that at the inner radius, nonlinear electromagnetic stabilization is dominant, and is critical for achieving simulated heat fluxes in agreement with the experiment. The enhancement of this effect by suprathermal pressure also remains significant. It is also found that flow shear stabilization is not effective at the inner radii. However, at outer radii the situation is reversed. Electromagnetic stabilization is negligible while the flow shear stabilization is significant. These results constitute the high-$\beta$ generalization of comparable observations found at low-$\beta$ at JET. This is encouraging for the extrapolation of electromagnetic ITG stabilization to future devices. An estimation of the impact of this effect on the ITER hybrid scenario leads to a 20% fusion power improvement.

68 citations


Journal ArticleDOI
TL;DR: In this paper, the electromagnetic energy equation is analyzed term by term in a 3D simulation of kinetic reconnection and the evolution presents the usual 2D-like topological structures caused by an initial perturbation independent of the third dimension.
Abstract: The electromagnetic energy equation is analyzed term by term in a 3D simulation of kinetic reconnection previously reported by Vapirev et al. [1]. The evolution presents the usual 2D-like topological structures caused by an initial perturbation independent of the third dimension. However, downstream of the reconnection site, where the jetting plasma encounters the yet unperturbed pre-existing plasma, a downstream front (DF) is formed and made unstable by the strong density gradient and the unfavorable local acceleration eld. The energy exchange between plasma and elds

61 citations


Journal ArticleDOI
TL;DR: In this paper, the authors examined electron bulk heating during magnetic reconnection with symmetric inflow conditions using kinetic particle-in-cell (PIC) simulations and found that the degree of electron heating is well correlated with the inflowing Alfv\'en speed.
Abstract: Electron bulk heating during magnetic reconnection with symmetric inflow conditions is examined using kinetic particle-in-cell (PIC) simulations. The degree of electron heating is well correlated with the inflowing Alfv\'en speed $c_{Ar}$ based on the reconnecting magnetic field through the relation $\Delta T_e = 0.033 \,m_i\,c_{Ar}^2$, where $\Delta T_{e}$ is the increase in electron temperature. For the range of simulations performed, the heating shows almost no correlation with inflow total temperature $T_{tot} = T_i + T_e$ or plasma $\beta$. An out-of-plane (guide) magnetic field of similar magnitude to the reconnecting field does not affect the total heating, but it does quench perpendicular heating, with almost all heating being in the parallel direction. These results are qualitatively consistent with a recent statistical survey of electron heating in the dayside magnetopause, which also found that $\Delta T_e$ was proportional to the inflowing Alfv\'en speed. The net electron heating varies very little with distance downstream of the x-line. The simulations show at most a very weak dependence of electron heating on the ion to electron mass ratio. In the antiparallel reconnection case, the largely parallel heating is eventually isotropized downstream due a scattering mechanism such as stochastic particle motion or instabilities. The study highlights key properties that must be satisfied by an electron heating mechanism: (1) Preferential heating in the parallel direction; (2) Heating proportional to $m_i\,c_{Ar}^2$; (3) At most a weak dependence on electron mass; and (4) An exhaust electron temperature that varies little with distance from the x-line.

61 citations


Journal ArticleDOI
TL;DR: In this article, the authors present an experimental-computational study of the parameters necessary to optimize a tunable, quasi-monoenergetic, efficient, low-background Compton backscattering (CBS) x-ray source that is based on the self-aligned combination of a laser-plasma accelerator (LPA) and a plasma mirror (PM).
Abstract: We present an in-depth experimental-computational study of the parameters necessary to optimize a tunable, quasi-monoenergetic, efficient, low-background Compton backscattering (CBS) x-ray source that is based on the self-aligned combination of a laser-plasma accelerator (LPA) and a plasma mirror (PM). The main findings are: (1) an LPA driven in the blowout regime by 30 TW, 30 fs laser pulses producesnot only a high-quality, tunable, quasi-monoenergetic electron beam, but also a high-quality, relativistically intense (a0~1) spent drive pulse that remains stable in profile and intensity over the LPA tuning range. (2) A thin plastic film near the gas jet exit retro-reflects the spent drive pulse efficiently into oncoming electrons to produce CBS x-rays without detectable bremsstrahlung background. Meanwhile anomalous far-field divergence of the retro-reflected light demonstrates relativistic "denting" of the PM. Exploiting these optimized LPA and PM conditions, we demonstrate quasi-monoenergetic (50% FWHM energy spread), tunable (75 to 200 KeV) CBS x-rays, characteristics previously achieved only on more powerful laser systems by CBS of a split-off, counter-propagating pulse. Moreover, laser-to-x-ray photon conversion efficiency ~6e12 exceeds that of any previous LPA-based quasi-monoenergetic Compton source. Particle-in-cell simulations agree well with the measurements.

Posted Content
TL;DR: In this paper, the progress made towards laser driven hadron cancer therapy and what has still to be accomplished to realize its inherent enormous potential is reviewed and compared with electron and photon therapy.
Abstract: It has been known for about sixty years that proton and heavy ion therapy is a very powerful radiation procedure for treating tumours. It has an innate ability to irradiate tumours with greater doses and spatial selectivity compared with electron and photon therapy and hence is a tissue sparing procedure. For more than twenty years powerful lasers have generated high energy beams of protons and heavy ions and hence it has been frequently speculated that lasers could be used as an alternative to RF accelerators to produce the particle beams necessary for cancer therapy. The present paper reviews the progress made towards laser driven hadron cancer therapy and what has still to be accomplished to realise its inherent enormous potential.

Journal ArticleDOI
TL;DR: In this paper, the authors summarized the recent advances and present conditions of the nanomaterials synthesis from the plasma-liquid interactions and presented a theoretical analysis for the synthesis process by analyzing the experimental data.
Abstract: In this review, we have summarized the recent advances and present conditions of the nanomaterials synthesis from the plasma-liquid interactions. A theoretical analysis for the nanomaterials synthesis process is presented by analyzing the experimental data. Besides the theoretical analysis, the practical applications in several nanomaterials syntheses of the the plasma-liquid interactions are also presented.

Journal ArticleDOI
TL;DR: In this article, a subproton scale magnetic hole was found in two-dimensional PIC simulations of decaying turbulence with a guide field, and the properties of these structures were performed with a realistic value for ion to electron mass ratio.
Abstract: We report the properties of a novel type of sub-proton scale magnetic hole found in two dimensional PIC simulations of decaying turbulence with a guide field. The simulations were performed with a realistic value for ion to electron mass ratio. These structures, electron vortex magnetic holes (EVMHs), have circular cross-section. The magnetic field depression is associated with a diamagnetic azimuthal current provided by a population of trapped electrons in petal-like orbits. The trapped electron population provides a mean azimuthal velocity and since trapping preferentially selects high pitch angles, a perpendicular temperature anisotropy. The structures arise out of initial perturbations in the course of the turbulent evolution of the plasma, and are stable over at least 100 electron gyroperiods. We have verified the model for the EVMH by carrying out test particle and PIC simulations of isolated structures in a uniform plasma. It is found that (quasi-)stable structures can be formed provided that there is some initial perpendicular temperature anisotropy at the structure location. The properties of these structures (scale size, trapped population etc.) are able to explain the observed properties of magnetic holes in the terrestrial plasma sheet. EVMHs may also contribute to turbulence properties, such as intermittency, at short scale lengths in other astrophysical plasmas.

Journal ArticleDOI
TL;DR: In this paper, a two-step Hartree-Fock-Slater model is proposed to describe the ionization potential depression (IPD) effect in dense plasmas.
Abstract: The charged environment within a dense plasma leads to the phenomenon of ionization potential depression (IPD) for ions embedded in the plasma. Accurate predictions of the IPD effect are of crucial importance for modeling atomic processes occurring within dense plasmas. Several theoretical models have been developed to describe the IPD effect, with frequently discrepant predictions. Only recently, first experiments on IPD in Al plasma have been performed with an x-ray free-electron laser (XFEL), where their results were found to be in disagreement with the widely-used IPD model by Stewart and Pyatt. Another experiment on Al, at the Orion laser, showed disagreement with the model by Ecker and Kroll. This controversy shows a strong need for a rigorous and consistent theoretical approach to calculate the IPD effect. Here we propose such an approach: a two-step Hartree-Fock-Slater model. With this parameter-free model we can accurately and efficiently describe the experimental Al data and validate the accuracy of standard IPD models. Our model can be a useful tool for calculating atomic properties within dense plasmas with wide-ranging applications to studies on warm dense matter, shock experiments, planetary science, inertial confinement fusion and studies of non-equilibrium plasmas created with XFELs.

Journal ArticleDOI
TL;DR: In this article, it was shown that an electromagnetic Weibel shock always forms when two relativistic collisionless, initially unmagnetized, plasma shells encounter, and the predicted shock formation time is in good agreement with 2D and 3D particle-in-cell simulations of counterstreaming pair plasmas.
Abstract: Collisionless shocks in plasmas play an important role in space physics (Earth's bow shock) and astrophysics (supernova remnants, relativistic jets, gamma-ray bursts, high energy cosmic rays). While the formation of a fluid shock through the steepening of a large amplitude sound wave has been understood for long, there is currently no detailed picture of the mechanism responsible for the formation of a collisionless shock. We unravel the physical mechanism at work and show that an electromagnetic Weibel shock always forms when two relativistic collisionless, initially unmagnetized, plasma shells encounter. The predicted shock formation time is in good agreement with 2D and 3D particle-in-cell simulations of counterstreaming pair plasmas. By predicting the shock formation time, experimental setups aiming at producing such shocks can be optimised to favourable conditions.

Journal ArticleDOI
TL;DR: In this article, the general two-fluid action functional is used to ensure that the reduced models are Hamiltonian, and the resulting equations are then mapped to the Eulerian fluid variables using a novel nonlocal Lagrange-Euler map.
Abstract: The general, non-dissipative, two-fluid model in plasma physics is Hamiltonian, but this property is sometimes lost or obscured in the process of deriving simplified (or reduced) two-fluid or one-fluid models from the two-fluid equations of motion. To ensure that the reduced models are Hamiltonian, we start with the general two-fluid action functional, and make all the approximations, changes of variables, and expansions directly within the action context. The resulting equations are then mapped to the Eulerian fluid variables using a novel nonlocal Lagrange-Euler map. Using this method, we recover L\"{u}st's general two-fluid model, extended MHD, Hall MHD, and electron MHD from a unified framework. The variational formulation allows us to use Noether's theorem to derive conserved quantities for each symmetry of the action.

Journal ArticleDOI
TL;DR: In this paper, the Coulomb exchange interaction in three and two-dimensional quantum plasmas was studied and the force densities were derived for all polarizations of the exchange interaction.
Abstract: We describe quantum hydrodynamic equations with the Coulomb exchange interaction for three and two dimensional plasmas. Explicit form of the force densities are derived. We present non-linear Schrodinger equations (NLSEs) for the Coulomb quantum plasmas with the exchange interaction. We show contribution of the exchange interaction in the dispersion of the Langmuir, and ion-acoustic waves. We consider influence of the spin polarization ratio on strength of the Coulomb exchange interaction. This is important since exchange interaction between particles with same spin direction and particles with opposite spin directions are different. At small particle concentrations $n_{0} >10^{25}cm^{-3}$ the Fermi pressure prevails over the exchange interaction for all polarizations. Similar picture we obtain for two dimensional quantum plasmas.

Journal ArticleDOI
Abstract: With the recent advances in laser technology, experimental investigation of radiation reaction phenomena is at last becoming a realistic prospect. A pedagogical introduction to electromagnetic radiation reaction is given with the emphasis on matter driven by ultra-intense lasers. Single-particle, multi-particle, classical and quantum aspects are all addressed.

Journal ArticleDOI
TL;DR: In this article, the potential energies curves for the N$_2$ electronic ground state, taken from literature, and the corresponding rate coefficients for electron-N$-2$ collisions occurring through the resonant state, obtained from $R$-matrix calculations are reviewed.
Abstract: Resonant vibrational excitation cross sections and the corresponding rate coefficients for electron-N$_2$ collisions occurring through the N$_2^-(\textrm{X}\ ^2\Pi_g)$ resonant state are reviewed. New calculations are performed using accurate potential energies curves for the N$_2$ electronic ground state, taken from literature, and for the N$_2^-$ resonant state, obtained from $R$-matrix calculations. The calculations are extended also to the resonant excitation processes involving the N$_2$ ground state vibrational continuum, leading to dissociation. Electron impact dissociation is found to be significant from higher vibrational levels. Accurate analytical fits for the complete set of the rate coefficients are provided. The behavior of the dissociative cross sections is investigated for rotationally excited N$_2$ molecules, with $J=50,100$ and 150 and for different vibrational levels.

Journal ArticleDOI
TL;DR: In this paper, a comparison of three intervals of Ulysses data with synthetic time series generated from a physically motivated turbulence model is made, and it is hypothesized that the observed signatures result from a perpendicular cascade of Alvenic fluctuations and a local, non-turbulent population of ion cyclotron or whistler waves generated by temperature anisotropy instabilities.
Abstract: Motivated by recent observations of distinct parallel and perpendicular signatures in magnetic helicity measurements segregated by wave period and angle between the local magnetic field and the solar wind velocity, this paper undertakes a comparison of three intervals of \emph{Ulysses} data with synthetic time series generated from a physically motivated turbulence model From these comparisons, it is hypothesized that the observed signatures result from a perpendicular cascade of \Alfvenic fluctuations and a local, non-turbulent population of ion cyclotron or whistler waves generated by temperature anisotropy instabilities By constraining the model's free parameters through comparison to \emph{in situ} data, it is found that, on average, $ \sim 95\%$ of the power near dissipative scales is contained in a perpendicular \Alfvenic cascade and that the parallel fluctuations are propagating nearly unidirectionally The effects of aliasing on magnetic helicity measurements are considered and shown to be significant near the Nyquist frequency

Journal ArticleDOI
TL;DR: In this paper, a general, fast, and effective approach for numerical calculation of kinetic plasma dispersion relations is developed, which is based on a standard matrix eigenvalue problem of an equivalent linear system.
Abstract: A general, fast, and effective approach is developed for numerical calculation of kinetic plasma dispersion relations. The plasma dispersion function is approximated by $J$-pole expansion. Subsequently, the dispersion relation is transformed to a standard matrix eigenvalue problem of an equivalent linear system. The result is accurate for $J=8$ except the solutions that are the little interesting heavily damped modes. In contrast to conventional approaches, such as Newton's iterative method, this approach can give either all the solutions in the system or a few solutions around the initial guess. It is also free from convergent problems. The approach is demonstrated from electrostatic one-dimensional and three-dimensional dispersion relations, to electromagnetic kinetic magnetized plasma dispersion relation for bi-Maxwellian distribution with parallel velocity drift.

Journal ArticleDOI
TL;DR: In this paper, the authors describe a realistic experimental configuration to measure collective effects that couple strong field quantum electrodynamics to plasma kinetics, and show that when quantum effects are included, a foil that classically would effectively transmit the laser pulse becomes opaque.
Abstract: With the aid of large-scale three-dimensional QED-PIC simulations, we describe a realistic experimental configuration to measure collective effects that couple strong field quantum electrodynamics to plasma kinetics. For two counter propagating lasers interacting with a foil at intensities exceeding $10^{22}$ Wcm$^{-2}$, a binary result occurs; when quantum effects are included, a foil that classically would effectively transmit the laser pulse becomes opaque. This is a dramatic change in plasma behavior, directly as a consequence of the coupling of radiation reaction and pair production to plasma dynamics.

Journal ArticleDOI
TL;DR: BOUT++ as mentioned in this paper is a 3D nonlinear finite-difference plasma simulation code, capable of solving quite general systems of PDEs, but targeted particularly on studies of the edge region of tokamak plasmas.
Abstract: BOUT++ is a 3D nonlinear finite-difference plasma simulation code, capable of solving quite general systems of PDEs, but targeted particularly on studies of the edge region of tokamak plasmas. BOUT++ is publicly available, and has been adopted by a growing number of researchers worldwide. Here we present improvements which have been made to the code since its original release, both in terms of structure and its capabilities. Some recent applications of these methods are reviewed, and areas of active development are discussed. We also present algorithms and tools which have been developed to enable creation of inputs from analytic expressions and experimental data, and for processing and visualisation of output results. This includes a new tool Hypnotoad for the creation of meshes from experimental equilibria. Algorithms have been implemented in BOUT++ to solve a range of linear algebraic problems encountered in the simulation of reduced MHD and gyro-fluid models: A preconditioning scheme is presented which enables the plasma potential to be calculated efficiently using iterative methods supplied by the PETSc library, without invoking the Boussinesq approximation. Scaling studies are also performed of a linear solver used as part of physics-based preconditioning to accelerate the convergence of implicit time-integration schemes.

Journal ArticleDOI
TL;DR: In this article, a 1D Vlasov-Maxwell code was used to evaluate the effect of the pump intensity on the energy transfer from the pump to the seed pulse in backward Raman amplifiers.
Abstract: In plasma-based backward Raman amplifiers, the output pulse intensity increases with the input pump pulse intensity, as long as the Langmuir wave mediating energy transfer from the pump to the seed pulse remains intact. However, at high pump intensity, the Langmuir wave breaks, at which point the amplification efficiency may no longer increase with the pump intensity. Numerical simulations presented here, employing a 1D Vlasov-Maxwell code, show that, although the amplification efficiency remains high when the pump only mildly exceeds the wavebreaking threshold, the efficiency drops precipitously at larger pump intensities.

Journal ArticleDOI
TL;DR: In this article, a semi-empirical approximation was constructed based on the results of the \textit{ab initio} caculation, which yields the dependences of the quantum molecular dynamics, density functional theory and the Kubo-Greenwood formula.
Abstract: This work is devoted to the investigation of transport and optical properties of liquid aluminum in the two-temperature case. At first optical properties, static electrical and thermal conductivities were obtained in the \textit{ab initio} calculation. The \textit{ab initio} calculation is based on the quantum molecular dynamics, density functional theory and the Kubo-Greenwood formula. The semiempirical approximation was constructed based on the results of the \textit{ab initio} caculation. The approximation yields the dependences $\sigma_{1_\mathrm{DC}}\propto1/T_i^{0.25}$ and $K\propto T_e/T_i^{0.25}$ for the static electrical conductivity and thermal conductivity, respectively. The approximation is valid for liquid aluminum at $\rho=2.70$~g/cm$^3$, 3~kK~$\leq T_i\leq T_e\leq20$~kK. Our results are well described by the Drude model with the effective relaxation time $\tau\propto T_i^{-0.25}$. We have compared our results with a number of other models. They are all reduced in the low-temperature limit to the Drude model with different expressions for the relaxation time $\tau$. Our results are not consistent with the models in which $\tau\propto T_i^{-1}$ and support the models which use the expressions with the slower decrease of the relaxation time.

Journal ArticleDOI
TL;DR: In this article, the fast electron beam guiding using externally applied longitudinal magnetic fields was proposed for enhancing the core heating efficiency in electron-driven fast ignition, which demonstrated the sufficient beam guiding performance in the collisional dense plasma by kT-class external magnetic fields for the case with moderate mirror ratio.
Abstract: For enhancing the core heating efficiency in electron-driven fast ignition, we proposed the fast electron beam guiding using externally applied longitudinal magnetic fields. Based on the PIC simulations for the FIREX-class experiments, we demonstrated the sufficient beam guiding performance in the collisional dense plasma by kT-class external magnetic fields for the case with moderate mirror ratio (~<10 ). Boring of the mirror field was found through the formation of magnetic pipe structure due to the resistive effects, which indicates a possibility of beam guiding in high mirror field for higher laser intensity and/or longer pulse duration.

Journal ArticleDOI
TL;DR: Parker et al. as discussed by the authors studied Landau damping in the 1+1D Vlasov-Poisson system using a Fourier-Hermite spectral representation, where the change in the electric field corresponds to the net Hermite flux via a free energy evolution equation.
Abstract: We study Landau damping in the 1+1D Vlasov-Poisson system using a Fourier-Hermite spectral representation. We describe the propagation of free energy in phase space using forwards and backwards propagating Hermite modes recently developed for gyrokinetics [Schekochihin et al. (2014)]. The change in the electric field corresponds to the net Hermite flux via a free energy evolution equation. In linear Landau damping, decay in the electric field corresponds to forward propagating Hermite modes; in nonlinear damping, the initial decay is followed by a growth phase characterised by the generation of backwards propagating Hermite modes by the nonlinear term. The free energy content of the backwards propagating modes increases exponentially until balancing that of the forward propagating modes. Thereafter there is no systematic net Hermite flux, so the electric field cannot decay and the nonlinearity effectively suppresses Landau damping. These simulations are performed using the fully-spectral 5D gyrokinetics code SpectroGK [Parker et al. 2014], modified to solve the 1+1D Vlasov-Poisson system. This captures Landau damping via an iterated L\'enard-Bernstein collision operator or via Hou-Li filtering in velocity space. Therefore the code is applicable even in regimes where phase-mixing and filamentation are dominant.

Journal ArticleDOI
TL;DR: In this paper, the authors unearth a basic Lie algebra that generates the Poisson bracket and a class of similar Poisson algebra may be generated by the same Lie algebra, which encompasses the Hall MHD system and inertial MHD systems.
Abstract: The extended magnetohydrodynamics (MHD) system, including the Hall effect and the electron inertia effect, has a Hamiltonian structure embodied by a noncanonical Poisson algebra on an infinite-dimensional phase space. A nontrivial part of the formulation is the proof of Jacobi's identity for the Poisson bracket. We unearth a basic Lie algebra that generates the Poisson bracket. A class of similar Poisson algebra may be generated by the same Lie algebra, which encompasses the Hall MHD system and inertial MHD system.

Journal ArticleDOI
TL;DR: The former all-carbon wall on JET has been replaced with beryllium in the main torus and tungsten in the divertor to mimic the surface materials envisaged for ITER as mentioned in this paper.
Abstract: The former all-carbon wall on JET has been replaced with beryllium in the main torus and tungsten in the divertor to mimic the surface materials envisaged for ITER. Comparisons are presented between Type I H-mode characteristics in each design by examining respective scans over deuterium fuelling and impurity seeding, required to ameliorate exhaust loads both in JET at full capability and in ITER.