Equations for nonlinear MHD convection in shearless magnetic systems

A closed set of reduced dynamic equations is derived that describe nonlinear low-frequency flute MHD convection and resulting nondiffusive transport processes in weakly dissipative plasmas with closed or open magnetic field lines. The equations obtained make it possible to self-consistently simulate transport processes and the establishment of the self-consistent plasma temperature and density profiles for a large class of axisymmetric nonparaxial shearless magnetic devices: levitated dipole configurations, mirror systems, compact tori, etc. Reduced equations that are suitable for modeling the long-term evolution of the plasma on time scales comparable to the plasma lifetime are derived by the method of the adiabatic separation of fast and slow motions.     

Adiabatically reduced magnetohydrodynamic equations for a cylindrical plasma with an anisotropic pressure

A closed set of reduced equations describing low-frequency nonlinear flute magnetohydrodynamic (MHD) convection and the resulting nondiffusive processes of particle and energy transport in a weakly collisional cylindrical plasma with an anisotropic pressure is derived. The Chew-Goldberger-Low anisotropic magnetohydrodynamics is used as the basic dynamic model, because this model is applicable to describing flute convection in a cylindrical plasma column even in the low-frequency limit. The reduced set of equations was derived using the method of adiabatic separation of fast and slow motions. It is shown that the structure of the adiabatic transformation and the corresponding velocity field are identical to those obtained earlier in the isotropic MHD model. However, the derived heat transfer equations differ drastically from the isotropic pressure model. In particular, they indicate a tendency toward maintaining different radial profiles of the longitudinal and transverse pressures.     

Numerical modeling of the dynamics of a slow Z-pinch

A study is made of the method for numerical modeling of pulsed plasma systems by simultaneously solving two-temperature MHD equations and the equations of ionization kinetics. As an example, the method is applied to simulate a relatively slow moderate-density Z-pinch, whose dynamics is well studied experimentally. A specially devised two-dimensional computer code makes use of a promising technique of parallel modeling.     

Reduced MHD equations with allowance for slow magnetosonic waves

The so-called reduced magnetohydrodynamics, which deals with the motion of incompressible fluids and is usually applied to describe plasma flows in a strong toroidal magnetic field, has a number of drawbacks and, in some cases, fails to produce correct results. The equations proposed here are simpler than the original MHD equations and are free of these drawbacks. These equations, like reduced MHD equations, make it possible to remove from consideration fast magnetosonic waves and to introduce the vector potential for the poloidal magnetic field. However, our equations differ from the reduced MHD equations in that they completely incorporate slow magnetosonic waves, the specific features of the toroidal geometry, and the effects of the toroidal velocity.     

Electric Field and Poloidal Rotation in the Turbulent Edge Plasma of the T-10 Tokamak

Plasma Physics Reports - The dynamics of turbulent edge plasma in the T-10 tokamak is simulated numerically by solving reduced nonlinear MHD equations of Braginskii’s two-fluid hydrodynamics....     

Grand Molecular Dynamics: A Method for Open Systems

We present a new molecular dynamics method for studying the dynamics of open systems. The method couples a classical system to a chemical potential reservior. In the formulation, following the extended system dynamics approach, we introduce a variable, v to represent the coupling to the chemical potential reservoir. The new variable governs the dynamics of the variation of number of particles in the system. The number of particles is determined by taking the integer part of v. The fractional part of the new variable is used to scale the potential energy and the kinetic energy of an additional particle: i.e., we introduce a fractional particle. We give the ansatz Lagrangians and equations of motion for both the isothermal and the adiabatic forms of grand molecular dynamics. The averages calculated over the trajectories generated by these equations of motion represent the classical grand canonical ensemble (μVT) and the constant chemical potential adiabatic ensemble (μVL) averages, respectively. The microcanonical phase space densities of the adiabatic and isothermal forms the molecular dynamics method are shown to be equivalent to adiabatic constant chemical potential ensemble, and grand canonical ensemble partition functions. We also discuss the extension to multi-component systems, molecular fluids, ionic solutions and the problems and solutions associated with the implementation of the method. The statistical expressions for thermodynamic functions such as specific heat; adiabatic bulk modulus, Grüneissen parameter and number fluctuations are derived. These expressions are used to analyse trajectories of constant chemical potential systems.     

Low-frequency instabilities of collisionless plasma and the 16-moment approximation

Using the 16-moment equations that take into account heat fluxes in anisotropic collisionless plasma, the properties of magnetohydrodynamic (MHD) instabilities are investigated. For all instabilities occurring in the MHD approach (the normal incompressible firehose instability, the second compressible almost longitudinal firehose instability, and the almost transverse mirror instability of slow magnetosonic modes, as well as thermal instability caused by the heat flux directed along the magnetic field), their kinetic analogs are considered. The kinetic dispersion relation in the low-frequency range in the vicinity of the ion thermal velocity is analyzed. The flow of plasma ions along the magnetic field is taken into account. The thresholds and instability growth rates obtained in the MHD and kinetic approaches are found to be in good agreement. This indicates that the 16-moment MHD equations adequately describe the dynamics of collisionless plasma.     

Magnetohydrodynamics of a weakly ionized plasma: Ambipolar magnetic diffusion and shock front structure

Kinetic equations with the BGK collision integral are used to derive MHD equations for a weakly ionized plasma that are applicable over a broad range of magnetic field strengths. In strong magnetic fields, a substantial contribution to the transverse diffusion of the magnetic field comes from the ambipolar magnetic diffusion, which is associated with the motion of both the charged component and the magnetic field against the background of the neutral plasma component. The problems of the magnetic field diffusion in a weakly ionized plasma and the shock wave structure are solved.     

Neural rate equations for bursting dynamics derived from conductance-based equations

A method of obtaining rate equations from conductance-based equations is developed and applied to fast-spiking and bursting neocortical neurons. It involves splitting systems of conductance-based equations into fast and slow subsystems, and averaging the effects of fast terms that drive the slowly varying quantities by showing that their average is closely proportional to the firing rate. The dependence of the firing rate on the injected current is then approximated in the analysis. The resulting behavior of the slow variables is then substituted back into the fast equations, with the further approximation of replacing the fast voltages in these terms by effective values. For bursting neurons the method yields two coupled limit-cycle oscillators: a self-exciting oscillator for the slow variables that commences limit-cycle oscillations at a critical current and modulates a fast spike-generating oscillator, thereby leading to slowly modulated bursts with a group of spikes in each burst. The dynamics of these coupled oscillators are then verified against those of the conductance-based equations. Finally, it is shown how to place the results in a form suitable for use in mean-field equations for neural population dynamics.     

Influence of the Hall effect on the plasma dynamics in an MTF/MAGO chamber

The role of the Hall effect in experiments on the magnetic implosion of a D-T plasma in a cylindrical MTF/MAGO chamber fed from a helical explosive magnetic generator is investigated. The plasma dynamics is simulated numerically by a 2D code developed for solving the set of MHD equations with account of the Hall effect. In simulations, the generator, the break switch, and other units were replaced with LR circuits. It is shown that taking into account the Hall effect provides better agreement between numerical simulations and experimental data.     

MHD processes during the cascade development of the neck and hot spot in an X-pinch

Results are presented from two-dimensional MHD simulations of X-pinch implosion. The simulations were performed in the (r, z) and (x, y) geometries for homogeneous (dense plasma) and heterogeneous (core-corona) loads. The formation of a minidiode, the development of a neck and an X-radiating hot spot, and the influence of the plasma corona on the implosion dynamics of the dense X-pinch plasma were investigated. For through simulations, the conical neck model was used, whereas a detailed analysis of the X-ray burst was performed in the parabolic neck model. The MHD processes occurring during the implosion of oblique shock waves and the onset of instability of the plasma column were examined. It is found that, due to the quasi-periodic character of these processes, the neck compression proceeds in a cascade fashion. The plasma state in a hot spot just before the break of the neck is analyzed, and the possibility of generating fast particle beams is considered.     

Reduced models of networks of coupled enzymatic reactions

The Michaelis-Menten equation has played a central role in our understanding of biochemical processes. It has long been understood how this equation approximates the dynamics of irreversible enzymatic reactions. However, a similar approximation in the case of networks, where the product of one reaction can act as an enzyme in another, has not been fully developed. Here we rigorously derive such an approximation in a class of coupled enzymatic networks where the individual interactions are of Michaelis-Menten type. We show that the sufficient conditions for the validity of the total quasi-steady state assumption (tQSSA), obtained in a single protein case by Borghans, de Boer and Segel can be extended to sufficient conditions for the validity of the tQSSA in a large class of enzymatic networks. Secondly, we derive reduced equations that approximate the network's dynamics and involve only protein concentrations. This significantly reduces the number of equations necessary to model such systems. We prove the validity of this approximation using geometric singular perturbation theory and results about matrix differentiation. The ideas used in deriving the approximating equations are quite general, and can be used to systematize other model reductions.     

Effect of helical magnetic fields on plasma stability in tokamaks

Physical mechanisms for destabilization of MHD perturbations by external quasistatic magnetic fields and rotating helical magnetic fields in a tokamak plasma are identified using a numerical model of tearing modes in a viscous high-temperature plasma. The critical conditions for the onset of MHD perturbations and their dynamic model are compared with the experimental results from the JET tokamak. The model is used to predict how the stray magnetic fields will influence plasma stability in a tokamak reactor (ITER). __________ Translated from Fizika Plazmy, Vol. 26, No. 8, 2000, pp. 675–682. Original Russian Text Copyright ￠ 2000 by Savrukhin.     

生物协同学,Lorenz模型和种群动力学

由协同学方程出发,可以描述种群的大迁徙,由此又能够得到Lorenz模型,它可以描述两种种群的变化关系．当取绝热近似时,还可以导致种群动力学的不同模型．因此,生物协同学能够深刻揭示不同物种之间,既竞争又协同的复杂的非线性关系．     

Skin Effects in a Dusty Plasma

A multifluid MHD model is applied to study the magnetic field dynamics in a dusty plasma. The motion of plasma electrons and ions is treated against the background of arbitrarily charged, immobile dust grains. When the dust density gradient is nonzero and when the inertia of the ions and electrons and the dissipation from their collisions with dust grains are neglected, we are dealing with a nonlinear convective penetration of the magnetic field into the plasma. When the dust density is uniform, the magnetic field dynamics is described by the nonlinear diffusion equations. The limiting cases of diffusion equations are analyzed for different parameter values of the problem (i.e., different rates of the collisions of ions and electrons with the dust grains and different ratios between the concentrations of the plasma components), and some of their solutions (including self-similar ones) are found. The results obtained can also be useful for research in solid-state physics, in which case the electrons and holes in a semiconductor may be analogues of plasma electrons and ions and the role of dust grains may be played by the crystal lattice and impurity atoms.     

Determination of the plasma velocity in an imploding wire array from magnetic field measurements by a gradient probe

A method for measuring the gradient of the magnetic field in the plasma of an imploding wire array is described. Results from measurements of the magnitude and gradient of the magnetic field in a tungsten wire array on the Angara-5-1 facility at currents of ∼3 MA are presented. A novel method for calculating the velocity of the current-carrying plasma in the framework of MHD equations from data on the magnitude and gradient of the magnetic field at a certain point inside the array is proposed. It is demonstrated that a gradient magnetic probe can be used to investigate the plasma current sheath in plasma focus facilities.     

Steps toward numerical mode analysis of organizing systems

A well established method to analyze dynamical systems described by coupled nonlinear differential equations is to determine their normal modes and reduce the dynamics, by adiabatic elimination of stable modes, to a much smaller system for the amplitudes of unstable modes and their nonlinear interactions. So far, this analysis is possible only for idealized symmetric model systems. We aim to build a framework in which realistic systems with less symmetry can be analyzed automatically. In this paper we present a first example of mode analysis with the assistance of numerical computation. Our method is illustrated using a model system for the ontogenesis of retinotopy, and the results reproduce those from theoretical analysis precisely. Aspects of organization generalized from this model system are discussed. This research was supported by EU projects Daisy and SECO, and the Hertie Foundation.