Effects of Stochastic Noises on Limit-Cycle Oscillations and Power Losses in Fusion Plasmas and Information Geometry.

We investigate the effects of different stochastic noises on the dynamics of the edge-localised modes (ELMs) in magnetically confined fusion plasmas by using a time-dependent PDF method, path-dependent information geometry (information rate, information length), and entropy-related measures (entropy production, mutual information). The oscillation quenching occurs due to either stochastic particle or magnetic perturbations, although particle perturbation is more effective in this amplitude diminishment compared with magnetic perturbations. On the other hand, magnetic perturbations are more effective at altering the oscillation period; the stochastic noise acts to increase the frequency of explosive oscillations (large ELMs) while decreasing the frequency of more regular oscillations (small ELMs). These stochastic noises significantly reduce power and energy losses caused by ELMs and play a key role in reproducing the observed experimental scaling relation of the ELM power loss with the input power. Furthermore, the maximum power loss is closely linked to the maximum entropy production rate, involving irreversible energy dissipation in non-equilibrium. Notably, over one ELM cycle, the information rate appears to keep almost a constant value, indicative of a geodesic. The information rate is also shown to be useful for characterising the statistical properties of ELMs, such as distinguishing between explosive and regular oscillations and the regulation between the pressure gradient and magnetic fluctuations.

Taylor-Couette and related flows on the centennial of Taylor's seminal Philosophical Transactions paper: part 2.

In 1923, the Philosophical Transactions published G. I. Taylor's seminal paper on the stability of what we now call Taylor-Couette flow. In the century since the paper was published, Taylor's ground-breaking linear stability analysis of fluid flow between two rotating cylinders has had an enormous impact on the field of fluid mechanics. The paper's influence has extended to general rotating flows, geophysical flows and astrophysical flows, not to mention its significance in firmly establishing several foundational concepts in fluid mechanics that are now broadly accepted. This two-part issue includes review articles and research articles spanning a broad range of contemporary research areas, all rooted in Taylor's landmark paper. This article is part of the theme issue 'Taylor-Couette and related flows on the centennial of Taylor's seminal Philosophical Transactions paper (Part 2)'.

Taylor-Couette and related flows on the centennial of Taylor's seminal Philosophical Transactions paper: part 1.

In 1923, the Philosophical Transactions published G. I. Taylor's seminal paper on the stability of what we now call Taylor-Couette flow. In the century since the paper was published, Taylor's ground-breaking linear stability analysis of fluid flow between two rotating cylinders has had an enormous impact on the field of fluid mechanics. The paper's influence has extended to general rotating flows, geophysical flows and astrophysical flows, not to mention its significance in firmly establishing several foundational concepts in fluid mechanics that are now broadly accepted. This two-part issue includes review articles and research articles spanning a broad range of contemporary research areas, all rooted in Taylor's landmark paper. This article is part of the theme issue 'Taylor-Couette and related flows on the centennial of Taylor's seminal Philosophical Transactions paper (part 1)'.

A stochastic model of edge-localized modes in magnetically confined plasmas.

Magnetically confined plasmas are far from equilibrium and pose considerable challenges in statistical analysis. We discuss a non-perturbative statistical method, namely a time-dependent probability density function (PDF) approach that is potentially useful for analysing time-varying, large, or non-Gaussian fluctuations and bursty events associated with instabilities in the low-to-high confinement transition and the H-mode. Specifically, we present a stochastic Langevin model of edge-localized modes (ELMs) by including stochastic noise terms in a previous ODE ELM model. We calculate exact time-dependent PDFs by numerically solving the Fokker-Planck equation and characterize time-varying statistical properties of ELMs for different energy fluxes and noise amplitudes. The stochastic noise is shown to introduce phase-mixing and plays a significant role in mitigating extreme bursts of large ELMs. Furthermore, based on time-dependent PDFs, we provide a path-dependent information geometric theory of the ELM dynamics and demonstrate its utility in capturing self-regulatory relaxation oscillations, bursts and a sudden change in the system. This article is part of a discussion meeting issue 'H-mode transition and pedestal studies in fusion plasmas'.

Information Geometry of Spatially Periodic Stochastic Systems.

We explore the effect of different spatially periodic, deterministic forces on the information geometry of stochastic processes. The three forces considered are f 0 = sin ( π x ) / π and f ± = sin ( π x ) / π ± sin ( 2 π x ) / 2 π , with f - chosen to be particularly flat (locally cubic) at the equilibrium point x = 0 , and f + particularly flat at the unstable fixed point x = 1 . We numerically solve the Fokker-Planck equation with an initial condition consisting of a periodically repeated Gaussian peak centred at x = µ , with µ in the range [ 0 , 1 ] . The strength D of the stochastic noise is in the range 10 - 4 - 10 - 6 . We study the details of how these initial conditions evolve toward the final equilibrium solutions and elucidate the important consequences of the interplay between an initial PDF and a force. For initial positions close to the equilibrium point x = 0 , the peaks largely maintain their shape while moving. In contrast, for initial positions sufficiently close to the unstable point x = 1 , there is a tendency for the peak to slump in place and broaden considerably before reconstituting itself at the equilibrium point. A consequence of this is that the information length L ∞ , the total number of statistically distinguishable states that the system evolves through, is smaller for initial positions closer to the unstable point than for more intermediate values. We find that L ∞ as a function of initial position µ is qualitatively similar to the force, including the differences between f 0 = sin ( π x ) / π and f ± = sin ( π x ) / π ± sin ( 2 π x ) / 2 π , illustrating the value of information length as a useful diagnostic of the underlying force in the system.

Information Geometry of Nonlinear Stochastic Systems.

We elucidate the effect of different deterministic nonlinear forces on geometric structure of stochastic processes by investigating the transient relaxation of initial PDFs of a stochastic variable x under forces proportional to -xn (n=3,5,7) and different strength D of δ-correlated stochastic noise. We identify the three main stages consisting of nondiffusive evolution, quasi-linear Gaussian evolution and settling into stationary PDFs. The strength of stochastic noise is shown to play a crucial role in determining these timescales as well as the peak amplitude and width of PDFs. From time-evolution of PDFs, we compute the rate of information change for a given initial PDF and uniquely determine the information length L(t) as a function of time that represents the number of different statistical states that a system evolves through in time. We identify a robust geodesic (where the information changes at a constant rate) in the initial stage, and map out geometric structure of an attractor as L(t→∞)∝µm, where µ is the position of an initial Gaussian PDF. The scaling exponent m increases with n, and also varies with D (although to a lesser extent). Our results highlight ubiquitous power-laws and multi-scalings of information geometry due to nonlinear interaction.

Time-Dependent Probability Density Functions and Attractor Structure in Self-Organised Shear Flows.

We report the time-evolution of Probability Density Functions (PDFs) in a toy model of self-organised shear flows, where the formation of shear flows is induced by a finite memory time of a stochastic forcing, manifested by the emergence of a bimodal PDF with the two peaks representing non-zero mean values of a shear flow. Using theoretical analyses of limiting cases, as well as numerical solutions of the full Fokker-Planck equation, we present a thorough parameter study of PDFs for different values of the correlation time and amplitude of stochastic forcing. From time-dependent PDFs, we calculate the information length ( L ), which is the total number of statistically different states that a system passes through in time and utilise it to understand the information geometry associated with the formation of bimodal or unimodal PDFs. We identify the difference between the relaxation and build-up of the shear gradient in view of information change and discuss the total information length ( L ∞ = L ( t → ∞ ) ) which maps out the underlying attractor structures, highlighting a unique property of L ∞ which depends on the trajectory/history of a PDF's evolution.

Dynamo Action in a Quasi-Keplerian Taylor-Couette Flow.

We numerically compute the flow of an electrically conducting fluid in a Taylor-Couette geometry where the rotation rates of the inner and outer cylinders satisfy Ω_{o}/Ω_{i}=(r_{o}/r_{i})^{-3/2}. In this quasi-Keplerian regime, a nonmagnetic system would be Rayleigh stable for all Reynolds numbers Re, and the resulting purely azimuthal flow incapable of kinematic dynamo action for all magnetic Reynolds numbers Rm. For Re = 10^{4} and Rm=10^{5}, we demonstrate the existence of a finite-amplitude dynamo, whereby a suitable initial condition yields mutually sustaining turbulence and magnetic fields, even though neither could exist without the other. This dynamo solution results in significantly increased outward angular momentum transport, with the bulk of the transport being by Maxwell rather than Reynolds stresses.

Geometric structure and information change in phase transitions.

We propose a toy model for a cyclic order-disorder transition and introduce a geometric methodology to understand stochastic processes involved in transitions. Specifically, our model consists of a pair of forward and backward processes (FPs and BPs) for the emergence and disappearance of a structure in a stochastic environment. We calculate time-dependent probability density functions (PDFs) and the information length L, which is the total number of different states that a system undergoes during the transition. Time-dependent PDFs during transient relaxation exhibit strikingly different behavior in FPs and BPs. In particular, FPs driven by instability undergo the broadening of the PDF with a large increase in fluctuations before the transition to the ordered state accompanied by narrowing the PDF width. During this stage, we identify an interesting geodesic solution accompanied by the self-regulation between the growth and nonlinear damping where the time scale τ of information change is constant in time, independent of the strength of the stochastic noise. In comparison, BPs are mainly driven by the macroscopic motion due to the movement of the PDF peak. The total information length L between initial and final states is much larger in BPs than in FPs, increasing linearly with the deviation γ of a control parameter from the critical state in BPs while increasing logarithmically with γ in FPs. L scales as |lnD| and D^{-1/2} in FPs and BPs, respectively, where D measures the strength of the stochastic forcing. These differing scalings with γ and D suggest a great utility of L in capturing different underlying processes, specifically, diffusion vs advection in phase transition by geometry. We discuss physical origins of these scalings and comment on implications of our results for bistable systems undergoing repeated order-disorder transitions (e.g., fitness).

Signature of nonlinear damping in geometric structure of a nonequilibrium process.

We investigate the effect of nonlinear interaction on the geometric structure of a nonequilibrium process. Specifically, by considering a driven-dissipative system where a stochastic variable x is damped either linearly (∝x) or nonlinearly (∝x^{3}) while driven by a white noise, we compute the time-dependent probability density functions (PDFs) during the relaxation towards equilibrium from an initial nonequilibrium state. From these PDFs, we quantify the information change by the information length L, which is the total number of statistically distinguishable states which the system passes through from the initial state to the final state. By exploiting different initial PDFs and the strength D of the white-noise forcing, we show that for a linear system, L increases essentially linearly with an initial mean value y_{0} of x as L∝y_{0}, demonstrating the preservation of a linear geometry. In comparison, in the case of a cubic damping, L has a power-law scaling as L∝y_{0}^{m}, with the exponent m depending on D and the width of the initial PDF. The rate at which information changes also exhibits a robust power-law scaling with time for the cubic damping.

Time-dependent probability density function in cubic stochastic processes.

We report time-dependent probability density functions (PDFs) for a nonlinear stochastic process with a cubic force using analytical and computational studies. Analytically, a transition probability is formulated by using a path integral and is computed by the saddle-point solution (instanton method) and a new nonlinear transformation of time. The predicted PDF p(x,t) in general involves a time integral, and useful PDFs with explicit dependence on x and t are presented in certain limits (e.g., in the short and long time limits). Numerical simulations of the Fokker-Planck equation provide exact time evolution of the PDFs and confirm analytical predictions in the limit of weak noise. In particular, we show that transient PDFs behave drastically differently from the stationary PDFs in regard to the asymmetry (skewness) and kurtosis. Specifically, while stationary PDFs are symmetric with the kurtosis smaller than 3, transient PDFs are skewed with the kurtosis larger than 3; transient PDFs are much broader than stationary PDFs. We elucidate the effect of nonlinear interaction on the strong fluctuations and intermittency in the relaxation process.

Geometric structure and geodesic in a solvable model of nonequilibrium process.

We investigate the geometric structure of a nonequilibrium process and its geodesic solutions. By employing an exactly solvable model of a driven dissipative system (generalized nonautonomous Ornstein-Uhlenbeck process), we compute the time-dependent probability density functions (PDFs) and investigate the evolution of this system in a statistical metric space where the distance between two points (the so-called information length) quantifies the change in information along a trajectory of the PDFs. In this metric space, we find a geodesic for which the information propagates at constant speed, and demonstrate its utility as an optimal path to reduce the total time and total dissipated energy. In particular, through examples of physical realizations of such geodesic solutions satisfying boundary conditions, we present a resonance phenomenon in the geodesic solution and the discretization into cyclic geodesic solutions. Implications for controlling population growth are further discussed in a stochastic logistic model, where a periodic modulation of the diffusion coefficient and the deterministic force by a small amount is shown to have a significant controlling effect.

Magnetic field evolution in magnetar crusts through three-dimensional simulations.

Current models of magnetars require extremely strong magnetic fields to explain their observed quiescent and bursting emission, implying that the field strength within the star's outer crust is orders of magnitude larger than the dipole component inferred from spin-down measurements. This presents a serious challenge to theories of magnetic field generation in a proto-neutron star. Here, we present detailed modeling of the evolution of the magnetic field in the crust of a neutron star through 3D simulations. We find that, in the plausible scenario of equipartition of energy between global-scale poloidal and toroidal magnetic components, magnetic instabilities transfer energy to nonaxisymmetric, kilometer-sized magnetic features, in which the local field strength can greatly exceed that of the global-scale field. These intense small-scale magnetic features can induce high-energy bursts through local crust yielding, and the localized enhancement of Ohmic heating can power the star's persistent emission. Thus, the observed diversity in magnetar behavior can be explained with mixed poloidal-toroidal fields of comparable energies.

A comparison of no-slip, stress-free and inviscid models of rapidly rotating fluid in a spherical shell.

We investigate how the choice of either no-slip or stress-free boundary conditions affects numerical models of rapidly rotating flow in Earth's core by computing solutions of the weakly-viscous magnetostrophic equations within a spherical shell, driven by a prescribed body force. For non-axisymmetric solutions, we show that models with either choice of boundary condition have thin boundary layers of depth E(1/2), where E is the Ekman number, and a free-stream flow that converges to the formally inviscid solution. At Earth-like values of viscosity, the boundary layer thickness is approximately 1 m, for either choice of condition. In contrast, the axisymmetric flows depend crucially on the choice of boundary condition, in both their structure and magnitude (either E(-1/2) or E(-1)). These very large zonal flows arise from requiring viscosity to balance residual axisymmetric torques. We demonstrate that switching the mechanical boundary conditions can cause a distinct change of structure of the flow, including a sign-change close to the equator, even at asymptotically low viscosity. Thus implementation of stress-free boundary conditions, compared with no-slip conditions, may yield qualitatively different dynamics in weakly-viscous magnetostrophic models of Earth's core. We further show that convergence of the free-stream flow to its asymptotic structure requires E ≤ 10(-5).

Nonaxisymmetric linear instability of cylindrical magnetohydrodynamic Taylor-Couette flow.

We consider the nonaxisymmetric modes of instability present in Taylor-Couette flow under the application of helical magnetic fields, mainly for magnetic Prandtl numbers close to the inductionless limit, and conduct a full examination of marginal stability in the resulting parameter space. We allow for the azimuthal magnetic field to be generated by a combination of currents in the inner cylinder and fluid itself and introduce a parameter governing the relation between the strength of these currents. A set of governing eigenvalue equations for the nonaxisymmetric modes of instability are derived and solved by spectral collocation with Chebyshev polynomials over the relevant parameter space, with the resulting instabilities examined in detail. We find that by altering the azimuthal magnetic field profiles the azimuthal magnetorotational instability, nonaxisymmetric helical magnetorotational instability, and Tayler instability yield interesting dynamics, such as different preferred mode types and modes with azimuthal wave number m>1. Finally, a comparison is given to the recent WKB analysis performed by Kirillov et al. [Kirillov, Stefani, and Fukumoto, J. Fluid Mech. 760, 591 (2014)JFLSA70022-112010.1017/jfm.2014.614] and its validity in the linear regime.

Experimental evidence for nonaxisymmetric magnetorotational instability in a rotating liquid metal exposed to an azimuthal magnetic field.

The azimuthal version of the magnetorotational instability (MRI) is a nonaxisymmetric instability of a hydrodynamically stable differentially rotating flow under the influence of a purely or predominantly azimuthal magnetic field. It may be of considerable importance for destabilizing accretion disks, and plays a central role in the concept of the MRI dynamo. We report the results of a liquid metal Taylor-Couette experiment that shows the occurrence of an azimuthal MRI in the expected range of Hartmann numbers.

Electromagnetically driven westward drift and inner-core superrotation in Earth's core.

A 3D numerical model of the earth's core with a viscosity two orders of magnitude lower than the state of the art suggests a link between the observed westward drift of the magnetic field and superrotation of the inner core. In our model, the axial electromagnetic torque has a dominant influence only at the surface and in the deepest reaches of the core, where it respectively drives a broad westward flow rising to an axisymmetric equatorial jet and imparts an eastward-directed torque on the solid inner core. Subtle changes in the structure of the internal magnetic field may alter not just the magnitude but the direction of these torques. This not only suggests that the quasi-oscillatory nature of inner-core superrotation [Tkalcic H, Young M, Bodin T, Ngo S, Sambridge M (2013) The shuffling rotation of the earth's inner core revealed by earthquake doublets. Nat Geosci 6:497-502.] may be driven by decadal changes in the magnetic field, but further that historical periods in which the field exhibited eastward drift were contemporaneous with a westward inner-core rotation. The model further indicates a strong internal shear layer on the tangent cylinder that may be a source of torsional waves inside the core.

Nonaxisymmetric magnetorotational instabilities in cylindrical Taylor-Couette flow.

We study the stability of cylindrical Taylor-Couette flow in the presence of azimuthal magnetic fields, and show that one obtains nonaxisymmetric magnetorotational instabilities, having azimuthal wave number m=1. For Omega{o}/Omega{i} only slightly greater than the Rayleigh value (r{i}/r{o}){2}, the critical Reynolds and Hartmann numbers are Re{c} approximately 10{3} and Ha{c} approximately 10{2}, independent of the magnetic Prandtl number Pm. These values are sufficiently small that it should be possible to obtain these instabilities in the Potsdam Rossendorf Magnetic Instability Experiment.

Helical magnetorotational instability in a Taylor-Couette flow with strongly reduced Ekman pumping.

The magnetorotational instability (MRI) is thought to play a key role in the formation of stars and black holes by sustaining the turbulence in hydrodynamically stable Keplerian accretion disks. In previous experiments the MRI was observed in a liquid metal Taylor-Couette flow at moderate Reynolds numbers by applying a helical magnetic field. The observation of this helical MRI (HMRI) was interfered with a significant Ekman pumping driven by solid end caps that confined the instability only to a part of the Taylor-Couette cell. This paper describes the observation of the HMRI in an improved Taylor-Couette setup with the Ekman pumping significantly reduced by using split end caps. The HMRI, which now spreads over the whole height of the cell, appears much sharper and in better agreement with numerical predictions. By analyzing various parameter dependencies we conclude that the observed HMRI represents a self-sustained global instability rather than a noise-sustained convective one.

Theory of current-driven instability experiments in magnetic Taylor-Couette flows.

We consider the linear stability of dissipative magnetic Taylor-Couette flow with imposed toroidal magnetic fields. The inner and outer cylinders can be either insulating or conducting; the inner one rotates, the outer one is stationary. The magnetic Prandtl number can be as small as 10(-5) , approaching realistic liquid-metal values. The magnetic field destabilizes the flow, except for radial profiles of B(phi)(R) close to the current-free solution. The profile with B(in)=B(out) (the most uniform field) is considered in detail. For weak fields the Taylor-Couette flow is stabilized, until for moderately strong fields the m=1 azimuthal mode dramatically destabilizes the flow again so that a maximum value for the critical Reynolds number exists. For sufficiently strong fields (as measured by the Hartmann number) the toroidal field is always unstable, even for the nonrotating case with Re=0 . The electric currents needed to generate the required toroidal fields in laboratory experiments are a few kA if liquid sodium is used, somewhat more if gallium is used. Weaker currents are needed for wider gaps, so a wide-gap apparatus could succeed even with gallium. The critical Reynolds numbers are only somewhat larger than the nonmagnetic values; hence such experiments would work with only modest rotation rates.