ALBERT

Description: This investigation analyzes the effect of vortex wakes on the Lagrangian displacement of particles induced by the passage of an obstacle in a two-dimensional incompressible and inviscid fluid. In addition to the trajectories of individual particles, we also study their drift and the corresponding total drift areas in the Föppl and Kirchhoff potential flow models. Our findings, which are obtained numerically and in some regimes are also supported by asymptotic analysis, are compared to the wakeless potential flow which serves as a reference. We show that in the presence of the Föppl vortex wake, some of the particles follow more complicated trajectories featuring a second loop. The appearance of an additional stagnation point in the Föppl flow is identified as a source of this effect. It is also demonstrated that, while the total drift area increases with the size of the wake for large vortex strengths, it is actually decreased for small circulation values. On the other hand, the Kirchhoff flow model is shown to have an unbounded total drift area. By providing a systematic account of the wake effects on the drift, the results of this study will allow for more accurate modeling of hydrodynamic stirring.

Description: The dynamics of the magnetic relaxation process during the sustainment of spheromak configurations at different helicity injection rates is studied. The three-dimensional activity is recovered using time-dependent resistive magnetohydrodynamic simulations. A cylindrical flux conserver with concentric electrodes is used to model configurations driven by a magnetized coaxial gun. Magnetic helicity is injected by tangential boundary flows. Different regimes of sustainment are identified and characterized in terms of the safety factor profile. The spatial and temporal behavior of fluctuations is described. The dynamo action is shown to be in close agreement with existing experimental data. These results are relevant to the design and operation of helicity injected devices, as well as to basic understanding of the plasma relaxation mechanism in quasi-steady state.

Description: The present experiments were performed on the PF-1000 plasma focus device at a current of 2 MA with the deuterium injected from the gas-puff placed in the axis of the anode face. The XUV frames showed, in contrast with the interferograms, the fine structure: filaments and spots up to 1 mm diameter. In the deuterium filling, the short filaments are registered mainly in the region of the internal plasmoidal structures and their number correlates with the intensity of neutron production. The longer filamentary structure was recorded close to the anode after the constriction decay. The long curve-like filaments with spots were registered in the big bubble formed after the pinch phase in the head of the umbrella shape of the plasma sheath. Filaments can indicate the filamentary structure of the current in the pinch. Together with the filaments, small compact balls a few mm in diameter were registered by both interferometry and XUV frame pictures. They emerge out of the dense column and their life-time can be greater than hundreds of ns.

Description: The radiation pressure acceleration of protons in the interaction of Gaussian laser pulses and surface modulated targets is examined by multi-dimensional particle-in-cell simulations. It is shown that strong longitudinal quasi-static magnetic field is generated on the modulated surface of the target, which significantly enhances the transverse diffusion of electrons. This is beneficial for suppressing the transverse Rayleigh-Taylor instability. Finally, the surface of the accelerated proton beams becomes smoother than that in the case of the planar target, and a final mono-energetic proton beam is obtained by using the surface modulated target.

Description: The dependence of divertor asymmetry and scrape-off layer (SOL) flow on heating power has been investigated in the Experimental Advanced Superconducting Tokamak (EAST). Divertor plasma exhibits an outboard-enhanced in-out asymmetry in heat flux in lower single null configuration for in reversed (ion ∇ B drift direction toward the upper X-point) field directions. Upper single null exhibits an inboard-favored asymmetry in low heating power condition, while exhibits an outboard-favored asymmetry when increasing the heating power. Double null has the strongest in-out asymmetry in heat flux, favoring the outer divertor. The in-out asymmetry ratios of q t , out / q t , in and P out / P total increase with the power across the separatrix P loss , which is probably induced by the enhanced radial particle transport due to a large pressure gradient. The characteristics of the measured SOL parallel flow under various discharge conditions are consistent with the Pfirsch-Schlüter (PS) flow with the parallel Mach number M ∥ decreasing with the line averaged density but increasing with P loss , in the same direction as the PS flow. The contributions of both poloidal E × B drift and parallel flow on poloidal particle transport in SOL on EAST are also assessed.

Description: Magnetic islands have been implemented in the gyrokinetic toroidal code to study the effects of the islands on microturbulence. The pressure profile flattening is verified in the simulation with the islands. Simulations of ion temperature gradient instability find that different toroidal modes are linearly coupled together and that toroidal spectra become broader when the island width increases. The real frequencies and growth rates of different toroidal modes approach each other with the averaged value independent of the island width. The linear mode structures are enhanced at the island separatrices and weakened at the island centers, consistent with the flattening of the pressure profile inside the islands.

Description: The slow motion of a circular cylinder in a plane Poiseuille flow in a microchannel is analyzed for a wide range of cylinder radii and positions across the channel. The cylinder translates parallel to the channel walls and rotates about its axis. The Stokes approximation is used and the problem is solved analytically using the Papkovich-Fadle eigenfunction expansion and the least-squares method. The stream function and the pressure distribution of the flow field are obtained as results. The force and moment exerted on the cylinder, and the pressure change far from the cylinder, are calculated and shown as functions of the size and location of the cylinder. The results confirm some reciprocal relations exactly. In particular, the translational and rotational velocities of the drifting cylinder in the existing Poiseuille flow are determined. The induced pressure change, when the cylinder drifts in the Poiseuille flow, is also calculated. Some typical streamline patterns, depending on the size and location of the cylinder, are shown and discussed. When the cylinder translates and/or rotates in the channel blocked at infinity, a series of Moffatt eddies appears far from the cylinder in the channel, as expected.

Description: Interactions between capillary and elastic effects are relevant to a variety of applications from micro- and nano-scale manufacturing to biological systems. In this work, we investigate capillary flows in flexible, millimeter-scale cylindrical elastic tubes. We demonstrate that surface tension can cause sufficiently flexible tubes to collapse and coalesce spontaneously through non-axisymmetric buckling, and develop criteria for the initial deformation and complete collapse of a circular tube. Experimental results for capillary rise and evaporation of a liquid in a flexible tube are presented, and several regimes are seen for the equilibrium state of a flexible tube deforming under capillary pressure. Deformations of the tube walls are measured in different regimes and compared with a shell theory model. Analysis and experimental results show that despite the complex and non-axisymmetric deformed shapes of cylindrical structures, the elastocapillary length used in previous literature for flat plates and sheets can be used to predict the behavior of flexible tubes.

Description: A study is presented for the oblique propagation of ion acoustic cnoidal waves in a magnetized plasma consisting of cold ions and two temperature superthermal electrons modelled by kappa-type distributions. Using the reductive perturbation method, the nonlinear Korteweg de-Vries equation is derived, which further gives the solutions with a special type of cnoidal elliptical functions. Both compressive and rarefactive structures are found for these cnoidal waves. Nonlinear periodic cnoidal waves are explained in terms of plasma parameters depicting the Sagdeev potential and the phase curves. It is found that the density ratio of hot electrons to ions μ significantly modifies compressive/refractive wave structures. Furthermore, the combined effects of superthermality of cold and hot electrons κ c , κ h , cold to hot electron temperature ratio σ , angle of propagation and ion cyclotron frequency ω ci have been studied in detail to analyze the height and width of compressive/refractive cnoidal waves. The findings in the present study could have important implications in understanding the physics of electrostatic wave structures in the Saturn's magnetosphere where two temperature superthermal electrons are present.

Description: We present laboratory experimental results demonstrating that librational forcing of an ellipsoidal container of water can produce intense motions through the mechanism of a libration driven elliptical instability (LDEI). These libration studies are conducted using an ellipsoidal acrylic container filled with water. A particle image velocimetry method is used to measure the 2D velocity field in the equatorial plane over hundreds libration cycles for a fixed Ekman number, E = 2 × 10 −5 . In doing so, we recover the libration induced base flow and a time averaged zonal flow. Further, we show that LDEI in non-axisymmetric container geometries is capable of driving both intermittent and saturated turbulent motions in the bulk fluid. Additionally, we measure the growth rate and amplitude of the LDEI induced excited flow in a fully ellipsoidal container at more extreme parameters than previously studied [Noir et al. , “Experimental study of libration-driven flows in nonaxisymmetric containers,” Phys. Earth Planet. Inter. 204-205 , 1 (2012); Cébron et al. , Phys. Fluids 24 , 061703, “Libration driven elliptical instability,” (2012)]. Excitation of bulk filling turbulence by librational forcing provides a mechanism for transferring rotational energy into turbulent fluid motion and thus can play an important role in the thermal evolution, interior dynamics, and magneto-hydrodynamics of librating bodies, as appear to be common in solar system settings [e.g., Comstock and Bills, “A solar system survey of forced librations in longitude,” J. Geophys. Res. Planets 108 , 1 (2003)].

Description: We numerically study the displacement flow of two iso-viscous Newtonian fluids in an inclined two-dimensional channel, formed by two parallel plates. The results are complementary to our previous studies on displacement flows in pipes and channels. The heavier displacing fluid moves the lighter displaced fluid in the downward direction. Three dimensionless groups largely describe these flows: the densimetric Froude number ( Fr ), the Reynolds number ( Re ), and the duct inclination (β). As a first order approximation, we are able to classify different flow regimes phenomenologically in a two-dimensional ( Fr ; Re cosβ/ Fr )-plane and provide leading order expressions for the transitions between different regimes. The stabilizing and/or de-stabilizing effects of the imposed mean flow on buoyant exchange flows (zero imposed velocity) are described for a broad range of dimensionless parameters.

Description: Compressible granular materials are involved in many applications, some of them being related to energetic porous media. Gas permeation effects are important during their compaction stage, as well as their eventual chemical decomposition. Also, many situations involve porous media separated from pure fluids through two-phase interfaces. It is thus important to develop theoretical and numerical formulations to deal with granular materials in the presence of both two-phase interfaces and gas permeation effects. Similar topic was addressed for fluid mixtures and interfaces with the Discrete Equations Method (DEM) [R. Abgrall and R. Saurel, “Discrete equations for physical and numerical compressible multiphase mixtures,” J. Comput. Phys. 186 (2), 361-396 (2003)] but it seemed impossible to extend this approach to granular media as intergranular stress [K. K. Kuo, V. Yang, and B. B. Moore, “Intragranular stress, particle-wall friction and speed of sound in granular propellant beds,” J. Ballist. 4 (1), 697-730 (1980)] and associated configuration energy [J. B. Bdzil, R. Menikoff, S. F. Son, A. K. Kapila, and D. S. Stewart, “Two-phase modeling of deflagration-to-detonation transition in granular materials: A critical examination of modeling issues,” Phys. Fluids 11 , 378 (1999)] were present with significant effects. An approach to deal with fluid-porous media interfaces was derived in Saurel et al. [“Modelling dynamic and irreversible powder compaction,” J. Fluid Mech. 664 , 348-396 (2010)] but its validity was restricted to weak velocity disequilibrium only. Thanks to a deeper analysis, the DEM is successfully extended to granular media modelling in the present paper. It results in an enhanced version of the Baer and Nunziato [“A two-phase mixture theory for the deflagration-to-detonation transition (DDT) in reactive granular materials,” Int. J. Multiphase Flow 12 (6), 861-889 (1986)] model as symmetry of the formulation is now preserved. Several computational examples are shown to validate and illustrate method’s capabilities.

Description: We perform a theoretical and numerical study of the Coulomb-driven electroconvection flow of a dielectric liquid between two coaxial cylinders. The specific case, where the inner to outer diameter ratio is 0.5, is analyzed. A strong unipolar injection of ions either from the inner or outer cylinder is considered to introduce free charge carriers into the system. A finite volume method is used to solve all governing equations including Navier-Stokes equations and a simplified set of Maxwell’s equations. The flow is characterized by a subcritical bifurcation in the finite amplitude regime. A linear stability criterion and a nonlinear one that correspond to the onset and stop of the flow motion, respectively, are linked with a hysteresis loop. In addition, we also explore the behavior of the system for higher values of the stability parameter. For inner injection, we observe a transition between the patterns made of 7 and 8 cells, before an oscillatory regime is attained. Such a transition leads to a second finite amplitude stability criterion. A simple modal analysis reveals that the competition of different modes is at the origin of this behavior. The charge density, as well as velocity field distributions is provided to help understand the bifurcation behavior.

Description: We report on the production of metal ions of magnesium and zinc in the beam plasma formed by a forevacuum-pressure electron source. Magnesium and zinc vapor were generated by electron beam evaporation from a crucible and subsequently ionized by electron impact from the e-beam itself. Both gaseous and metallic plasmas were separately produced and characterized using a modified RGA-100 quadrupole mass-spectrometer. The fractional composition of metal isotopes in the plasma corresponds to their fractional natural abundance.

Description: To better understand the role of residual electrons in the repeatability of an atmospheric pressure plasma plume, the characteristics of a helium plasma jet from the 1st, 2nd,… until the repeatable discharge pulse are investigated for the first time. It's found that the longest plasma plume is achieved in the 1st discharge pulse. The length of the plasma plume becomes shorter and shorter and reaches a constant value in the 3rd discharge pulse and keeps the same for the following pulses. The dynamics of the 1st discharge pulse show that the plasma bullet appears random in nature. Two photomultiplier tubes are used to distinguish the two potential factors that could result in the stochastic dynamics of the plasma bullet, i.e., stochastic ignition of the plasma and the stochastic propagation velocity. The results show that the stochastic propagation velocity occurs only in the 1st and the 2nd discharge pulses, while the stochastic ignition of the plasma presents until the 100th pulse. The dynamics of the plasma propagation become repeatable after about 100 pulses. Detail analysis shows that the repeatability of plasma bullet is due to the residual electrons density. The residual electron density of 10 9  cm −3 or higher is needed for repeatable discharges mode.

Description: We present a full two-fluid magnetohydrodynamic (MHD) description for a completely ionized hydrogen plasma, retaining the effects of the Hall current, electron pressure, and electron inertia. According to this description, each plasma species introduces a new spatial scale: the ion inertial length λ i and the electron inertial length λ e , which are not present in the traditional MHD description. In the present paper, we seek for possible changes in the energy power spectrum in fully developed turbulent regimes, using numerical simulations of the two-fluid equations in two-and-a-half dimensions. We have been able to reproduce different scaling laws in different spectral ranges, as it has been observed in the solar wind for the magnetic energy spectrum. At the smallest wavenumbers where plain MHD is valid, we obtain an inertial range following a Kolmogorov k −5∕3 law. For intermediate wavenumbers such that λ i − 1 ≪ k ≪ λ e − 1 , the spectrum is modified to a k −7∕3 power-law, as has also been obtained for Hall-MHD neglecting electron inertia terms. When electron inertia is retained, a new spectral region given by k 〉 λ e − 1 arises. The power spectrum for magnetic energy in this region is given by a k −11∕3 power law. Finally, when the terms of electron inertia are retained, we study the self-consistent electric field. Our results are discussed and compared with those obtained in the solar wind observations and previous simulations.

Description: A novel mechanism of collisionless heating in large planar arrays of small inductive coils operated at radio frequencies is presented. In contrast to the well-known case of non-local heating related to the transversal conductivity, when the electrons move perpendicular to the planar coil, we investigate the problem of electrons moving in a plane parallel to the coils. Two types of periodic structures are studied. Resonance velocities where heating is efficient are calculated analytically by solving the Vlasov equation. Certain scaling parameters are identified. The concept is further investigated by a single particle simulation based on the ergodic principle and combined with a Monte Carlo code allowing for collisions with Argon atoms. Resonances, energy exchange, and distribution functions are obtained. The analytical results are confirmed by the numerical simulation. Pressure and electric field dependences are studied. Stochastic heating is found to be most efficient when the electron mean free path exceeds the size of a single coil cell. Then the mean energy increases approximately exponentially with the electric field amplitude.

Description: The collision-less transfer of momentum and energy from explosive debris plasma to magnetized background plasma is a salient feature of various astrophysical and space environments. While much theoretical and computational work has investigated collision-less coupling mechanisms and relevant parameters, an experimental validation of the results demands the measurement of the complex, collective electric fields associated with debris-background plasma interaction. Emission spectroscopy offers a non-interfering diagnostic of electric fields via the Stark effect. A unique experiment at the University of California, Los Angeles, that combines the Large Plasma Device (LAPD) and the Phoenix laser facility has investigated the marginally super-Alfvénic, quasi-perpendicular expansion of a laser-produced carbon (C) debris plasma through a preformed, magnetized helium (He) background plasma via emission spectroscopy. Spectral profiles of the He II 468.6 nm line measured at the maximum extent of the diamagnetic cavity are observed to intensify, broaden, and develop equally spaced modulations in response to the explosive C debris, indicative of an energetic electron population and strong oscillatory electric fields. The profiles are analyzed via time-dependent Stark effect models corresponding to single-mode and multi-mode monochromatic (single frequency) electric fields, yielding temporally resolved magnitudes and frequencies. The proximity of the measured frequencies to the expected electron plasma frequency suggests the development of the electron beam-plasma instability, and a simple saturation model demonstrates that the measured magnitudes are feasible provided that a sufficiently fast electron population is generated during C debris–He background interaction. Potential sources of the fast electrons, which likely correspond to collision-less coupling mechanisms, are briefly considered.

Description: It is shown that co-linear injection of electrons or positrons into the wakefield of the self-modulating particle beam is possible and ensures high energy gain. The witness beam must co-propagate with the tail part of the driver, since the plasma wave phase velocity there can exceed the light velocity, which is necessary for efficient acceleration. If the witness beam is many wakefield periods long, then the trapped charge is limited by beam loading effects. The initial trapping is better for positrons, but at the acceleration stage a considerable fraction of positrons is lost from the wave. For efficient trapping of electrons, the plasma boundary must be sharp, with the density transition region shorter than several centimeters. Positrons are not susceptible to the initial plasma density gradient.

Description: The foundational theory for dusty plasmas is the dust charging theory that provides the dust potential and charge arising from the dust interaction with a plasma. The most widely used dust charging theory for negatively charged dust particles is the so-called orbital motion limited (OML) theory, which predicts the dust potential and heat collection accurately for a variety of applications, but was previously found to be incapable of evaluating the dust charge and plasma response in any situation. Here, we report a revised OML formulation that is able to predict the plasma response and hence the dust charge. Numerical solutions of the new OML model show that the widely used Whipple approximation of dust charge-potential relationship agrees with OML theory in the limit of small dust radius compared with plasma Debye length, but incurs large (order-unity) deviation from the OML prediction when the dust size becomes comparable with or larger than plasma Debye length. This latter case is expected for the important application of dust particles in a tokamak plasma.

Description: Using the ideal magnetohydrodynamic model, we calculate the temporal evolution of initial ripples on the boundaries of a planar plasma slab that is subjected to the magneto-Rayleigh-Taylor instability. The plasma slab consists of three regions. We assume that in each region the plasma density is constant with an arbitrary value and the magnetic field is also constant with an arbitrary magnitude and an arbitrary direction parallel to the interfaces. Thus, the instability may be driven by a combination of magnetic pressure and kinetic pressure. The general dispersion relation is derived, together with the feedthrough factor between the two interfaces. The temporal evolution is constructed from the superposition of the eigenmodes. Previously established results are recovered in the various limits. Numerical examples are given on the temporal evolution of ripples on the interfaces of the finite plasma slab.

Description: Phytoplankton patchiness, namely the heterogeneous distribution of microalgae over multiple spatial scales, dramatically impacts marine ecology. A spectacular example of such heterogeneity occurs in thin phytoplankton layers (TPLs), where large numbers of photosynthetic microorganisms are found within a small depth interval. Some species of motile phytoplankton can form TPLs by gyrotactic trapping due to the interplay of their particular swimming style (directed motion biased against gravity) and the transport by a flow with shear along the direction of gravity. Here we consider gyrotactic swimmers in numerical simulations of the Kolmogorov shear flow, both in laminar and turbulent regimes. In the laminar case, we show that the swimmer motion is integrable and the formation of TPLs can be fully characterized by means of dynamical systems tools. We then study the effects of rotational Brownian motion or turbulent fluctuations (appearing when the Reynolds number is large enough) on TPLs. In both cases, we show that TPLs become transient, and we characterize their persistence.

Description: Gas-puff Z-pinch experiments were conducted on the 1 MA, 200 ns pulse duration Cornell Beam Research Accelerator (COBRA) pulsed power generator in order to achieve an understanding of the dynamics and instability development in the imploding and stagnating plasma. The triple-nozzle gas-puff valve, pre-ionizer, and load hardware are described. Specific diagnostics for the gas-puff experiments, including a Planar Laser Induced Fluorescence system for measuring the radial neutral density profiles along with a Laser Shearing Interferometer and Laser Wavefront Analyzer for electron density measurements, are also described. The results of a series of experiments using two annular argon (Ar) and/or neon (Ne) gas shells (puff-on-puff) with or without an on- (or near-) axis wire are presented. For all of these experiments, plenum pressures were adjusted to hold the radial mass density profile as similar as possible. Initial implosion stability studies were performed using various combinations of the heavier (Ar) and lighter (Ne) gasses. Implosions with Ne in the outer shell and Ar in the inner were more stable than the opposite arrangement. Current waveforms can be adjusted on COBRA and it was found that the particular shape of the 200 ns current pulse affected on the duration and diameter of the stagnated pinched column and the x-ray yield.

Description: During July 2012, 150 almost identical H-mode plasmas were consecutively created in the Joint European Torus, providing a combined total of approximately 8 minutes of steady-state plasma with 15 000 Edge Localised Modes (ELMs). In principle, each of those 15 000 ELMs are statistically equivalent. Here, the changes in edge density and plasma energy associated with those ELMs are explored, using the spikes in Beryllium II (527 nm) radiation as an indicator for the onset of an ELM. Clearly different timescales are observed during the ELM process. Edge temperature falls over a 2 ms timescale, edge density and pressure fall over a 5 ms timescale, and there is an additional 10 ms timescale that is consistent with a resistive relaxation of the plasma's edge. The statistical properties of the energy and density losses due to the ELMs are explored. For these plasmas the ELM energy (δ E ) is found to be approximately independent of the time between ELMs, despite the average ELM energy ( 〈 δ E 〉 ) and average ELM frequency ( f ) being consistent with the scaling of 〈 δ E 〉 ∝ 1 / f . Instead, beyond the first 0.02 s of waiting time between ELMs, the energy losses due to individual ELMs are found to be statistically the same. Surprisingly no correlation is found between the energies of consecutive ELMs either. A weak link is found between the density drop and the ELM waiting time. Consequences of these results for ELM control and modelling are discussed.

Description: The compression of a relativistic Gaussian laser pulse in a magnetized plasma is investigated. By considering relativistic nonlinearity and using non-linear Schrödinger equation with paraxial approximation, a second-order differential equation is obtained for the pulse width parameter (in time) to demonstrate the longitudinal pulse compression. The compression of laser pulse in a magnetized plasma can be observed by the numerical solution of the equation for the pulse width parameter. The effects of magnetic field and chirping are investigated. It is shown that in the presence of magnetic field and negative initial chirp, compression of pulse is significantly enhanced.

Description: Nonlinear saturation amplitudes (NSAs) of the first two harmonics in classical Rayleigh-Taylor instability (RTI) in cylindrical geometry for arbitrary Atwood numbers have been analytically investigated considering nonlinear corrections up to the fourth-order. The NSA of the fundamental mode is defined as the linear (purely exponential) growth amplitude of the fundamental mode at the saturation time when the growth of the fundamental mode (first harmonic) is reduced by 10% in comparison to its corresponding linear growth, and the NSA of the second harmonic can be obtained in the same way. The analytic results indicate that the effects of the initial radius of the interface ( r 0 ) and the Atwood number ( A ) play an important role in the NSAs of the first two harmonics in cylindrical RTI. On the one hand, the NSA of the fundamental mode first increases slightly and then decreases quickly with increasing A . For given A , the smaller the r 0 / λ (with λ perturbation wavelength) is, the larger the NSA of the fundamental mode is. When r 0 / λ is large enough ( r 0 ≫ λ ), the NSA of the fundamental mode is reduced to the prediction of previous literatures within the framework of third-order perturbation theory [J. W. Jacobs and I. Catton, J. Fluid Mech. 187 , 329 (1988); S. W. Haan, Phys. Fluids B 3 , 2349 (1991)]. On the other hand, the NSA of the second harmonic first decreases quickly with increasing A , reaching a minimum, and then increases slowly. Furthermore, the r 0 can reduce the NSA of the second harmonic for arbitrary A at r 0 ≲ 2 λ while increase it for A   ≲   0.6 at r 0 ≳ 2 λ . Thus, it should be included in applications where the NSA has a role, such as inertial confinement fusion ignition target design.

Description: The δ   f particle-in-cell code GEM is used to study the transport “shortfall” problem of gyrokinetic simulations. In local simulations, the GEM results confirm the previously reported simulation results of DIII-D [Holland et al. , Phys. Plasmas 16 , 052301 (2009)] and Alcator C-Mod [Howard et al. , Nucl. Fusion 53 , 123011 (2013)] tokamaks with the continuum code GYRO. Namely, for DIII-D the simulations closely predict the ion heat flux at the core, while substantially underpredict transport towards the edge; while for Alcator C-Mod, the simulations show agreement with the experimental values of ion heat flux, at least within the range of experimental error. Global simulations are carried out for DIII-D L-mode plasmas to study the effect of edge turbulence on the outer core ion heat transport. The edge turbulence enhances the outer core ion heat transport through turbulence spreading. However, this edge turbulence spreading effect is not enough to explain the transport underprediction.

Description: The model of gas bubble growth in high-viscous gas-saturated magmatic melt, subjected to rapid decompression, is presented in the current study. It is shown that consideration of unsteady character of the process is extremely important in a wide range of supersaturation. The analytical solution is found for the profile of dissolved gas concentration and the rate of bubble growth. The model of kinetics of overall degassing is developed. This model is based on distinguishing the so-called “forbidden” zone in the melt volume with suppressed formation of the new nucleation sites. The simple analytical dependences of the number of nucleating bubbles and typical nucleation time on the value of initial decompression were derived together with time dependence of volumetric concentration of the gas phase. Our results match the available experimental data.

Description: In the present work, we extend results of a previous paper [Peter et al. , Phys. Plasmas 20 , 12 3104 (2013)] and develop a semi-analytical model to account for thermal effects on the nonlinear dynamics of the electron beam in free-electron lasers. We relax the condition of a cold electron beam but still use the concept of compressibility, now associated with a warm beam model, to evaluate the time scale for saturation and the peak laser intensity in high-gain regimes. Although vanishing compressibilites and the associated divergent densities are absent in warm models, a series of discontinuities in the electron density precede the saturation process. We show that full wave-particle simulations agree well with the predictions of the model.

Description: We investigate the effect of viscosity contrast on the stability of gravitationally unstable, diffusive layers in porous media. Our analysis helps evaluate experimental observations of various diffusive (boundary) layer models that are commonly used to study the sequestration of CO 2 in brine aquifers. We evaluate the effect of viscosity contrast for two basic models that are characterized with respect to whether or not the interface between CO 2 and brine is allowed to move. We find that diffusive layers are in general more unstable when viscosity decreases with depth within the layer compared to when viscosity increases with depth. This behavior is in contrast to the one associated with the classical displacement problem of gravitationally unstable diffusive layers that are subject to mean flow. For the classical problem, a greater instability is associated with the displacement of a more viscous, lighter fluid along the direction of gravity by a less viscous, heavier fluid. We show that the contrasting behavior highlighted in this study is a special case of the classical displacement problem that depends on the relative strength of the displacement and buoyancy velocities. We demonstrate the existence of a critical viscosity ratio that determines whether the flow is buoyancy dominated or displacement dominated. We explain the new behaviors in terms of the interaction of vorticity components related to gravitational and viscous effects.

Description: Particle pairing in a complex plasma was experimentally studied with the emphasis on pair spatial extent and stability. Micron-size particles were suspended in the (pre)sheath area above the lower electrode in a capacitively coupled radio-frequency discharge in argon. They formed vertical pairs due to the ion wakes created by the flow of ions past particles. We discuss the confinement mechanism for the lower particle, resulting from a combination of the wake field and the field of non-uniform sheath. A model of particle pairs is proposed, which provides good description for the dependence of pair size and stability on experimental parameters.

Description: Significantly collimated fast electron beam with a divergence angle 10° (FWHM) is observed when an ultra-intense laser pulse ( I  = 10 14  W/cm 2 , 300 fs) irradiates a uniform critical density plasma. The uniform plasma is created through the ionization of an ultra-low density (5 mg/c.c.) plastic foam by X-ray burst from the interaction of intense laser ( I  = 10 14  W/cm 2 , 600 ps) with a thin Cu foil. 2D Particle-In-Cell (PIC) simulation well reproduces the collimated electron beam with a strong magnetic field in the region of the laser pulse propagation. To understand the physical mechanism of the collimation, we calculate energetic electron motion in the magnetic field obtained from the 2D PIC simulation. As the results, the strong magnetic field (300 MG) collimates electrons with energy over a few MeV. This collimation mechanism may attract attention in many applications such as electron acceleration, electron microscope and fast ignition of laser fusion.

Description: In this work, stochastic behavior of atmospheric pressure dielectric barrier discharge (DBD) has been investigated. The experiment is performed in a DBD reactor consisting of a pair of stainless steel parallel plate electrodes powered by a 50 Hz ac high voltage source. Current pulse amplitude distributions for different space gaps and the time separation between consecutive current pulses are studied. A probability distribution function is proposed to predict the experimental distribution function for the current pulse amplitudes and the occurrence of the transition regime of the pulse distribution. Breakdown voltage at different positions on the dielectric surface is suggested to be stochastic in nature. The simulated results based on the proposed distribution function agreed well with the experimental results and able to predict the regime of transition voltage. This model would be useful for the understanding of stochastic behaviors of DBD and the design of DBD device for effective operation and applications.

Description: We investigate the dispersion relation for low frequency electromagnetic waves propagating parallel to the ambient magnetic field, considering that the velocity distributions of ions and electrons can be either bi-Maxwellian of product bi-kappa distributions. The effect of the anisotropy and non-thermal features associated to the product-bi-kappa distributions on the firehose instability are numerically investigated. The general conclusion to be drawn from the results obtained is that the increase in non-thermal features which is consequence of the decrease of the κ indexes in the ion distribution contributes to increase the instability in magnitude and wave number range, in comparison with bi-Maxwellian distributions with similar temperature anisotropy, and that the increase of non-thermal features in the electron distribution contributes to the quenching of the instability, which is nevertheless driven by the anisotropy in the ion distribution. Significant differences between results obtained either considering product-bi-kappa distributions or bi-kappa distributions are also reported.

Description: The effects of shock dynamics on compressibility of indirect-drive ignition-scale surrogate implosions, CH shells filled with D 3 He gas, have been studied using charged-particle spectroscopy. Spectral measurements of D 3 He protons produced at the shock-bang time probe the shock dynamics and in-flight characteristics of an implosion. The proton shock yield is found to vary by over an order of magnitude. A simple model relates the observed yield to incipient hot-spot adiabat, suggesting that implosions with rapid radiation-power increase during the main drive pulse may have a 2× higher hot-spot adiabat, potentially reducing compressibility. A self-consistent 1-D implosion model was used to infer the areal density ( ρR ) and the shell center-of-mass radius ( R cm ) from the downshift of the shock-produced D 3 He protons. The observed ρR at shock-bang time is substantially higher for implosions, where the laser drive is on until near the compression bang time (“short-coast”), while longer-coasting implosions have lower ρR . This corresponds to a much larger temporal difference between the shock- and compression-bang time in the long-coast implosions (∼800 ps) than in the short-coast (∼400 ps); this will be verified with a future direct bang-time diagnostic. This model-inferred differential bang time contradicts radiation-hydrodynamic simulations, which predict constant 700–800 ps differential independent of coasting time; this result is potentially explained by uncertainties in modeling late-time ablation drive on the capsule. In an ignition experiment, an earlier shock-bang time resulting in an earlier onset of shell deceleration, potentially reducing compression and, thus, fuel ρR .

Description: In this work, the low frequency waves and gravitational (Jeans) instability of a homogeneous gyrotropic, magnetized, quantum plasma is investigated using the quantum magnetohydrodynamic and Chew-Goldberger-Low fluid models. An analytical dispersion relation for the considered system is obtained solving the linearized perturbations equations employing the Fourier transformation. The onset criterion of the “firehose” instability is retained in parallel propagation, which is unaffected due to the presence of quantum corrections. The gravitational mode modified by the quantum corrections is obtained separately along with the “firehose” mode. In perpendicular propagation, the quantum diffraction term is coupled with the Jeans and Alfven modes whereas in parallel propagation, the Alfven mode does not contribute to the dispersion characteristics as it leads to the “firehose” instability criterion in terms of quantum pressure anisotropy. The stabilizing influences of the quantum diffraction parameter and magnetic field on the growth rates of Jeans instability are examined. It is observed that the growth rate stabilizes much faster in transverse mode due to Alfven stabilization as compared to the longitudinal mode of propagation.

Description: Despite its affordability, the krypton Hall-effect thruster in applications always had problems in regard to performance. The reason for this degradation is studied from the perspective of the near-wall conductivity of electrons. Using the particle-in-cell method, the sheath oscillation characteristics and its effect on near-wall conduction are compared in the krypton and xenon Hall-effect thrusters both with wall material composed of BNSiO 2 . Comparing these two thrusters, the sheath in the krypton-plasma thruster will oscillate at low electron temperatures. The near-wall conduction current is only produced by collisions between electrons and wall, thereby causing a deficiency in the channel current. The sheath displays spatial oscillations only at high electron temperature; electrons are then reflected to produce the non-oscillation conduction current needed for the krypton-plasma thruster. However, it is accompanied with intensified oscillations.

Description: In toroidally confined plasmas, the Grad-Shafranov equation, in general, a non-linear partial differential equation, describes the hydromagnetic equilibrium of the system. This equation becomes linear when the kinetic pressure is proportional to the poloidal magnetic flux and the squared poloidal current is a quadratic function of it. In this work, the eigenvalue of the associated homogeneous equation is related with the safety factor on the magnetic axis, the plasma beta, and the Shafranov shift; then, the adjustable parameters of the particular solution are bounded through physical constrains. The poloidal magnetic flux becomes a linear superposition of independent solutions and its parameters are adjusted with a non-linear fitting algorithm. This method is used to find hydromagnetic equilibria with normal and reversed magnetic shear and defined values of the elongation, triangularity, aspect-ratio, and X-point(s). The resultant toroidal and poloidal beta, the safety factor at the 95% flux surface, and the plasma current are in agreement with usual experimental values for high beta discharges and the model can be used locally to describe reversed magnetic shear equilibria.

Description: The development of a round liquid jet under the influence of a confined coaxial flow of an immiscible liquid of comparable density (central to annular flow density ratio of 8:10) was investigated in the vicinity of the nozzle exit. Two flow regimes were considered; one where the annular flow is faster than the central jet, so the central liquid jet is accelerated and one where the annular flow is slower, so the central liquid jet is decelerated. The central jet was visualised by high speed photography. Three modes of jet development were identified and classified in terms of the Reynolds number, Re, of the central jet which was in the range of 525 〈 Re 〈 2725, a modified definition of the Weber number, We, which allows the distinction between accelerating and deceleration flows and was in the range of −22 〈 We 〈 67 and the annular to central Momentum Ratio, MR, of the two streams which was in the range of 3.6 〈 MR 〈 91. By processing the time resolved jet images using Proper Orthogonal Decomposition (POD), it was possible to reduce the description of jet morphology to a small number of spatial modes, which isolated the most significant morphologies of the jet development. In this way, the temporal and spatial characteristics of the instabilities on the interface were clearly identified which highlights the advantages of POD over direct observation of the images. Relationships between the flow parameters and the interfacial waves were established. The wavelength of the interfacial instability was found to depend on the velocity of the fastest moving stream, which is contrary to findings for fluids with large density differences.

Description: Interaction of a vortex ring impinging on multiple permeable screens orthogonal to the ring axis was studied to experimentally investigate the persistence and decay of vortical structures inside the screen array using digital particle image velocimetry in a refractive index matched environment. The permeable screens had porosities (open area ratios) of 83.8%, 69.0%, and 55.7% and were held by a transparent frame that allowed the screen spacing to be changed. Vortex rings were generated using a piston-cylinder mechanism at nominal jet Reynolds numbers of 1000, 2000, and 3000 with piston stroke length-to-diameter ratios of 2 and 3. The interaction of vortex rings with the porous medium showed a strong dependence of the overall flow evolution on the screen porosity, with a central flow being preserved and vortex ring-like structures (with smaller diameter than the primary vortex ring) being generated near the centerline. Due to the large rod size used in the screens, immediate reformation of the transmitted vortex ring with size comparable to the primary ring (as has been observed with thin screens) was not observed in most cases. Since the screens have lower complexity and high open area ratios, centerline vortex ring-like flow structures formed with comparable size to the screen pore size and penetrated through the screens. In the case of low porosity screens (55.7%) with large screen spacing, re-emergence of large scale (large separation), weak vortical structures/pairs (analogous to a transmitted vortex ring) was observed downstream of the first screen. Additional smaller scale vortical structures were generated by the interaction of the vortex ring with subsequent screens. The size distribution of the generated vortical structures were shown to be strongly affected by porosity, with smaller vortical structures playing a stronger role as porosity decreased. Finally, porosity significantly affected the decay of total energy, but the effect of screen spacing decreased as porosity decreased.

Description: Vortex cavitation forming in the leading-edge vortices of a delta wing was examined to determine how the individual cavitation bubbles incepted, grew, interacted with the underlying vortical flow and produced acoustic tones. The non-cavitating vortical flow over the delta wing was chosen to be similar to those previously reported in the literature. It was found that vortex breakdown was unaffected by the presence of incipient and developed vortex cavitation bubbles in the vortex core. While some cavitation bubbles incepted, grew, and collapsed relatively quickly, others reached an equilibrium position wherein the bubble tip was stationary in the laboratory frame at a particular location along the vortex axis. For a given attack angle, the equilibrium location moved upstream with a reduction in free stream cavitation number. It is shown that the existence of these stationary vortex bubbles is possible when there is a balance between the axial growth of the bubble along the vortex axis and the opposite motion of the axial jetting flow in the vortex core, and only a single equilibrium position is possible along the axially evolving vortex for a given free stream cavitation number. These transient and stationary vortex bubbles emit significant cavitation noise upon inception, growth, and collapse. The spectral content of the noise produced was expected to be related to the interaction of the bubble with the surrounding vortical flow in a manner similar to that reported in previous studies, where sustained tones were similar to the underlying vortex frequency. However, in the present study, the dominant frequency and higher harmonics of the tones occur at a higher frequency than that of the underlying vortex. Hence, it is likely that the highly elongated stationary bubbles have higher-order volume oscillations compared to the two-dimensional radial mode of the vortex cores of vortex cavitation bubbles with much smaller diameter-to-length ratios.

Description: Geodesic acoustic mode in anisotropic tokamak plasmas is theoretically analyzed by using double adiabatic model and gyro-kinetic equation. The bi-Maxwellian distribution function for guiding-center ions is assumed to obtain a self-consistent form, yielding pressures satisfying the magnetohydrodynamic (MHD) anisotropic equilibrium condition. The double adiabatic model gives the dispersion relation of geodesic acoustic mode (GAM), which agrees well with the one derived from gyro-kinetic equation. The GAM frequency increases with the ratio of pressures, p ⊥ / p ∥ , and the Landau damping rate is dramatically decreased by p ⊥ / p ∥ . MHD result shows a low-frequency zonal flow existing for all p ⊥ / p ∥ , while according to the kinetic dispersion relation, no low-frequency branch exists for p ⊥ / p ∥ ≳ 2.

Description: Nanoparticles are grown from the sputtering of a tungsten cathode in a direct current argon glow discharge. Laser light scattering of a vertical laser sheet going through the plasma reveals that the dust particle cloud is compressed and pushed towards the anode during the discharge. Scanning electron microscopy images of substrates exposed to the plasma for given durations show that dust particles are continuously falling down on the anode during the discharge. These observations are explained by the fact that the electrostatic force at the negative glow-anode sheath boundary cannot balance the ion drag, gravity, and thermophoresis forces for particles of more than a few tens of nanometres in diameter.

Description: Laminar flow over a periodic array of cylindrical surface roughness elements is simulated with an immersed boundary spectral method both to validate the method for subsequent studies and to examine how persistent streamwise vortices are introduced by a low Reynolds number roughness element. Direct comparisons are made with prior studies at a roughness-based Reynolds number Re k (= U ( k ) k / ν ) of 205 and a diameter to spanwise spacing ratio d / λ of 1/3. Downstream velocity contours match present and past experiments very well. The shear layer developed over the top of the roughness element produces the downstream velocity deficit. Upstream of the roughness element, the vortex topology is found to be consistent with juncture flow experiments, creating three cores along the recirculation line. Streamtraces stemming from these upstream cores, however, have unexpectedly little effect on the downstream flowfield as lateral divergence of the boundary layer quickly dissipates their vorticity. Long physical relaxation time of the recirculating wake behind the roughness remains a prominent issue for simulating this type of flowfield.

Description: We investigate the generation of intense pulsed focused ion beams at the 6 MeV level using an inductive voltage adder (IVA) pulsed-power generator, which employs a magnetically insulated transmission line (MITL). Such IVA machines typical run at an impedance of few tens of Ohms. Previous successful intense ion beam generation experiments have often featured an “axial” pinch-reflex ion diode (i.e., with an axial anode-cathode gap) and operated on a conventional Marx generator/water line driver with an impedance of a few Ohms and no need for an MITL. The goals of these experiments are to develop a pinch-reflex ion diode geometry that has an impedance to efficiently match to an IVA, produces a reasonably high ion current fraction, captures the vacuum electron current flowing forward in the MITL, and focuses the resulting ion beam to small spot size. A new “radial” pinch-reflex ion diode (i.e., with a radial anode-cathode gap) is found to best demonstrate these properties. Operation in both positive and negative polarities was undertaken, although the negative polarity experiments are emphasized. Particle-in-cell (PIC) simulations are consistent with experimental results indicating that, for diode impedances less than the self-limited impedance of the MITL, almost all of the forward-going IVA vacuum electron flow current is incorporated into the diode current. PIC results also provide understanding of the diode-impedance and ion-focusing properties of the diode. In addition, a substantial high-energy ion population is also identified propagating in the “reverse” direction, i.e., from the back side of the anode foil in the electron beam dump.

Description: Electron bulk heating during magnetic reconnection with symmetric inflow conditions is examined using kinetic particle-in-cell simulations. Inflowing plasma parameters are varied over a wide range of conditions, and the increase in electron temperature is measured in the exhaust well downstream of the x-line. The degree of electron heating is well correlated with the inflowing Alfvén speed c Ar based on the reconnecting magnetic field through the relation Δ T e = 0.033   m i   c A r 2 , where Δ 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 t o t = T i + T e or plasma β . 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 (Phan et al ., Geophys. Res. Lett. 40 , 4475, 2013), which also found that Δ T e was proportional to the inflowing Alfvén 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 simulation size is large enough to be directly relevant to reconnection in the Earth's magnetosphere, and the present findings may prove to be universal in nature with applications to the solar wind, the solar corona, and other astrophysical plasmas. 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 A r 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.

Description: Previous nonlinear gyrokinetic simulations of specific DIII-D L-mode cases have been found to significantly underpredict the ion heat transport and associated density and temperature fluctuation levels by up to almost one of order of magnitude in the outer-core domain, i.e., roughly in the last third of the minor radius. Since then, this so-called shortfall issue has been subject to various speculations on possible reasons and furthermore motivation for a number of dedicated comparisons for L-mode plasmas in comparable machines. However, only a rather limited number of simulations and gyrokinetic codes has been applied to the original scenario, thus calling for further dedicated investigations in order to broaden the scientific basis. The present work contributes along these lines by employing another well-established gyrokinetic code in a numerically and physically comprehensive manner. Contrary to the previous studies, only a mild underprediction is observed at the outer radial positions which can furthermore be overcome by varying the ion temperature gradient within the error bars associated with the experimental measurement. The significance and reliability of these simulations are demonstrated by benchmarks, numerical convergence tests, and furthermore by extensive validation studies. The latter involve cross-phase and cross-power spectra analyses of various fluctuating quantities and confirm a high degree of realism. The code discrepancies come as a surprise since the involved software packages had been benchmarked repeatedly and very successfully in the past. Further collaborative effort in identifying the underlying difference is hence required.

Description: We examined the plasma flow response to meso- and microscale magnetic dipoles by performing three-dimensional full particle-in-cell simulations. We particularly focused on the formation of a magnetosphere and its dependence on the intensity of the magnetic moment. The size of a magnetic dipole immersed in a plasma flow can be characterized by a distance L from the dipole center to the position where the pressure of the local magnetic field becomes equal to the dynamic pressure of the plasma flow under the magnetohydrodynamics (MHD) approximation. In this study, we are interested in a magnetic dipole whose L is smaller than the Larmor radius of ions r iL calculated with the unperturbed dipole field at the distance L from the center. In the simulation results, we confirmed the clear formation of a magnetosphere consisting of a magnetopause and a tail region in the density profile, although the spatial scale is much smaller than the MHD scale. One of the important findings in this study is that the spatial profiles of the plasma density as well as the current flows are remarkably affected by the finite Larmor radius effect of the plasma flow, which is different from the Earth's magnetosphere. The magnetopause found in the upstream region is located at a position much closer to the dipole center than L . In the equatorial plane, we also found an asymmetric density profile with respect to the plasma flow direction, which is caused by plasma gyration in the dipole field region. The ion current layers are created in the inner region of the dipole field, and the electron current also flows in the region beyond the ion current layer because ions with a large inertia can closely approach the dipole center. Unlike the ring current structure of the Earth's magnetosphere, the current layers in the microscale dipole fields are not circularly closed around the dipole center. Since the major current is caused by the particle gyrations, the current is independently determined to be in the direction of the electron and ion gyrations, which are the same in both the upstream and downstream regions. The present analysis on the formation of a magnetosphere in the regime of a microscale magnetic dipole is significant for understanding the solar wind response to the crustal magnetic anomalies on the Moon surface, such as were recently observed by spacecraft.

Description: The electrostatic potential profile and the particle densities of a simplified auroral double layer are found using a relaxation method to solve Poisson's equation in one dimension. The electron and ion distribution functions for the ionosphere and magnetosphere are specified at the boundaries, and the particle densities are found from a collisionless kinetic model. The ion distribution function includes the gravitational potential energy; hence, the unperturbed ionospheric plasma has a density gradient. The plasma potential at the upper boundary is given a large negative value to accelerate electrons downward. The solutions for a wide range of dimensionless parameters show that the double layer forms just above a critical altitude that occurs approximately where the ionospheric density has fallen to the magnetospheric density. Below this altitude, the ionospheric ions are gravitationally confined and have the expected scale height for quasineutral plasma in gravity.

Description: Solutal Marangoni instability (SMI) is investigated both in 2D and 3D using a combined Cahn-Hilliard and Navier-Stokes model in a finite system. Fe-Sn is chosen as a representative alloy system since the phase diagram reveals a region with a miscibility gap, where two liquid phases, namely, the Fe-rich phase L 1 and the Sn-rich phase L 2 , are in chemical equilibrium. In 3D, considering a perturbed liquid cylinder ( L 2 phase) with a length of λ and a radius of R 0 embedded in the middle of a simulation box of λ × H × H (length × width × height) surrounded by the phase L 1 , we find that the perturbation induced Marangoni flow is either clockwise or anti-clockwise depending on the mean curvature difference between the convex and concave regions which is affected by the ratio of λ/ R 0 . The critical ratio of λ/ R 0 for SMI is shown to be invariant for different Marangoni numbers as well as independent of the geometrical properties of the L 1 phase. In 2D, a perturbed liquid pipe with a length of λ and a radius of R 0 embedded in the middle of a simulation box of λ × H (length × height) is taken into account. Due to different curvature constitution, the critical ratio of λ/ R 0 for SMI depends on the height of the L 1 phase.

Description: Harmonic multiplying gyrotron traveling-wave amplifiers (gyro-TWAs) permit for magnetic field reduction and frequency multiplication. A high-order-mode harmonic multiplying gyro-TWA with large circuit dimensions and low ohmic loss can achieve a high average power. By amplifying a fundamental harmonic TE 01 drive wave, the second harmonic component of the beam current initiates a TE 02 wave to be amplified. Wall losses can suppress some competing modes because they act as an effective sink of the energy of the modes. However, such wall losses do not suppress all competing modes as the fields are contracted in the copper section in the gyro-TWA. An improved mode-selective circuit, using circular waveguides with the specified radii, can provide the rejection points within the frequency range to suppress the competing modes. The simulated results reveal that the mode-selective circuit can provide an attenuation of more than 10 dB to suppress the competing modes ( TE 21 , TE 51 , TE 22 , and TE 03 ). A G-band second harmonic multiplying gyro-TWA with the mode-selective circuit is predicted to yield a peak output power of 50 kW at 198.8 GHz, corresponding to a saturated gain of 55 dB at an interaction efficiency of 10%. The full width at half maximum bandwidth is 5 GHz.

Description: For several years, a promising Plasma Synthetic Jet actuator for high-speed flow control has been under development at ONERA. So far, its confined geometry and small space-time scales at play have prevented its full experimental characterization. Complementary accurate numerical simulations are then considered in this study in order to provide a complete aerothermodynamic description of the actuator. Two major obstacles have to be overcome with this approach: the modeling of the energy deposited by the electric arc and the accurate computation of the transient response of the cavity generating the pulsed jet. To solve the first problem, an Euler solver coupled with an electric circuit model was used to evaluate the energy deposition in the cavity. Such a coupling is performed by considering the electric field between the two electrodes. The second issue was then addressed by injecting these source terms in large Eddy simulations of the entire actuator. Aerodynamic results were finally compared with Schlieren visualizations. Using the proposed methodology, the temporal evolution of the jet front is remarkably well predicted.

Description: Hysteresis in amplitudes of the self-excited oscillations of the floating potential and discharge current is observed in an unmagnetized co-axial electrode-geometry DC glow discharge plasma system. The nonlinearities of these oscillations are studied using standard dynamical analysis tools. The characterization revealed the transition of low-amplitude high-frequency period-n oscillations to a large amplitude low frequency period-1 oscillations through a chaotic intermediate route. The transition of the low amplitude, high frequency period-n oscillations to chaotic type is observed to be linked to the dynamical change in the plasma system, i.e., after a negative differential resistance (NDR) region, whereas the transition from chaotic to period-1 is observed to be linked to a discharge current threshold.

Description: We have simulated the L-H transition on the EAST tokamak [Baonian Wan, EAST and HT-7 Teams, and International Collaborators, “Recent experiments in the EAST and HT-7 superconducting tokamaks,” Nucl. Fusion 49 , 104011 (2009)] using a predictive transport code where ion and electron temperatures, electron density, and poloidal and toroidal momenta are simulated self consistently. This is, as far as we know, the first theory based simulation of an L-H transition including the whole radius and not making any assumptions about where the barrier should be formed. Another remarkable feature is that we get H-mode gradients in agreement with the α – α d diagram of Rogers et al . [Phys. Rev. Lett. 81 , 4396 (1998)]. Then, the feedback loop emerging from the simulations means that the L-H power threshold increases with the temperature at the separatrix. This is a main feature of the C-mod experiments [Hubbard et al ., Phys. Plasmas 14 , 056109 (2007)]. This is also why the power threshold depends on the direction of the grad B drift in the scrape off layer and also why the power threshold increases with the magnetic field. A further significant general H-mode feature is that the density is much flatter in H-mode than in L-mode.

Description: Space-charge-limited (SCL) current can always be obtained from the blade surface of annular cathodes in foil-less diodes which are widely used in O-type relativistic high power microwave generators. However, there is little theoretical analysis regarding it due to the mathematical complexity, and almost all formulas about the SCL current in foil-less diodes are based on numerical simulation results. This paper performs an initial trial in calculation of the SCL current from annular cathodes theoretically under the ultra-relativistic assumption and the condition of infinitely large guiding magnetic field. The numerical calculation based on the theoretical research is coherent with the particle-in-cell (PIC) simulation result to some extent under a diode voltage of 850 kV. Despite that the theoretical research gives a much larger current than the PIC simulation (41.3 kA for the former and 9.7 kA for the latter), which is induced by the ultra-relativistic assumption in the theoretical research, they both show the basic characteristic of emission from annular cathodes in foil-less diodes, i.e., the emission enhancement at the cathode blade edges, especially at the outer edge. This characteristic is confirmed to some extent in our experimental research of cathode plasma photographing under the same diode voltage and a guiding magnetic field of 4 T.

Description: Snapshot and classical proper orthogonal decomposition (POD) are used to examine the large-scale, energetic motions in fully developed turbulent pipe flow at Re D = 47,000 and 93,000. The snapshot POD modes come in pairs, representing the same azimuthal mode number but with a simple phase shift. The first 10 snapshot POD modes, associated with the very large scale motions (VLSMs), contribute 43% of the average Reynolds shear stress, and for first 80 modes u ′ and v ′ are anti-correlated so that they all contribute to positive shear stress events. The attached motions are contained in the lower order modes, and detached motions do not appear until snapshot POD mode numbers ≥15. We find that snapshot POD can introduce mode mixing, which is avoided in classical POD. Classical POD also gives frequency information, confirming that the low order modes capture well the behavior of the very large scale motions.

Description: We derive two equations describing the coupling between electromagnetic and electrostatic oscillations in one-dimensional geometry in a magnetized cold and non-relativistic plasma. The nonlinear interaction between the wave modes is studied numerically. The effects of the external magnetic field strength and the initial electromagnetic polarization are of particular interest here. New results can, thus, be identified.

Description: We investigate changes in the intermittent sequence of bursts in the electrostatic turbulence due to imposed positive bias voltage applied to control the plasma radial electric field in Texas Helimak [K. W. Gentle and H. He, Plasma Sci. Technol. 10 , 284 (2008)]—a toroidal plasma device with a one-dimensional equilibrium, magnetic curvature, and shear. We identify the burst characteristics by analyzing ion saturation current fluctuations collected in a large set of Langmuir probes. The number of bursts increase with positive biasing, giving rise to a long tailed skewed turbulence probability distribution function. The burst shape does not change much with the applied bias voltage, while their vertical velocity increases monotonically. For high values of bias voltage, the bursts propagate mainly in the vertical direction which is perpendicular to the radial density gradient and the toroidal magnetic field. Moreover, in contrast with the bursts in tokamaks, the burst velocity agrees with the phase velocity of the overall turbulence in both vertical and radial directions. For a fixed bias voltage, the time interval between bursts and their amplitudes follows exponential distributions. Altogether, these burst characteristics indicate that their production can be modelled by a stochastic process.

Description: Simulations of doublet III-D, the joint European tokamak, and the tokamak fusion test reactor L-mode tokamak plasmas are carried out using the PTRANSP predictive integrated modeling code. The simulation and experimental temperature profiles are compared. The time evolved temperature profiles are computed utilizing the Multi-Mode anomalous transport model version 7.1 (MMM7.1) which includes transport associated with drift-resistive-inertial ballooning modes (the DRIBM model [T. Rafiq et al ., Phys. Plasmas 17 , 082511 (2010)]). The tokamak discharges considered involved a broad range of conditions including scans over gyroradius, ITER like current ramp-up, with and without neon impurity injection, collisionality, and low and high plasma current. The comparison of simulation and experimental temperature profiles for the discharges considered is shown for the radial range from the magnetic axis to the last closed flux surface. The regions where various modes in the Multi-Mode model contribute to transport are illustrated. In the simulations carried out using the MMM7.1 model it is found that: The drift-resistive-inertial ballooning modes contribute to the anomalous transport primarily near the edge of the plasma; transport associated with the ion temperature gradient and trapped electron modes contribute in the core region but decrease in the region of the plasma boundary; and neoclassical ion thermal transport contributes mainly near the center of the discharge.

Description: The determination of plasma equilibria from diagnostic information is a fundamental issue. V3FIT is a fully three-dimensional reconstruction code capable of solving the inverse problem using both magnetic and kinetic measurements. It uses VMEC as core equilibrium solver and supports both free- and fixed-boundary reconstruction approaches. In fixed-boundary mode VMEC does not use explicit information about currents in external coils, even though it has important effects on the shape of the safety factor profile. Indeed, the edge safety factor influences the reversal position in RFP plasmas, which then determines the position of the m =  0 island chain and the edge transport properties. In order to exploit such information a new virtual diagnostic has been developed, that thanks to Ampère's law relates the external current through the center of the torus to the circulation of the toroidal magnetic field on the outermost flux surface. The reconstructions that exploit the new diagnostic are indeed found to better interpret the experimental data with respect to edge physics.

Description: The structural stability of the laser/plasma interaction is discussed, for the case of a linearly polarized laser beam interacting with a solid at normal incidence. Using a semi-analytical cold fluid model, the dynamics of the electron plasma boundary (EPB), usually related to the high-order harmonic generation and laser absorption, are presented. While the well-known J  ×  B plasma oscillations at two times the laser frequency are recovered by the model, several other periodic in time stable solutions exist for exactly the same value of the physical parameters. This novel behavior highlights the importance of the laser pulse history among other factors. Some important features, such as the synchronization between the incident laser and the EPB oscillation, depend on the solution under consideration. A description of the possible types of stable oscillations in a parametric plane involving plasma density and laser amplitude is presented. The semi-analytical model is compared with particle-in-cell and semi-Lagrangian Vlasov simulations. They show that, among all the stable solutions, the plasma preferentially evolves to a state with the EPB oscillating twice faster than the laser. The effect of the plasma temperature and the existence of a ramp in the ion density profile are also discussed.

Description: Existence of the density limit in the helicon plasma sources is revisited. The low- and high-frequency regimes of a helicon plasma source operation are distinguished. In the low-frequency regime with ω 〈 ω c i ω c e , the density limit is deduced from the Golant-Stix criterion of the accessibility of the lower hybrid resonance. In the high-frequency case, ω 〉 ω c i ω c e , an appropriate limit is given by the Shamrai-Taranov criterion. Both these criteria are closely related to the phenomenon of the coalescence of the helicon wave with the Trivelpiece-Gould mode. We draw a conclusion that the derived density limits are not currently achieved in existing devices, perhaps, because of high energy cost of gas ionization.

Description: Energy transfer to damped modes in gyrokinetic ion temperature gradient driven turbulence is studied to understand the transfer dynamics and find scaling representations for the heavily populated mode space. Proper orthogonal and linear eigenmode decompositions are introduced and examined to assess whether modes are well-resolved and what scales they encompass. It is observed that damped modes across a range of inhomogeneous scales receive energy simultaneously and directly from the unstable mode, constituting a form of parallel transfer, distinct from the serial mode-to-mode transfer of the wavenumber cascade of hydrodynamic turbulence. Controlling for modes that are well resolved and labeling the modes of the linear decomposition in order of damping rate, energy transfer in the mode space satisfies an equipartition of the energy dissipation rate, leading to a simple rule for the distribution of energy in the space of damped modes. Energy dissipation rate equipartition is the form that the canonical nonlinear invariance of energy transfer assumes in a dissipation range with parallel rather than serial transfer.

Description: By combining the Cerenkov-type generator with the cascaded resonators, this paper proposes a Ka-band relativistic backward wave oscillator operating under the guide magnetic field 1.0 T with high power handling capability and high conversion efficiency. It is found that TM 02 can be selected as the operation mode in order to increase the power handling capability and provide sufficient coupling with the electron beam. In slow wave structure (SWS), ripples composed of semicircle on top of the rectangle enhance the wave-beam interaction and decrease the intensity of the electric field on the metallic surface. Taking advantage of the resonator cascades, the output power and the conversion efficiency are promoted greatly. The front cascaded resonators efficiently prevent the power generated in SWS from leaking into the diode region, and quicken the startup of the oscillation due to the premodulation of the beam. However, the post cascade slightly postpones the startup because of the further energy extraction from the electron beam. The numerical simulation shows that generation with power 514 MW and efficiency 41% is obtained under the diode voltage 520 kV and current 2.4 kA. And the microwave with the pure frequency spectrum of 29.35 GHz radiates in the pure TM 01 mode.

Description: The electrons in capacitively coupled plasmas (CCPs) absorb energy via ohmic heating due to electron-neutral collisions and stochastic heating due to momentum transfer from high voltage moving sheaths. We use Particle-in-Cell (PIC) simulations to explore these heating mechanisms and to compare the PIC results with available theories on ohmic and stochastic heating. The PIC results for ohmic heating show good agreement with the ohmic heating calculation of Lafleur et al . [Phys. Plasmas 20 , 124503 (2013)]. The PIC results for stochastic heating in low pressure CCPs with collisionless sheaths show good agreement with the stochastic heating model of Kaganovich et al . [IEEE Trans. Plasma Sci. 34 , 696 (2006)], which revises the hard wall asymptotic model of Lieberman [IEEE Trans. Plasma Sci. 16 , 638 (1988)] by taking current continuity and bulk oscillation into account.

Description: A 3-dimensional particle-in-cell/Monte Carlo collision simulation that is fully implemented on a graphics processing unit (GPU) is described and used to determine low-temperature plasma characteristics at high reduced electric field, E/n, in nitrogen gas. Details of implementation on the GPU using the NVIDIA Compute Unified Device Architecture framework are discussed with respect to efficient code execution. The software is capable of tracking around 10 × 10 6 particles with dynamic weighting and a total mesh size larger than 10 8 cells. Verification of the simulation is performed by comparing the electron energy distribution function and plasma transport parameters to known Boltzmann Equation (BE) solvers. Under the assumption of a uniform electric field and neglecting the build-up of positive ion space charge, the simulation agrees well with the BE solvers. The model is utilized to calculate plasma characteristics of a pulsed, parallel plate discharge. A photoionization model provides the simulation with additional electrons after the initial seeded electron density has drifted towards the anode. Comparison of the performance benefits between the GPU-implementation versus a CPU-implementation is considered, and a speed-up factor of 13 for a 3D relaxation Poisson solver is obtained. Furthermore, a factor 60 speed-up is realized for parallelization of the electron processes.

Description: The characterization of plasma bursts produced after the pinch phase in a plasma focus of hundreds of joules, using pulsed optical refractive techniques, is presented. A pulsed Nd-YAG laser at 532 nm and 8 ns FWHM pulse duration was used to obtain Schlieren images at different times of the plasma dynamics. The energy, interaction time with a target, and power flux of the plasma burst were assessed, providing useful information for the application of plasma focus devices for studying the effects of fusion-relevant pulses on material targets. In particular, it was found that damage factors on targets of the order of 10 4 (W/cm 2 )s 1/2 can be obtained with a small plasma focus operating at hundred joules.

Description: We explore the instabilities developed in a fluid in which viscosity depends on temperature. In particular, we consider a dependency that models a very viscous (and thus rather rigid) lithosphere over a convecting mantle. To this end, we study a 2D convection problem in which viscosity depends on temperature by abruptly changing its value by a factor of 400 within a narrow temperature gap. We conduct a study which combines bifurcation analysis and time-dependent simulations. Solutions such as limit cycles are found that are fundamentally related to the presence of symmetry. Spontaneous plate-like behaviors that rapidly evolve towards a stagnant lid regime emerge sporadically through abrupt bursts during these cycles. The plate-like evolution alternates motions towards either the right or the left, thereby introducing temporary asymmetries on the convecting styles. Further time-dependent regimes with stagnant and plate-like lids are found and described.

Description: We study the shock wave structure in a rarefied polyatomic gas based on a simplified model of extended thermodynamics in which the dissipation is due only to the dynamic pressure. In this case the differential system is very simple because it is a variant of Euler system with a new scalar equation for the dynamic pressure [T. Arima, S. Taniguchi, T. Ruggeri, and M. Sugiyama, Phys. Lett. A376, 2799–2803 (2012)]. It is shown that this theory is able to describe the three types of the shock wave structure observed in experiments: the nearly symmetric shock wave structure (Type A, small Mach number), the asymmetric structure (Type B, moderate Mach number), and the structure composed of thin and thick layers (Type C, large Mach number).

Description: In this paper, the scaling property of the inverse energy cascade and forward enstrophy cascade of the vorticity filed ω( x , y ) in two-dimensional (2D) turbulence is analyzed. This is accomplished by applying a Hilbert-based technique, namely Hilbert-Huang transform, to a vorticity field obtained from a 8192 2 grid-points direct numerical simulation of the 2D turbulence with a forcing scale k f = 100 and an Ekman friction. The measured joint probability density function p ( C , k ) of mode C i ( x ) of the vorticity ω and instantaneous wavenumber k ( x ) is separated by the forcing scale k f into two parts, which correspond to the inverse energy cascade and the forward enstrophy cascade. It is found that all conditional probability density function p ( C | k ) at given wavenumber k has an exponential tail. In the inverse energy cascade, the shape of p ( C | k ) does collapse with each other, indicating a nonintermittent cascade. The measured scaling exponent ζ ω I ( q ) is linear with the statistical order q , i.e., ζ ω I ( q ) = − q / 3 , confirming the nonintermittent cascade process. In the forward enstrophy cascade, the core part of p ( C | k ) is changing with wavenumber k , indicating an intermittent forward cascade. The measured scaling exponent ζ ω F ( q ) is nonlinear with q and can be described very well by a log-Poisson fitting: ζ ω F ( q ) = 1 3 q + 0.45 1 − 0 . 43 q . However, the extracted vorticity scaling exponents ζ ω ( q ) for both inverse energy cascade and forward enstrophy cascade are not consistent with Kraichnan's theory prediction. New theory for the vorticity field in 2D turbulence is required to interpret the observed scaling behavior.

Description: We investigate theoretically the influence of static plasma screening on relativistic spin-orbit interaction-induced fine-structure splitting of the D-line doublet arising from the transitions 3p 1/2 –3s 1/2 and 3p 3/2 –3s 1/2 of the valence electron of a sodium atom embedded in a model plasma environment. The many-electron atomic problem is formulated first as an effective one-electron problem in which the interaction between the optically active valence electron and the atomic ion core is represented by an accurate parametric model potential including core-polarization correction, and then the plasma effect on the atomic system is simulated by the Debye-screening model for the valence-core interaction. It is observed that the magnitude of spin-orbit energy shift reduces for both the upper component 3p 3/2 and the lower component 3p 1/2 with increasing plasma screening strength, thereby reducing the spin-orbit energy separation between these two components as the screening becomes stronger. As a consequence, the magnitude of fine-structure splitting between the D 1 and D 2 line energies of sodium drops significantly with stronger plasma screening. The optical (absorption) oscillator strength for 3s → 3p transition is seen to reduce with stronger screening and this leads to a screening-induced gradual suppression of the 3p → 3s spontaneous decay rate.

Description: Here, we show that defocusing of the laser in the interaction of a picosecond duration, 1.053  μ m wavelength, high energy pulse with a cone-wire target does not significantly affect the laser energy coupling efficiency, but does result in a drop in the fast electron effective temperature. This may be beneficial for fast ignition, since not only were more electrons with lower energies seen in the experiment but also the lower prepulse intensity will reduce the amount of preplasma present on arrival of the main pulse, reducing the distance the hot electrons have to travel. We used the Vulcan Petawatt Laser at the Rutherford Appleton Laboratory and gold cone targets with approximately 1 mm long, 40  μ m diameter copper wires attached to their tip. Diagnostics included a quartz crystal imager, a pair of highly oriented pyrolytic graphite crystal spectrometers and a calibrated CCD operating in the single photon counting regime, all of which looked at the copper K α emission from the wire. A short pulse optical probe, delayed 400 ps relative to the main pulse was employed to diagnose the extent of plasma expansion around the wire. A ray-tracing code modeled the change in intensity on the interior surface of the cone with laser defocusing. Using a model for the wire copper K α emission coupled to a hybrid Vlasov-Fokker-Planck code, we ran a series of simulations, holding the total energy in electrons constant whilst varying the electron temperature, which support the experimental conclusions.

Description: In this paper, a computational simulation of the sheath formation during the streamer-surface interaction at atmospheric pressure is presented. A two-dimensional fluid model of a point-to-plane configuration is applied to investigate the evolution of the discharge in the vicinity of cathode plane. The effects of the surfaces on the properties of streamer have been studied for three cases, i.e., conductive surface with secondary electron emission (SEE), conductive surface without SEE, and dielectric surface. In all cases, we found that the axial propagation velocity of the streamer front decreases as the streamer arrives at the boundary of the cathode sheath. And the simulation results showed that the properties of the surface have a significant effect on the streamer. Besides the influences, the secondary emission coefficient and the relative permittivity on the streamer-surface interactions are also studied.

Description: It was recently shown that a promising way to accelerate protons in the forward direction to high energies is to use under-dense or near-critical density targets instead of solids. Simulations have revealed that the acceleration process depends on the density gradients of the plasma target. Indeed, under certain conditions, the most energetic protons are predicted to be accelerated by a collisionless shock mechanism that significantly increases their energy. We report here the results of a recent experiment dedicated to the study of longitudinal ion acceleration in partially exploded foils using a high intensity (∼5 × 10 18  W/cm 2 ) picosecond laser pulse. We show that protons accelerated using targets having moderate front and rear plasma gradients (up to ∼8  μ m gradient length) exhibit similar maximum proton energy and number compared to proton beams that are produced, in similar laser conditions, from solid targets, in the well-known target normal sheath acceleration regime. Particle-In-Cell simulations, performed in the same conditions as the experiment and consistent with the measurements, allow laying a path for further improvement of this acceleration scheme.

Description: We report experimental observations of the controlled deformation of a dielectric liquid jet subjected to a local high-voltage electrostatic field in the direction normal to the jet. The jet deforms to the shape of an elliptic cylinder upon application of a normal electrostatic field. As the applied electric field strength is increased, the elliptic cylindrical jet deforms permanently into a flat sheet, and eventually breaks-up into droplets. We interpret this observation—the stretch of the jet is in the normal direction to the applied electric field—qualitatively using the Taylor-Melcher leaky dielectric theory, and develop a simple scaling model that predicts the critical electric field strength for the jet-to-sheet transition. Our model shows a good agreement with experimental results, and has a form that is consistent with the classical drop deformation criterion in the Taylor-Melcher theory. Finally, we statistically analyze the resultant droplets from sheet breakup, and find that increasing the applied electric field strength improves droplet uniformity and reduces droplet size.

Description: Normal impingement of a single droplet on a thin liquid film is investigated numerically solving the axisymmetric Navier-Stokes equations. Gravity and viscosity are taken into account whereas compressibility effects are neglected. Two phases are tracked by means of volume of fluid method and adaptive mesh refinement is used to increase accuracy of the interface. Numerical results are validated both qualitatively and quantitatively using experimental measurements. Effects of gas density, gas viscosity, and film thickness on the crown behavior are studied. Influence of droplet deviation from spherical shape on the crown behavior is investigated. It is shown that increasing the gas density leads to reduction of crown radius evolution rate, while gas viscosity does not affect the rate of crown radius evolution. Development rate of crown height decreases by increasing the gas density. Reynolds number and splashing regime can change the effect of gas viscosity on the crown height evolution. Deviation of droplet from sphere can change behavior of crown completely as result of change in droplet mass center position. Difference between numerical results and experimental ones is justified using different droplet shapes.

Description: The heat flux patterns measured in low-collisionality DIII-D H-mode plasmas strongly deviate from simultaneously measured CII emission patterns, used as indicator of particle flux, during applied resonant magnetic perturbations. While the CII emission clearly shows typical striations, which are similar to magnetic footprint patterns obtained from vacuum field line tracing, the heat flux is usually dominated by one large peak at the strike point position. The vacuum approximation, which only considers applied magnetic fields and neglects plasma response and plasma effects, cannot explain the shape of the observed heat flux pattern. One possible explanation is the effect of particle drifts. This is included in the field line equations and the results are discussed with reference to the measurement. Electrons and ions show different drift motions at thermal energy levels in a guiding center approximation. While electrons hardly deviate from the field lines, ions can drift several centimetres away from field line flux surfaces. A model is presented in which an ion heat flux, based on the ion drift motion from various kinetic energies as they contribute to a thermal Maxwellian distribution, is calculated. The simulated heat flux is directly compared to measurements with a varying edge safety factor q 95 . This analysis provides evidence for the dominate effect of high-energy ions in carrying heat from the plasma inside the separatrix to the target. High-energy ions are deposited close to the unperturbed strike line, while low-energy ions can travel into the striated magnetic topology.

Description: Two fundamental challenging problems of laboratory and astrophysical plasmas are the understanding of the relaxation of a collisionless plasmas with nearly isotropic velocity distribution functions and the resultant state of nearly equipartition energy density with electromagnetic plasma turbulence. Here, we present the results of a study which shows the role that higher-order-modes play in limiting the electromagnetic whistler-like fluctuations in a thermal and non-thermal plasma. Our main results show that for a thermal plasma the magnetic fluctuations are confined by regions that are bounded by the least-damped higher order modes. We further show that the zone where the whistler-cyclotron normal modes merges the electromagnetic fluctuations shifts to longer wavelengths as the β e increases. This merging zone has been interpreted as the beginning of the region where the whistler-cyclotron waves losses their identity and become heavily damped while merging with the fluctuations. Our results further indicate that in the case of nonthermal plasmas, the higher-order modes do not confine the fluctuations due to the effective higher-temperature effects and the excess of suprathermal plasma particles. The analysis presented here considers the second-order theory of fluctuations and the dispersion relation of weakly transverse fluctuations, with wave vectors parallel to the uniform background magnetic field, in a finite temperature isotropic bi-Maxwellian and Tsallis-kappa-like magnetized electron–proton plasma. Our results indicate that the spontaneously emitted electromagnetic fluctuations are in fact enhanced over these quasi modes suggesting that such modes play an important role in the emission and absorption of electromagnetic fluctuations in thermal or quasi-thermal plasmas.

Description: We present an analytical study of peak mode isotachophoresis (ITP), and provide closed form solutions for sample distribution and electric field, as well as for leading-, trailing-, and counter-ion concentration profiles. Importantly, the solution we present is valid not only for the case of fully ionized species, but also for systems of weak electrolytes which better represent real buffer systems and for multivalent analytes such as proteins and DNA. The model reveals two major scales which govern the electric field and buffer distributions, and an additional length scale governing analyte distribution. Using well-controlled experiments, and numerical simulations, we verify and validate the model and highlight its key merits as well as its limitations. We demonstrate the use of the model for determining the peak concentration of focused sample based on known buffer and analyte properties, and show it differs significantly from commonly used approximations based on the interface width alone. We further apply our model for studying reactions between multiple species having different effective mobilities yet co-focused at a single ITP interface. We find a closed form expression for an effective-on rate which depends on reactants distributions, and derive the conditions for optimizing such reactions. Interestingly, the model reveals that maximum reaction rate is not necessarily obtained when the concentration profiles of the reacting species perfectly overlap. In addition to the exact solutions, we derive throughout several closed form engineering approximations which are based on elementary functions and are simple to implement, yet maintain the interplay between the important scales. Both the exact and approximate solutions provide insight into sample focusing and can be used to design and optimize ITP-based assays.

Description: Theoretical refinements to an existing model for the loss of ions by drifting across the last closed flux surface are presented.

Description: A novel waveform modified from the standard-sinusoidal function is adopted to enhance the virtual aeroshaping effect of the synthetic jets positioned at the front stagnation point of a circular cylinder. The waveform is characterized by a control parameter, namely, the suction duty cycle factor k , which is the ratio of the time duration of the suction cycle to that of the blowing cycle. The strength of the synthetic jet vortex pair is enhanced by increasing the suction duty cycle factor. The periodic closed envelope forms upstream of the circular cylinder for k ≤ 1.00, while the quasi-steady open envelope forms for k ≥ 2.00, acting the virtual aeroshaping effect. As a result, both the statistical characteristics and the vortex dynamics of the near-wake flow field change with the suction duty cycle factor. The recirculation region downstream of the circular cylinder becomes smaller or even disappears, and thus, the drag coefficient over the circular cylinder is reduced by increasing the suction duty cycle factor to k ≥ 1.00. The statistical mean and fluctuating velocities show corresponding changes in the near wake with the different wake patterns. For k ≤ 0.50, the wake vortex shows the antisymmetric shedding mode which is similar with the natural case. For 1.00 ≤ k ≤ 2.00, the wake vortex shows the bistable state mode, where vortex sheds with symmetric or antisymmetric mode; the antisymmetric shedding mode dominates the global flow field for k = 1.00, while it is the symmetric shedding mode that dominates the flow field for k = 2.00. For k = 4.00, it shows the antisymmetric shedding mode with a shorter vortex formation length than the natural case. The above findings indicate that the virtual aeroshaping effect of the synthetic jets can be enhanced by increasing the suction duty cycle factor so as to increase the momentum coefficient while keeping other control parameters unchanged, providing us another way for effective flow control.

Description: The dynamics of vapor-liquid interface are important because interfacial instability determines bubble growth, detachment frequency, waiting time, shape of bubbles, and the interrelationship between bubble formation sites. In this study, a detailed numerical simulation has been performed to understand the transition in bubble release pattern and multimode bubble formation in saturated pool boiling. The interfaces drop down alternatively at the nodes and antinodes of the wavelengths dictated by Rayleigh-Taylor instability and Taylor-Helmholtz instability. Due to higher degrees of superheat, vapor jets emanate from nodes and antinodes. An attempt has been made to predict the maximum and minimum heat fluxes during saturated pool boiling.

Description: We present a combination of experiment, theory, and modelling on laminar mixing at large Péclet number. The flow is produced by oscillating electromagnetic forces in a thin electrolytic fluid layer, leading to oscillating dipoles, quadrupoles, octopoles, and disordered flows. The numerical simulations are based on the Diffusive Strip Method (DSM) which was recently introduced (P. Meunier and E. Villermaux, “The diffusive strip method for scalar mixing in two-dimensions,” J. Fluid Mech.662, 134–172 (2010)) to solve the advection-diffusion problem by combining Lagrangian techniques and theoretical modelling of the diffusion. Numerical simulations obtained with the DSM are in reasonable agreement with quantitative dye visualization experiments of the scalar fields. A theoretical model based on log-normal Probability Density Functions (PDFs) of stretching factors, characteristic of homogeneous turbulence in the Batchelor regime, allows to predict the PDFs of scalar in agreement with numerical and experimental results. This model also indicates that the PDFs of scalar are asymptotically close to log-normal at late stages, except for the large concentration levels which correspond to low stretching factors.

Description: Janus droplets are compound droplets that consist of two adhering drops of different fluids that are suspended in a third fluid. We use the Shan-Chen lattice Boltzmann method for multicomponent mixtures to simulate Janus droplets at rest and in shear. In this simulation model, interfacial tensions are not known a priori from the model parameters and must be determined using numerical experiments. We show that interfacial tensions obtained with the Young-Laplace law are consistent with those measured from the equilibrium geometry. The regimes of adhering, separated, and engulfing droplets were explored. Two different adhesion geometries were considered for two-dimensional simulations of Janus droplets in shear. The first geometry resembles two adhering circles with small overlap. In the second geometry, the two halves are semicircular. For both geometries, the rotation rate of the droplet depends on its orientation. The width of the periodic simulation domain also affects the rotation rate of both droplet types up to an aspect ratio of 6:1 (width:height). While the droplets with the first geometry oscillated about the middle of the domain, the droplets of the second geometry did not translate while rotating. A four-pole vortex structure inside droplets of the second geometry was found. These simulations of single Janus droplets reveal complex behaviour that implies a rich range of possibilities for the rheology of Janus emulsions.

Description: Two issues in the cavity design for a Megawatt-class, 240 GHz gyrotron are addressed. Those are first, the effect of a misaligned electron beam on the gyrotron efficiency and second, a possible azimuthal instability of the gyrotron. The aforementioned effects are important for any gyrotron operation, but could be more critical in the operation of Megawatt-class gyrotrons at frequencies above 200 GHz, which will be the anticipated requirement of DEMO. The target is to provide some basic trends to be considered during the refinement and optimization of the design. Self-consistent calculations are the base for simulations wherever possible. However, in cases for which self-consistent models were not available, fixed-field results are presented. In those cases, the conservative nature of the results should be kept in mind.

Description: The connection between impurity ion heating and other physical processes in the plasma is evaluated by studying variations in the amount of ion heating at reconnection events in the Madison Symmetric Torus (MST). Correlation of the change in ion temperature with individual tearing mode amplitudes indicates that the edge-resonant modes are better predictors for the amount of global ion heating than the core-resonant modes. There is also a strong correlation between ion heating and current profile relaxation. Simultaneous measurements of the ion temperature at different toroidal locations reveal, for the first time, a toroidal asymmetry to the ion heating in MST. These results present challenges for existing heating theories and suggest a stronger connection between edge-resonant tearing modes, current profile relaxation, and ion heating than has been previously thought.

Description: This paper provides 2D comparative study of results obtained using laminar and turbulent flow model for RF (radio frequency) Inductively Coupled Plasma (ICP) torch. The study was done for the RF-ICP torch operating at 50 kW DC power and 3 MHz frequency located at BARC. The numerical modeling for this RF-ICP torch is done using ANSYS software with the developed User Defined Function. A comparative study is done between laminar and turbulent flow model to investigate how temperature and flow fields change when using different operating conditions such as (a) swirl and no swirl velocity for sheath gas flow rate, (b) variation in sheath gas flow rate, and (c) variation in plasma gas flow rate. These studies will be useful for different material processing applications.

Description: Effects of pressure anisotropy on magnetospheric magnetohydrodynamics equilibrium are studied analytically and numerically, where the plasma is confined by only poloidal magnetic field generated by an internal ring current. The plasma current due to finite pressure can be divided into two components; one remains at isotropic pressure and the other arises from pressure anisotropy. When p ⊥ , the pressure perpendicular to the magnetic field, is larger than p ∥ , the pressure parallel to the magnetic field, those two components of plasma current tend to cancel each other to reduce the total amount of plasma current. Equilibrium beta limit is also examined, where the beta is a ratio of the plasma pressure to the magnetic pressure. The equilibrium beta limit decreases as the pressure anisotropy becomes strong. The beta value is strictly limited by ellipticity of the equilibrium equation when p ∥ 〉 p ⊥ . On the other hand, when p ⊥ 〉 p ∥ , although the tendency of the beta limit agrees with the ellipticity condition of the equilibrium equation, equilibria with a hyperbolic region can be obtained by iterative procedure with practically reasonable convergence criteria.

Description: Understanding the physics of water evaporation from saline porous media is important in many natural and engineering applications such as durability of building materials and preservation of monuments, water quality, and mineral-fluid interactions. We applied synchrotron x-ray micro-tomography to investigate the pore-scale dynamics of dissolved salt distribution in a three dimensional drying saline porous media using a cylindrical plastic column (15 mm in height and 8 mm in diameter) packed with sand particles saturated with CaI 2 solution (5% concentration by mass) with a spatial and temporal resolution of 12 μ m and 30 min, respectively. Every time the drying sand column was set to be imaged, two different images were recorded using distinct synchrotron x-rays energies immediately above and below the K-edge value of Iodine. Taking the difference between pixel gray values enabled us to delineate the spatial and temporal distribution of CaI 2 concentration at pore scale. Results indicate that during early stages of evaporation, air preferentially invades large pores at the surface while finer pores remain saturated and connected to the wet zone at bottom via capillary-induced liquid flow acting as evaporating spots. Consequently, the salt concentration increases preferentially in finer pores where evaporation occurs. Higher salt concentration was observed close to the evaporating surface indicating a convection-driven process. The obtained salt profiles were used to evaluate the numerical solution of the convection-diffusion equation (CDE). Results show that the macro-scale CDE could capture the overall trend of the measured salt profiles but fail to produce the exact slope of the profiles. Our results shed new insight on the physics of salt transport and its complex dynamics in drying porous media and establish synchrotron x-ray tomography as an effective tool to investigate the dynamics of salt transport in porous media at high spatial and temporal resolution.

Description: Experimental studies by Poehlmann et al. [Phys. Plasmas 17 (12), 123508 (2010)] on a coaxial electrode magnetohydrodynamic (MHD) plasma accelerator have revealed two modes of operation. A deflagration or stationary mode is observed for lower power settings, while higher input power leads to a detonation or snowplow mode. A numerical modeling study of a coaxial plasma accelerator using the non-ideal MHD equations is presented. The effect of plasma conductivity on the axial distribution of radial current is studied and found to agree well with experiments. Lower conductivities lead to the formation of a high current density, stationary region close to the inlet/breech, which is a characteristic of the deflagration mode, while a propagating current sheet like feature is observed at higher conductivities, similar to the detonation mode. Results confirm that plasma resistivity, which determines magnetic field diffusion effects, is fundamentally responsible for the two modes.

Description: In this work, we re-examine collisional absorption of 800 nm wavelength laser pulses in under-dense plasma. For a given temperature and density of the plasma, most of the conventional models of the electron-ion collision frequency ν ei , with a Coulomb logarithm independent of the electron-ponderomotive velocity, show that ν ei and the corresponding fractional laser absorption α remain almost constant (or decrease slowly) up to a value I c of the peak intensity I 0 of the laser pulse, and then ν ei and α decrease as ≈ I 0 − 3 / 2 when I 0 is increased beyond I c . On the contrary, below some temperature ( ≲ 10   eV ) and density, with a total-velocity (thermal velocity plus the ponderomotive velocity) dependent Coulomb logarithm, we find that ν ei and α grow hand in hand up to a maximum value around I c followed by the conventional I 0 − 3 / 2 decrease when I 0 〉 I c . Such a non-conventional anomalous variation of α with I 0 was observed in some earlier experiments, but no explanation has been given so far. The modified Coulomb logarithm considered in this work may be responsible for those experimental observations. With increasing temperature and density, the anomalous behavior is found to disappear even with the modified Coulomb logarithm, and the variation of ν ei and α with I 0 approach to the conventional scenario.

Description: Since ignition target design with layered deuterium and triterium ice had been proposed several decades ago, much effort was devoted to fabricate and implode cryogenic targets. Until recently, direct-drive cryogenic target implosion experiment was carried out on SGIII prototype laser facility. The target consisted of a plastic capsule supported by fill tube. Cryogenic helium gas was used to cool the capsule to a few degrees below the deuterium triple point. The resulting deuterium ice layer was characterized by optical shadowgraph and smoothed by applied temperature gradient. Eight laser beams with total energy of 7 kJ were used to directly drive the implosion. On the path of laser light to the capsule, there were 500 nm sealing film and helium gas of mm length. X-ray pinhole images were analyzed to confirm that the sealing film, and helium gas had little effect on aiming accuracy but caused some loss of laser energy especially when condensation on the sealing film was observed.

Description: The magnetoacoustic cyclotron instability (MCI) probably underlies observations of ion cyclotron emission (ICE) from energetic ion populations in tokamak plasmas, including fusion-born alpha-particles in JET and TFTR [Dendy et al. , Nucl. Fusion 35 , 1733 (1995)]. ICE is a potential diagnostic for lost alpha-particles in ITER; furthermore, the MCI is representative of a class of collective instabilities, which may result in the partial channelling of the free energy of energetic ions into radiation, and away from collisional heating of the plasma. Deep understanding of the MCI is thus of substantial practical interest for fusion, and the hybrid approximation for the plasma, where ions are treated as particles and electrons as a neutralising massless fluid, offers an attractive way forward. The hybrid simulations presented here access MCI physics that arises on timescales longer than can be addressed by fully kinetic particle-in-cell simulations and by analytical linear theory, which the present simulations largely corroborate. Our results go further than previous studies by entering into the nonlinear stage of the MCI, which shows novel features. These include stronger drive at low cyclotron harmonics, the re-energisation of the alpha-particle population, self-modulation of the phase shift between the electrostatic and electromagnetic components, and coupling between low and high frequency modes of the excited electromagnetic field.

Description: The behaviour of low Reynolds number, non-Boussinesq fountains from four different nozzle geometries (one circular and three rectangular nozzles) are studied. High speed laser schlieren imaging is used to study the fountain behaviour (frequency and penetration height). Bi-orthogonal decomposition and dynamic mode decomposition (DMD) are used to understand the unsteady characteristics of fountains. The flow regimes of fountains are classified as steady, flapping, and flapping-bobbing type. The DMD technique successfully separates the bobbing oscillation from the combined flapping-bobbing oscillation of the fountain. The frequency of the flapping oscillation, and the frequency of the bobbing oscillation in the flapping-bobbing regime scales as St h Fr h = C 1 and S t h F r h 2 = C 2 , respectively, where the characteristic length scale is the smallest dimension ( h ) of the nozzle. The mean steady state penetration heights ( Z s / h ) of “forced” low Reynolds number non-Boussinesq fountains are independent of nozzle shape (circular and rectangular), and scales linearly with the Froude number.

Description: The large-scale properties of self-similar unstably stratified homogeneous (USH) turbulence are investigated using an eddy-damped quasi-normal markovianized approximation of the nonlinear term. This analysis shows that a special role is played by the wave vectors contained in the equatorial plane, i.e., the plane perpendicular to gravity. It is indeed in this plane that turbulent spectra reach their maxima and evolve linearly from their initial condition when their initial infrared exponent is smaller than 4. At other angles, this property is not satisfied and turbulent spectra eventually undergo an evolution dominated by nonlinear backscattering processes. The self-similar evolution of USH turbulence is also shown to be related to the properties of large scales. In particular, the asymptotic growth rate of the mixing length depends on the initial infrared exponent in the equatorial plane. Besides, the self-similar asymptotic values of the concentration and velocity correlations also depend on the properties of large scales. This allows to derive relations between the correlations and the growth rate parameter.

Description: Secondary flow cells are commonly observed in straight laboratory channels, where they are often associated with duct corners. Here, we present velocity measurements acquired with an acoustic Doppler current profiler in a straight reach of the Seine river (France). We show that a remarkably regular series of stationary flow cells spans across the entire channel. They are arranged in pairs of counter-rotating vortices aligned with the primary flow. Their existence away from the river banks contradicts the usual interpretation of these secondary flow structures, which invokes the influence of boundaries. Based on these measurements, we use a depth-averaged model to evaluate the momentum transfer by these structures, and find that it is comparable with the classical turbulent transfer.

Description: We propose a spherical hohlraum with octahedral six laser entrance holes at a specific hohlraum-to-capsule radius ratio of 5.14 for inertial fusion study, which has robust high symmetry during the capsule implosion and low backscatter without supplementary technology. To produce an ignition radiation pulse of 300 eV, it needs 1.5 MJ absorbed laser energy in such a golden octahedral hohlraum, about 30% more than a traditional cylinder. Nevertheless, it is worth for a high symmetry and low backscatter. The proposed octahedral hohlraum is also flexible and can be applicable to diverse inertial fusion drive approaches.

Description: Analysis of fluxes across the turbulent/non-turbulent interface (TNTI) of turbulent boundary layers is performed using data from two-dimensional particle image velocimetry (PIV) obtained at high Reynolds numbers. The interface is identified with an iso-surface of kinetic energy, and the rate of change of total kinetic energy ( K ) inside a control volume with the TNTI as a bounding surface is investigated. Features of the growth of the turbulent region into the non-turbulent region by molecular diffusion of K , viscous nibbling, are examined in detail, focussing on correlations between interface orientation, viscous stress tensor elements, and local fluid velocity. At the level of the ensemble (Reynolds) averaged Navier-Stokes equations (RANS), the total kinetic energy K is shown to evolve predominantly due to the turbulent advective fluxes occurring through an average surface which differs considerably from the local, corrugated, sharp interface. The analysis is generalized to a hierarchy of length-scales by spatial filtering of the data as used commonly in Large-Eddy-Simulation (LES) analysis. For the same overall entrainment rate of total kinetic energy, the theoretical analysis shows that the sum of resolved viscous and subgrid-scale advective flux must be independent of scale. Within the experimental limitations of the PIV data, the results agree with these trends, namely that as the filter scale increases, the viscous resolved fluxes decrease while the subgrid-scale advective fluxes increase and tend towards the RANS values at large filter sizes. However, a definitive conclusion can only be made with fully resolved three-dimensional data, over and beyond the large dynamic spatial range presented here. The qualitative trends from the measurement results provide evidence that large-scale transport due to the energy-containing eddies determines the overall rate of entrainment, while viscous effects at the smallest scales provide the physical mechanism ultimately responsible for entrainment. Data spanning over a decade in Reynolds number suggest that the fluxes (or the entrainment velocity) scale with the friction velocity (or equivalently the local turbulent fluctuating velocity), whereas Taylor microscale and boundary-layer thickness are the appropriate length scales at small and large filter sizes, respectively.

Description: We present an analytic study of momentum transport of tokamak plasmas in the vicinity of minimum safety factor ( q ) position in reversed magnetic shear configuration. Slab ion temperature gradient modes with an equilibrium flow profile are considered in this study. Quasi-linear calculations of momentum flux clearly show the novel effects of q -curvature on the generation of intrinsic rotation and mean poloidal flow without invoking reflectional symmetry breaking of parallel wavenumber ( k ∥ ). This q -curvature effect originates from the inherent asymmetry in k ∥ populations with respect to a rational surface due to the quadratic proportionality of k ∥ when q -curvature is taken into account. Discussions are made of possible implications of q -curvature induced plasma flows on internal transport barrier formation in reversed shear tokamaks.

Description: In a recent direct numerical simulation (DNS) study [P. K. Yeung and K. R. Sreenivasan, “ Spectrum of passive scalars of high molecular diffusivity in turbulent mixing,” J. Fluid Mech.716, R14 (2013)] with Schmidt number as low as 1/2048, we verified the essential physical content of the theory of Batchelor, Howells, and Townsend [“Small-scale variation of convected quantities like temperature in turbulent fluid. 2. The case of large conductivity,” J. Fluid Mech.5, 134 (1959)] for turbulent passive scalar fields with very strong diffusivity, decaying in the absence of any production mechanism. In particular, we confirmed the existence of the −17/3 power of the scalar spectral density in the so-called inertial-diffusive range. In the present paper, we consider the DNS of the same problem, but in the presence of a uniform mean gradient, which leads to the production of scalar fluctuations at (primarily) the large scales. For the parameters of the simulations, the presence of the mean gradient alters the physics of mixing fundamentally at low Peclet numbers. While the spectrum still follows a −17/3 power law in the inertial-diffusive range, the pre-factor is non-universal and depends on the magnitude of the mean scalar gradient. Spectral transfer is greatly reduced in comparison with those for moderately and weakly diffusive scalars, leading to several distinctive features such as the absence of dissipative anomaly and a new balance of terms in the spectral transfer equation for the scalar variance, differing from the case of zero gradient. We use the DNS results to present an alternative explanation for the observed scaling behavior, and discuss a few spectral characteristics in detail.