Nonlinearity in a dynamo

Using a rotating flat layer heated from below as an example, we consider effects which lead to stabilizing an exponentially growing magnetic field in magnetostrophic convection in transition from the kinematic dynamo to the full non-linear dynamo. We present estimates of the energy redistribution over the spectrum and helicity quenching by the magnetic field. We also study the alignment of the velocity and magnetic fields. These regimes are similar to those in planetary dynamo simulations.Comment: Accepted to Geophys. Astrophys. Fluid Dyna

Convective plan-form two-scale dynamos in a plane layer

We study generation of magnetic fields, involving large spatial scales, by convective plan-forms in a horizontal layer. Magnetic modes and their growth rates are expanded in power series in the scale ratio, and the magnetic eddy diffusivity (MED) tensor is derived for flows, symmetric about the vertical axis in a layer. For convective rolls magnetic eddy correction is demonstrated to be always positive. For rectangular cell patterns, the region in the parameter space of negative MED coincides with that of small-scale magnetic field generation. No instances of negative MED in hexagonal cells are found. A family of plan-forms with a smaller symmetry group than that of rectangular cell patterns has been found numerically, where MED is negative for molecular magnetic diffusivity over the threshold for the onset of small-scale magnetic field generation.Comment: Latex. 24 pages with 3 Postscript figures, 19 references. Final version (expanded Appendix 2, 4 references added, notation changed to a more "user-friendly"), accepted in Geophysical and Astrophysical Fluid Dynamic

On the effects of turbulence on a screw dynamo

In an experiment in the Institute of Continuous Media Mechanics in Perm (Russia) an non--stationary screw dynamo is intended to be realized with a helical flow of liquid sodium in a torus. The flow is necessarily turbulent, that is, may be considered as a mean flow and a superimposed turbulence. In this paper the induction processes of the turbulence are investigated within the framework of mean--field electrodynamics. They imply of course a part which leads to an enhanced dissipation of the mean magnetic field. As a consequence of the helical mean flow there are also helical structures in the turbulence. They lead to some kind of $\alpha$--effect, which might basically support the screw dynamo. The peculiarity of this $\alpha$--effect explains measurements made at a smaller version of the device envisaged for the dynamo experiment. The helical structures of the turbulence lead also to other effects, which in combination with a rotational shear are potentially capable of dynamo action. A part of them can basically support the screw dynamo. Under the conditions of the experiment all induction effects of the turbulence prove to be rather weak in comparison to that of the main flow. Numerical solutions of the mean--field induction equation show that all the induction effects of the turbulence together let the screw dynamo threshold slightly, at most by one per cent, rise. The numerical results give also some insights into the action of the individual induction effects of the turbulence.Comment: 15 pages, 7 figures, in GAFD prin

Towards a finite-time singularity of the Navier-Stokes equations Part 1. Derivation and analysis of dynamical system

© 2018 Cambridge University Press. The evolution towards a finite-time singularity of the Navier-Stokes equations for flow of an incompressible fluid of kinematic viscosity is studied, starting from a finite-energy configuration of two vortex rings of circulation and radius , symmetrically placed on two planes at angles to a plane of symmetry . The minimum separation of the vortices, , and the scale of the core cross-section, , are supposed to satisfy the initial inequalities , and the vortex Reynolds number is supposed very large. It is argued that in the subsequent evolution, the behaviour near the points of closest approach of the vortices (the 'tipping points') is determined solely by the curvature at the tipping points and by and , where is a dimensionless time variable. The Biot-Savart law is used to obtain analytical expressions for the rate of change of these three variables, and a nonlinear dynamical system relating them is thereby obtained. The solution shows a finite-time singularity, but the Biot-Savart law breaks down just before this singularity is realised, when and become of order unity. The dynamical system admits 'partial Leray scaling' of just and , and ultimately full Leray scaling of and , conditions for which are obtained. The tipping point trajectories are determined; these meet at the singularity point at a finite angle. An alternative model is briefly considered, in which the initial vortices are ovoidal in shape, approximately hyperbolic near the tipping points, for which there is no restriction on the initial value of the parameter ; however, it is still the circles of curvature at the tipping points that determine the local evolution, so the same dynamical system is obtained, with breakdown again of the Biot-Savart approach just before the incipient singularity is realised. The Euler flow situation is considered, and it is conjectured on the basis of the above dynamical system that a finite-time singularity can indeed occur in this case

Helicity within the vortex filament model

Kinetic helicity is one of the invariants of the Euler equations that is associated with the topology of vortex lines within the fluid. In superfluids, the vorticity is concentrated along vortex filaments. In this setting, helicity would be expected to acquire its simplest form. However, the lack of a core structure for vortex filaments appears to result in a helicity that does not retain its key attribute as a quadratic invariant. By defining a spanwise vector to the vortex through the use of a Seifert framing, we are able to introduce twist and henceforth recover the key properties of helicity. We present several examples for calculating internal twist to illustrate why the centreline helicity alone will lead to ambiguous results if a twist contribution is not introduced. Our choice of the spanwise vector can be expressed in terms of the tangential component of velocity along the filament. Since the tangential velocity does not alter the configuration of the vortex at later times, we are able to recover a similar equation for the internal twist angle to that of classical vortex tubes. Our results allow us to explain how a quasi-classical limit of helicity emerges from helicity considerations for individual superfluid vortex filaments

The gradient of potential vorticity, quaternions and an orthonormal frame for fluid particles

The gradient of potential vorticity (PV) is an important quantity because of the way PV (denoted as $q$) tends to accumulate locally in the oceans and atmospheres. Recent analysis by the authors has shown that the vector quantity \bdB = \bnabla q\times \bnabla\theta for the three-dimensional incompressible rotating Euler equations evolves according to the same stretching equation as for \bom the vorticity and \bB, the magnetic field in magnetohydrodynamics (MHD). The \bdB-vector therefore acts like the vorticity \bom in Euler's equations and the \bB-field in MHD. For example, it allows various analogies, such as stretching dynamics, helicity, superhelicity and cross helicity. In addition, using quaternionic analysis, the dynamics of the \bdB-vector naturally allow the construction of an orthonormal frame attached to fluid particles\,; this is designated as a quaternion frame. The alignment dynamics of this frame are particularly relevant to the three-axis rotations that particles undergo as they traverse regions of a flow when the PV gradient \bnabla q is large.Comment: Dedicated to Raymond Hide on the occasion of his 80th birthda

Dynamics of the Tippe Top -- properties of numerical solutions versus the dynamical equations

We study the relationship between numerical solutions for inverting Tippe Top and the structure of the dynamical equations. The numerical solutions confirm oscillatory behaviour of the inclination angle $\theta(t)$ for the symmetry axis of the Tippe Top. They also reveal further fine features of the dynamics of inverting solutions defining the time of inversion. These features are partially understood on the basis of the underlying dynamical equations

Writhe in the Stretch-Twist-Fold Dynamo

This is an Author's Original Manuscript of an article whose final and definitive form, the Version of Record, has been published in Geophysical and Astrophysical Fluid Dynamics (2008) Copyright © 2008 Taylor & Francis, available online at: http://www.tandfonline.com/10.1080/03091920802531791This article looks at the influence of writhe in the stretch-twist-fold dynamo. We consider a thin flux tube distorted by simple stretch, twist, and fold motions and calculate the helicity and energy spectra. The writhe number assists in the calculations, as it tells us how much the internal twist changes as the tube is distorted. In addition it provides a valuable diagnostic for the degree of distortion. Non mirror-symmetric dynamos typically generate magnetic helicity of one sign on large-scales and the opposite sign on small scales. The calculations presented here confirm the hypothesis that the large-scale helicity corresponds to writhe and the small scale corresponds to twist. In addition, the writhe helicity spectrum exhibits an interesting oscillatory behavior. The technique of calculating Fourier spectra for the writhe helicity may be useful in other areas of research, for example, the study of highly coiled molecules

Similarity solutions for unsteady shear-stress-driven flow of Newtonian and power-law fluids : slender rivulets and dry patches

Unsteady flow of a thin film of a Newtonian fluid or a non-Newtonian power-law fluid with power-law index N driven by a constant shear stress applied at the free surface, on a plane inclined at an angle α to the horizontal, is considered. Unsteady similarity solutions representing flow of slender rivulets and flow around slender dry patches are obtained. Specifically, solutions are obtained for converging sessile rivulets (0 < α < π/2) and converging dry patches in a pendent film (π/2 < α < π), as well as for diverging pendent rivulets and diverging dry patches in a sessile film. These solutions predict that at any time t, the rivulet and dry patch widen or narrow according to |x|3/2, and the film thickens or thins according to |x|, where x denotes distance down the plane, and that at any station x, the rivulet and dry patch widen or narrow like |t|−1, and the film thickens or thins like |t|−1, independent of N