Recent analyses of the helioseismic data have produced evidence for a variety of interesting dynamical behaviour associated with torsional oscillations. What is not so far clear is whether these oscillations extend all the way to the bottom of the convection zone and, if so, whether the oscillatory behaviour at the top and the bottom of the convection zone is different. Attempts have been made to understand such modes of behaviour within the framework of nonlinear dynamo models which include the nonlinear action of the Lorentz force of the dynamo generated magnetic field on the solar angular velocity. One aspect of these models that remains uncertain is the nature of the boundary conditions on the magnetic field. Here by employing a range of physically plausible boundary conditions, we show that for near-critical and moderately supercritical dynamo regimes, the oscillations extend all the way down to the bottom of the convection zone. Thus, such penetration is an extremely robust feature of the models considered. We also find parameter ranges for which the supercritical models show spatiotemporal fragmentation for a range of choices of boundary conditions. Given their observational importance, we also make a comparative study of the amplitude of torsional oscillations as a function of the boundary conditions.

We find concrete evidence for a recently discovered form of intermittency, referred to as in--out intermittency, in both PDE and ODE models of mean field dynamos. This type of intermittency (introduced in Ashwin et al 1999) occurs in systems with invariant submanifolds and, as opposed to on--off intermittency which can also occur in skew product systems, it requires an absence of skew product structure. By this we mean that the dynamics on the attractor intermittent to the invariant manifold cannot be expressed simply as the dynamics on the invariant subspace forcing the transverse dynamics; the transverse dynamics will alter that tangential to the invariant subspace when one is far enough away from the invariant manifold. Since general systems with invariant submanifolds are not likely to have skew product structure, this type of behaviour may be of physical relevance in a variety of dynamical settings. The models employed here to demonstrate in--out intermittency are axisymmetric mean--field dynamo models which are often used to study the observed large scale magnetic variability in the Sun and solar-type stars. The occurrence of this type of intermittency in such models may be of interest in understanding some aspects of such variabilities.

Recent analysis of the helioseismic observations indicate that the previously observed surface torsional oscillations extend significantly downwards into the solar convection zone. In an attempt to understand these oscillations, we study the nonlinear coupling between the magnetic field and the solar differential rotation in the context of a mean field dynamo model, in which the nonlinearity is due to the action of the azimuthal component of the Lorentz force of the dynamo generated magnetic field on the solar angular velocity. The underlying zero order angular velocity is chosen to be consistent with the most recent helioseismic data. The model produces butterfly diagrams which are in qualitative agreement with the observations. It displays torsional oscillations that penetrate into the convection zone, and which with time migrate towards the equator. The period of these oscillations is found to be half that of the period of the global magnetic fields. This is compatible with the observed period of the surface torsional oscillations. Inside the convection zone, this is a testable prediction that is not ruled out by the observations so far available.

We study the changes in the dynamical behaviour of axisymmetric spherical mean-field dynamo models produced by changes in their geometry and topology, by considering a two parameter family of models, ranging from a full sphere to spherical shell, torus and disc-like configurations, within a unified framework. We find that the two parameter space of the family of models considered here separates into at least three different regions with distinct characteristics for the onset of dynamo action. In two of these regions, the most easily excited fields are oscillatory, in one case with dipolar symmetry, and in the other with quadrupolar, whereas in the third region the most easily excited field is steady and quadrupolar. In the nonlinear regime, we find that topological changes can alter significantly the dynamical behaviour, whilst modest changes in geometry can produce qualitative changes, particularly for thin disc-like configurations. This is of potential importance, since the exact shapes of astrophysical bodies, especially accretion discs and galaxies, are usually not precisely known.

Calculations for mean field dynamo models (in both full spheres and spherical shells), with both algebraic and dynamic $\alpha$--quenchings, show qualitative as well as quantitative differences and similarities in the dynamical behaviour of these models. We summarise and enhance recent results with extra examples. Overall, the effect of using a dynamic $\alpha$ appears to be complicated and is affected by the region of parameter space examined.

We study axisymmetric mean field spherical and spherical shell dynamo models, with both dynamic and algebraic $\alpha$--quenchings. Our results show that there are qualitative as well as quantitative differences and similarities between these models. Regarding similarities, both groups of models exhibit symmetric, antisymmetric and mixed modes of behaviour. As regards differences, the important feature in the full sphere models is the occurrence of chaotic behaviour in the algebraic $\alpha$--quenching models. For the spherical shell models with dynamic $\alpha$ the main features include the possibility of multi-attractor regimes with final state sensitivity with respect to small changes in the magnitude of $\alpha$ and the initial parity. We find the effect of introducing a dynamic $\alpha$ is likely to be complicated and depend on the region of the parameter space considered, rather than a uniform change towards simplicity or complexity.

We investigate the behaviour of $\alpha\Omega$ dynamos with a dynamic $\alpha$, whose evolution is governed by the imbalance between a driving and a damping term. We focus on truncated versions of such dynamo models which are often studied in connection with solar and stellar variability. Given the approximate nature of such models, it is important to study how robust they are with respect to reasonable changes in the formulation of the driving and damping terms. For each case, we also study the effects of changes of the dynamo number and its sign, the truncation order and initial conditions. Our results show that changes in the formulation of the driving term have important consequences for the dynamical behaviour of such systems, with the detailed nature of these effects depending crucially on the form of the driving term assumed, the value and the sign of the dynamo number and the initial conditions. On the other hand, the change in the damping term considered here seems to produce little qualitative effect.

In a recent work (Covas et al. 1997), the behaviour and the robustness of truncated $\alpha\Omega$ dynamos with a dynamic $\alpha$ were studied with respect to a number of changes in the driving term of the dynamic $\alpha$ equation, which was considered previously by Schmalz and Stix (1991) to be of the form $\sim A_{\phi}B_{\phi}$. Here we review and extend our previous work and consider the effect of adding a quadratic quenching term of the form $\alpha|B|^2$. We find that, as before, such a change can have significant effects on the dynamics of the related truncated systems. We also find intervals of (negative) dynamo numbers, in the system considered by Schmalz and Stix (1991), for which there is sensitivity with respect to small changes in the dynamo number and the initial conditions, similar to what was found in our previous work. This latter behaviour may be of importance in producing the intermittent type of behaviour observed in the sun.