The possible role of meridional flows in suppressing magnetic buoyancy

The equation of motion for a toroidal flux ring in a stellar convective envelope is derived, and the equilibrium of such a ring is considered. Necessary conditions for the stability of toroidal flux rings are derived, and results of stability calculations for a particular model of the meridional flow are presented. The motions of the flux rings when the rings are far from their equilibrium position or when equilibrium does not exist are considered. The results confirm the linear stability analysis, and show that in the absence of stable equilibrium, the rings move toward the solar surface along a trajectory which is parallel to the rotation axis. It is expected that viscosity will tend to reduce the rotational velocity difference between the flux ring and its surroundings, thus reducing the Coriolis force and altering the equilibrium. The storage time of toroidal flux rings is estimated, and some implications for the sun are discussed.

On the Large-Scale Diffuse Magnetic Field of the Sun

Although the sunspots migrate towards the equator, the large-scale weak diffuse magnetic fields of the Sun migrate poleward with the solar cycle, the polar field reversing at the time of the sunspot maxima. We apply the vector model of Dikpati and Choudhuri (1994, Paper I) to fit these observations. The dynamo layer at the base of the convection zone is taken to be the source of the diffuse field, which is then evolved in the convection zone subject to meridional circulation and turbulent diffusion. We find that the longitudinally averaged observational data can be fitted reasonably well both for positive and negative values of the alpha-effect by adjusting the subsurface meridional flow suitably. The model will be extended in a future paper to include the decay of active regions as an extra source of the diffuse field, which may be necessary to explain the probable phase lag between B-tau and B-phi at lower latitudes.

Interpreting the upper level structure of the Madden-Julian oscillation

The nonlinear response of a spherical shallow water model to an imposed heat source in the presence of realistic zonal mean zonal winds is investigated numerically. The solutions exhibit elongated, meridionally tilted ridges and troughs indicative of a poleward dispersion of wave activity. As the speed of the jets is increased, the equatorial Kelvin wave is unaffected but the global Rossby wave train coalesces to form a compact, amplified quadrupole structure that bears a striking resemblance to the observed upper level structure of the Madden-Julian oscillation. In the presence of strong subtropical westerly jets, the advection of planetary vorticity by the meridional flow and relative vorticity by the zonally averaged background flow conspire to create the distinctive quadrupole configuration of flanking Rossby waves.

Development of a turning mid turbine frame with embedded design-part II: Unsteady measurements

© 2014 by ASME. This paper, the second of two parts, presents a new setup for the two-stage two-spool facility located at the Institute for Thermal Turbomachinery and Machine Dynamics (ITTM) of Graz University of Technology. The rig was designed to reproduce the flow behavior of a transonic turbine followed by a counter-rotating low pressure stage such as those in high bypass aero-engines. The meridional flow path of the machine is characterized by a diffusing S-shaped duct between the two rotors. The role of wide chord vanes placed into the mid turbine frame is to lead the flow towards the low pressure (LP) rotor with appropriate swirl. Experimental and numerical investigations performed on this setup showed that the wide chord struts induce large wakes and extended secondary flows at the LP inlet flow. Moreover, large deterministic fluctuations of pressure, which may cause noise and blade vibrations, were observed downstream of the LP rotor. In order to minimize secondary vortices and to damp the unsteady interactions, the mid turbine frame was redesigned to locate two zero-lift splitters into each vane passage. While in the first part of the paper the design process of the splitters and the time-averaged flow field were presented, in this second part the measurements performed by means of a fast response probe will support the explanation of the time-resolved field. The discussion will focus on the comparison between the baseline case (without splitters) and the embedded design.

Development of a turning mid turbine frame with embedded design-part I: Design and steady measurements

© 2014 by ASME. The paper presents a new setup for the two-stage two-spool facility located at the Institute for Thermal Turbomachinery and Machine Dynamics (ITTM) of Graz University of Technology. The rig was designed in order to simulate the flow behavior of a transonic turbine followed by a counter-rotating low pressure (LP) stage like the spools of a modern high bypass aeroengine. The meridional flow path of the machine is characterized by a diffusing S-shaped duct between the two rotors. The role of turning struts placed into the mid turbine frame is to lead the flow towards the LP rotor with appropriate swirl. Experimental and numerical investigations performed on the setup over the last years, which were used as baseline for this paper, showed that wide chord vanes induce large wakes and extended secondary flows at the LP rotor inlet flow. Moreover, unsteady interactions between the two turbines were observed downstream of the LP rotor. In order to increase the uniformity and to decrease the unsteady content of the flow at the inlet of the LP rotor, the mid turbine frame was redesigned with two zero-lifting splitters embedded into the strut passage. In this first part of the paper the design process of the splitters and its critical points are presented, while the time-averaged field is discussed by means of five-hole probe measurements and oil flow visualizations. The comparison between the baseline case and the embedded design configuration shows that the new design is able to reduce the flow gradients downstream of the turning struts, providing a more suitable inlet condition for the low pressure rotor. The improvement in the flow field uniformity is also observed downstream of the turbine and it is, consequently, reflected in an enhancement of the LP turbine performance. In the second part of this paper the influence of the embedded design on the time-resolved field is investigated.

Les oscillations torsionnelles dans la zone de convection solaire

Nous analysons les oscillations torsionnelles se développant dans une simulation magnétohydrodynamique de la zone de convection solaire produisant des champs magnétiques de type solaire (champs axisymétriques subissant des inversions de polarités régulières sur des échelles temporelles décadaires). Puisque ces oscillations sont également similaires à celles observées dans le Soleil, nous analysons les dynamiques zonales aux grandes échelles. Nous séparons donc les termes aux grandes échelles (force de Coriolis exercée sur la circulation méridienne et les champs magnétiques aux grandes échelles) de ceux aux petites échelles (les stress de Reynolds et de Maxwell). En comparant les flux de moments cinétiques entre chacune des composantes, nous nous apercevons que les oscillations torsionnelles sont maintenues par l’écoulement méridien aux grandes échelles, lui même modulé par les champs magnétiques. Une analyse d’échange d’énergie confirme ce résultat, puisqu’elle montre que seul le terme comprenant la force de Coriolis injecte de l’énergie dans l’écoulement. Une analyse de la dynamique rotationnelle ayant lieu à la limite de la zone stable et de la zone de convection démontre que celle-ci est fortement modifiée lors du passage de la base des couches convectives à la base de la fine tachocline s’y formant juste en-dessous. Nous concluons par une discussion au niveau du mécanisme de saturation en amplitude dans la dynamo s’opérant dans la simulation ainsi que de la possibilité d’utiliser les oscillations torsionnelles comme précurseurs aux cycles solaires à venir.

A Coupled 2 × 2D Babcock–Leighton Solar Dynamo Model. I. Surface Magnetic Flux Evolution

The need for reliable predictions of the solar activity cycle motivates the development of dynamo models incorporating a representation of surface processes sufficiently detailed to allow assimilation of magnetographic data. In this series of papers we present one such dynamo model, and document its behavior and properties. This first paper focuses on one of the model's key components, namely surface magnetic flux evolution. Using a genetic algorithm, we obtain best-fit parameters of the transport model by least-squares minimization of the differences between the associated synthetic synoptic magnetogram and real magnetographic data for activity cycle 21. Our fitting procedure also returns Monte Carlo-like error estimates. We show that the range of acceptable surface meridional flow profiles is in good agreement with Doppler measurements, even though the latter are not used in the fitting process. Using a synthetic database of bipolar magnetic region (BMR) emergences reproducing the statistical properties of observed emergences, we also ascertain the sensitivity of global cycle properties, such as the strength of the dipole moment and timing of polarity reversal, to distinct realizations of BMR emergence, and on this basis argue that this stochasticity represents a primary source of uncertainty for predicting solar cycle characteristics.

The evolution of the sun’s open magnetic flux: II. full solar cycle simulations

In this paper the origin and evolution of the Sun’s open magnetic flux is considered by conducting magnetic flux transport simulations over many solar cycles. The simulations include the effects of differential rotation, meridional flow and supergranular diffusion on the radial magnetic field at the surface of the Sun as new magnetic bipoles emerge and are transported poleward. In each cycle the emergence of roughly 2100 bipoles is considered. The net open flux produced by the surface distribution is calculated by constructing potential coronal fields with a source surface from the surface distribution at regular intervals. In the simulations the net open magnetic flux closely follows the total dipole component at the source surface and evolves independently from the surface flux. The behaviour of the open flux is highly dependent on meridional flow and many observed features are reproduced by the model. However, when meridional flow is present at observed values the maximum value of the open flux occurs at cycle minimum when the polar caps it helps produce are the strongest. This is inconsistent with observations by Lockwood, Stamper and Wild (1999) and Wang, Sheeley, and Lean (2000) who find the open flux peaking 1–2 years after cycle maximum. Only in unrealistic simulations where meridional flow is much smaller than diffusion does a maximum in open flux consistent with observations occur. It is therefore deduced that there is no realistic parameter range of the flux transport variables that can produce the correct magnitude variation in open flux under the present approximations. As a result the present standard model does not contain the correct physics to describe the evolution of the Sun’s open magnetic flux over an entire solar cycle. Future possible improvements in modeling are suggested.

THE ROLE OF DIFFUSIVITY QUENCHING IN FLUX-TRANSPORT DYNAMO MODELS

In the nonlinear phase of a dynamo process, the back-reaction of the magnetic field upon the turbulent motion results in a decrease of the turbulence level and therefore in a suppression of both the magnetic field amplification (the alpha-quenching effect) and the turbulent magnetic diffusivity (the eta-quenching effect). While the former has been widely explored, the effects of eta-quenching in the magnetic field evolution have rarely been considered. In this work, we investigate the role of the suppression of diffusivity in a flux-transport solar dynamo model that also includes a nonlinear alpha-quenching term. Our results indicate that, although for alpha-quenching the dependence of the magnetic field amplification with the quenching factor is nearly linear, the magnetic field response to eta-quenching is nonlinear and spatially nonuniform. We have found that the magnetic field can be locally amplified in this case, forming long-lived structures whose maximum amplitude can be up to similar to 2.5 times larger at the tachocline and up to similar to 2 times larger at the center of the convection zone than in models without quenching. However, this amplification leads to unobservable effects and to a worse distribution of the magnetic field in the butterfly diagram. Since the dynamo cycle period increases when the efficiency of the quenching increases, we have also explored whether the eta-quenching can cause a diffusion-dominated model to drift into an advection-dominated regime. We have found that models undergoing a large suppression in eta produce a strong segregation of magnetic fields that may lead to unsteady dynamo-oscillations. On the other hand, an initially diffusion-dominated model undergoing a small suppression in eta remains in the diffusion-dominated regime.

Upwelling in the Eastern Subtropical North Atlantic Ocean

[EN]Coastal upwelling in the eastern margin and offshore curl-driven upwelling in the southeastern margin, make the subtropical Northeast Atlantic a region of major primary productivity. When examining a broad zonal area, from the coast to 40_W, we find that the upward transport of nutrients due to offshore curl-driven upwelling becomes the main control on productivity. Nevertheless, despite its relatively small zonal extension of about 100 km, coastal upwelling extends its impact towards the open ocean through offshore Ekman transport and convergence of the meridional flow at Cape Blanc (21_N).

The solar dynamo with meridional circulation

We show that meridional circulation can have a profound influence on dynamo models for the solar cycle. Motivated by the observed tilt angles of sunspot groups we assume that the generation of the poloidal field takes place near the surface, while a shear layer of radial differential rotation produces the toroidal field at the bottom of the convection zone. Both layers are coupled by a circulation with a poleward directed flow in the upper part and an equatorward flow in the deep layers of the convection zone. The circulation forces the toroidal field belts (which are responsible for the surface activity) to move equatorward. This leads to butterfly diagrams in qualitative agreement with the observations, even if the dynamo wave would propagate poleward in the absence of circulation. This result opens the possibility to construct models for the solar cycle which are based on observational data (tilt angles, differential rotation, and meridional circulation).

Unsteady MUD rotating flow over a rotating sphere near the equator

The unsteady rotating flow of a laminar incompressible viscous electrically conducting fluid over a rotating sphere in the vicinity of the equator has been studied. The fluid and the body rotate either in the same direction or in opposite directions. The effects of surface suction and magnetic field have been included in the analysis. There is an initial steady state that is perturbed by a sudden change in the rotational velocity of the sphere, and this causes unsteadiness in the flow field. The nonlinear coupled parabolic partial differential equations governing the boundary-layer flow have been solved numerically by using an implicit finite-difference scheme. For large suction or magnetic field, analytical solutions have also been obtained. The magnitude of the radial, meridional and rotational velocity components is found to be higher when the fluid and the body rotate in opposite directions than when they rotate in the same direction. The surface shear stresses in the meridional and rotational directions change sign when the ratio of the angular velocities of the sphere and the fluid lambda greater than or equal to lambda(0). The final (new) steady state is reached rather quickly which implies that the spin-up time is small. The magnetic field and surface suction reduce the meridional shear stress, but increase the surface shear stress in the rotational direction.

IS A DEEP ONE-CELL MERIDIONAL CIRCULATION ESSENTIAL FOR THE FLUX TRANSPORT SOLAR DYNAMO?

The solar activity cycle is successfully modeled by the flux transport dynamo, in which the meridional circulation of the Sun plays an important role. Most of the kinematic dynamo simulations assume a one-cell structure of the meridional circulation within the convection zone, with the equatorward return flow at its bottom. In view of the recent claims that the return flow occurs at a much shallower depth, we explore whether a meridional circulation with such a shallow return flow can still retain the attractive features of the flux transport dynamo (such as a proper butterfly diagram, the proper phase relation between the toroidal and poloidal fields). We consider additional cells of the meridional circulation below the shallow return flow-both the case of multiple cells radially stacked above one another and the case of more complicated cell patterns. As long as there is an equatorward flow in low latitudes at the bottom of the convection zone, we find that the solar behavior is approximately reproduced. However, if there is either no flow or a poleward flow at the bottom of the convection zone, then we cannot reproduce solar behavior. On making the turbulent diffusivity low, we still find periodic behavior, although the period of the cycle becomes unrealistically large. In addition, with a low diffusivity, we do not get the observed correlation between the polar field at the sunspot minimum and the strength of the next cycle, which is reproduced when diffusivity is high. On introducing radially downward pumping, we get a more reasonable period and more solar-like behavior even with low diffusivity.

Numerical study of bifurcation solutions of spherical Taylor-Couette flow

The steady bifurcation flows in a spherical gap (gap ratio sigma=0.18) with rotating inner and stationary outer spheres are simulated numerically for Re(c1)less than or equal to Re less than or equal to 1 500 by solving steady axisymmetric incompressible Navier-Stokes equations using a finite difference method. The simulation shows that there exist two steady stable flows with 1 or 2 vortices per hemisphere for 775 less than or equal to Re less than or equal to 1 220 and three steady stable flows with 0, 1, or 2 vortices for 1 220flows at the same Reynolds number is related with different initial conditions which on be generated by different accelerations of the inner sphere. Generation of zero-or two-vortex flow depends mainly on the acceleratio n, but that of one-vortex flow also depends on the perturbation breaking the equatorial symmetry. The mechanism of development of a saddle point in the meridional plane at higher Re number and its role in the formation of two-vortex flow are analyzed.