Energy approach to stability analysis of the locked and rotating resistive wall modes in tokamaks

A method is proposed for stability analysis of the locked and rotating resistive wall modes (RWMs) in tokamaks. The method is based on the relations describing the balance of energy permeating through the vessel wall. This is a natural extension of the traditional energy approach to the plasma stability tasks which allows incorporation of the energy outflow (absent in the classical energy principle) and its dissipation in the wall. The proposed method covers the locked and rotating modes with a complex growth rate. Its efficiency is proved by derivation of a general dispersion relation for such modes with further reduction to particular consequences for slow and fast RWMs. It is shown that in the latter case, when the skin depth becomes smaller than the wall thickness, the mode rotation essentially amplifies its damping, weakening and even suppressing the instability. This effect was earlier found in the frame of the slab model [V. D. Pustovitov, Phys. Plasmas 19, 062503 (2012)]. Here, it is confirmed with equations valid for toroidal geometry, which are obtained as a supplement to the standard energy principle. The presented results predict strong rotational stabilization of the fast RWMs, which occurs at the mode rotation frequency above a critical level. The estimates are given to allow comparison of these predictions with experimental results.     

Modeling of the feedback stabilization of the resistive wall modes in a tokamak

The problem of feedback stabilization of the resistive wall modes (RWMs) in a tokamak is discussed. An equilibrium configuration with the parameters accepted for the stationary ITER scenario 4A is considered as the main scenario. The effect of the vacuum chamber's shape on the plasma stability is studied. Ideal MHD stability is analyzed numerically by using the KINX code. It is shown that, in a tokamak with the parameters of the designed T-15M tokamak, RWMs can be stabilized by a conventional stabilizing system made of framelike coils. However, the maximum possible gain in β in such a tokamak is found to be smaller than that in ITER. It is shown that, in this case, a reduction in the plasma—wall gap width by 10 cm allows one to substantially increase the β limit, provided that RWMs are stabilized by active feedback.     

The linear model and experimentally observed resonant field amplification in tokamaks and reversed field pinches

A review is given of the experimentally observed effects related to the resonant field amplification (RFA) and the Resistive Wall Mode (RWM) instability in tokamaks and reversed field pinches (RFPs). This includes the feedback rotation of RWM in RFX-mod RFP, dependence of the RWM growth rate on the plasma-wall separation observed in JT-60U, appearance of the slowly growing RWM precursors in JT-60U and similar phenomena in other devices. The experimental results are compared with theoretical predictions based on the model comprising the Maxwell equations, Ohm’s law for the conducting wall, the boundary conditions and assumption of linear plasma response to the external magnetic perturbations. The model describes the plasma reaction to the error field as essentially depending on two factors: the plasma proximity to the RWM stability threshold and the natural rotation frequency of the plasma mode. The linear response means that these characteristics are determined by the plasma equilibrium parameters only. It is shown that the mentioned effects in different devices under different conditions can be described on a common basis with only assumption that the plasma behaves as a linear system. To extend the range of the model validation, some predictions are derived with proposals for experimental studies of the RFA dynamics.     

Error field amplification near the stability boundary of the modes interacting with a conducting wall

The effect is considered of the amplification of an external resonant error field near the stability boundary of the so-called resistive wall modes observed in the DIII-D tokamak. The analysis is performed in a standard cylindrical approximation. The model is based on Maxwell’s equations and Ohm’s law; therefore, the results of the analysis are valid for any large-scale modes interacting with a conducting wall. In contrast to earlier approaches, the model considers the resonant field amplification as a dynamic effect. It is shown that, when the effect is the strongest, the stationary solutions proposed earlier are inapplicable. The problem of plasma response to a probing pulse of the resonant field of a given amplitude and duration is analyzed. The relationships obtained explain the basic features of the observed phenomena in the DIII-D tokamak and allow direct experimental verification.     

Drift stabilization of internal resistive-wall modes in tokamaks

The problem of drift stabilization of the internal resistive-wall modes (RWMs) in tokamaks is theoretically investigated. The basic assumption of the model is that, when the drift effects are neglected, these modes are unstable in the absence of a conducting wall and stable in the presence of a close-fitting perfectly conducting wall. In the former case, the instability condition is expressed as Δ′_∞>0, where Δ′_∞ is the matching parameter calculated under the assumption that the wall is removed to infinity. In the latter case, one has Δ _W ^′ <0, where Δ _W ^′ is the external matching parameter of tearing modes calculated assuming a perfectly conducting wall at the plasma boundary. In the case with a resistive wall, the relevant parameter can be either Δ′_∞ or Δ _W ^′ , depending on whether the value of the dimensionless parameter ωτ_s/2m is small or large, respectively (here ω is the mode frequency, τ_s is the resistive time constant of the wall, and m is the poloidal mode number). In the presence of drift effects, the mode frequency ω is approximately equal to the electron drift frequency, ω≈ω_*e. The value of the parameter ω_*eτ_s/2m, which therefore determines the behavior of internal RWMs, is estimated for several existing tokamaks, namely, AUG (ASDEX-Upgrade), DIII-D, JET, TFTR, and JT-60U, as well as for the projected ITER-FEAT. It is shown that, although drift effects do not stabilize internal RWMs in current devices, they should be efficient in suppressing these modes in reactor-grade tokamaks.     

Long-wavelength instability of periodic flows and helicon waves in electron magnetohydrodynamics

The stability of periodic flows and helicon waves against large-scale perturbations is investigated analytically in resistive electron magnetohydrodynamics by the method of two-scale expansions. It is shown that long-wavelength perturbations of a Kolmogorov-type flow are destabilized by the effect of negative resistivity. The destabilization of long-wavelength perturbations of a Beltrami-type helical flow and helicon waves is related to the microhelicity of the primary flow (wave). The instability of long-wavelength perturbations of an anisotropic helical flow is found to result from both the effect of negative resistivity and the effect associated with the microhelical nature of the flow. The criteria for the onset of the corresponding instabilities are derived. Numerical simulations are carried out based on nonlinear electron magnetohydrodynamic equations with initial conditions corresponding to the analytic formulation of the problem. The results of simulations on the whole confirm analytical results in the parameter range in which the latter are applicable and, in addition, extend the stability analysis to the parameter ranges that are beyond the scope of analytic approximations.     

Optimal Foraging and Predator–Prey Dynamics

A system consisting of a population of predators and two types of prey is considered. The dynamics of the system is described by differential equations with controls. The controls model how predators forage on each of the two types of prey. The choice of these controls is based on the standard assumption in the theory of optimal foraging which requires that each predator maximizes the net rate of energy intake during foraging. Since this choice depends on the densities of populations involved, this allows us to link the optimal behavior of an individual with the dynamics of the whole system. Simple qualitative analysis and some simulations show the qualitative behavior of such a system. The effect of the optimal diet choice on the stability of the system is discussed.     

Inner ear gene transfection in neonatal mice using adeno-associated viral vector: a comparison of two approaches

Local gene transfection is a promising technique for the prevention and/or correction of inner ear diseases, particularly those resulting from genetic defects. Adeno-associated virus (AAV) is an ideal viral vector for inner ear gene transfection because of its safety, stability, long-lasting expression, and its high tropism for many different cell types. Recently, a new generation of AAV vectors with a tyrosine mutation (mut-AAV) has demonstrated significant improvement in transfection efficiency. A method for inner ear gene transfection via the intact round window membrane (RWM) has been developed in our laboratory. This method has not been tested in neonatal mice, an important species for the study of inherited hearing loss. Following a preliminary study to optimize the experimental protocol in order to reduce mortality, the present study investigated inner ear gene transfection in mice at postnatal day 7. We compared transfection efficiency, the safety of the scala tympani injection via RWM puncture, and the trans-RWM diffusion following partial digestion with an enzyme technique. The results revealed that approximately 47% of inner hair cells (IHCs) and 17% of outer hair cells (OHCs) were transfected via the trans-RWM approach. Transfection efficiency via RWM puncture (58% and 19% for IHCs and OHCs, respectively) was slightly higher, but the difference was not significant.     

Pulse wave propagation in a model human arterial network: assessment of 1-D numerical simulations against in vitro measurements

A numerical model based on the nonlinear, one-dimensional (1-D) equations of pressure and flow wave propagation in conduit arteries is tested against a well-defined experimental 1:1 replica of the human arterial tree. The tree consists of 37 silicone branches representing the largest central systemic arteries in the human, including the aorta, carotid arteries and arteries that perfuse the upper and lower limbs and the main abdominal organs. The set-up is mounted horizontally and connected to a pulsatile pump delivering a periodic output similar to the aortic flow. Terminal branches end in simple resistance models, consisting of stiff capillary tubes leading to an overflow reservoir that reflects a constant venous pressure. The parameters required by the numerical algorithm are directly measured in the in vitro set-up and no data fitting is involved. Comparison of experimental and numerical pressure and flow waveforms shows the ability of the 1-D time-domain formulation to capture the main features of pulse wave propagation measured throughout the system test. As a consequence of the simple resistive boundary conditions used to reduce the uncertainty of the parameters involved in the simulation, the experimental set-up generates waveforms at terminal branches with additional non-physiological oscillations. The frequencies of these oscillations are well captured by the 1-D model, even though amplitudes are overestimated. Adding energy losses in bifurcations and including fluid inertia and compliance to the purely resistive terminal models does not reduce the underdamped effect, suggesting that wall visco-elasticity might play an important role in the experimental results. Nevertheless, average relative root-mean-square errors between simulations and experimental waveforms are smaller than 4% for pressure and 19% for the flow at all 70 locations studied.     

Nonlinear modulation of an extraordinary wave under the conditions of parametric decay

A self-consistent set of Hamilton equations describing nonlinear saturation of the amplitude of oscillations excited under the conditions of parametric decay of an elliptically polarized extraordinary wave in cold plasma is solved analytically and numerically. It is shown that the exponential increase in the amplitude of the secondary wave excited at the half-frequency of the primary wave changes into a reverse process in which energy is returned to the primary wave and nonlinear oscillations propagating across the external magnetic field are generated. The system of ??slow?? equations for the amplitudes, obtained by averaging the initial equations over the high-frequency period, is used to describe steady-state nonlinear oscillations in plasma.     

Microplane constitutive model and computational framework for blood vessel tissue

This paper presents a nonlinearly elastic anisotropic microplane formulation in 3D for computational constitutive modeling of arterial soft tissue in the passive regime. The constitutive modeling of arterial (and other biological) soft tissue is crucial for accurate finite element calculations, which in turn are essential for design of implants, surgical procedures, bioartificial tissue, as well as determination of effect of progressive diseases on tissues and implants. The model presented is defined at a lower scale (mesoscale) than the conventional macroscale and it incorporates the effect of all the (collagen) fibers which are anisotropic structural components distributed in all directions within the tissue material in addition to that of isotropic bulk tissue. It is shown that the proposed model not only reproduces Holzapfel's recent model but also improves on it by accounting for the actual three-dimensional distribution of fiber orientation in the arterial wall, which endows the model with advanced capabilities in simulation of remodeling of soft tissue. The formulation is flexible so that its parameters could be adjusted to represent the arterial wall either as a single material or a material composed of several layers in finite element analyses of arteries. Explicit algorithms for both the material subroutine and the explicit integration with dynamic relaxation of equations of motion using finite element method are given. To circumvent the slow convergence of the standard dynamic relaxation and small time steps dictated by the stability of the explicit integrator, an adaptive dynamic relaxation technique that ensures stability and fastest possible convergence rates is developed. Incompressibility is enforced using penalty method with an updated penalty parameter. The model is used to simulate experimental data from the literature demonstrating that the model response is in excellent agreement with the data. An experimental procedure to determine the distribution of fiber directions in 3D for biological soft tissue is suggested in accordance with the microplane concept. It is also argued that this microplane formulation could be modified or extended to model many other phenomena of interest in biomechanics.     

The permeability of SPION over an artificial three-layer membrane is enhanced by external magnetic field

Background  

Sensorineural hearing loss, a subset of all clinical hearing loss, may be correctable through the use of gene therapy. We are testing a delivery system of therapeutics through a 3 cell-layer round window membrane model (RWM model) that may provide an entry of drugs or genes to the inner ear. We designed an in vitro RWM model similar to the RWM (will be referred to throughout the paper as RWM model) to determine the feasibility of using superparamagnetic iron oxide (Fe₃O₄) nanoparticles (SPION) for targeted delivery of therapeutics to the inner ear.     

Muscle pathway geometry in the heart wall

Muscle fiber pathways in the heart wall are described. Procedures are introduced which permit data to be standardized from cadaver and animal hearts fixed at different points of the cardiac cycle or obtained in vivo from patients with different ejection fractions and heart masses. Design criteria are also developed here to construct a hypothetical standard left ventricle to compare the data from different hearts. The equations allow the nested set of toroidal fiber-shells to be depicted with typical muscle fiberpaths. With this formulation the heart wall and typical elements in it can be shown computergraphically as they move from the contracted state to the distended. Man-made fiber structures that simulate the fail-safe shockload absorbing features of the heart can now be designed and tested computergraphically by use of the mathematical procedures described here.     

Permselectivity of the glomerular capillary wall to macromolecules. I. Theoretical considerations.

The transport of macromolecules across the renal glomerular capillary wall has been described theoretically using flux equations based on (a) restricted transport through small pores, and (b) the Kedem-Katchalsky formulation. The various assumptions and limitations inherent in these two approaches are discussed. To examine the coupling between macromolecular solute transport and the determinants of glomerular filtration rate, these flux equations were combined with mass balance relations which allow for variations in the transmembrane driving forces along a glomerular capillary. It was predicted, using both pore theory and the Kedem-Katchalsky equations, that fractional solute clearance should be strongly dependent on the determinants of glomerular filtration rate when convection and diffusion both contribute to solute transport. When convection becomes the sole mechanism for transcapillary solute transport, however, fractional solute clearance is essentially independent of changes in the determinants of glomerular filtration rate. Consequently, unless diffusion is absent, fractional solute clearances alone are insufficient to characterize the permselective properties of the glomerular capillary wall, since these values may be altered by changes in glomerular pressures and flows as well as changes in the properties of the capillary wall per se.     

A patient-specific aortic valve model based on moving resistive immersed implicit surfaces

In this paper, we propose a full computational framework to simulate the hemodynamics in the aorta including the valve. Closed and open valve surfaces, as well as the lumen aorta, are reconstructed directly from medical images using new ad hoc algorithms, allowing a patient-specific simulation. The fluid dynamics problem that accounts from the movement of the valve is solved by a new 3D–0D fluid–structure interaction model in which the valve surface is implicitly represented through level set functions, yielding, in the Navier–Stokes equations, a resistive penalization term enforcing the blood to adhere to the valve leaflets. The dynamics of the valve between its closed and open position is modeled using a reduced geometric 0D model. At the discrete level, a finite element formulation is used and the SUPG stabilization is extended to include the resistive term in the Navier–Stokes equations. Then, after time discretization, the 3D fluid and 0D valve models are coupled through a staggered approach. This computational framework, applied to a patient-specific geometry and data, allows to simulate the movement of the valve, the sharp pressure jump occurring across the leaflets, and the blood flow pattern inside the aorta.     

Numerical simulation of blood flow in the left ventricle and aortic sinus using magnetic resonance imaging and computational fluid dynamics

Understanding cardiac blood flow patterns has many applications in analysing haemodynamics and for the clinical assessment of heart function. In this study, numerical simulations of blood flow in a patient-specific anatomical model of the left ventricle (LV) and the aortic sinus are presented. The realistic 3D geometry of both LV and aortic sinus is extracted from the processing of magnetic resonance imaging (MRI). Furthermore, motion of inner walls of LV and aortic sinus is obtained from cine-MR image analysis and is used as a constraint to a numerical computational fluid dynamics (CFD) model based on the moving boundary approach. Arbitrary Lagrangian–Eulerian finite element method formulation is used for the numerical solution of the transient dynamic equations of the fluid domain. Simulation results include detailed flow characteristics such as velocity, pressure and wall shear stress for the whole domain. The aortic outflow is compared with data obtained by phase-contrast MRI. Good agreement was found between simulation results and these measurements.     

Self-consistent modeling of radio-frequency plasma generation in stellarators

A self-consistent model of radio-frequency (RF) plasma generation in stellarators in the ion cyclotron frequency range is described. The model includes equations for the particle and energy balance and boundary conditions for Maxwell’s equations. The equation of charged particle balance takes into account the influx of particles due to ionization and their loss via diffusion and convection. The equation of electron energy balance takes into account the RF heating power source, as well as energy losses due to the excitation and electron-impact ionization of gas atoms, energy exchange via Coulomb collisions, and plasma heat conduction. The deposited RF power is calculated by solving the boundary problem for Maxwell’s equations. When describing the dissipation of the energy of the RF field, collisional absorption and Landau damping are taken into account. At each time step, Maxwell’s equations are solved for the current profiles of the plasma density and plasma temperature. The calculations are performed for a cylindrical plasma. The plasma is assumed to be axisymmetric and homogeneous along the plasma column. The system of balance equations is solved using the Crank-Nicholson scheme. Maxwell’s equations are solved in a one-dimensional approximation by using the Fourier transformation along the azimuthal and longitudinal coordinates. Results of simulations of RF plasma generation in the Uragan-2M stellarator by using a frame antenna operating at frequencies lower than the ion cyclotron frequency are presented. The calculations show that the slow wave generated by the antenna is efficiently absorbed at the periphery of the plasma column, due to which only a small fraction of the input power reaches the confinement region. As a result, the temperature on the axis of the plasma column remains low, whereas at the periphery it is substantially higher. This leads to strong absorption of the RF field at the periphery via the Landau mechanism.