МЕХАНИЗЪМ НА 160-мин

Ivan T. Ivanov⁽¹⁾, Valentina T. Taneva⁽²⁾ and Boris Komitov⁽³⁾

(1) Department of Physics and Biophysics, Medical Institute, Thracian University, 6000 Stara Zagora;

(3) Institute of Astronomy, Bulgarian Academy of Sciences, 6000, Stara Zagora, P.O. Box 179

Scientific conference with international participation. Stara Zagora. 2005. Proceedings, vol. 1, Technical sciencies, Chemistry and Physics. pp. 300 - 304

The visible diameter of Sun oscilates with a period of 160 min. The same type of periodicity is also found in a huge number of solar radiation parameters. To elucidate the origin of these longitudinal radial pulsations we have used the equation for the equilibrium of inner layers which, after linearization, turned into the harmonic oscilator equation. The latter equation allows radial pulsations whose period and wave length were calculated using regresion expresions for the gas presure and density in various layers. The radial pulsations originate at the surface of active zone and propagate til the litosphere, where they undergo full inner reflection producing undersurface stationary waves with a period of 150 min.

INTRODUCTION. Ever since the discovery of the 160-min radial pulsations of the Sun ample data have been collected suggesting their impact on the sollar physics [1,2]. 160-min oscilations are observed in the intensity of light, radio and infrared radiation emitted from the resting regions of Sun [3]. The pulsations observed in the radiobrightness lag behind the cyclic radial movements by about 12 min [3]. A 0.02% variation in the intensity of the 1.65 mm infrared sollar radiation is detected with a 160-min period [4]. Oscilations with the same period are established in the cyrcular polarized radio emition of Sun at 13.5 mm, which lag behind the oscilations of the radiant velocity by 34 minutes [5].

These observations could not be explained as pure atmospheric effects [6]. The radiant pulsations of Sun have been unsuccefully related to the particularities of the radiant energy transfer [7], to the possible interferense between the Sun gravity modes [8] or even to the possible existence of a particular object orbiting the Sun center at appr. 20 000 km under the surface of the Sun [9]. Thus at present there is no commonly accepted theory that could explain the radiant pulsations of the Sun. In this paper we offer a model based on the linearization of the Sun status equations. This model reprisents the radiant oscilations as staying waves generated on the surface of the active zone and propagating to the surface of the Sun where they have a period of about 160 minutes.

DESCRIPTION OF THE MODEL. The equilibrium status of the star core is usually found as a solution of the system of basic equations describing the star physics. For the Sun, several solutions are obtained which give the equlibrium values of pressure P, temperature T etc as functions of the distance r from solar center. These functions are averaged in the summarised model [10] which will be used futher.

The status of the Sun intestine is a stable equilibrium. If a thin concentric layer inside the Sun is displaced from its equilibrium distance r_o from the center of gravity, it will commence to oscilate. What will be the frequency of these oscilations? The equation that describes the movement of the layer is

where m is the mass of the gas placed under the layer. Let this layer be shifted from a place with distance r_o to a place with a distance r_o + Dr. Since Dr <<r_o, the pressure gradient will change as given by

The gravity force attracting the layer to the center of Sun will also change according to the variations in the distance r

Let us combine the equations (2) and (3) with the equation (1). Inserting a new variable

¶ ²x/¶ t² + ( 1/r_o.¶ ²P/¶ r² – 2.G.m / r_o³).x = 1/ r_o. r_o.(¶P/¶r) - G.m / r_o³ (4)

The right side of the latter equation is equal to zero since a hydrostatic equilibrium holds at r = r_o. Equalizing the left side of the equation to zero a new equation will be obtained which resembles the harmonic oscilator relation. Thus, the respective frequency f of the adiabatic radiant oscilations will be given by the formula

The period of oscilation of the layer will be П = 2 п /f. Once generated these oscilations will apparently subside down provided there is no source of energy for their rejuvanation.

The above linearizations appear fairly correct since the visible diameter of the Sun is found to oscilate with an amplitude of only a few km and a radiant velocity of about several m/s. Similar pulsations are found for the variable stars however the corresponding parameters of pulsations are three orders of magnitude greater in respect to these of Sun. Nevertheless, similar linearization theory has been also applied for the variable stars giving constant amplitude, frequency and phase at different distances from the star center [11]. In contrast to the oscilations of variable stars, the oscilations that could occur within the Sun are very small and should be consequently allowed to have different amplitude, frequency and phase at different distances from the center of Sun as given by the formula (5).

RESULRS. According to formula (5), the period П of oscilation of different layers could be calculated using proper expresions for the dependences of P, r and m on the relative distance x = r/R_¤.from the Sun center. Using the model [10] and 16 cited values for each parameter, these expressions were obtained as regression formulae and are shown in Table I. For each parameter the correlation coefficient K_r, reliability factor F and probability weight factor P_w between the cited data and the calculated values demonstrate that the obtained expressions are fairly satisfactory.

Table I. Regression formulae for the density p and gas presure P within the Sun as functions of the relative distance

y = r/R from its center. m is the mass of solar gas beneath a concentric layer with a radius r. R and M are the Sun radius and mass correspondingly.

Regression formula	Kr	F	Pw
m = M. (2.84 y - 1.70 y - 0.06)	0.991	52.1	4.10
P = exp (38.94 - 24.68 y )	0.92	6.08	0.0008
r = exp ( 3.85 -17.36 y )	0.942	8.28	0.00015

Allowing x to vary with a step of Dx = 0.01, a set of discrete values of П were calculated which are shown on plot (Fig.1). Generally, the intestine of Sun is sharply divided into three zones; active core, intermidiant zone and convective zone, according to the nature of the physical processes in them. Each of these zones is clearly distinguished on the given plot.

Fig.1. Period of radial pulsations within the Sun as a function of the distance from its center.

Oscilations should be apparently imposible throughout the active zone (0 <x< 0.38) as f² < 0. It could be assumed that the period of oscillations there will be close to infinity. In the vast convective zone of radiative energy transfer, oscilations might occur with a nearly constant period of about 1000 s. However, at the boundary between this zone and the active zone the period of pulsations sharply increased inclining to infinity (Fig.1) in accordance with the result that ascilations are not allowed in the active core. This result indicates that any slow mechanical disturbance that might occur on the surface of the active zone could be transfered into the interior of the middle zone as periodic pulsations. Hence, the radial oscilations that might occur within the Sun could originate from this boundary. Such a conclusion is apprehensible since the surface of the active zone must be mechanically unstable.

According to Fig.1, the period of pulsations slowly increases near the outside boundary of the convective zone reaching the the value of 8900 s (148 min) just on the Sun surface. This value differs by about 8 % from the experimentally obtained value of 160 min. The striking coincidence between the calculated and measured values of solar pulsations supports the proposed model.

DISCUSSIONS. According to the proposed model, the pulsations of the Sun radius are brought about by a slight disturbance of the dynamic equilibrium of forces balancing the inner layers and could be decsribed by the mechanism of the charmonic oscilator. The energy source generating the oscilations and maintaining their amplitude constant is possibly found at the boundary between the active and the middle zones. Once generated, the oscilations should further propagate towards the Sun surface with the velocity "C" of the sound. This velosity could be calculated using the formula C = (P/r)^1/2 and the expressions for the P, and р as given in Table I. We have calculated the velocity C and the wave length l = П.C of these mechanical pulsations as a function of the distance x form the center of Sun. Suprisingly l was nearly constant for both the middle and convective zones (data not shown).

When the radial oscilations reach the Solar atmosphere (the chromosphere), they should sustain a complete backward reflection without change in phase as suggested by the strong fulfilment of the respective mandatory condition l > 4пH. Here H is the depth of the solar atmosphere. As the wave length is practically constant within a broad segment under the visible surface of Sun, the formation of staying mechanical wave shoud be allowed deep under the chromosphere. These oscillations have the capacity to impact a huge number of physical phenomena close to the surface of Sun introducing a variable component with the same 160 min periodicity in the observed parameters of Sun.

1. Severny A. B., V. A. Kotov, T. T. Tsap. Nature (1976) Observations of solar pulsations. № 5539, v. 259, pp. 87-89

2. Scherrer P. H., J. M. Wilson (1983) Structure of the solar oscilations with period near 160 - minutes. Solar Phys., v.82, pp. 37-42

3. Kotov V. A., A. B. Severny, T. T. Tsap, I. G. Moiseev, V. A. Efanov, H. S. Nesterov. (1983) Manifestation of the 160-min solar oscilations in velocity and brightness (optical and radio observations). Solar Physics, v. 82 pp. 9-19

4. Kotov V.A., S. Kouchmy, G. Kouchmy (1983) Observation of global 160 -min infrared (differential) intensity variation of the Sun. Solar Physics, v. 82, pp. 21-35

5. Efanov V. A., V. A. Kotov, I. G. Moiseev, N. S. Nesterov, A. B. Severny (1983) Reports of Krim astrophysical observatory, v. 67, pp. 111-118

6. Hill H. A., T. J. Boss, T. P. Candel (1983) On the origin of oscilations in a solar diameter observed through the Barth`s atmosphere: a terrestrial atmospheric or a solar phemomenon. Solar Phisics, v. 82, pp. 129 - 138

7. Christensen - Dalsgaard J., S. Frandsen (1983) Radiative transfer and oscilations of the Sun. Solar Physics, v. 82, pp. 165-204

8. Dsiombovsky W. (1983) Resonant coupling between solar gravity modes. Solar Physics, v. 82, pp. 253-266

9. Blinnikov S. I., M. Yu. Khlopov (1983) Excitation of the solar oscilations by objects consisting of y-matter. Solar Physics, v. 82, pp. 383-385

10. Allen C. W. (1973) Astrophysical Quantities. University of London. The Athlon Press.

11. Cox J. P. In: Theory of stellar pulsations. Princeton, New Jersy, Princeton University Press, 1980

МЕХАНИЗЪМ НА 160-мин. РАДИАЛНИ ПУЛСАЦИИ НА СЛЪНЦЕТО

Иван Т. Иванов⁽¹⁾, Валентина Т. Танева ⁽²⁾, Борис Комитов ⁽³⁾

⁽¹⁾Катедра Физика и Биофизика, Медицински факултет на Тракийски университет, Стара Загора

⁽²⁾Институт по приложна физика, Технически Университет-София

⁽³⁾ Институт за космични изследвания, ЦЛКИ към БАН, Стара Загора

Научна конференция с международно участие. Стара Загора. 2005. Том 1, технически науки, Химия и Физика. с. 300 - 304.

Видимият диаметър на Слънцето трепти с период 160 минути. Същата периодичност е установена в много параметри на слънчевото излъчване. За обяснение на тези надлъжни радиални пулсации, използвахме уравнението за състояние на слънчевите недра, което след линеаризиране се превръща в уравнение на хармоничния осцилатор. Това уравнение позволява радиални пулсации, чийто период на трептене и дължина на вълната бяха пресметнати, използвайки регресионни изрази за газовото налягане и плътност за различните слоеве. Пулсациите започват от повърхността на активната зона и се разпространяват до литосферата, където търпят пълно вътрешно отражение, образувайки стоящи вълни под слънчевата повърхност с период около 150 мин.

Ключови думи: Слънце; 160-мин радиални пулсации; механични вълни; уравнение за състояние.

УВОД. След откриването на 160-мин. радиални колебания на Слънцето до днес са установени многочислени данни за тяхното влияние върху слънчевата физика [1,2]. 160 - мин. пулсации се наблюдават в амплитудата на радио- и инфрачервеното излъчване от спокойните участъци на Слънцето, а също и в излъчването в оптичната област [3]. Наблюдаваните колебания в радиояркостта на Слънцето с период 160 мин закъсняват спрямо 160 - мин радиални движения с около 12 мин [3]. В [4] се съобщава за наблюдение на 160-мин вариации в диференциалната интензивност на ИЧ - лъчението при 1.65 nm с амплитуда около 0.02 % от средната интензивност. Установени са и 160 -мин вариации в кръговополяризираното радиоизлъчване на Слънцето при 13.5 mm с амплитуда 0.003 % от средната, които закъсняват с около 34 мин спрямо колебанията на лъчевата скорост [5].

Тези наблюдения не могат да се обяснят с ефекти, свързани със земната атмосфера [6]. Все още липсва общоприета теория за произхода на радиалните колебания на Слънцето. Правени са опити те да се обяснят чрез особеностите на радиалния пренос на енергия [7], на основата на нелинейни ефекти - резонансна връзка между слънчевите гравитационни модове [8] и даже чрез допускане съществуването на специална у-планета, обикаляща на 20 000 km под Слънчевата повърхност [9]. По-долу ние предлагаме модел, основан на линеаризиране на уравненията за състоянието на Слънцето и представящ радиалните колебания като стоящи вълни, генерирани от повърхността на активната зона и отразяващи се от слънчевата повърхност, където те имат период 160 - мин.

ОПИСАНИЕ НА МЕХАНИЗМА. Равновестното състояние на звездните недра се намира като решение на система от основни уравнения, описващи физичните процеси и състоянието на звездната плазма. Известни са няколко решения на тази система, които дават равновестните стойности на налягането Р, плътността r, температурата Т и др., като функция на отстоянието от центъра на Слънцето r. Тези решения са осреднени в обобщения модел [10], който ще използваме по-нататък.

Състоянието на недрата на Слънцето представлява устойчиво равновесие. Ако тънък сферичен слой, намиращ се на разтояние r от центъра на Слънцето бъде изместен от своето равновесно положение, ще възникне затихващо трептене. Нека пресметнем честотата на този колебателен процес. Уравнението за движение на този слой ще бъде

където m е масата на веществото, намиращо се под слоя. Нека този слой е изместен адиабатно от равновесното си положение r. Ако отместването Dr е пренебрежимо малко спрямо r, изменението на градиента на налягането ще бъде

Поради промяната на отстоянието от центъра, силата на тежестта също ще се измени съгласно уравнението

Заместваме уравненията (2) и (3) в (1) и след въвеждане на нова променлива x = Dr/r_o получаваме:

¶ ²x/¶ t² + ( 1/r_o. ¶ ²P/¶ r² – 2.G.m / r_o³).x = 1/ r_o. r_o.(¶P/¶r) - G.m / r_o³ (4)

Дясната страна на това уравнение е равна на нула поради наличието на хидростатично равновесие при r₀ и приравнявайки на нула лявата страна, получаваме известното уравнение на хармоничния осцилатор. За честотата на радиалните адиабатични пулсации получаваме израза

а периодът на пулсациите ще бъде T = 2п / f . Несъмнено, тези пулсации ще бъдат затихващи, ако липсва източник на енергия за тяхното непрекъснато подновяване.

Подобна теория е развита за пулсирането на променливите звезди. При тях, амплитудата на трептене и скоростта на трептене са изключително големи [11], поради което е абсолютно задължително трептенията да имат еднаква амплитуда, честота и фаза в целия обем на звездата. При Слънцето, амплитудата на колебание на слънчевия диаметър е само няколко километра, а измененията на лъчевата скорост - няколко m/s, което е с около три порядъка по-малко спрямо пулсиращите звезди. Това позволява колебанията на отделните слоеве във вътрешността на Слънцето да се извършват с различна амплитуда, честота и фаза, както се дава от формула (5).

РЕЗУЛТАТИ. Периодът Т на радиалните пулсации за отделните слоеве на Слънцето може да се изчисли по формула (5). За целта предварително трябва да се намерят Р, r и m като функция на относителното разтояние у = r / R от центъра на Слънцето. Използвайки обобщения модел [10] и на основата на 16 цитирани стойности, за всяка една от тези величини намерихме регресионни изрази, представени в Таблица 1. Точността на тези изрази е добра, което се оценява по получените стойности на коефициента на корелация K_r, фактора на достоверност F и вероятностното тегло P_w.

Табл.1. Регресионни формули за плътността r и газовото налягане Р във вътрешността на Слънцето като функции на относителното отстояние y = r/ R_Оот центъра на Слънцето. С m е означена масата на тази част от Слънцето, която се намира под концентричен слой с радиус r. С R_О и M_О са означени радиусът и общата маса на Слънцето. Точността на регресионните формули се потвърждава от стойностите на коефициента на корелация Kr, фактора на достоверност F и вероятностното тегло Pw.

Регресионна формула	Kr	F	Pw
m = M. (2.84 y - 1.70 y - 0.06)	0.991	52.1	4.10
P = exp (38.94 - 24.68 y )	0.92	6.08	0.0008
r = exp ( 3.85 -17.36 y )	0.942	8.28	0.00015

Фиг. 1. Период T на радиалните пулсации във вътрешността на Слънцето като функция на относителното отстояние r/R_О от центъра. R_О е радиусът на Слънцето.

Периодът на радиалните пулсации Т е показан на Фиг. 1 като функция на отстоянието от центъра на Слънцето. Както е известно, вътрешността на Слънцето се разделя на три зони; активно ядро, междинна зона и конвективна зона, съобразно характера на процесите протичащи в тях. Тези три зони ясно се открояват и на показаната фигура. Навсякъде в обширната междинна зона могат да възникнат пулсации с почти еднакъв период от около 1000 s. На границата между активната и междинната зони периодът на пулсациите разте и клони към безкрайност. Вътре в активната зона (y ≤ 0.32) пулсации не могат да възникнат, тъй като се получава, че f ² < 0. Може да се приеме, че вътре в активната зона пулсациите ще имат период равен на безкрайност. На границата между активната зона и междинната зона се получават пулсации, при това с рязко нарастваща честота при слабо отдалечаване от междинната граница. Това може би показва, че пулсациите се пораждат, т.е., енергийният източник на пулсациите се намира на границата между активната и междинната зона. Физически това е допустимо, тъй като повърхността на активната зона трябва да се характеризира с висока нестабилност. В зоната на конвективен пренос, периодът на пулсациите бавно разте и на самата повърхност на Слънцето достига до 8900 s. Тази стойност се отличава с около 8 % от експериментално установения период от 160 мин, което подкрепя справедливостта на представения модел.

ОБСЪЖДАНЕ НА РЕЗУЛТАТИТЕ. Според този модел, радиалните пулсации на Слънцето възникват като хармонични колебания на вътрешните слоеве, дължащи се на нарушаване на динамичното равновесие на силите. Енергийният източник на тези нарушения вероятно се намира на повърхността на активната зона, откъдето те се разпространяват навън със скоростта на звука. Скоростта на разпространение на звука С във вътрешните слоеве може да се пресметне по формулата C = (P/r)^1/2, изпозвайки получените регресионни зависимости за налягането Р и плътността r на плазмата в отделните слоеве от Табл. 1 (не е дадено). Очевидно, произведението Т.С = l ще бъде дължината на получената механична вълна l. Замествайки получените стойности на Т и С се получава неочакваният и много интересен резултат, че дължината на вълната на тези механични пулсации ще бъде практически една и съща от мястото, където те възникват (повърхността на активната зона) до крайната точка на своето разпространение (повърхността на Слънцето), т.е., l = Т . С @ const.

Достигайки до повърхността на Слънцето, механичните пулсации със сигурност търпят пълно отражение от слънчевата атмосфера (хромосферата) без смяна на фазата и образуват стоящи вълни, тъй като условието за отражение и образуване на стоящи вълни (l > 4пH ) при тези условия се изпълнява извънредно силно. В тази формула H е височината на еднородния слой (хромосферата), от който настъпва отражението. Тъй като l е практически една и съща за отделните слоеве, стоящата вълна може да се образува не само непосредствено под хромосферата, но и на значителна дълбочина във вътрешните слоеве на Слънцето. Следователно, отразените от повърхността механични вълни няма да гасят идващите отдолу, а ще ги усилват. Образуваната по този начин стояща вълна може да засегне в силна степен всички физични процеси в близост до повърхността на Слънцето и да индуцира появата на една променлива компонента с период около 130 мин във техните параметри. Това може да обясни и посочената по-горе променлива компонента в параметрите на голям брой наблюдавани процеси, както и на самото пулсиране на видимия диаметър на Слънцето.