绝热非径向太阳振动的数值计算方法

本文采用改进后的平行Shooting方法，即连续正交归一化方法，通过增加系统的秩和把微分方程从线性变成非线性的方法，使得病态方程变成非病态方程，克服了平行Shooting方法（不连续正交归一化方法）的主要困难：一是选择径向网格点数目的问题，如果被检验振动模的径向节点数大或小，所需要的网格点数则多或少；二是在远离阈值的时候，尚未有合适的方法提高计算精度的问题．数值计算出了标准太阳模型的绝热非径向太阳振动的理论本征频率谱．从数值计算的结果发现改进后的平行Shooting方法比线性化的Henyey方法和平行Shooting方法收敛性更好且精度更高．其误差可以达到1％的精度．     

绝热非径向太阳振动的数值计算方法

本采用改进后的平行Ｓｈｏｏｔｉｎｇ方法，即连续正交归一化方法，通过增加系统的秩和把微分方程从线性变成非线性的方法，使得病态方程变成非病态方程，克服了平行Ｓｈｏｏｔｉｎｇ方法（不连续正交归一化方法）的主要困难：一是选择径向网格点数目的问题，如果被检验振动模的径向节点数大或小，所需要的网格点数则多或少；二是在远离阈值的时候，尚未有合适的方法提高计算精度的问题。数值计算出了标准太阳模型的绝热非径向太阳     

太阳P模振动理论的研究

日震学是太阳物理的一个前沿分支学科,是根据太阳振动的观测来研究太阳的内部结构与运动的一种方法学。太阳５ｍｉｎ振动频率的理论计算和实测之间存在显著偏差和振动模的激发问题一直是困扰日震学的两大难题。经过多年的研究仍然没有解决。然而太阳表面层内绝热假设条件与真实情况有很大的偏差,我们认为绝大多数标准太阳模型的Ｐ模频率计算忽略了非绝热效应对频率的影响,忽略了振动的激发和衰减机制以及缺乏振动与对流湍流相互作     

太阳表面震波能量与磁场强度的关系

过去的一系列研究表明在太阳光球表面,随着纵向磁场的增强,对应区域所观测到的震波能量随之降低.为了进一步研究震波和磁场的关系,使用日震和磁成像仪(Helioseismic and Magnetic Imager,HMI)的数据,通过对视向磁场进行势场外推得到磁场总强度,与震波能量进行关联.对一定总磁场区域的p模式震波能量拟合可以得到其随总磁场强度而变化的关系.研究证实了震波能量与磁场总强度有更强关联.不同活动区的结果表明p模式震波能量随着磁场总强度增加而下降的趋势在500 Gs以下非常接近,在500 Gs至1 000 Gs范围内有略微差别.并且p模式震波能量的对数与总磁场强度有较强的线性相关性.对比迈克尔逊多普勒成像仪(Michelson Doppler Imager,MDI)和HMI得到的结果时发现,虽然MDI观测的磁场结果整体会比HMI观测结果偏大,但归一化p模式震波能量与磁场强度的关系分析对比结果相差不大.     

磁对流湍流对太阳P模振动的影响(英文)

为了进一步探索太阳对流区的物理性质 ,我们利用高精度N的日震观测数据来研究太阳内部扰动磁场对低阶太阳P模振动的影响。对于一个时间相关的MHD湍流源 ,我们给出了导致频率变化的各种可能性。如果只考虑磁扰动的贡献 ,不同值的振动模的频率变化仅只是涨落磁场能谱的函数。我们发现频率的变化随着太阳内部磁场强度增加而变大 ,并且和太阳活动周期密切相关。我们的计算表明太阳磁活动导致的频率变化可达 0 .3μΗz。     

非局部对流理论及其在太阳模型中的应用(英文)

在Li&Yang (2 0 0 1 )所给出的局部对流理论的基础上 ,我们进一步采用梯度型方案给出了非局部对流理论 ,并将它用于太阳模型中。这一理论考虑了恒星对流区内的非局部效应 ,它得到了一个与原来用混合长理论或局部理论给出的结果有所区别的对流区 ,扩散效应很明显。但是 ,目前我们的理论还不能处理时间相关的对流以及对流超射等问题。这些问题将在后续工作中加以考虑。当把这一理论应用于太阳模型中时 ,我们发现它对标准太阳模型的改正非常微小。我们讨论了这一现象 ,并对其加以解释     

太阳p模振动的激发机制

本估计了激发单个太阳ｐ模振动需要的能量输入率Ｅ≡ＥＴ,其中谱线宽度Г假设是振动模能量衰减率的观测值,在改进后的对流湍流三维时空分离描述的基础上,利用太阳对流区的混合长流模型计算了振动模的激发率,认为太阳Ｐ模振动主要是由雷诺应力扰动激发的,即Ｐ模激发对应着地流湍流运动引起声波发射,对于频率ｖ〈３ｍＨｚ,振动模的激发率为Ｅ∝ｖ＾７频率ｖ〉３ｍＨｚ,振动模的激发率为Ｅ∝ｖ＾－５,理论计算的振动能谱与观     

非绝热过程对太阳P模振动的影响

本讨论了与非绝热性有关的辐射损失和对流转移对太阳ｐ模振动的影响。在非绝热情况下，Ｐ模的本征频率增加虚σｉ和σ＾２（）ｉ。本试图探讨一种渐进方法研究非绝热产应对太阳Ｐ模振动的影响。在渐近近拟效的太阳外大气层，利用表面相移的相关关系给出了非绝热振动方程的严格解。     

非绝热过程对太阳p模振动的影响

本文讨论了与非绝热性有关的辐射损失和对流转移对太阳ｐ模振动的影响．在非绝热情况下，ｐ模的本征频率增加了虚部σ（１）ｉ和σ（２）ｉ．本文试图探讨一种渐进方法研究非绝热效应对太阳ｐ模振动的影响．在渐进近似失效的太阳外大气层，利用表面相移的相关关系给出了非绝热振动方程的严格解．对低、中间频率的振动模，通过渐进解和表面解在外大气层的拟合，得到表面相移只是频率的函数．与绝热振动相比，考虑非绝热效应有可能改善太阳５分钟振动的理论频率和观测频率之间存在的偏差．     

Modeling helioseismic modes in the magnetic atmosphere of the Sun

The pulsation of the solar surface is caused by acoustic waves traveling in the solar interior. Thorough analyses of observational data indicate that these f and p helioseismic oscillation modes are not bounced back completely at the surface but they partially penetrate into the atmosphere. Atmospheric effects and their possible observational application are investigated in one‐dimensional magnetohydrodynamic models. It is found that f and p mode frequencies are shifted of the order of μHz due to the presence of an atmospheric magnetic field. This shift varies with the direction of the wave propagation.Resonant coupling of global helioseismic modes to local Alfvén and slow waves reduce the life time of the global modes. The resulting line width of the frequency line is of the order of nHz, and it also varies with propagation angle. These features enable us to use helioseismic observations in magnetic diagnostics of the lower atmosphere. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Comparison of Large-Scale Flows on the Sun Measured by Time-Distance Helioseismology and Local Correlation Tracking

We present a direct comparison between two different techniques: time-distance helioseismology and a local correlation tracking method for measuring mass flows in the solar photosphere and in a near-surface layer. We applied both methods to the same dataset (MDI high-cadence Dopplergrams covering almost the entire Carrington rotation 1974) and compared the results. We found that, after necessary corrections, the vector flow fields obtained by these techniques are very similar. The median difference between directions of corresponding vectors is 24°, and the correlation coefficients of the results for mean zonal and meridional flows are 0.98 and 0.88, respectively. The largest discrepancies are found in areas of small velocities where the inaccuracies of the computed vectors play a significant role. The good agreement of these two methods increases confidence in the reliability of large-scale synoptic maps obtained by them.     

The Accuracy of Atmospheric Models in Solar Maximum Years Checked with the Data of CHAMP Accelerometer

A method of processing the data obtained with the CHAMP satelliteborne accelerometer is introduced. Through an analysis of the data the fault in the x-axis of the instrument is verified. We examined the accuracy of the atmospheric models in solar maximum years with the accelerometer data obtained from 2001 to 2002 and the inversion of the atmospheric density from the tangential non-gravitational acceleration. Then from a statistical analysis of the errors of the models, DTM94 and MSIS90, we obtained the following quantitative conclusions. During the solar maximum years the error of DTM94 is 30% to 35%, and that of MSIS90 is smaller at 25% to 30%. The greatest difference between the two models is about 10%.     

p‐mode power variation with solar atmosphere as observed in the Na D1 and K spectral lines

In this work we investigate p‐mode power variation with solar atmosphere. To this aim, we use THÉMIS observations of the Na D1 (λ 5896 Å) and K (λ 7699 Å) spectral lines. While the formation heights of the K spectral line are essentially located in the photospheric layer, the formation heights of the Na D1 line span a much wider region: from photosphere up to chromosphere. Hence, we had the opportunity to infer p‐mode power variation up to the chromospheric layer. By analyzing power spectra obtained by temporal series at different points of the Na D1 and K spectral lines, we confirm and quantify the increase in p‐mode power towards higher atmospheric layers. Furthermore, the large span in formation heights of the Na D1 line induces a larger enhancement of p‐mode power with solar atmosphere compared to the K spectral line. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Brightness fluctuations of the Sun and p‐mode oscillations: The inverse problem and nonadiabatic waves in the photosphere

Using photometric observations of the Sun as a star (DIFOS, SoHO) we were able to solve the inverse heloiseismic problem and determine the global time‐dependent relative temperature fluctuations as functions of the geometric height. This was done under the adiabatic assumption. A mathematical tool was developed to solve the inverse problem, which is ill‐posed. The calculations were done using the numerical software Matlab 7. The adiabatic solution shows signs of temperature waves in the lower photosphere, which agrees with calculations done by Rodríguez Hidalgo et al. (2001) and Stodilka (2011). (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

日冕中无衰减振荡的研究进展

日冕中冕环的无衰减横向振荡(简称无衰减振荡)自2012年被发现以来,受到了广泛的关注.无衰减振荡具有在日冕中广泛存在以及振幅无明显变小的特性,使之在解释日冕加热和冕震学诊断上都具有相当的潜力.总结了日冕中无衰减振荡的研究进展,包括观测研究得到的一系列结果、提出的理论和数值模型以及基于无衰减振荡进行冕震学诊断的一些尝试,并且展望了未来可以进一步开展的研究.     

Radiative Cooling in MHD Models of the Quiet Sun Convection Zone and Corona

We present a series of numerical simulations of the quiet-Sun plasma threaded by magnetic fields that extend from the upper convection zone into the low corona. We discuss an efficient, simplified approximation to the physics of optically thick radiative transport through the surface layers, and investigate the effects of convective turbulence on the magnetic structure of the Sun’s atmosphere in an initially unipolar (open field) region. We find that the net Poynting flux below the surface is on average directed toward the interior, while in the photosphere and chromosphere the net flow of electromagnetic energy is outward into the solar corona. Overturning convective motions between these layers driven by rapid radiative cooling appears to be the source of energy for the oppositely directed fluxes of electromagnetic energy.