Search

Based on the magnetic field configuration of the Chinese First Quasi-axisymmetric Stellarator (CFQS) device, three types of orbits such as the passing orbit, blocked trapped orbit, and localized trapped orbit are simulated using the Boris algorithm. Also, the orbital topology, orbit transition, and loss characteristics of these particles under different initial conditions are studied in this paper. It is found that there exists a transition from blocked trapped orbits to localized trapped orbits due to small continuous helical ripples. This phenomenon is analyzed as follows: we define the angle between the particle drift velocity ( v → D ) and the radial direction ( ρ ^ ) as β. If the blocked-localized transition appears in the region where β > 9 0 ○ , the localized particles will return to the blocked particles and be constrained. However, if the blocked-localized transition happens in the region where β < 9 0 ○ , localized particles will drift out of the last closed flux surface and be lost. The simulation results show that the loss of localized particles in the CFQS is mainly located in some specific regions within the vicinity of the toroidal angle φ = 0 ○ , 18 0 ○ and the poloidal angle θ = 9 0 ○ . [ABSTRACT FROM AUTHOR]

Proton imaging is a powerful tool for probing electromagnetic fields in a plasma, providing a path-integrated map of the field topology. However, in cases where the field structure is highly inhomogeneous, inferring spatial properties of the underlying field from proton images can be difficult. This problem is exemplified by recent experiments, which used proton imaging to probe the filamentary magnetic field structures produced by the Weibel instability in collisionless counterstreaming plasmas. In this paper, we perform analytical and numerical analyses of proton images of systems containing many magnetic filaments. We find that, in general, the features observed on proton images do not directly correspond to the spacing between magnetic filaments (the magnetic wavelength) as has previously been assumed and that they instead correspond to the filament size. We demonstrate this result by Fourier analysis of synthetic proton images for many randomized configurations of magnetic filaments. Our results help guide the interpretation of experimental proton images of filamentary magnetic structures in plasmas. [ABSTRACT FROM AUTHOR]

In this paper, a coaxial gyrotron with a double electron beam is investigated on the ohmic loss, starting current, mode competition, and beam wave interaction. Its beam-wave interaction equation and dependence of the starting current on the static magnetic field and parameter K (the ratio of the current of an electron beam to the total current of a double electron beam) are given. On this basis, a 4 MW 170 GHz coaxial gyrotron with a double electron beam is studied. By the numerical calculation, the influence of K on the starting current is analyzed under a given static magnetic field and at different guiding center radii of the double electron beam. Studies show that the ability of the gyrotron to suppress the mode competition has a relation with the selection of the guiding center radii of the two electron beams. Then, the investigation of a time-dependent multi-mode competition of the gyrotron is performed, which shows that the gyrotron can stably operate in the T E 38 , 18 − mode. The results show that by considering the ohmic loss, the gyrotron can operate at 170.4 GHz with the output power of 4.04 MW and the beam-wave interaction efficiency of 35.47% when B = 6.82 T, U = 88 KV, I = 129 A, and K = 0.4. Each of the two electron beams in the gyrotron has different beam-wave interaction efficiencies. Compared to the coaxial gyrotron with one electron beam, the beam-wave interaction efficiency of each electron beam in the coaxial gyrotron with a double electron beam has a different increase. [ABSTRACT FROM AUTHOR]

With the rapid development of automation techniques, the claw-poled stepping motors are widely used, especially in optical scanners. Most research on stepping motors has concentrated on electrical properties, but a few papers have reported on the performance in relation to the magnetic-field distribution. The purpose of this study is to analyze the effect of the shape of the stator tooth and to search for the optimal parameters using a magnetic-circuit simulation technique. Through numerical calculations, the optimal performance conditions can be obtained by varying the ratio between the two bases of the trapezoid teeth between 1/5 and 1/3, and varying the ratio between the heights of individual teeth in relation to the separation between the two yokes from 11/15 to 13/15. [ABSTRACT FROM AUTHOR]

This paper describes a large-scale and highly accurate magnetostatic field analysis targeting a magnetic shield. The hybrid finite element–boundary element method and the magnetic moment method do not require mesh division for the free space and can easily treat the nonlinearity of magnetic property. Therefore, these methods are considered very effective for the analysis of a magnetic shield which has a high aspect ratio of the size scale to the thickness. However, large memory and computational time have been required due to the dense matrices generated by those integral-based formulations. To overcome the difficulties, we introduce the fast multipole method into the hybrid method and the magnetic moment method. Furthermore, to achieve more reduction of CPU time, we propose an effective preconditioning technique suitable for the hybrid method. Some numerical results that demonstrate the effectiveness of these approaches are also presented. [ABSTRACT FROM AUTHOR]

In this paper, we study the swimming of a model organism, the so-called Taylor's swimming sheet, in a viscoelastic fluid close to a solid boundary. This situation comprises natural habitats of many swimming microorganisms, and while previous investigations have considered the effects of both swimming next to a boundary and swimming in a viscoelastic fluid, seldom have both effects been considered simultaneously.We re-visit the small wave amplitude result obtained by Elfring and Lauga ["Theory of locomotion through complex fluids," in Complex Fluids in Biological Systems, Biological and Medical Physics, Biomedical Engineering, edited by S. E. Spagnolie (Springer New York, New York, NY, 2015), pp. 283-317] and give a mechanistic explanation to the decoupling of the effects of viscoelasticity, which tend to slow the sheet, and the presence of the boundary, which tends to speed up the sheet. We also develop a numerical spectral method capable of finding the swimming speed of a waving sheet with an arbitrary amplitude and waveform. We use it to show that the decoupling mentioned earlier does not hold at finite wave amplitudes and that for some parameters the presence of a boundary can cause the viscoelastic effects to increase the swimming speed of microorganisms. [ABSTRACT FROM AUTHOR]

This paper numerically investigates the influence of a uniform magnetic field on the droplet formation process at a microfluidic flow focusing configuration. The mathematical model was formulated by considering the balance of forces such as interfacial tension, magnetic force, and viscous stress across the liquid/liquid interface. A linearly magnetizable fluid was assumed. The magnetic force acts as a body force where the magnetic permeability jumps across the interface. The governing equations were solved with finite volume method on a Cartesian fixed staggered grid. The evolution of the interface was captured by the particle level set method. The code was validated with the equilibrium steady state of a ferrofluid droplet exposed to a uniform magnetic field. The evolution of the droplet formation in a flow focusing configuration was discussed. The paper mainly analyzes the effects of magnetic Bond number and the susceptibility on the velocity field and the droplet size. The droplet size increased with increasing magnetic strength and susceptibility. [ABSTRACT FROM AUTHOR]

Lower-hybrid drift and Buneman instabilities operate in current sheets with or without the guide field. The lower-hybrid drift instability is a universal instability in that it operates for all parameters. In contrast, the excitation of Buneman instability requires sufficiently thin current sheet. That is, the relative electron-ion drift speed must exceed the threshold in order for Buneman instability to operate. Traditionally, the two instabilities were treated separately with different mathematical formalisms. In a recent paper, an improved electrostatic dispersion relation was derived that is valid for both unstable modes [P. H. Yoon and A. T. Y. Lui, Phys. Plasmas 15, 072101 (2008)]. However, the actual numerical analysis was restricted to a one-dimensional situation. The present paper generalizes the previous analysis and investigates the two-dimensional nature of both instabilities. It is found that the lower-hybrid drift instability is a flute mode satisfying k·B=0 and k·∇n=0, where k represents the wave number for the most unstable mode, B stands for the total local magnetic field, and ∇n is the density gradient. This finding is not totally unexpected. However, a somewhat surprising finding is that the Buneman instability is a field-aligned mode characterized by k×B=0 and k·∇n=0, rather than being a beam-aligned instability. [ABSTRACT FROM AUTHOR]

In this paper, the generation of high-quality and high-energy micro electron beam in vacuum by a chirped Gaussian laser pulse in the presence of an axial magnetic field is numerically investigated. The features of energy and angular spectra, emittances, and position distribution of electron beam are compared in two cases, i.e., in the presence and absence of an external magnetic field. The electron beam is accelerated with higher energy and qualified in spatial distribution in the presence of the magnetic field. The presence of an axial magnetic field improves electron beam spatial quality as well as its gained energy through keeping the electron motion parallel to the direction of propagation for longer distances. It has been found that a 64 μm electron bunch with about MeV initial energy becomes a 20 μm electron beam with high energy of the order of GeV, after interacting with a laser pulse in the presence of an external magnetic field. [ABSTRACT FROM AUTHOR]

Experimental data compiled over five decades of dense plasma focus research are consistent with the snowplow model of sheath propagation, based on the hypothetical balance between magnetic pressure driving the plasma into neutral gas ahead and "wind pressure" resisting its motion. The resulting sheath velocity, or the numerically proportional "drive parameter," is known to be approximately constant for devices optimized for neutron production over 8 decades of capacitor bank energy. This paper shows that the validity of the snowplow hypothesis, with some correction, as well as the non-dependence of sheath velocity on device parameters, have their roots in local conservation laws for mass, momentum, and energy coupled with the ionization stability condition. Both upper and lower bounds on sheath velocity are shown to be related to material constants of the working gas and independent of the device geometry and capacitor bank impedance. [ABSTRACT FROM AUTHOR]

The plasma scrape-off layer (SOL) in a tokamak is characterized by ion flow down a long narrow flux tube terminating on a solid surface. The ion flow velocity along a magnetic field line can be equal to or greater than sonic at the entrance of a Debye sheath or upstream in the presheath. This paper presents a numerical study of the transition between subsonic and supersonics flows. A quasineutral one-dimensional (1D) fluid code has been used for modeling of plasma transport in the SOL along magnetic field lines, both in steady state and under transient conditions. The model uses coupled equations for continuity, momentum, and energy balance with ionization, radiation, charge exchange, and recombination processes. The recycled neutrals are described in the diffusion approximation. Standard Bohm sheath criterion is used as boundary conditions at the material surface. Three conditions conducive for the generation of supersonic flows in SOL plasmas have been explored. It is found that in steady state high (attached) and low (detached) divertor temperatures cases, the role of particle, momentum, and energy loss is critical. For attached case, the appearance of shock waves in the divertor region if the incoming plasma flow is supersonic and its effect on impurity retention is presented. In the third case, plasma expansion along the magnetic field can yield time-dependent supersonic solutions in the quasineutral rarefaction wave. Such situations can arise in the parallel transport of intermittent structures such as blobs and edge localized mode filaments along field lines. [ABSTRACT FROM AUTHOR]

Magnetic microbubbles are a relatively recent development with the potential to greatly improve the efficacy of the minimally invasive drug-delivery procedure sonoporation. However, very little is known about the dynamics of magnetic microbubbles, in general. In this paper, a novel mathematical model and numerical method are developed to simulate the dynamics of non-spherical magnetic microbubbles in vitro. The ambient fluid is assumed to be inviscid and the flow irrotational, enabling a generalized Bernoulli equation to be derived that includes surface tension effects and the effect of the applied magnetic field. The governing equations are solved using the boundary element method in which both the bubble surface and the velocity potential are represented by cubic splines. Results show that magnetic microbubble dynamics are highly dependent on the magnetic susceptibility difference, ∆X, between the bubble and the ambient fluid, with the sign and magnitude of ∆X dictating the direction and velocity of any formed liquid jets. Importantly, it is shown that the magnetic field can provide an additional means of flow control to the experimental investigator: in the presence of surface tension, weak magnetic fields do not generate jets. However, increasing the magnitude of the magnetic field can instigate jet formation, and increase the maximum and time-averaged jet velocities. Experimentally relevant parameter values are also considered, and results suggest that a combined application of magnetic and ultrasound fields is required to generate the high-speed bubble collapse events most likely to maximise cell poration and drug delivery. [ABSTRACT FROM AUTHOR]

A magnetic field is defined by the property that its divergence is zero in a threedimensional oriented Riemannian manifold. Each magnetic field generates amagnetic flow whose trajectories are curves called as magnetic curves. In this paper, we give a new variational approach to study the magnetic flow associated with the Killing magnetic field in a three-dimensional oriented Riemann manifold, (M³, g). And then, we investigate the trajectories of the magnetic fields called as N-magnetic and Bmagnetic curves. [ABSTRACT FROM AUTHOR]

Implosions of (quasi-)spherical loads with mass redistribution and displacement modification are investigated numerically. Both methods can theoretically counterbalance the nonuniformity of magnetic pressure along the load surface and realize quasi-spherical Z-pinch implosions. Mass redistribution is feasible for spherical loads with large radius and weight, while the displacement modification is more suitable for light loads, such as those composed of wire arrays. Simulation results suggest that, for mass redistributed spherical loads, wall instabilities induced by polar mass flows will deform the imploding shell. For prolate spherical loads, in which the wall instability cannot develop, the kinetic energy distribution is disturbed at high latitude. These passive behaviors and their possible mitigation methods, such as reshaping the electrode, are investigated numerically in this paper. [ABSTRACT FROM AUTHOR]

In this paper we extend the new techniques of W. Tierens and D. D. Zutter, J. Comput. Phys. 231, 5144 (2012) to include finite Larmor radius effects up to second order in the Larmor radius. We limit ourselves to the case of propagation perpendicular to the background magnetic field Bvector 0. We show that our time-domain technique is able to produce the lowest-order Bernstein wave (a wave believed to be useful for heating fusion devices [H. P. Laqua, Plasma Phys. Controlled Fusion 49, R1 (2007)]). The discrete equations retain many of the favourable properties described in W. Tierens and D. D. Zutter, J. Comput. Phys. 231, 5144 (2012), i.e., unconditional stability and a straightforward relation between the second-order accurate continuous dispersion relation and the dispersion relation of the discretized problem. The theory is illustrated by a place-independent and a place-dependent temperature numerical example. [ABSTRACT FROM AUTHOR]

The wandering of magnetic field lines is an important topic in theoretical plasma physics and astrophysics. Previous analytical work, as well as computer simulations, is based on magnetostatic models to warrant mathematical and numerical tractability. Recently, we have studied the first time field line random walk in dynamical turbulence by using a field line tracing method. These calculations were performed for a slab model of the turbulence. It is the purpose of the present paper to use the latter method to compute the field line diffusion coefficient for dynamical two-dimensional turbulence. Two models for the dynamical correlation function are used, namely the damping model of dynamical turbulence and the nonlinear anisotropic dynamical turbulence model. It is shown that the largest scales of the turbulence and the choice of the dynamical turbulence model have a strong influence on the diffusivity of the field lines and the absolute value of the diffusion coefficient. [ABSTRACT FROM AUTHOR]

A novel large-orbit electron gun with a gradual reversal of the magnetic field for the 8 mm third-harmonic peniotron is numerically studied in this paper. In the permanent magnetic condition, an electron beam with axial velocity spread of 4.78%, guiding center deviation of 7.18%, and velocity ratio of 2.2 is generated in the gun, satisfying the special requirements of beam-wave interaction in the third-harmonic peniotron. Compared with the traditional three stage method, it is unnecessary here to generate a thin tubular electron beam before the reversal point and the mutational distribution of reversal magnetic field. In addition, the emission band can be placed in the area of axial magnetic field where its magnitude decreases gradually, which reduces the structure complexity and the difficulties in tube-making process greatly. Finally, the peniotron is predicted to yield an output power of 31.9 kW at 30 GHz, with its microwave conversion efficiency of up to 49.4%, which shows the perfect matching between its high frequency system and its electron-optical system. According to the numerical calculations, the device performance is very sensitive to the relative axial position between the magnetic system and the electrode system. If the magnetic system shifts 0.8 mm to the right, the device efficiency would decrease from 49.4% to 31.7%. This phenomenon provides meaningful information in the device hot-test. [ABSTRACT FROM AUTHOR]

In this paper, we discuss the formation of ring-shaped and gyro-phase restricted velocity distribution functions (VDFs) at the edges of a cloud of protons injected into non-uniform distributions of the electromagnetic field. The velocity distribution function is reconstructed using the forward test-kinetic method. We consider two profiles of the electric field: (1) a non-uniform E-field obtained by solving the Laplace equation consistent with the conservation of the electric drift and (2) a constant and uniform E-field. In both cases, the magnetic field is similar to the solutions obtained for tangential discontinuities. The initial velocity distribution function is Liouville mapped along numerically integrated trajectories. The numerical results show the formation of an energy-dispersed structure due to the energy-dependent displacement of protons towards the edges of the cloud by the gradient-B drift. Another direct effect of the gradient-B drift is the formation of ring-shaped velocity distribution functions within the velocity-dispersed structure. Higher energy particles populate the edges of the proton beam, while smaller energies are located in the core. Non-gyrotropic velocity distribution functions form on the front-side and trailing edge of the cloud; this effect is due to remote sensing of energetic particles with guiding centers inside the beam. The kinetic features revealed by the test-kinetic solutions have features similar to in-situ velocity distribution functions observed by Cluster satellites in the magnetotail, close to the neutral sheet. [ABSTRACT FROM AUTHOR]

In this paper, the effects of magnetic field gradient (i.e., the magnetic field transition layer effects) on the Rayleigh-Taylor instability (RTI) with continuous magnetic field and density profiles are investigated analytically. The transition layers of magnetic field and density with two different typical profiles are studied and the analytic expressions of the linear growth rate of the RTI are obtained. It is found that the magnetic field effects strongly reduce the linear growth rate of the RTI, especially when the perturbation wavelength is short. The linear growth rate of the RTI increases with the thickness of the magnetic field transition layer, especially for the case of small thickness of the magnetic field transition layer. When the magnetic field transition layer width is long enough, the linear growth rate of the RTI can be saturated. Thus when one increases the width of the magnetic field transition layer, the linear growth rate of the RTI increases only in a certain range, which depends on the magnetic field strength. The numerical results are compared with the analytic linear growth rates and they agree well with each other. [ABSTRACT FROM AUTHOR]

The Rayleigh-Taylor instability arises whenever two fluids with different densities are arranged such that the heavier fluid sits above the lighter fluid, with a sharp interface in between. The magnetic Rayleigh-Taylor instability has the further complication due to the presence of a magnetic field throughout both media. The two fluids in question may also have differing magnetic properties, such as the magnetic permeability. When the fluids in consideration are in fact plasmas comprised of charged particles, induced currents, magnetic fields, and Lorentz forces can all act in ways that will affect the stability of the system. Stable base flows exist for the 2D case, and small sinusoidal disturbances to the base flow will grow in the unstable scenario. The numerical method described in this paper calculates the growth of the interface in the nonlinear regime, since closed form solutions are obtained only in the linear approximation. Through the analysis of both the fluid and magnetic vorticities and streamfunctions, the simulated results can be explained from the principles of magnetohydrodynamics. A range of simulations is presented, looking at cases with different initial conditions, cases with strong and weak magnetic fields, and cases with magnetic fields oriented at different angles relative to the interface of the two fluids. It is shown in particular how different initial conditions give rise to outcomes that are very different in terms of the geometry of the interface between the two fluids, primarily the differences between a single mode disturbance and a multimode disturbance to the interface at time t = 0. [ABSTRACT FROM AUTHOR]

In recent work, a technique was developed to extract high quality beams from single-component plasmas confined in a Penning-Malmberg trap in a 4.8 T magnetic field. In this paper, a procedure is developed to extract these beams from the confining magnetic field and then focus them to create especially tailored electrostatic beams. Electron beams are extracted from the field in two stages: they are first transported to a region of reduced field (1 mT), and then taken to zero field with a nonadiabatic, fast extraction. Once in the field-free region, the beams are focused using an Einzel lens. Experimental results and numerical simulations are presented to illustrate the extraction and focusing process. Theoretical expressions are developed to describe the modifications to the relevant beam energy and spatial distributions. Where possible, analytic expressions are presented for the case relevant here of beams with Gaussian radial profiles. Beam emittance considerations are discussed as well as prospects for further development of these techniques. Application of these techniques to provide high-quality positron beams is also discussed. [ABSTRACT FROM AUTHOR]

Magnetic fields are an essential component of a plasma. In many astrophysical, solar, magnetospheric, and laboratory situations the magnetic field in the plasma can be very dynamic and form highly complex structures. One approach to unraveling these structures is to determine the magnetic skeleton of the field, a set of topological features that divide the magnetic field into topologically distinct domains. In general, the features of the magnetic skeleton are difficult to locate, in particular those given by numerical experiments. In this paper, we propose a new set of tools to find the skeleton of general magnetic fields including null points, spines, separatrix surfaces, and separators. This set of tools is found to be considerably better at finding the skeleton than the currently favored methods used in magnetohydrodynamics. [ABSTRACT FROM AUTHOR]

The paper investigates the role of relative motion between the fluid components of a plasma model, which is simulated by concatenation of two anisotropic magnetohydrodynamic (MHD) fluids, on the propagation of low-frequency waves and instabilities. The gyrotropic pressure of both the MHD components is given by generalized polytropic laws that allow the system to reduce to a variety of states. The linearized analysis is carried out and dispersion relation is derived using the normal mode technique. The dispersion relation, which gives numerous earlier results as special cases, is discussed both analytically as well as numerically for parameters appropriate for the space plasma. It is found that the relative motion between the two components, besides modifying the condition for fire-hose instability, causes the concatenated system to exhibit different kinds of relative ordering of the magnitudes of phase speeds of the various MHD modes in the directions parallel and antiparallel to the magnetic field. The relative motion between the components also influences the phase relationship between density and magnetic field fluctuations for the various compressive modes in the background ambient plasma. [ABSTRACT FROM AUTHOR]

Previous work on asymmetry-induced transport in a modified Malmberg–Penning trap showed that the radial particle flux was empirically constrained to be of the form Γ(ε)=-(B0/B)1.33D(ε)[∇n0+f(ε)], where ε=ω-lωR, ωR(r)=vθ/r is the column rotation frequency, ω and l are the asymmetry frequency and azimuthal mode number, ∇n0 is the radial density gradient, B is the magnetic field, B0 is an empirical constant, and D(ε) and f(ε) are unknown functions. In this paper, it is shown that further constraints can be placed upon D(ε) and f(ε) by comparing data near the ε=0 points to a first order expansion of Γ(ε). It is shown that dD/dε(0)≠0, in contradiction to resonant particle theory, and that f(ε) can only be a fraction of the size predicted by that theory. Finally, it is shown that dD/dε(0) exhibits a power-law scaling with radius, magnetic field, and the bias of the center conductor of the trap. [ABSTRACT FROM AUTHOR]

In this paper, memory and nonlocal effects on fluctuating mass diffusion are addressed in the context of fusion plasmas. Nonlocal effects are included by considering a diffusivity coefficient depending on the size of the container in the transverse direction to the applied magnetic field. It is obtained by resorting to the general formulation of the extended version of irreversible thermodynamics in terms of the higher order dissipative fluxes. The developed model describes two different types of the particle density time correlation function. Both have been observed in tokamak and nontokamak devices. These two kinds of time correlation function characterize the wave and the diffusive transport mechanisms of particle density perturbations. A transition between them is found, which is controlled by the size of the container. A phase diagram in the {L,2π/k} space describes the relation between the dynamics of particle density fluctuations and the size L of the system together with the oscillating mode k of the correlation function. [ABSTRACT FROM AUTHOR]

Plasma from the Earth’s magnetosheath has previously been observed inside the magnetosphere. Inhomogeneities in the magnetosheath plasma, here called plasmoids, can impact the magnetopause and doing so set up a polarizing field that allows it to penetrate the magnetopause and enter the magnetosphere. A set of simulations of plasmoids with different dimensions is presented in this paper. For plasmoids that are longer than those previously published, waves propagating upstream from the barrier are found. It is also found that the penetration process causes the part of the plasmoid that is upstream of the barrier to rotate. The role of plasmoid width and cross sectional shape in penetration is studied, and for plasmoids that are less than half an ion gyroradius wide, the plasmoid is compressed to obtain a vertically oriented elliptical cross section, regardless of the initial shape. When the initial plasmoid width exceeds the ion gyroradius, the plasmoid still penetrates through a mechanism involving a potential that propagates upstream from the magnetic barrier. [ABSTRACT FROM AUTHOR]

Drift wave turbulence is known to self-organize to form axisymmetric macroscopic flows. The basic mechanism for macroscopic flow generation is called inverse energy cascade. Essentially, it is an energy transfer from the short wavelengths to the long wavelengths in the turbulent spectrum due to nonlinear interactions. A class of macroscopic flows, the poloidally symmetric zonal flows, is widely recognized as a key constituent in nearly all cases and regimes of microturbulence, also because of the realization that zonal flows are a critical agent of self-regulation for turbulent transport. In tokamaks and other toroidal magnetic confinement systems, axisymmetric flows exist in two branches, a zero frequency branch and a finite frequency branch, named Geodesic Acoustic Modes (GAMs). The finite frequency is due to the geodesic curvature of the magnetic field. There is a growing body of evidence that suggests strong GAM activity in most devices. Theoretical investigation of the GAMs is still an open field of research. Part of the difficulty of modelling the GAMs stems from the requirement of running global codes. Another issue is that one cannot determine a simple one to one relation between turbulence stabilization and GAM activity. This paper focuses on the study of ion temperature gradient turbulence in realistic tokamak magnetohydrodynamic equilibria. Analytical and numerical analyses are applied to the study of geometrical effects on zonal flows oscillations. Results are shown on the effects of the plasma elongation on the GAM amplitude and frequency and on the zonal flow residual amplitude. [ABSTRACT FROM AUTHOR]

The moment equation approach to neoclassical transport is used to calculate neoclassical particle and heat fluxes, impurity transport, the ambipolar electric field, and momentum damping rates. These equations are often written in Hamada coordinates which makes it easier to obtain analytic solutions. However, previous simplifying assumptions used to evaluate the basis vectors analytically are often invalid for advanced stellarator configurations. In this paper, a numerical method is presented by which the Hamada basis set can be determined for an arbitrary three dimensional toroidal confinement device by integrating along a magnetic field line. The method is applied to the magnetic configuration in the Helically Symmetric Experiment [F. S. B. Anderson, A. F. Almagri, D. T. Anderson, P. G. Matthews, J. N. Talmadge, and J. L. Shohet, Fusion Technol. 27, 273 (1995)] and compared to the large-aspect-ratio tokamak approximation to the basis set. The results indicate that the numerical technique is a more accurate method to specify the basis vectors, especially in a device with negligible toroidal curvature. [ABSTRACT FROM AUTHOR]

This paper presents a combined numerical and analytical study of the impulsive axisymmetric spin-up from rest of an isothermal liquid metal in a closed cylinder. The motion of the liquid is caused by the action of a low-frequency, low-induction rotating magnetic field, whose magnetic Taylor number is in the range (0.01–0.9) Tacr 3D with Tacr 3D given by Grants and Gerbeth. The computations were performed for cylindrical containers of aspect ratios (diameter/ height) R equal to 0.5, 1, and 2. The numerical simulations are compared with the predictions of an analytical model, valid for small Ekman numbers E extending a former work by Ungarish. The first phase of the motion from rest is an initial adjustment: the inviscid fluid begins to rotate due to the externally forced azimuthal acceleration, and concomitantly a meridional velocity field is induced by the unbalanced centrifugal effects; this occurs during, approximately, the first revolution about the axis. Subsequently, the spin-up flow is dominated by an axially independent swirl in the core (which can be regarded as a geostrophic mode) and Bödewadt (Ekman) layers on the top-bottom boundaries. This is accompanied by inertial oscillations. Both the development of the swirl motion and the decay of the inertial mode occur on the spin-up time scale predicted by the asymptotic model. The investigation also points out the important role of the aspect-ratio parameter in the analysis of the magnetohydrodynamics-driven spin-up. The efficiency of the driving decays strongly for R>0.6 and becomes close to zero at R ≈3.3. The numerical results confirm the analytical inference that the thickness of the sidewall layer for the angular velocity adjustment is E¼ ∕ √R, and when this parameter is larger than 0.3 the sidewall influence gains dominance over the Ekman-layer effects. The frequency of the inertial oscillations increases with R. [ABSTRACT FROM AUTHOR]

In this paper, a new magnetically levitated synchronous permanent magnet planar motor (MLSPMPM) driven by composite-current is proposed, of which the mover is made of a copper coil array and the stator are magnets and magnetic conductor. The coil pitch τt and permanent magnet pole pitch τp satisfy the following relationship 3nτt = (3n ± 1)τp. Firstly, an analytical model of the planar motor is established, flux density distribution of the two-dimensional magnet array is obtained by solving the equations of the scalar magnetic potential. Secondly, the expressions of the electromagnetic forces induced by magnetic field and composite current are derived. To verify the analytical model and the electromagnetic forces, finite element method (FEM) is used for calculating the flux density and electromagnetic forces of the MLSPMPM. And the results from FEM are in good agreement with the results from the analytical equations. This indicates that the analytical model is reasonable. [ABSTRACT FROM AUTHOR]