Solar rotation rate and its gradients during cycle 23

Available helioseismic data now span almost the entire solar activity cycle 23 making it possible to study solar-cycle related changes of the solar rotation rate in detail. In this paper we study how the solar rotation rate, in particular, the zonal flows change with time. In addition to the zonal flows that show a well known pattern in the solar convection zone, we also study changes in the radial and latitudinal gradients of the rotation rate, particularly in the shear layer that is present in the immediate sub-surface layers of the Sun. In the case of the zonal-flow pattern, we find that the band indicating fast rotating region close to the equator seems to have bifurcated around 2005. Our investigation of the rotation-rate gradients show that the relative variation in the rotation-rate gradients is about 20% or more of their average values, which is much larger than the relative variation in the rotation rate itself. These results can be used to test predictions of various solar dynamo models.Comment: To appear in ApJ. Fig 5 has been corrected in this versio

Diffusive transport in networks built of containers and tubes

We developed analytical and numerical methods to study a transport of non-interacting particles in large networks consisting of M d-dimensional containers C_1,...,C_M with radii R_i linked together by tubes of length l_{ij} and radii a_{ij} where i,j=1,2,...,M. Tubes may join directly with each other forming junctions. It is possible that some links are absent. Instead of solving the diffusion equation for the full problem we formulated an approach that is computationally more efficient. We derived a set of rate equations that govern the time dependence of the number of particles in each container N_1(t),N_2(t),...,N_M(t). In such a way the complicated transport problem is reduced to a set of M first order integro-differential equations in time, which can be solved efficiently by the algorithm presented here. The workings of the method have been demonstrated on a couple of examples: networks involving three, four and seven containers, and one network with a three-point junction. Already simple networks with relatively few containers exhibit interesting transport behavior. For example, we showed that it is possible to adjust the geometry of the networks so that the particle concentration varies in time in a wave-like manner. Such behavior deviates from simple exponential growth and decay occurring in the two container system.Comment: 21 pages, 18 figures, REVTEX4; new figure added, reduced emphasis on graph theory, additional discussion added (computational cost, one dimensional tubes

Changes in Solar Dynamics from 1995 to 2002

Data obtained by the GONG and MDI instruments over the last seven years are used to study how solar dynamics -- both rotation and other large scale flows -- have changed with time. In addition to the well known phenomenon of bands of faster and slower rotation moving towards the equator and pole, we find that the zonal flow pattern rises upwards with time. Like the zonal flows, the meridional flows also show distinct solar activity related changes. In particular, the anti-symmetric component of the meridional flow shows a decrease in speed with activity. We do not see any significant temporal variations in the dynamics of the tachocline region where the solar dynamo is believed to be operating.Comment: To appear in ApJ, March 1 200

Exact, E=0, Solutions for General Power-Law Potentials. I. Classical Orbits

For zero energy, $E=0$, we derive exact, classical solutions for {\em all} power-law potentials, $V(r)=-\gamma/r^\nu$, with $\gamma>0$ and $-\infty <\nu<\infty$. When the angular momentum is non-zero, these solutions lead to the orbits $\r(t)= [\cos \mu (\th(t)-\th_0(t))]^{1/\mu}$, for all $\mu \equiv \nu/2-1 \ne 0$. When $\nu>2$, the orbits are bound and go through the origin. This leads to discrete discontinuities in the functional dependence of $\th(t)$ and $\th_0(t)$, as functions of $t$, as the orbits pass through the origin. We describe a procedure to connect different analytic solutions for successive orbits at the origin. We calculate the periods and precessions of these bound orbits, and graph a number of specific examples. Also, we explain why they all must violate the virial theorem. The unbound orbits are also discussed in detail. This includes the unusual orbits which have finite travel times to infinity and also the special $\nu = 2$ case.Comment: LaTeX, 27 pages with 12 figures available from the authors or can be generated from Mathematica instructions at end of the fil

Direct Singular Positions of the Parallel Manipulator Tricept

[[abstract]]In this article, the direct singular positions of the parallel manipulator Tricept are determined. An alternative 3 x 3 Jacobian matrix, simpler than the existing one, is obtained in this study. For a given moving platform's orientation, the determinant of this Jacobian matrix may be expressed as a cubic polynomial in moving platform's equation length. Direct singular positions may thus be obtained by solving cubic polynomial equations. For an arbitrarily chosen moving platform's orientation, there exists at least one moving platform's extension length that causes direct kinematic singularity. It is found that if moving platform's size is larger than a specific value, then within the moving platform's domain there exist two regions, in which direct kinematic singularities can only occur at positions impossible to reach.[[notice]]補正完畢[[booktype]]紙本[[countrycodes]]GB

Global-Scale Turbulent Convection and Magnetic Dynamo Action in the Solar Envelope

The operation of the solar global dynamo appears to involve many dynamical elements. Self-consistent MHD simulations which realistically incorporate all of these processes are not yet computationally feasible, though some elements can now be studied with reasonable fidelity. Here we consider the manner in which turbulent compressible convection within the bulk of the solar convection zone can generate large-scale magnetic fields through dynamo action. We accomplish this through a series of three-dimensional numerical simulations of MHD convection within rotating spherical shells using our ASH code on massively parallel supercomputers. Since differential rotation is a key ingredient in all dynamo models, we also examine here the nature of the rotation profiles that can be sustained within the deep convection zone as strong magnetic fields are built and maintained. We find that the convection is able to maintain a solar-like angular velocity profile despite the influence of Maxwell stresses which tend to oppose Reynolds stresses and thus reduce the latitudinal angular velocity contrast throughout the convection zone. The dynamo-generated magnetic fields exhibit a complex structure and evolution, with radial fields concentrated in downflow lanes and toroidal fields organized into twisted ribbons which are extended in longitude and which achieve field strengths of up to 5000 G. The flows and fields exhibit substantial kinetic and magnetic helicity although systematic hemispherical patterns are only apparent in the former. Fluctuating fields dominate the magnetic energy and account for most of the back-reaction on the flow via Lorentz forces. Mean fields are relatively weak and do not exhibit systematic latitudinal propagation or periodic polarity reversals as in the sun. This may be attributed to the absence of a tachocline.Comment: 55 pages (ApJ refereeing format), 15 figures (low res), published by ApJ on October 2004 (abstract slightly reduced in order to fit in 24 lines limit) see also Browning, Miesch, Brun & Toomre 2006, ApJL, 648, 157 (astro-ph/0609153) for the effect of a tachocline in organizing the mean field

In--out intermittency in PDE and ODE models

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.Comment: To be published in Chaos, June 2001, also available at http://www.eurico.web.co

Causal Relativistic Fluid Dynamics

We derive causal relativistic fluid dynamical equations from the relaxation model of kinetic theory as in a procedure previously applied in the case of non-relativistic rarefied gases. By treating space and time on an equal footing and avoiding the iterative steps of the conventional Chapman-Enskog --- CE---method, we are able to derive causal equations in the first order of the expansion in terms of the mean flight time of the particles. This is in contrast to what is found using the CE approach. We illustrate the general results with the example of a gas of identical ultrarelativistic particles such as photons under the assumptions of homogeneity and isotropy. When we couple the fluid dynamical equations to Einstein's equation we find, in addition to the geometry-driven expanding solution of the FRW model, a second, matter-driven nonequilibrium solution to the equations. In only the second solution, entropy is produced at a significant rate.Comment: 23 pages (CQG, in press

Phase transitions in self-gravitating systems and bacterial populations with a screened attractive potential

We consider a system of particles interacting via a screened Newtonian potential and study phase transitions between homogeneous and inhomogeneous states in the microcanonical and canonical ensembles. Like for other systems with long-range interactions, we obtain a great diversity of microcanonical and canonical phase transitions depending on the dimension of space and on the importance of the screening length. We also consider a system of particles in Newtonian interaction in the presence of a ``neutralizing background''. By a proper interpretation of the parameters, our study describes (i) self-gravitating systems in a cosmological setting, and (ii) chemotaxis of bacterial populations in the original Keller-Segel model