SOHO contribution to the understanding of mass supply and flows in the solar corona

We expect a variety of dynamic phenomena in the quiescent non-flaring corona. Plasma flows, such as siphon flows or convective flows of chromospheric material evaporating into the corona, are expected whenever a pressure differences is established either between the footpoints or between the coronal and chromospheric segments of a coronal loop. Such flows can induce phenomena of spatial and temporal brightness variability of the corona. In particular, evaporation induces a net mass input into the corona and consequently coronal density enhancements. Flows are also expected in the regions where energy is released during magnetic reconnection. From the observational point of view the dynamics of the solar atmosphere has been investigated in great detail mostly in the lower transition region with the HRTS, and during flares with theSolar Maximum Mission andYohkoh. The high spectral, temporal and spatial resolution of theSOHO ultraviolet spectrometers should enable us in the near future to fill the gap providing a continuous coverage from the chromosphere to the corona, in the 10⁴–10⁶ K domain, and therefore to best study the dynamics throughout the solar atmosphere.     

High resolution UV solar spectroscopy

The advantages of high resolution UV spectroscopy for the investigation of the solar atmosphere are stressed while the limitations in the areas of instrumentation and diagnosis are discussed. The recent achievements (made essentially by Skylab, OSO-8 and rocket instruments) are reviewed and discussed.It is shown that high resolution UV solar spectroscopy has improved our knowledge of the dynamics of the upper layers of the solar atmosphere. Within the present instrument capabilities the birth of coronal expansion is shown to take place at the top of the transition region. The existence of downward flows over the bright regions of the network is evidenced from redshifts or transition region and chromospheric optically thin lines: velocities as large as 22 km s^-1 have been measured in O vi. Short period waves (95 s) have been detected in lines of Si ii at chromospheric levels in addition to the well known 300s and 180 s photospheric and chromospheric oscillations. There is strong evidence that optically thin chromospheric and transition region lines are broadened by a nonthermal velocity component which is maximum at 1.3 × 10⁵ K and decreases at higher temperatures. This may indicate the presence of unresolved acoustic or magnetohydrodynamic waves so oftenly set fourth as the source of chromospheric and coronal heating.Contradictions between the various results are pointed out and discussed. They might be attributed to the different angular resolution of the instruments, a key parameter for future space observations. It is suggested that the Solar Optical Telescope (SOT) and the Grazing Incidence Solar Telescope (GRIST) which are presently under phase A studies at NASA and ESA be considered as a tandem of instruments to fly on Spacelab in the 1980's. Both their angular and spectral resolution appear sufficient to resolve most of the problems under discussion today.Review presented at the Vth Conference on UV and X-ray Spectroscopy of Astrophysical and Laboratory Plasmas, London, July 4–7, 1977.     

Self-consistent Simulations of Alfv��n Wave Driven Winds from the Sun and Stars

We review our recent results of Alfvén wave-driven winds. First, we present the result of self-consistent 1D MHD simulations for solar winds from the photosphere to interplanetary region. Here, we emphasize the importance of the reflection of Alfvén waves in the density stratified corona and solar winds. We also introduce the recent Hinode observation that might detect the reflection signature of transverse (Alfvénic) waves by Fujimura and Tsuneta (Astrophys. J. 702:1443, 2009). Then, we show the results of Alfvén wave-driven winds from red giant stars. As a star evolves to the red giant branch, the properties of stellar winds drastically change from steady coronal winds to intermittent chromospheric winds. We also discuss how the stellar evolution affects the wave reflection in the stellar atmosphere and similarities and differences of accretion disk winds by MHD turbulence.     

Chromospheric and coronal heating mechanisms II

We review the mechanisms which are thought to provide steady heating of chromospheres and coronae. It appears now fairly well established that nonmagnetic chromospheric regions of latetype stars are heated by shock dissipation of acoustic waves which are generated in the stellar surface convection zones. In the case of late-type giants there is additional heating by shocks from pulsational waves. For slowly rotating stars, which have weak or no magnetic fields, these two are the dominant chromospheric heating mechanisms.Except for F-stars, the chromospheric heating of rapidly rotating late-type stars is dominated by magnetic heating either through MHD wave dissipation (AC mechanisms) or through magnetic field dissipation (DC mechanisms). The MHD wave and magnetic field energy comes from fluid motions in the stellar convection zones. Waves are also generated by reconnective events at chromospheric and coronal heights. The high-frequency part of the motion spectrum leads to AC heating, the low frequency part to DC heating. The coronae are almost exclusively heated by magnetic mechanisms. It is not possible to say at the moment whether AC or DC mechanisms are dominant, although presently the DC mechanisms (e.g., nanoflares) appear to be the more important. Only a more detailed study of the formation of and the dissipation in small-scale structures can answer this question.The X-ray emission in early-type stars shows the presence of coronal structures which are very different from those in late-type stars. This emission apparently arises in the hot post-shock regions of gas blobs which are accelerated in the stellar wind by the intense radiation field of these stars.     

A CDS Observation of the Relationship Between a Coronal Hole and Chromospheric Structure

This paper studies the relationship between chromospheric structure and the location of the coronal hole boundary, using a SOHO-CDS observation made on 13 June 1996. The cell structure visible in the oxygen and helium data is smoothed and a basin algorithm is used to find the individual cells. As the smoothing scale increases, cells merge to form larger ones. A trade off between merging and smoothing was made, a scale of 5,000 km being used. There is good agreement between the edges of cells and the coronal hole boundary, as determined in Mg IX 368 Å, although it is not perfect. This preliminary analysis shows that the coronal hole boundary is located along chromospheric network cell edges in most places, as would be expected if the chromosphere is the source of the coronal structure and solar wind.     

The dynamics of chromospheric spicules

We present the observational results on chromospheric spicules obtained at the Sayan observatory 50 cm coronograph. To investigate the evolution of chromospheric spicules, we analysed spicule spectra of strong chromospheric lines measured simultaneously at three altitudes above the solar limb during 5–60 min with a time resolution of 10 to 20 s. The spatial resolution was better than 1, and the spectral resolution was 0.03Å in 6563Å. The appearance of a spicule at a given altitude is preceded by an sharp increase in line-of-sight velocity and/or in line half-width at a lower level. Generally, the evolution has a non-monotonous impulsive character. Changes of line-of-sight velocities and other parameters of the line profile can be represented as the superposition of slow, evolutionary changes and fluctuations with periods of about 80 to 120 s. The amplitude of line-of-sight velocity fluctuations is 2–3 km/sec and tends to increase with height. By studying the phase delays of the fluctuations at different heights, we found that the propagation velocity exceeds 300 km s^–1, and that the disturbances do not necessarily propagate upwards.     

Solar flares and particle acceleration

Energy release in solar flares occurs during the impulsive phase, which is a period of a few to about ten minutes, during which energy is injected into the flare region in bursts with durations of various time scales, from a few tens of seconds down to 0.1 s or even shorter. Non-thermal heating is observed during a short period, not longer than a few minutes, in the very first part of the impulsive phase; in average flares, with ambient particle densities not larger than a few times 10¹⁰ cm^–3 it is due to thick-target electron beam injection, causing chromospheric ablation followed by convection. In flares with larger densities the heating is due to thermal fronts (Section 1). The average energy released in chromospheric regions is a few times 10³⁰ erg, and an average number of 10³⁸ electrons with E 15 keV is accelerated. In subsecond pulses these values are about 10³⁵ electrons and about 10²⁷ erg per subsecond pulse. The total energy released in flares is larger than these values (Section 2). Energization occurs gradually, in a series of fast non-explosive flux-thread interactions, on the average at levels about 10⁴ km above the solar photosphere, a region permeated by a large number ( 10) of fluxthreads, each carrying electric currents of 10¹⁰–10¹¹ A. The energy is fed into the flare by differential motions of magnetic fields driven by photospheric-chromospheric movements (Section 3). In contrast to these are the high-energy flares, characterized by the emission of gamma-radiation and/or very high-frequency (millimeter) radiobursts. Observations of such flares, of the flare neutron emission, as well as the observation of ³He-rich interplanetary plasma clouds from flares all point to a common source, identified with shortlived ( 0.1 s) superhot ( 10⁸ K) flare knots, situated in chromospheric levels (Section 4). Pre-flare phenomena and the existence of homologous flares prove that flare energization can occur repeatedly in the same part of an active region: the consequent conclusions are that only seldom the full energy of an active region is exhausted in one flare, or that the flare energy is generated anew between homologous flares; this latter case looks more probable (Section 5). Flare energization requires the formation of direct electric fields, in value comparable with, or somewhat smaller than the Dreicer field (Section 6). Such fields originate by current-thread reconnection in a regime in which the current sheet is thin enough to let resistive instability originate (Section 7). Particle acceleration occurs by fast reconnection in magnetic fields 100 G and electric fields exceeding about 0.3 times the Dreicer field at fairly low particle densities ( 10¹⁰ cm^–3); for larger densities plasma heating is expected to occur (Section 8). Transport of accelerated particles towards interplanetary space demands a field-line configuration open to space. Such a configuration originates mainly after the gradual gamma-ray/proton flares, and particularly after two-ribbon flares; these flares belong to the dynamic flares in Sturrock and vestka's flare classification. Acceleration to GeV energies occurs subsequently in shock waves, probably by first-order Fermi acceleration (Section 9).     

Hydrodynamic shock wave formation in the solar chromosphere and corona during flares

Basic mechanisms of the hydrodynamic shock wave formation in the solar atmosphere during flares are considered. Hydrodynamic plasma flows during flares arise due to fast energy release which is accumulated in the magnetic field of currents in the solar atmosphere. Shock waves arise as a result of rapid heating of the chromospheric upper layers from accelerated particles or heat fluxes. Powerful hydrodynamic phenomena can also arise due to explosive current sheet disruption in the region of strong magnetic field reconnection. Fundamental questions of shock wave formation and propagation in a non-homogeneous emitting solar atmosphere are discussed.An invited paper presented at STIP Workshop on Shock Waves in the Solar Corona and Interplanetary Space, 15–19 June, 1980, Smolenice, Czechoslovakia.     

Structure of reconnection layers in the magnetosphere

Magnetic reconnection can lead to the formation of observed boundary layers at the dayside magnetopause and in the nightside plasma sheet of the earth's magnetosphere. In this paper, the structure of these reconnection layers is studied by solving the one-dimensional Riemann problem for the evolution of a current sheet. Analytical method, resistive MHD simulations, and hybrid simulations are used. Based on the ideal MHD formulation, rotational discontinuities, slow shocks, slow expansion waves, and contact discontinuity are present in the dayside reconnection layer. Fast expansion waves are also present in the solution of the Riemann problem, but they quickly propagate out of the reconnection layer. Our study provides a coherent picture for the transition from the reconnection layer with two slow shocks in Petschek's model to the reconnection layer with a rotational discontinuity and a slow expansion wave in Levy et al's model. In the resistive MHD simulations, the rotational discontinuities are replaced by intermediate shocks or time-dependent intermediate shocks. In the hybrid simulations, the time-dependent intermediate shock quickly evolves to a steady rotational discontinuity, and the contact discontinuity does not exist. The magnetotail reconnection layer consists of two slow shocks. Hybrid simulations of slow shocks indicate that there exists a critical number,M _c, such that for slow shocks with an intermediate Mach numberM _IM _c, a large-amplitude rotational wavetrain is present in the downstream region. For slow shocks withM _I<M _c, the downstream wavetrain does not exist. Chaotic ion orbits in the downstream wave provide an efficient mechanism for ion heating and wave damping and explain the existence of the critical numberM _c in slow shocks.     

高压内冷涡轮的气热耦合计算(英文)

本文目的在于验证所开发的适用于涡轮气热耦合模拟的计算程序,并研究转捩对气热耦合计算结果影响。首先在流场求解程序HIT-3D中加入考虑转捩影响的q-ω低雷诺数二方程模型以及AGS代数转捩模型模块,然后采用直接耦合方法关联HIT-3D与自主开发的温度场求解程序,使新程序具备气热耦合求解功能。选取MarkⅡ叶片的三个不同试验工况作为验证算例,在计算中考虑管内流动影响,并对流道内流动分别采用了B-L代数模型、q-ω二方程模型以及B-L&AGS模型,而对管内流动则分别采用了B-L模型与B-L&AGS模型。计算表明采用各模型预测的压力分布与试验吻合较好,而在层流转捩区域采用B-L&AGS模型预测的温度分布与实验吻合最好,而在湍流流动区域,各模型预测的温度分布接近。这一方面表明所开发程序具备较准确进行内冷涡轮气热耦合计算的能力,另一方面也证明了考虑转捩影响对提高气热耦合计算精度的重要性。     

Magnetohydrodynamic simulation of a streamer beside a realistic coronal hole

Existing models of coronal streamers establish their credibility and act as the initial state for transients. The models have produced satisfactory streamer simulations, but unsatisfactory coronal hole simulations. This is a consequence of the character of the models and the boundary conditions. The models all have higher densities in the magnetically open regions than occur in coronal holes (Noci,et al., 1993).     

Thermodynamical Properties of the ICM from Hydrodynamical Simulations

Modern hydrodynamical simulations offer nowadays a powerful means to trace the evolution of the X-ray properties of the intra-cluster medium (ICM) during the cosmological history of the hierarchical build up of galaxy clusters. In this paper we review the current status of these simulations and how their predictions fare in reproducing the most recent X-ray observations of clusters. After briefly discussing the shortcomings of the self-similar model, based on assuming that gravity only drives the evolution of the ICM, we discuss how the processes of gas cooling and non-gravitational heating are expected to bring model predictions into better agreement with observational data. We then present results from the hydrodynamical simulations, performed by different groups, and how they compare with observational data. As terms of comparison, we use X-ray scaling relations between mass, luminosity, temperature and pressure, as well as the profiles of temperature and entropy. The results of this comparison can be summarised as follows: (a) simulations, which include gas cooling, star formation and supernova feedback, are generally successful in reproducing the X-ray properties of the ICM outside the core regions; (b) simulations generally fail in reproducing the observed “cool core” structure, in that they have serious difficulties in regulating overcooling, thereby producing steep negative central temperature profiles. This discrepancy calls for the need of introducing other physical processes, such as energy feedback from active galactic nuclei, which should compensate the radiative losses of the gas with high density, low entropy and short cooling time, which is observed to reside in the innermost regions of galaxy clusters.     

Alfvén Waves in the Solar Atmosphere

Alfvén waves are considered to be viable transporters of the non-thermal energy required to heat the Sun’s quiescent atmosphere. An abundance of recent observations, from state-of-the-art facilities, have reported the existence of Alfvén waves in a range of chromospheric and coronal structures. Here, we review the progress made in disentangling the characteristics of transverse kink and torsional linear magnetohydrodynamic (MHD) waves. We outline the simple, yet powerful theory describing their basic properties in (non-)uniform magnetic structures, which closely resemble the building blocks of the real solar atmosphere.     

Simulation Techniques for Cosmological Simulations

Modern cosmological observations allow us to study in great detail the evolution and history of the large scale structure hierarchy. The fundamental problem of accurate constraints on the cosmological parameters, within a given cosmological model, requires precise modelling of the observed structure. In this paper we briefly review the current most effective techniques of large scale structure simulations, emphasising both their advantages and shortcomings. Starting with basics of the direct N-body simulations appropriate to modelling cold dark matter evolution, we then discuss the direct-sum technique GRAPE, particle-mesh (PM) and hybrid methods, combining the PM and the tree algorithms. Simulations of baryonic matter in the Universe often use hydrodynamic codes based on both particle methods that discretise mass, and grid-based methods. We briefly describe Eulerian grid methods, and also some variants of Lagrangian smoothed particle hydrodynamics (SPH) methods.     

Solar mass flow in fine-scale structures

Observations of transient and steady velocities at chromospheric, transition region and coronal temperatures in the quiet Sun and coronal holes are reviewed. The relevance of fine-scale structures in governing the mass balance of the solar atmosphere is stressed. At present, a coherent picture of these mass flows does not exist. However, the current observational base of transition region and coronal velocity information is limited but should greatly improve with measurements from the SOHO satellite.     

The Structure of the Chromosphere Properties Pertaining to Element Fractionation

We review the structure and dynamics of the solar chromosphere with emphasis on the quiet Sun and properties that are relevant to element fractionation mechanisms. Attention is given to the chromospheric magnetic field, its connections to the photosphere, and to the dynamical evolution of the chromosphere. While some profound advances have been made in the “unmagnetized” chromosphere, our knowledge of the magnetically controlled chromosphere, more relevant for the discussion of element fractionation, is limited. Given the dynamic nature of the chromosphere and the poorly understood magnetic linkage to the corona, it is unlikely that we will soon know the detailed processes leading to FIP fractionation. This revised version was published online in June 2006 with corrections to the Cover Date.     

Cosmic Ray Transport in a Heliospheric Magnetic Field with Non-Polar Coronal Holes

The simple tilted dipole picture of Corotating Interaction Regions which prevailed during the first polar pass of Ulysses no longer applies since the Sun entered a more active phase. Recent observations show that CIRs still persist, though the large polar coronal holes of solar minimum shrink to smaller areas and move to lower latitudes. We present 3-D simulations for the cosmic-ray intensity variations in a model with non-polar high speed streams. Latitudinal and recurrent time-variations are discussed, but more detailed and realistic simulations are required before quantitative comparisons with observations can be made. This revised version was published online in August 2006 with corrections to the Cover Date.     

RANS analysis of the low-Reynolds number flow around the SD7003 airfoil

This paper presents a numerical analysis of the incompressible flow at Reynolds number 6.0×¹⁰⁴ around the Selig–Donovan 7003 airfoil. The airfoil performances have been computed by the Reynolds averaged Navier–Stokes equations and large eddy simulations. The airfoil stall and preliminary post-stall have been obtained by both the methods. Some limitations of the RANS turbulence models for low-Reynolds number flows have been overcome by the κ–ω SST-LR model, a recent modification of the well-known SST model. Large-eddy simulations have also been performed for a more detailed analysis of the results. The relevance in the stall mechanism of the laminar separation bubble arising on the airfoil is highlighted. The stall occurs when the laminar bubble present in the leading edge zone and a separated region forming on the central part of the airfoil join together. The κ–ω SST-LR model returns the same stall mechanism as the large eddy simulation. Flows at low-Reynolds numbers can be simulated by the RANS methods, but the choice of the turbulence model is crucial. The κ–ω SST-LR model has provided results in good agreement with the large eddy simulation and the available experimental data.