A magnetohydrodynamic model applied to the lower convective region in the Sun including the radial components of magnetic field and flow

In this work, some numerical solutions of magnetohydrodynamic equations are investigated in the presence of radial and azimuthal components of magnetic field with the use of previously developed algorithm. In this algorithm, the thin shell approximation and a special separation of variables is used to obtain the radial and latitudinal variations of physical parameters in spherical coordinates. The solutions are obtained via this separation of variables in the components of momentum transfer equation. The analysis yields three important parameters which are the sphericity, density and radial components shape parameters in the latitudinal distributions of physical variables. The magnetic field profile, used here, produces comparable magnetic fluxes found in previous works. There is a considerable change in density with respect to reference model. Other physical parameters also reveal important physical results. It is as well shown that the spherical symmetric distributions of physical parameters are broken for the region of study.     

Modelling the magnetic field and differential rotation effects for the vicinity of the base of convective zone in the Sun

In this work, some numerical solutions of magnetohydrodynamics equations are investigated in the presence of differential rotation with the use of previously developed algorithm. This algorithm includes the thin shell approximation and a special separation of variables which were used to obtain the radial and latitudinal variations of physical parameters in spherical coordinates. The magnetic field profile is chosen to produce comparable magnetic fluxes found in previous works. The sphericity and density shape parameters relevant to model is determined by using two different known differential rotation profiles. It is found that the shape of variations in physical parameters is strongly dependent to magnetic field profile and there is a considerable change in density with respect to reference model. It is as well shown that the spherical symmetric distributions of physical parameters are broken for the region of study.     

Self‐similar evolutionary solutions for an accreting magneto‐fluid around a compact object with finite electrical conductivity

In this paper, we investigate the time evolution of an accreting magneto‐fluid with finite conductivity. For the case of a thin disk, the fluid equations along with Maxwell's equations are derived in a simplified, one‐dimensional model that neglects the latitudinal dependence of the flow. The finite electrical conductivity of the plasma is taken into account by Ohm's law; however, the shear viscous stress is neglected, as well as the self‐gravity of the disk. In order to solve the integrated equations that govern the dynamical behaviour of the magneto‐fluid, we have used a self‐similar solution. We introduce two dimensionless variables, S₀ and ε_ϱ, which represent the size of the electrical conductivity and the density behaviour with time, respectively. The effect of each of these on the structure of the disk is studied. While the pressure is obtained simply by solving an ordinary differential equation, the density, the magnetic field, the radial velocity, and the rotational velocity are presented analytically. The solutions show that the S0 and ε_ϱ parameters affect the radial thickness of the disk. Also, radial velocity and gas pressure are more sensitive to the electrical conductivity in the inner regions of disk. Moreover, the parameter ε_ϱ has a more significant effect on the physical quantities for small radii. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modelling RR Lyrae pulsation: mission (im)possible?

In order to make an in-depth comparison between theory and observations, we analyse the light and velocity curves of various hydrodynamical models simulating RRab stars. The observations are represented by empirical formulae, derived in this and our earlier papers. It is shown that the overwhelming majority of the models tested do not follow the empirical relations regarding the shape of the light curves and the physical parameters. In almost all cases the luminosities predicted from the model light curves are significantly lower than the corresponding model values. The overall discrepancy of the models is an important indication of the limitation of the applicability of the present theoretical light and velocity curves in the determination of the physical parameters of these stars. In transforming the theoretical data to the observed light curves in V colour and in computing the observed radial velocities, it is shown that both bolometric correction and tracing the line-forming regions have considerable effects on the evaluation of the observed quantities. In an effort to resolve the discrepancy between theory and observations, it is suggested that a proper evaluation of the bolometric correction and radial velocity based on complete dynamical atmosphere models may be a useful step in this direction.     

Evolution of the Radial Abundance Gradient and Cold Gas of the Galactic Disk

In order to understand the forming mechanism of the radial abun- dance gradient of the Galactic disk and the evolution of cold gas, we have con- structed a chemical evolution model of the Galactic disk, in which the star for- mation law concerned with molecular hydrogens is adopted, and the evolution of mass surface density is calculated for the molecular and atomic hydrogens separately, then the model predictions and the observed radial distributions of some physical quantities are compared. The result indicates that the model prediction is sensitive to the adopted infall timescale, the model which adopts the star formation law concerned with the molecular hydrogens can agree well with the major observed properties of the Galactic disk, especially can obtain naturally the radial oxygen abundance gradient of the Galactic disk, and the radial surface density profile of cold gas. The assumption of instantaneous or non-instantaneous recycling approximation has a small effect on the evolution of cold gas, especially in the case of rather low gas density.     

Dynamics of the physical state during two-current-loop collisions

A model of two-current-loop collisions is presented to explain the impulsive nature of solar flares. From MHD equations considering the gravity and resistivity effects we find self-consistent expressions and a set of equations governing the behavior of all physical quantities just after magnetic reconnection has taken place. Numerical simulations have revealed that the most important parameters of the problem are the plasma and the ratio of initial values of pressure gradient in the longitudinal and radial directions. Thus, the low plasma case during aY-type interaction (initial longitudinal pressure gradient is comparable with initial radial pressure gradient) shows a rapid pinch and simultaneous enhancement of all physical quantities, including the electric field components, which are important for high-energy particle acceleration. However, an increase of the plasma causes a weakening of the pinch effect and a decrease of extreme values of all physical quantities. On the other hand, for anX-type collision (initial longitudinal pressure gradient is much greater than initial radial pressure gradient), which is able to provide a jet, the increase of the plasma causes a high velocity jet. As for aI-type collision (initial longitudinal pressure gradient is much less than initial radial pressure gradient) it shows neither jet production nor very strong enhancement of physical quantities. We also consider direct and oblique collisions, taking into account both cases of partial and complete reconnection.     

两维银河系旋臂对比度

为改善文献上惯用的表现银河系分子谱线巡视结果的完全平滑了方位信息的径向分布方法,我们发展了原子气体或分子云参量的分环银经分布图,(X—l)_i图,它在某种程度上给出了方位信息。用现存旋臂模型结合这种图我们得到的银道面旋臂区和臂间区的E(HI),E(CO),E(~(13)CO)和N/S(~(13)CO)的两维对比度约为1—2。     

First order latitude effects in the solar wind

The Weber-Davis model of the solar wind is generalized to include the effects of latitude. The principal assumptions of perfect electrical conductivity, rotational symmetry, a polytropic relation between pressure and density, and a flow aligned magnetic field in a system rotating with the Sun, are retained. A flow aligned magnetic field in the rotating system may be expressed in terms of the flow velocity and density. Rotational symmetry fixes the longitudinal flow velocity V_φ in terms of the flow in the r?θ plane. Thus, the original three dimensional magnetohydrodynamic flow problem is reduced to a two dimensional hydrodynamic flow problem in the r?θ plane.There are three critical surfaces associated with the equations which supply conditions to determine three of six required boundary conditions. The specified boundary conditions at the base of the corona are the temperature, density, and magnitude of the magnetic field. The equations are then expanded about the radial, nonrotating Parker solution and an analytic solution is obtained for the resulting first order equations. The results show that for constant coronal boundary conditions there is a latitudinal flow toward the solar poles, as a result of magnetic stresses, which persists out to large distances for the Sun. Associated with this flow is a latitudinal component of the magnetic field. The radial flow parameters are, to within small first order differences, in agreement with those of the Parker and the Weber-Davis models of the solar wind.The equations are further generalized to permit first order latitudinal variations in the specified coronal boundary conditions. Results at 1 a.u. are presented for 5 per cent latitudinal differences between the equatorial and polar values. These results show that the solution at 1 a.u. is most sensitive to a latitudinal dependence in the boundary temperature and least sensitive to a latitudinal dependence in the magnetic field magnitude.A solution is then obtained for an approximate dipolar variation in the coronal magnetic field magnitude. This solution predicts that the latitudinal flow is initially toward the Equator due to magnetic channeling; however, this effect is rapidly overcome and the latitudinal flow at 1 a.u. is toward the pole and not significantly different from the solution for constant boundary conditions.     

Differential rotation of solar filaments

The latitudinal component of solar differential rotation and the possibility of a radial component are discussed and compared to the observed rotational velocities of solar filaments. Our values of rotational rate versus heliographic latitude for 100 points in the solar atmosphere derived from 17 quiescent filaments are comparable to the rates found by d'Azambuja and d'Azambuja (1948). The filament rate is significantly greater than the spot rate (Newton and Nunn, 1951); the difference cannot be accounted for by the poleward migration of filaments and seems to reflect a true radial gradient of rotational velocity in the Sun. We show that filaments in closer proximity to active regions usually exhibit no differential rotation, while those far from active regions generally show it clearly. Comparison with Mt. Wilson photospheric Doppler measurements shows that filaments rotate faster than the general photosphere and that, as is well known, the spot rate exceeds that for the general photosphere.     

Corrugation Instability of a Coronal Arcade

We analyse the behaviour of linear magnetohydrodynamic perturbations of a coronal arcade modelled by a half-cylinder with an azimuthal magnetic field and non-uniform radial profiles of the plasma pressure, temperature, and the field. Attention is paid to the perturbations with short longitudinal (in the direction along the arcade) wavelengths. The radial structure of the perturbations, either oscillatory or evanescent, is prescribed by the radial profiles of the equilibrium quantities. Conditions for the corrugation instability of the arcade are determined. It is established that the instability growth rate increases with decreases in the longitudinal wavelength and the radial wave number. In the unstable mode, the radial perturbations of the magnetic field are stronger than the longitudinal perturbations, creating an almost circularly corrugated rippling of the arcade in the longitudinal direction. For coronal conditions, the growth time of the instability is shorter than one minute, decreasing with an increase in the temperature. Implications of the developed theory for the dynamics of coronal active regions are discussed.     

Photometric and spectroscopic investigation of the oscillating Algol type binary: EW Boo

We obtained the physical and geometrical parameters of the EW Boo system, which exhibits short period and small amplitude pulsations as well as brightness variations due to orbital motion of components. Towards this end we carried out photometric observations at Ankara University Kreiken Observatory (AUKO) as well as spectroscopic observations at TUBITAK National Observatory (TNO). The light and radial velocity curves obtained from these observations have been simultaneously analyzed with PHOEBE and the absolute parameters of the system along with the geometric parameters of the components have been determined. Using model light curves of EW Boo, light curve regions in which the pulsations are active have been determined and as a result of analyses performed in the frequency region, characteristic parameters of pulsations have been obtained. We find that the results are compatible with current parameters of similar systems in the literature. The evolutionary status of the components is propounded and discussed.     

Solar internal rotation rate and the latitudinal variation of the tachocline

A new set of accurately measured frequencies of solar oscillations is used to infer the rotation rate inside the Sun, as a function of radial distance as well as latitude. We have adopted a regularized least-squares technique with iterative refinement for both 1.5D inversion, using the splitting coefficients, and 2D inversion using individual m splittings. The inferred rotation rate agrees well with earlier estimates showing a shear layer just below the surface and another one around the base of the convection zone. The tachocline or the transition layer where the rotation rate changes from differential rotation in the convection zone to an almost latitudinally independent rotation rate in the radiative interior is studied in detail. No compelling evidence for any latitudinal variation in the position and width of the tachocline is found, although it appears that the tachocline probably shifts to a slightly larger radial distance at higher latitudes and possibly also becomes thicker. However, these variations are within the estimated errors and more accurate data would be needed to make a definitive statement about latitudinal variations.     

Hα measurements using DEFPOS/RTT150 Telescope： Instrumentation and observations

To investigate the physical properties of HII regions and some PNe about 4' in size, a DEFPOS Fabry-Perot spectrometer has been redesigned and set up at the (Antalya/Bakirlitepe, Turkey), The spectrometer has a 4' circular field of view and a velocity resolution of 27.27 km s- 1 (a spectral resolving power of 11000) over a 200 km s-1 spectral window near Hα. This work presents the details of the newly redesigned inphysical results of our new Hα observations selected from the Reynolds et al. (2005) and Fich et al. (1990) papers. The DEFPOS system has been fully tested and the first observations of HII regions and PNe in the Galaxy are used to illustrate the power of the instrument. We feel that our first Fabry-Perot observations can provide a powerful tool for the study of objects with small angular size. In the future, we are planning to prepare a catalog including some physical properties such as radial velocity, line width, and intensity of some PNe and HII regions close to the 4' field of view.     

The Alpha Effect and the Observed Twist and Current Helicity of Solar Magnetic Fields

We present a straightforward comparison of model calculations for the α-effect, helicities, and magnetic field line twist in the solar convection zone with magnetic field observations at atmospheric levels. The model calculations are carried out in a mixing-length approximation for the turbulence with a profile of the solar internal rotation rate obtained from helioseismic inversions. The magnetic field data consist of photospheric vector magnetograms of 422 active regions for which spatially-averaged values of the force-free twist parameter and of the current helicity density are calculated, which are then used to determine latitudinal profiles of these quantities. The comparison of the model calculations with the observations suggests that the observed twist and helicity are generated in the bulk of the convection zone, rather than in a layer close to the bottom. This supports two-layer dynamo models where the large-scale toroidal field is generated by differential rotation in a thin layer at the bottom while the α-effect is operating in the bulk of the convection zone. Our previous observational finding was that the moduli of the twist factor and of the current helicity density increase rather steeply from zero at the equator towards higher latitudes and attain a certain saturation at about 12 – 15^∘. In our dynamo model with algebraic nonlinearity, the increase continues, however, to higher latitudes and is more gradual. This could be due to the neglect of the coupling between small-scale and large-scale current and magnetic helicities and of the latitudinal drift of the activity belts in the model.     

Motion of charged particles in the magnetosphere

The adiabatic motion of charged particles in the magnetosphere has been investigated using Mead-Fairfield magnetospheric field model (Mead and Fairfield, 1975). Since the motion of charged particles in a dipolar field geometry is well understood, we bring out in this paper some important features in characteristic motion due to non-dipolar distortions in the field geometry. We look at the tilt averaged picture of the field configuration and estimate theoretically the parameters like bounce period, longitudinal invariant and the bounce averaged drift velocities of the charged particle in the Mead-Fairfield field geometry. These parameters are evaluated as a function of pitch angle and azimuthal position in the region of ring current (5 to 7 Earth radii from the centre of the Earth) for four ranges of magnetic activity. At different longitudes the non-dipolar contribution as a percentage of dipole value in bounce period and longitudinal invariant show maximum variation for particles close to 90° pitch angles. For any low pitch angle, these effects maximize at the midnight meridian. The radial component of the bounce averaged drift velocity is found to be greatest at the dawn-dusk meridians and the contribution vanishes at the day and midnight meridians for all pitch angles. In the absence of tilt-dependent terms in the model, the latitudinal component of the drift velocity vanishes. On the other hand, the relative non-dipolar contribution to bounce averaged azimuthal drift velocity is very high as compared to similar contribution in other characteristic parameters of particle motion. It is also shown that non-dipolar contribution in bounce period, longitudinal invariant and bounce averaged drift velocities increases in magnitude with increase in distance and magnetic activity.     

Rotation of solar structures in the upper chromosphere. II. Time variations in the latitudinal distribution of the rotation of active regions and coronal holes

The time variations in the latitudinal distribution of the rotation of active regions and coronal holes are investigated. The synoptic maps obtained from observations in the He I 1083 nm line at Kitt Peak Observatory over almost three solar cycles are used as observational data. A Fourier analysis of the time series constructed from synoptic maps has yielded the following results. The rotation of active regions differs significantly from the rotation of coronal holes in all parameters: the set of the most significant rotation periods, their latitudinal distribution, and time variations. The rotation of active regions and coronal holes is characterized by variations from cycle to cycle, a time-varying north-south asymmetry. The power spectra for consecutive cycles of solar activity differ significantly for both epochs of high activity and minima. Analysis of the total power of the spectra within four selected intervals of periods from 21 to 33 days has shown that the total power is highest in the intervals of periods 24–27 and 27–30 days. This is valid for both active regions and coronal holes. The correlation between the total powers in the above intervals of periods changes noticeably with time. Long-lived or successively appearing active regions with rotation periods in the range 24–30 days are typical of the time of a sharp decrease in the total equivalent width of active regions. This includes not only the decline time of the 11-year cycles, but also the minima between recurrent activity maxima during one cycle. A predominance of long-lived coronal holes as their total equivalent width decreases is noticeable for coronal holes with rotation periods in the interval 30–33 days. All of the above results suggest that the rotation of solar structures is determined mainly by the subphotospheric sources of specific structures, not by the rotation of the main volumes of solar plasma of the quiet Sun.     

Extragalactic jets on subpc and large scales

Jets can be probed in their innermost regions (d≲0.1 pc) through the study of the relativistically boosted emission of blazars. On the other extreme of spatial scales, the study of structure and dynamics of extragalactic relativistic jets received renewed impulse after the discovery, made by Chandra, of bright X-ray emission from regions at distances larger than hundreds of kpc from the central engine. At both scales it is thus possible to infer some of the basic parameters of the flow (speed, density, magnetic field intensity, power). After a brief review of the available observational evidence, I discuss how the comparison between the physical quantities independently derived at the two scales can be used to shed light on the global dynamics of the jet, from the innermost regions to the hundreds of kpc scale.     

Dynamo-generated magnetic fields at the surface of a massive star

Spruit has shown that an astrophysical dynamo can operate in the non-convective material of a differentially rotating star as a result of a particular instability in the magnetic field (the Tayler instability). By assuming that the dynamo operates in a state of marginal instability, Spruit has obtained formulae which predict the equilibrium strengths of azimuthal and radial field components in terms of local physical quantities. Here, we apply Spruit's formulae to our previously published models of rotating massive stars in order to estimate Tayler dynamo field strengths. There are no free parameters in Spruit's formulae. In our models of 10- and  50-M_⊙  stars on the zero-age main sequence, we find internal azimuthal fields of up to 1 MG, and internal radial components of a few kG. Evolved models contain weaker fields. In order to obtain estimates of the field strength at the stellar surface, we examine the conditions under which the Tayler dynamo fields are subject to magnetic buoyancy. We find that conditions for Tayler instability overlap with those for buoyancy at intermediate to high magnetic latitudes. This suggests that fields emerge at the surface of a massive star between magnetic latitudes of about 45° and the poles. We attempt to estimate the strength of the field which emerges at the surface of a massive star. Although these estimates are very rough, we find that the surface field strengths overlap with values which have been reported recently for line-of-sight fields in several O and B stars.     

Meteor radiant activity mapping using single-station radar observations

In this paper, we present a study of the general physical and chemical properties and radial velocity monitoring of young active stars. We derive temperatures, log  g , [Fe/H], v sin  i and R _spec values for eight stars. The detailed analysis reveals that the stars are not homogeneous in their principal physical parameters or in the age distribution. In 4/5, we found a periodic radial velocity signal which originates in surface features; the fifth is surprisingly inactive and shows little variation.     

The SHASTA Code Modified by Self-adaptive Mesh and Numerical Experiment of Magnetic Reconnections

SHASTA (Sharp and Smooth Transport Algorithm) is a code with single mesh to solve the 2-dimensional magnetohydrodynamic (MHD) equations. When SHASTA is used to the numerical simulation of magnetic reconnection problem, it is modified to be the code which adopts the method of the selfadaptive mesh. The modified code can carry out refined calculations in diffusion regions. In the process of the self-adaptive calculations with SHASTA, a “plugand-play” strategy is adopted and the original algorithm to solve 2-dimensional MHD partial differential equations is treated as an independent cell. In addition, the hierarchical data structure is used in this modification and parameters in each refined level are described by a 2-dimensional variable array. The regions where the distributions of magnetic field and pressure exhibit steep variations are marked as the refined regions. Then, the distributions of physical quantities and the boundary conditions in the grid points of refined levels are deduced via interpolation method. Finally, the refined calculated results of refined regions are assigned to the previous level of mesh and the existing results are updated. The numerical experiment of magnetic reconnections which adopts refined calculations indicates that compared with the code with single mesh, the resolution of details is improved and the corresponding increment of computing time is related to the selection of parameters in the simulation. The calculation accuracy and effect on instability, which are caused by a part of the self-adaptive code, depend on the boundary settings, push strategy over each single step as well as the interpolation algorithm.