首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
We carry out the adaptive mesh refinement (AMR) implementation of our solar–interplanetary space-time conservation element and solution element (CESE) magnetohydrodynamic model (SIP–CESE MHD model) using a six-component grid system (Feng, Zhou, and Wu, Astrophys. J. 655, 1110, 2007; Feng et al., Astrophys. J. 723, 300, 2010). By transforming the governing MHD equations from the physical space (x,y,z) to the computational space (ξ,η,ζ) while retaining the form of conservation (Jiang et al., Solar Phys. 267, 463, 2010), the SIP–AMR–CESE MHD model is implemented in the reference coordinates with the aid of the parallel AMR package PARAMESH available at http://sourceforge.net/projects/paramesh/ . Meanwhile, the volumetric heating source terms derived from the topology of the magnetic-field expansion factor and the minimum angular separation (at the photosphere) between an open-field foot point and its nearest coronal-hole boundary are also included. We show the preliminary results of applying the SIP–AMR–CESE MHD model for simulating the solar-wind background of different solar-activity phases by comparison with SOHO observations and other spacecraft data from OMNI. Our numerical results show overall good agreements in the solar corona and in interplanetary space with these multiple-spacecraft observations.  相似文献   

2.
In this paper, the Space–Time Conservation Element and Solution Element (CESE) method is applied to 2.5-dimensional resistive magnetohydrodynamics (MHD) equations in Cartesian coordinates, with the purpose of modeling the magnetic reconnection study. To show the validity and capacity of its application to MHD reconnection problem, spontaneous fast reconnection and magnetic reconnection in multiple heliospheric current sheets are studied, which show good consistency with those obtained formerly by other authors. In order to assess the ∇ ⋅ B = 0 constraint numerically, the contours and evolution of ∇ ⋅ B are analyzed. The numerical results tell us that the CESE numerical scheme not only has good numerical resolution but also can keep the divergence-free condition for magnetic fields in the reconnection problems during the evolutionary process without any special treatment.  相似文献   

3.
R. P. Kane 《Solar physics》2007,243(2):205-217
For many purposes (e.g., satellite drag, operation of power grids on Earth, and satellite communication systems), predictions of the strength of a solar cycle are needed. Predictions are made by using different methods, depending upon the characteristics of sunspot cycles. However, the method most successful seems to be the precursor method by Ohl and his group, in which the geomagnetic activity in the declining phase of a sunspot cycle is found to be well correlated with the sunspot maximum of the next cycle. In the present communication, the method is illustrated by plotting the 12-month running means aa(min ) of the geomagnetic disturbance index aa near sunspot minimum versus the 12-month running means of the sunspot number Rz near sunspot maximum [aa(min ) versus Rz(max )], using data for sunspot cycles 9 – 18 to predict the Rz(max ) of cycle 19, using data for cycles 9 – 19 to predict Rz(max ) of cycle 20, and so on, and finally using data for cycles 9 – 23 to predict Rz(max ) of cycle 24, which is expected to occur in 2011 – 2012. The correlations were good (∼+0.90) and our preliminary predicted Rz(max ) for cycle 24 is 142±24, though this can be regarded as an upper limit, since there are indications that solar minimum may occur as late as March 2008. (Some workers have reported that the aa values before 1957 would have an error of 3 nT; if true, the revised estimate would be 124±26.) This result of the precursor method is compared with several other predictions of cycle 24, which are in a very wide range (50 – 200), so that whatever may be the final observed value, some method or other will be discredited, as happened in the case of cycle 23.  相似文献   

4.
The exact nonlinear cylindrical solution for incompressible Hall – magnetohydrodynamic (HMHD) waves, including dissipation, essentially from electron – neutral collisions, is obtained in a uniformly rotating, weakly ionized plasma such as exists in photospheric flux tubes. The ω – k relation of the waves, called here Hall – MHD waves, demonstrates the dispersive nature of the waves, introduced by the Hall effect, at large axial and radial wavenumbers. The Hall – MHD waves are in general elliptically polarized. The partially ionized plasma supports lower frequency modes, lowered by the factor δ≡ratio of the ion mass density to the neutral particle mass density, as compared to the fully ionized plasma (δ=1). The relation between the velocity and the magnetic field fluctuations departs significantly from the equipartition found in Alfvén waves. These short-wavelength and arbitrarily large amplitude waves could contribute toward the heating of the solar atmosphere.  相似文献   

5.
Some 8000 images obtained with the Solar Eclipse Coronal Imaging System (SECIS) fast-frame CCD camera instrument located at Lusaka, Zambia, during the total eclipse of 21 June 2001 have been analysed to search for short-period oscillations in intensity that could be a signature of solar coronal heating mechanisms by MHD wave dissipation. Images were taken in white-light and Fe xiv green-line (5303 ?) channels over 205 seconds (frame rate 39 s−1), approximately the length of eclipse totality at this location, with a pixel size of four arcseconds square. The data are of considerably better quality than those that we obtained during the 11 August 1999 total eclipse (Rudawy et al.: Astron. Astrophys. 416, 1179, 2004), in that the images are much better exposed and enhancements in the drive system of the heliostat used gave a much improved image stability. Classical Fourier and wavelet techniques have been used to analyse the emission at 29 518 locations, of which 10 714 had emission at reasonably high levels, searching for periodic fluctuations with periods in the range 0.1 – 17 seconds (frequencies 0.06 – 10 Hz). While a number of possible periodicities were apparent in the wavelet analysis, none of the spatially and time-limited periodicities in the local brightness curves was found to be physically important. This implies that the pervasive Alfvén wave-like phenomena (Tomczyk et al.: Science 317, 1192, 2007) using polarimetric observations with the Coronal Multi-Channel Polarimeter (CoMP) instrument do not give rise to significant oscillatory intensity fluctuations.  相似文献   

6.
The Fresnel Interferometric Imager has been proposed to the European Space Agency (ESA) Cosmic Vision plan as a class L mission. This mission addresses several themes of the CV Plan: Exoplanet study, Matter in extreme conditions, and The Universe taking shape. This paper is an abridged version of the original ESA proposal. We have removed most of the technical and financial issues, to concentrate on the instrumental design and astrophysical missions. The instrument proposed is an ultra-lightweight telescope, featuring a novel optical concept based on diffraction focussing. It yields high dynamic range images, while releasing constraints on positioning and manufacturing of the main optical elements. This concept should open the way to very large apertures in space. In this two spacecraft formation-flying instrument, one spacecraft holds the focussing element: the Fresnel interferometric array; the other spacecraft holds the field optics, focal instrumentation, and detectors. The Fresnel array proposed here is a 3.6 ×3.6 m square opaque foil punched with 105 to 106 void “subapertures”. Focusing is achieved with no other optical element: the shape and positioning of the subapertures (holes in the foil) is responsible for beam combining by diffraction, and 5% to 10% of the total incident light ends up into a sharp focus. The consequence of this high number of subapertures is high dynamic range images. In addition, as it uses only a combination of vacuum and opaque material, this focussing method is potentially efficient over a very broad wavelength domain. The focal length of such diffractive focussing devices is wavelength dependent. However, this can be corrected. We have tested optically the efficiency of the chromatism correction on artificial sources (500 < λ < 750 nm): the images are diffraction limited, and the dynamic range measured on an artificial double source reaches 6.2 10 − 6. We have also validated numerical simulation algorithms for larger Fresnel interferometric arrays. These simulations yield a dynamic range (rejection factor) close to 10 − 8 for arrays such as the 3.6 m one we propose. A dynamic range of 10 − 8 allows detection of objects at contrasts as high as than 10 − 9 in most of the field. The astrophysical applications cover many objects in the IR, visible an UV domains. Examples are presented, taking advantage of the high angular resolution and dynamic range capabilities of this concept.  相似文献   

7.
Quasi-periodic pulsations (QPP) of microwave emission generated in single flaring loops observed with the Nobeyama Radioheliograph (NoRH) and Nobeyama Radio Polarimeters (NoRP) are studied. Specific features of the time profiles, i.e. the visible presence or absence of QPPs, are not accounted for in the selection. The time evolution of the periods of the QPPs is examined using wavelet and correlation analyses. In ten out of twelve considered events, at least one or more significant spectral components with periods from 5 – 60 s have been found. The quality of the oscillations is rather low: Q=π N, where N is the number of cycles, mostly varies in the range 12 to 40, with an average of 25. We suggest that the detected QPPs can be classified into four types: i) those with stable mean periods (e.g. of 15 – 20 s or 8 – 9 s, the prevailing type); ii) those with spectral drift to shorter periods (mostly in the rise phase of the microwave emission); iii) those with drift to longer periods (mostly in the decay phase); iv)  those with multiple periods showing an X-shaped drift (e.g. in the range from 20 – 40 s in the rise phase).  相似文献   

8.
We performed a detailed analysis of 27 slow coronal mass ejections (CMEs) whose heights were measured in at least 30 coronagraphic images and were characterized by a high quality index (≥4). Our primary aim was to study the radial evolution of these CMEs and their properties in the range 2 – 30 solar radii. The instantaneous speeds of CMEs were calculated by using successive height – time data pairs. The obtained speed – distance profiles [v(R)] are fitted by a power law v = a(Rb) c . The power-law indices are found to be in the ranges a=30 – 386, b=1.95 – 3.92, and c=0.03 – 0.79. The power-law exponent c is found to be larger for slower and narrower CMEs. With the exception of two events that had approximately constant velocity, all events were accelerating. The majority of accelerating events shows a v(R) profile very similar to the solar-wind profile deduced by Sheeley et al. (Astrophys. J. 484, 472, 1997). This indicates that the dynamics of most slow CMEs are dominated by the solar wind drag.  相似文献   

9.
This contribution is a follow-up to the recent paper of Kuznetsov et al. (Contrib. Astron. Obs. Skalnaté Pleso 36, 85, 2006) on the ground level enhancement (GLE) on 20 January 2005. We focused on a study of Forbush decrease (FD) of 17 – 18 and 21 – 22 January 2005, respectively. The data from the neutron monitor at Lomnicky Štít (1 min counts) and from the Geomagnetic Observatory in Hurbanovo, both in Slovakia, were used as the basis for our investigation. The data on magnetic field and solar wind from GOES 10 and 12, SOHO-CELIAS, ACE and WIND satellites were used for better understanding of the global evolution of the event. The magnetic field is transformed to the RTN (Radial – Tangential – Normal) system where only the disturbed part of the field is compared, i.e., daily variations and a constant part are subtracted. The field reduction method is described. Our results are temporal vector diagrams of variation of all parameters at all positions from where we used the data. The amplitudes of |B| exceed 100 nT and variations during the arrival of the wavefront of CME take place simultaneously at the ground-based station and at GOES satellites. The character of the variations is as if there would be regions with the dominant electric charge of opposite signs, or electric currents with different orientations in the CME. On the basis of the values v p and n p and using certain assumptions we determined the mass of CME on 17 January and 21 January, respectively, of 1012 kg. A decrease of the cosmic ray level runs suddenly (during 10 minutes), starting, however, about two hours after a sudden change of the magnetic field.  相似文献   

10.
We compare the number and characteristics of interplanetary coronal mass ejections (ICMEs) to those of magnetic clouds (MCs) by using in-situ solar wind plasma and magnetic field observations made at 1 AU during solar cycle 23. We found that ≈ 28% of ICMEs appear to contain MCs, since 103 magnetic clouds (MCs) occurred during 1995  – 2006, and 307 ICMEs occurred during 1996 – 2006. For the period between 1996 and 2006, 85 MCs are identified as part of ICMEs, and six MCs are not associated with ICMEs, which conflicts with the idea that MCs are usually a subset of ICMEs. It was also found that solar wind conditions inside MCs and ICMEs are usually similar, but the linear correlation between geomagnetic storm intensity (Dst min ) and relevant solar wind parameters is better for MCs than for ICMEs. The differences between average event duration (Δt) and average proton plasma β (〈β〉) are two of the major differences between MCs and ICMEs: i) the average duration of ICMEs (29.6 h) is 44% longer than for MCs (20.6 hours), and ii) the average of 〈β〉 is 0.01 for MCs and 0.24 for ICMEs. The difference between the definition of a MC and that for an ICME is one of the major reasons for these average characteristics being different (i.e., listed above as items i) and ii)), and it is the reason for the frequency of their occurrences being different.  相似文献   

11.
Recent numerical investigations of wave propagation near coronal magnetic null points (McLaughlin and Hood: Astron. Astrophys. 459, 641, 2006) have indicated how a fast MHD wave partially converts into a slow MHD wave as the disturbance passes from a low-β plasma to a high-β plasma. This is a complex process and a clear understanding of the conversion mechanism requires the detailed investigation of a simpler model. An investigation of mode conversion in a stratified, isothermal atmosphere with a uniform, vertical magnetic field is carried out, both numerically and analytically. In contrast to previous investigations of upward-propagating waves (Zhugzhda and Dzhalilov: Astron. Astrophys. 112, 16, 1982a; Cally: Astrophys. J. 548, 473, 2001), this paper studies the downward propagation of waves from a low-β to high-β environment. A simple expression for the amplitude of the transmitted wave is compared with the numerical solution.  相似文献   

12.
G. Michalek 《Solar physics》2010,261(1):107-114
A set of 106 limb CMEs which are wide and could be possible halo events, when directed towards Earth, are used to check the accuracy of the asymmetric cone model. For this purpose characteristics of CMEs (widths and radial speeds) measured for the possible halo CMEs are compared with these obtained for halo CMEs using the asymmetric cone model (Michalek, Solar Phys. 237, 101, 2006). It was shown that the width and speed distributions for both datasets are very similar and with a probability of p>0.93 (using the Kolmogorov – Smirnov test) were drawn from the same distribution of events. We also determined the accurate relationship between radial (V rad) and expansion (V exp) speeds of halo CMEs. This relation for the halo CMEs is simply V rad=V exp and could be very useful for space weather application.  相似文献   

13.
We present the first in-depth statistical survey of flare source heights observed by RHESSI. Flares were found using a flare-finding algorithm designed to search the 6 – 10 keV count-rate when RHESSI’s full sensitivity was available in order to find the smallest events (Christe et al. in Astrophys. J. 677, 1385, 2008). Between March 2002 and March 2007, a total of 25 006 events were found. Source locations were determined in the 4 – 10 keV, 10 – 15 keV, and 15 – 30 keV energy ranges for each event. In order to extract the height distribution from the observed projected source positions, a forward-fit model was developed with an assumed source height distribution where height is measured from the photosphere. We find that the best flare height distribution is given by g(h)∝exp (−h/λ) where λ=6.1±0.3 Mm is the scale height. A power-law height distribution with a negative power-law index, γ=3.1±0.1 is also consistent with the data. Interpreted as thermal loop-top sources, these heights are compared to loops generated by a potential-field model (PFSS). The measured flare heights distribution are found to be much steeper than the potential-field loop height distribution, which may be a signature of the flare energization process.  相似文献   

14.
Zhao and Kosovichev (Astrophys. J. 591, 446, 2003) found two opposite sub-photospheric vortical flows in the depth range of 0 – 12 Mm around a fast rotating sunspot. So far there is no theoretical model explaining such flow motions. In this paper, we try to explain this phenomenon from the point of view of magnetic flux tubes interacting with large-scale vortical motions of plasma. In the deeper zone under the photosphere, the magnetic force may be less than the nonmagnetic force of plasma. The vortical flow located there twists the flux tube and magnetic free energy is built up in the tube. In the shallower zone under the photosphere, the magnetic force may be greater than the nonmagnetic force. Thus, part of the stored magnetic free energy is released to drive the plasma to rotate in two opposite directions, e.g., in the depth ranges of 0 – 3(5) and 9 – 12 Mm. In addition, we also define a vector of nonpotential magnetic stress τ, which can be related to flare occurrence. It is calculated for the active region NOAA 10930 on 11 December 2006. We find that: i) the integral of its line-of-sight (LOS) stress successively increases around the magnetic neutral line (MNL) prior to and during the flare and decreases to a minimum after the flare; ii) the integral of its transverse stress exceeds the integral of its LOS component by one order of magnitude over the whole field of view; iii) the transverse stress first points toward the MNL, then along it, and finally it points away from it. We need other data to verify whether or not the magnetic energy is transported in the horizontal direction to the neutral line, and then partly changes into the energy in LOS direction before and during the flare.  相似文献   

15.
R. P. Kane 《Solar physics》2008,249(2):369-380
The sunspot number series at the peak of sunspot activity often has two or three peaks (Gnevyshev peaks; Gnevyshev, Solar Phys. 1, 107, 1967; Solar Phys. 51, 175, 1977). The sunspot group number (SGN) data were examined for 1997 – 2003 (part of cycle 23) and compared with data for coronal mass ejection (CME) events. It was noticed that they exhibited mostly two Gnevyshev peaks in each of the four latitude belts 0° – 10°, 10° – 20°, 20 ° – 30°, and > 30°, in both N (northern) and S (southern) solar hemispheres. The SGN were confined to within latitudes ± 50° around the Equator, mostly around ± 35°, and seemed to occur later in lower latitudes, indicating possible latitudinal migration as in the Maunder butterfly diagrams. In CMEs, less energetic CMEs (of widths < 71°) showed prominent Gnevyshev peaks during sunspot maximum years in almost all latitude belts, including near the poles. The CME activity lasted longer than the SGN activity. However, the CME peaks did not match the SGN peaks and were almost simultaneous at different latitudes, indicating no latitudinal migration. In energetic CMEs including halo CMEs, the Gnevyshev peaks were obscure and ill-defined. The solar polar magnetic fields show polarity reversal during sunspot maximum years, first at the North Pole and, a few months later, at the South Pole. However, the CME peaks and gaps did not match with the magnetic field reversal times, preceding them by several months, rendering any cause – effect relationship doubtful.  相似文献   

16.
We report solar flare plasma to be multi-thermal in nature based on the theoretical model and study of the energy-dependent timing of thermal emission in ten M-class flares. We employ high-resolution X-ray spectra observed by the Si detector of the “Solar X-ray Spectrometer” (SOXS). The SOXS onboard the Indian GSAT-2 spacecraft was launched by the GSLV-D2 rocket on 8 May 2003. Firstly we model the spectral evolution of the X-ray line and continuum emission flux F(ε) from the flare by integrating a series of isothermal plasma flux. We find that the multi-temperature integrated flux F(ε) is a power-law function of ε with a spectral index (γ)≈−4.65. Next, based on spectral-temporal evolution of the flares we find that the emission in the energy range E=4 – 15 keV is dominated by temperatures of T=12 – 50 MK, while the multi-thermal power-law DEM index (δ) varies in the range of −4.4 and −5.7. The temporal evolution of the X-ray flux F(ε,t) assuming a multi-temperature plasma governed by thermal conduction cooling reveals that the temperature-dependent cooling time varies between 296 and 4640 s and the electron density (n e) varies in the range of n e=(1.77 – 29.3)×1010 cm−3. Employing temporal evolution technique in the current study as an alternative method for separating thermal from nonthermal components in the energy spectra, we measure the break-energy point, ranging between 14 and 21±1.0 keV.  相似文献   

17.
Using nine years of solar wind plasma and magnetic field data from the Wind mission, we investigated the characteristics of both magnetic clouds (MCs) and magnetic cloud-like structures (MCLs) during 1995 – 2003. A MCL structure is an event that is identified by an automatic scheme (Lepping, Wu, and Berdichevsky, Ann. Geophys. 23, 2687, 2005) with the same criteria as for a MC, but it is not usually identifiable as a flux rope by using the MC (Burlaga et al., J. Geophys. Res. 86, 6673, 1981) fitting model developed by Lepping, Jones, and Burlaga (Geophys. Res. Lett. 95(11), 957, 1990). The average occurrence rate is 9.5 for MCs and 13.6 for MCLs per year for the overall period of interest, and there were 82 MCs and 122 MCLs identified during this period. The characteristics of MCs and MCL structures are as follows: (1) The average duration, Δt, of MCs is 21.1 h, which is 40% longer than that for MCLs (Δt=15 h); (2) the average (minimum B z found in MC/MCL measured in geocentric solar ecliptic coordinates) is −10.2 nT for MCs and −6 nT for MCLs; (3) the average Dstmin  (minimum Dst caused by MCs/MCLs) is −82 nT for MCs and −37 nT for MCLs; (4) the average solar wind velocity is 453 km s−1 for MCs and 413 km s−1 for MCLs; (5) the average thermal speed is 24.6 km s−1 for MCs and 27.7 km s−1 for MCLs; (6) the average magnetic field intensity is 12.7 nT for MCs and 9.8 nT for MCLs; (7) the average solar wind density is 9.4 cm−3 for MCs and 6.3 cm−3 for MCLs; and (8) a MC is one of the most important interplanetary structures capable of causing severe geomagnetic storms. The longer duration, more intense magnetic field and higher solar wind speed of MCs, compared to those properties of the MCLs, are very likely the major reasons for MCs generally causing more severe geomagnetic storms than MCLs. But the fact that a MC is an important interplanetary structure with respect to geomagnetic storms is not new (e.g., Zhang and Burlaga, J. Geophys. Res. 93, 2511, 1988; Bothmer, ESA SP-535, 419, 2003).  相似文献   

18.
We summarize the response of the galactic cosmic ray (CGR) intensity to the passage of the more than 300 interplanetary coronal mass ejections (ICMEs) and their associated shocks that passed the Earth during 1995 – 2009, a period that encompasses the whole of Solar Cycle 23. In ∼ 80% of cases, the GCR intensity decreased during the passage of these structures, i.e., a “Forbush decrease” occurred, while in ∼ 10% there was no significant change. In the remaining cases, the GCR intensity increased. Where there was an intensity decrease, minimum intensity was observed inside the ICME in ∼ 90% of these events. The observations confirm the role of both post-shock regions and ICMEs in the generation of these decreases, consistent with many previous studies, but contrary to the conclusion of Reames, Kahler, and Tylka (Astrophys. J. Lett. 700, L199, 2009) who, from examining a subset of ICMEs with flux-rope-like magnetic fields (magnetic clouds) argued that these are “open structures” that allow free access of particles including GCRs to their interior. In fact, we find that magnetic clouds are more likely to participate in the deepest GCR decreases than ICMEs that are not magnetic clouds.  相似文献   

19.
The solar wind quasi-invariant (QI) has been defined by Osherovich, Fainberg, and Stone (Geophys. Res. Lett. 26, 2597, 1999) as the ratio of magnetic energy density and the energy density of the solar wind flow. In the regular solar wind QI is a rather small number, since the energy of the flow is almost two orders of magnitude greater than the magnetic energy. However, in magnetic clouds, QI is the order of unity (less than 1) and thus magnetic clouds can be viewed as a great anomaly in comparison with its value in the background solar wind. We study the duration, extent, and amplitude of this anomaly for two groups of isolated magnetic clouds: slow clouds (360<v<450 km s−1) and fast clouds (450≤v<720 km s−1). By applying the technique of superposition of epochs to 12 slow and 12 fast clouds from the catalog of Richardson and Cane (Solar Phys. 264, 189, 2010), we create an average slow cloud and an average fast cloud observed at 1 AU. From our analysis of these average clouds, we obtain cloud boundaries in both time and space as well as differences in QI amplitude and other parameters characterizing the solar wind state. Interplanetary magnetic clouds are known to cause major magnetic storms at the Earth, especially those clouds which travel from the sun to the Earth at high speeds. Characterizing each magnetic cloud by its QI value and extent may help in understanding the role of those disturbances in producing geomagnetic activity.  相似文献   

20.
Results are presented from a study of various sunspot contrast parameters in broadband red (672.3 nm) Cartesian full-disk digital images taken at the San Fernando Observatory (SFO) over eight years, 1997 – 2004, of the twenty-third sunspot cycle. A subset of over 2700 red sunspots was analyzed and values of average and maximum sunspot contrast as well as maximum umbral contrast were compared to various sunspot parameters. Average and maximum sunspot contrasts were found to be significantly correlated with sunspot area (r s=− 0.623 and r s=− 0.714, respectively). Maximum umbral contrast was found to be significantly correlated with umbral area (r s=− 0.535). These results are in agreement with the works of numerous other authors. No significant dependence was detected between average contrast, maximum contrast, or maximum umbral contrast during the rising phase of the solar cycle (r s=0.024, r s=0.033, and r s=0.064, respectively). During the decay phase, no significant correlation was found between average contrast or maximum contrast and time (r s=− 0.057 and r s=0.009, respectively), with a weak dependence seen between maximum umbral contrast and cycle (r s=0.102).  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号