Stellar radiative zones are typically assumed to be motionless in standard models of stellar structure but there is sound theoretical and observational evidence that this cannot be the case. We investigate by direct numerical simulations a three-dimensional and time-dependent model of stellar radiation zones consisting of an electrically-conductive and stably-stratified anelastic fluid confined to a rotating spherical shell and driven by a baroclinic torque. As the baroclinic driving is gradually increased a sequence of transitions from an axisymmetric and equatorially-symmetric time-independent flow to flows with a strong poloidal component and lesser symmetry are found. It is shown that all flow regimes characterised with significant non-axisymmetric components are capable of generating self-sustained magnetic field. As the value of the Prandtl number is decreased and the value of the Ekman number is decreased flows become strongly time-dependent with progressively complex spatial structure and dynamos can be generated at lower values of the magnetic Prandtl number.

The dynamics of stably stratified stellar radiative zones is of considerable interest due to the availability of increasingly detailed observations of Solar and stellar interiors. This article reports the first non-axisymmetric and time-dependent simulations of flows of anelastic fluids driven by baroclinic torques in stably stratified rotating spherical shells -- a system serving as an elemental model of a stellar radiative zone. With increasing baroclinicity a sequence of bifurcations from simpler to more complex flows is found in which some of the available symmetries of the problem are broken subsequently. The poloidal component of the flow grows relative to the dominant toroidal component with increasing baroclinicity. The possibility of magnetic field generation thus arises and this paper proceeds to provide some indications for self-sustained dynamo action in baroclinically-driven flows. We speculate that magnetic fields in stably stratified stellar interiors are thus not necessarily of fossil origin as it is often assumed.

Numerical MHD simulations play increasingly important role for understanding mechanisms of stellar magnetism. We present simulations of convection and dynamos in density-stratified rotating spherical fluid shells. We employ a new 3D simulation code for the solution of a physically consistent anelastic model of the process with a minimum number of parameters. The reported dynamo simulations extend into a "buoyancy-dominated" regime where the buoyancy forcing is dominant while the Coriolis force is no longer balanced by pressure gradients and strong anti-solar differential rotation develops as a result. We find that the self-generated magnetic fields, despite being relatively weak, are able to reverse the direction of differential rotation from anti-solar to solar-like. We also find that convection flows in this regime are significantly stronger in the polar regions than in the equatorial region, leading to non-oscillatory dipole-dominated dynamo solutions, and to concentration of magnetic field in the polar regions. We observe that convection has different morphology in the inner and at the outer part of the convection zone simultaneously such that organized geostrophic convection columns are hidden below a near-surface layer of well-mixed highly-chaotic convection. While we focus the attention on the buoyancy-dominated regime, we also demonstrate that conical differential rotation profiles and persistent regular dynamo oscillations can be obtained in the parameter space of the rotation-dominated regime even within this minimal model.

The onset of convection in a rotating cylindrical annulus with parallel ends filled with a compressible fluid is studied in the anelastic approximation. Thermal Rossby waves propagating in the azimuthal direction are found as solutions. The analogy to the case of Boussinesq convection in the presence of conical end surfaces of the annular region is emphasized. As in the latter case the results can be applied as an approximation for the description of the onset of anelastic convection in rotating spherical fluid shells. Reasonable agreement with three-dimensional numerical results published by Jones et al. (J. Fluid Mech., vol. 634, 2009, pp. 291-319) for the latter problem is found. As in those results the location of the onset of convection shifts outward from the tangent cylinder with increasing number $N_\rho$ of density scale heights until it reaches the equatorial boundary. A new result is that at a much higher number $N_\rho$ the onset location returns to the interior of the fluid shell.

A physically consistent model of magnetic field generation by convection in a rotating spherical shell with a minimum of parameters is applied to the Sun. Despite its unrealistic features the model exhibits a number of properties resembling those observed on the Sun. The model suggests that the large scale solar dynamo is dominated by a non-axisymmetric $m=1$ component of the magnetic field.

Simple models of magnetic field generation by convection in rotating spherical shells exhibit properties resembling those observed on the sun. The {assumption of the Boussinesq approximation made in these models} prevents a realistic description of the solar cycle, but through a physically motivated change in the boundary condition for the differential rotation the propagation of dynamo waves towards higher latitudes can be reversed at least at low latitudes.

We present results on modeling solar pores and sunspots using 2D axisymmetric magneto-hydrodynamic (MHD) simulations. These models are helpful for understanding the mechanisms of magnetic field concentration in sunspots, and the large-scale flow patterns associated with them. The simulations provide consistent, self-maintained, although not fully realistic, models of concentrated magnetic field near the solar surface. In this paper, we explore under which conditions the associated flows are converging or diverging near the surface. We find that in most cases in which a stable, pore-like concentration of magnetic field forms, a configuration with converging over diverging flow is established.

Possibilities and difficulties of applying the theory of magnetic field generation by convection flows in rotating spherical fluid shells to the Giant Planets are outlined. Recent progress in the understanding of the distribution of electrical conductivity in the Giant Planets suggests that the dynamo process occurs predominantly in regions of semiconductivity. In contrast to the geodynamo the magnetic field generation in the Giant Planets is thus characterized by strong radial conductivity variations. The importance of the constraint on the Ohmic dissipation provided by the planetary luminosity is emphasized. Planetary dynamos are likely to be of an oscillatory type, although these oscillations may not be evident from the exterior of the planets.

A numerical model of axisymmetric convection in the presence of a vertical magnetic flux bundle and rotation about the axis is presented. The model contains a compressible plasma described by the nonlinear MHD equations, with density and temperature gradients simulating the upper layer of the sun's convection zone. The solutions exhibit a central magnetic flux tube in a cylindrical numerical domain, with convection cells forming collar flows around the tube. When the numerical domain is rotated with a constant angular velocity, the plasma forms a Rankine vortex, with the plasma rotating as a rigid body where the magnetic field is strong, as in the flux tube, while experiencing sheared azimuthal flow in the surrounding convection cells, forming a free vortex. As a result, the azimuthal velocity component has its maximum value close to the outer edge of the flux tube. The azimuthal flow inside the magnetic flux tube and the vortex flow are prograde relative to the rotating cylindrical reference frame. A retrograde flow appears at the outer wall. The most significant convection cell outside the flux tube is the location for the maximum value of the azimuthal magnetic field component. The azimuthal flow and magnetic structure are not generated spontaneously, but decay exponentially in the absence of any imposed rotation of the cylindrical domain.

A convection-driven MHD dynamo in a rotating spherical shell, with clearly defined structural elements in the flow and magnetic field, is simulated numerically. Such dynamos can be called deterministic, in contrast to those explicitly dependent on the assumed properties of turbulence. The cases most interesting from the standpoint of studying the nature of stellar magnetism demonstrate the following features. On a global scale, the convective flows can maintain a ``general'' magnetic field with a sign-alternating dipolar component. Local (in many cases, bipolar) magnetic structures are associated with convection cells. Disintegrating local structures change into background fields, which drift toward the poles. From time to time, reversals of the magnetic fields in the polar regions occur, as ``new'' background fields expel the ``old'' fields.

It is shown that protostrange stars (PSSs) can be convective and that there are two possible scenarios describing their turbulence. Besides the local turbulence on the scale which is less than the mean free path of neutrinos, large-scale (~1 km) convection also may occur with properties that are similar to those of convection in protoneutron stars (PNSs). We thus suggest that strange stars can also create dynamo-originated magnetic fields during the deleptonization episode soon after a supernova explosion. Further detailed investigations are needed to see whether or not strange stars and neutron stars can be distinguished according to the differences in dynamo actions in strange quark matter and in neutron matter. The magnetic fields of strange stars and neutron stars may also behave very differently during the accretion-phase when the fields decay.